c c------------------------------------------------------------------ c Begin file /home/sting/foster/timegcm/mksrc/hdrsum.f c------------------------------------------------------------------ c subroutine hdrsum c c Write header and summary to unit lutime c include 'mksrc.h' include 'tgcmhdr.h' dimension ihead(512) c c Load header: c read(lu,end=900) iter,nday,nhr,nmin,label, c + date,output,start,stp,hist,sav,step,mag,difhor,iuivi, c + sdtide,ipower,aurora,dispos,data,source,sourct,dtide, c + dum,rdate,naurp,hp,cp,byimf c ihead(1) = iter x = iter*step + 1.e-4 ihead(2) = x/86400. x = x-float(ihead(2))*86400. ihead(3) = x/3600. x = x-float(ihead(3))*3600. ihead(4) = x/60. c do 50 i=1,ninpt c50 ihead(i+4) = inpt(i) write(6,"('hdrsum: ihead(1-4)=',4i5)") (ihead(i),i=1,4) c c Write header: buffer out(lutime,1) (ihead,ihead(512)) if (unit(lutime)) 100,105,105 c c Write summary: 100 continue buffer out(lutime,1) (f,f(100)) if (unit(lutime)) 110,115,115 c 110 continue return c c Errors in buffer outs: 105 write(6,"('hdrsum: error ',i3,' writing header')") unit(lutime) stop 'header' 115 write(6,"('hdrsum: error ',i3,' writing summary')") unit(lutime) stop 'summary' c end c