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,*)
call flush(34)
enddo