c program tiegcm_p7 include "tiegcm_p7.h" integer ihead(512) real summary(100) character*80 errmsg c c Copy a p7 dynamo volume (tiegcm with Zp -7 to +7) to a new c volume and write summary info so rdtgcm will recognize it c as a p7 volume. This code for tiegcm with ZST=+7 only. c (use job-stepping for multiple volumes) c call getinp c c Loop through all histories until EOF is encountered: c 10 continue read(lurd,end=900) ihead write(6,"('Read time ',i2,':',i2,':',i2,' from ', + a,' and copy to ',a)") (ihead(i),i=2,4), + rdmss(1:lenstr(rdmss)), wrmss(1:lenstr(wrmss)) write(luwr) ihead read(lurd) summary call wrsumm ! write summary info for new history c c Copy data for each lat slice of current history: c do j=1,jmx read(lurd) f write(luwr) f enddo goto 10 900 continue c c Dispose new mss volume: c write(6,"('EOF encountered on lurd=',i2)") lurd close(lurd) close(luwr) lwrmss = lenstr(wrmss) if (lwrmss.gt.0) then write(6,"('Disposing file ',a,' ...')") wrmss(1:lwrmss) call mswrite(ier,'wrfile',wrmss,',ECRIDLEY',367,' ') if (ier.eq.0) then write(6,"('Disposed mss file ',a)") wrmss(1:lwrmss) else call mserror(errmsg) write(6,"('>>> Error from mswrite of disk file wrfile ', + ' to mspath ',a)") wrmss(1:lwrmss) write(6,"(' Error message: ',a)") errmsg endif endif stop end