wham gfortran single chain corrections

[unres.git] / source / wham / src / wham_calc1.F

diff --git a/source/wham/src/wham_calc1.F b/source/wham/src/wham_calc1.F
index 41f42bc..a15186d 100644 (file)
--- a/source/wham/src/wham_calc1.F
+++ b/source/wham/src/wham_calc1.F
@@ -1444,9 +1444,16 @@ c        write (iout,*) "i",i," potEmin",potEmin
         write (iout,*) 
         write (34,'(f7.1,2f15.5,$)') startGridT+i*delta_T,
      &   sumW(i,iparm),sumE(i,iparm)
-        write (34,'(f10.5,$)') (sumQ(j,i,iparm),j=1,nQ+2)
-        write (34,'(e15.5,$)') sumEsq(i,iparm)-sumEbis(i,iparm),
-     &   (sumQsq(j,i,iparm),j=1,nQ+2),(sumEQ(j,i,iparm),j=1,nQ+2)
+        do j=1,nQ+2
+         write (34,'(f10.5,$)') sumQ(j,i,iparm)
+        enddo
+        write (34,'(e15.5,$)') sumEsq(i,iparm)-sumEbis(i,iparm)
+        do j=1,nQ+2
+          write (34,'(e15.5,$)') sumQsq(j,i,iparm)
+        enddo
+        do j=1,nQ+2
+          write (34,'(e15.5,$)') sumEQ(j,i,iparm)
+        enddo
         write (34,*) 
       enddo
       close(34)