program form c Reads unformatted output files y11res etc. and makes c formatted file itg.resf c include 'itg.par' c c ********* merge up to 7 output (-res) files ******** c c.... if newg=1, then introduce newgrid size and mode numbers. c Right now, only decrease mode numbers. c include 'post.cmn' namelist/datlis/ld,ne *,tim,shr,qsf,epsn,etai,rmu1,rmu2,rkap,aex,rmu2pa,rkpara *,dt0,dt,x0,y0,z0,nfile,kyspec,conplot,tkap,qkap,qkappar namelist/wmod1/newg,ld1,nmin1,nmax1,mmin1,mmax1 c initialization call needed to enable restarting on the NMFECC CRAY: call dropfile(0) !cray-only open(unit=7,file='itg.resf',status='new') open(unit=12,file='form.in',status='old') c c open all of the existing resuts files (y11res-y17res), otherwise c skip down to line 7. (At least itg.res must exist): c c bill this is awkward, but I just want to get it working right now c open(unit=5,file='itg.res',status='old',form='unformatted') open(unit=8,file='y12res',status='old',form='unformatted',err=7) open(unit=9,file='y13res',status='old',form='unformatted',err=7) open(unit=10,file='y14res',status='old',form='unformatted',err=7) open(unit=11,file='y15res',status='old',form='unformatted',err=7) open(unit=13,file='y16res',status='old',form='unformatted',err=7) open(unit=14,file='y17res',status='old',form='unformatted',err=7) 7 continue c input files: c c 5 itg.res a copy of itg.res, the results file from an ITG run. c 8-11,13-14 y12res-y17res additional results files to be concatenated c to the end of y11res. c c output files: c c 7 itg.resf a concatenation (unformatted) c of the results files itg.res, y12res-y17res c implement LABELS and and REPEAT options for namelists: call ddiopon('ibm') ! cray-only c allow for comments in namelists which begin with "!": call ddilist('!',22,1) ! cray-only c c read form.in c read(12,datlis) read(12,wmod1) c nd1=nmax1-nmin1+1 ldb1=ld1-1 c GWH: there was a (probably insignificant) bug in the original c version of the code, mmin(n) was nmin(n), and ndd1 was not c defined as an array!: c c ndd1(n)=mmax(n)-mmin(n)+1 c do 10 n=1,nd1 c ndd1(n)=mmax(n)-mmin(n)+1 c 10 continue c nff=0 net=0 ncnt=0 c nff=nff+1 if(nff.gt.nfile) goto 1000 netb=net+1 ncntb=ncnt+1 read(5)ld,ldb,nmin,nmax,nd,lin,ne,nfreq,ntotal,md write(*,*) ld,ldb,nmin,nmax,nd,lin,ne,nfreq,ntotal,md read(5)(mmin(n),n=1,nd) write(*,*) mmin(1) read(5)(mmax(n),n=1,nd) write(*,*) mmax(1) c read(5)tim,shr,etai,rmu1,rmu2,rkap,aex,rmu2pa,rkpara c mbeer changed this line for the better read(5)tim,shr,qsf,epsn,etai,rmu1,rmu2,rkap,tiovte write(*,*) tim,shr,qsf,epsn,etai c read(5)dt0,dt,x0,y0,z0,rkkmult c mbeer changed this line read(5) dt0,dt,x0,y0,z0 write(*,*)dt0,dt,x0,y0,z0 net=net+ne read(5)(((ws(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((wc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((us(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((uc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((tpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((tparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((prs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((prc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((vs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((vc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((ts(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((tc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((qs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((qc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((qpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(5)(((qparc(l,m,n),l=1,ld),m=1,md),n=1,nd) write(*,*) ok read(5)(timo(i),i=netb,net) read(5)(wpfx(i),i=netb,net) read(5)(gamx(i),i=netb,net) read(5)(wrhx(i),i=netb,net) read(5)(wrhy(i),i=netb,net) read(5)(wdux(i),i=netb,net) read(5)(wu2x(i),i=netb,net) read(5)(wp2x(i),i=netb,net) read(5)(wv2x(i),i=netb,net) read(5)(wlux(i),i=netb,net) read(5)(wdpx(i),i=netb,net) read(5)(wdvx(i),i=netb,net) read(5)(wenx(i),i=netb,net) read(5)(wdfx(i),i=netb,net) read(5)(pcerr(i),i=netb,net) read(5)(((wkif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(5)(((wpif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(5)(((wvif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(5)(((wplx(i,m,n),i=netb,net),m=1,md),n=1,nd) read(5)(wakx(i),i=netb,net) read(5)(waky(i),i=netb,net) read(5)(therm(i),i=netb,net) read(5)(wxsp(i),i=netb,net) read(5)(((grtmx(i,m,n),i=netb,net),m=1,md),n=1,nd) ncnt=ncnt+ntotal read(5)(((uctim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(5)(((ustim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(5)(timf(i),i=1,ntotal) write(*,*) ok read(5)avgflag write(*,*) avgflag netm=net-1 do 11 n=1,nd do 11 m=1,md grtmx(net,m,n)=grtmx(netm,m,n) 11 continue c nff=nff+1 if(nff.gt.nfile) goto 1000 netb=net+1 ncntb=ncnt+1 read(8)ld,ldb,nmin,nmax,nd,lin,ne,nfreq,ntotal,md read(8)(mmin(n),n=1,nd) read(8)(mmax(n),n=1,nd) read(8)tim,shr,qsf,epsn,etai,rmu1,rmu2,rkap,aex,rmu2pa,rkpara read(8)dt0,dt,x0,y0,z0,rkkmult net=net+ne read(8)(((ws(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((wc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((us(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((uc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((tpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((tparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((prs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((prc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((vs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((vc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((ts(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((tc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((qs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((qc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((qpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(((qparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(8)(timo(i),i=netb,net) read(8)(wpfx(i),i=netb,net) read(8)(gamx(i),i=netb,net) read(8)(wrhx(i),i=netb,net) read(8)(wrhy(i),i=netb,net) read(8)(wdux(i),i=netb,net) read(8)(wu2x(i),i=netb,net) read(8)(wp2x(i),i=netb,net) read(8)(wv2x(i),i=netb,net) read(8)(wlux(i),i=netb,net) read(8)(wdpx(i),i=netb,net) read(8)(wdvx(i),i=netb,net) read(8)(wenx(i),i=netb,net) read(8)(wdfx(i),i=netb,net) read(8)(pcerr(i),i=netb,net) read(8)(((wkif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(8)(((wpif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(8)(((wvif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(8)(((wplx(i,m,n),i=netb,net),m=1,md),n=1,nd) read(8)(wakx(i),i=netb,net) read(8)(waky(i),i=netb,net) read(8)(therm(i),i=netb,net) read(8)(wxsp(i),i=netb,net) read(8)(((grtmx(i,m,n),i=netb,net),m=1,md),n=1,nd) ncnt=ncnt+ntotal read(8)(((uctim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(8)(((ustim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(8)(timf(i),i=1,ntotal) read(8)x netm=net-1 do 12 n=1,nd do 12 m=1,md grtmx(net,m,n)=grtmx(netm,m,n) 12 continue c nff=nff+1 if(nff.gt.nfile) goto 1000 netb=net+1 ncntb=ncnt+1 read(9)ld,ldb,nmin,nmax,nd,lin,ne,nfreq,ntotal,md read(9)(mmin(n),n=1,nd) read(9)(mmax(n),n=1,nd) read(9)tim,shr,qsf,epsn,etai,rmu1,rmu2,rkap,aex,rmu2pa,rkpara read(9)dt0,dt,x0,y0,z0,rkkmult net=net+ne read(9)(((ws(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((wc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((us(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((uc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((tpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((tparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((prs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((prc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((vs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((vc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((ts(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((tc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((qs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((qc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((qpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(((qparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(9)(timo(i),i=netb,net) read(9)(wpfx(i),i=netb,net) read(9)(gamx(i),i=netb,net) read(9)(wrhx(i),i=netb,net) read(9)(wrhy(i),i=netb,net) read(9)(wdux(i),i=netb,net) read(9)(wu2x(i),i=netb,net) read(9)(wp2x(i),i=netb,net) read(9)(wv2x(i),i=netb,net) read(9)(wlux(i),i=netb,net) read(9)(wdpx(i),i=netb,net) read(9)(wdvx(i),i=netb,net) read(9)(wenx(i),i=netb,net) read(9)(wdfx(i),i=netb,net) read(9)(pcerr(i),i=netb,net) read(9)(((wkif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(9)(((wpif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(9)(((wvif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(9)(((wplx(i,m,n),i=netb,net),m=1,md),n=1,nd) read(9)(wakx(i),i=netb,net) read(9)(waky(i),i=netb,net) read(9)(therm(i),i=netb,net) read(9)(wxsp(i),i=netb,net) read(9)(((grtmx(i,m,n),i=netb,net),m=1,md),n=1,nd) ncnt=ncnt+ntotal read(9)(((uctim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(9)(((ustim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(9)(timf(i),i=1,ntotal) read(9)x netm=net-1 do 13 n=1,nd do 13 m=1,md grtmx(net,m,n)=grtmx(netm,m,n) 13 continue c nff=nff+1 if(nff.gt.nfile) goto 1000 netb=net+1 ncntb=ncnt+1 read(10)ld,ldb,nmin,nmax,nd,lin,ne,nfreq,ntotal,md read(10)(mmin(n),n=1,nd) read(10)(mmax(n),n=1,nd) read(10)tim,shr,qsf,epsn,etai,rmu1,rmu2,rkap,aex,rmu2pa,rkpara read(10)dt0,dt,x0,y0,z0,rkkmult net=net+ne read(10)(((ws(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((wc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((us(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((uc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((tpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((tparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((prs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((prc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((vs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((vc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((ts(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((tc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((qs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((qc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((qpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(((qparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(10)(timo(i),i=netb,net) read(10)(wpfx(i),i=netb,net) read(10)(gamx(i),i=netb,net) read(10)(wrhx(i),i=netb,net) read(10)(wrhy(i),i=netb,net) read(10)(wdux(i),i=netb,net) read(10)(wu2x(i),i=netb,net) read(10)(wp2x(i),i=netb,net) read(10)(wv2x(i),i=netb,net) read(10)(wlux(i),i=netb,net) read(10)(wdpx(i),i=netb,net) read(10)(wdvx(i),i=netb,net) read(10)(wenx(i),i=netb,net) read(10)(wdfx(i),i=netb,net) read(10)(pcerr(i),i=netb,net) read(10)(((wkif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(10)(((wpif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(10)(((wvif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(10)(((wplx(i,m,n),i=netb,net),m=1,md),n=1,nd) read(10)(wakx(i),i=netb,net) read(10)(waky(i),i=netb,net) read(10)(therm(i),i=netb,net) read(10)(wxsp(i),i=netb,net) read(10)(((grtmx(i,m,n),i=netb,net),m=1,md),n=1,nd) ncnt=ncnt+ntotal read(10)(((uctim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(10)(((ustim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(10)(timf(i),i=1,ntotal) read(10)x netm=net-1 do 14 n=1,nd do 14 m=1,md grtmx(net,m,n)=grtmx(netm,m,n) 14 continue c nff=nff+1 if(nff.gt.nfile) goto 1000 netb=net+1 ncntb=ncnt+1 read(11)ld,ldb,nmin,nmax,nd,lin,ne,nfreq,ntotal,md read(11)(mmin(n),n=1,nd) read(11)(mmax(n),n=1,nd) read(11)tim,shr,qsf,epsn,etai,rmu1,rmu2,rkap,aex,rmu2pa,rkpara read(11)dt0,dt,x0,y0,z0,rkkmult net=net+ne read(11)(((ws(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((wc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((us(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((uc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((tpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((tparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((prs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((prc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((vs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((vc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((ts(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((tc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((qs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((qc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((qpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(((qparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(11)(timo(i),i=netb,net) read(11)(wpfx(i),i=netb,net) read(11)(gamx(i),i=netb,net) read(11)(wrhx(i),i=netb,net) read(11)(wrhy(i),i=netb,net) read(11)(wdux(i),i=netb,net) read(11)(wu2x(i),i=netb,net) read(11)(wp2x(i),i=netb,net) read(11)(wv2x(i),i=netb,net) read(11)(wlux(i),i=netb,net) read(11)(wdpx(i),i=netb,net) read(11)(wdvx(i),i=netb,net) read(11)(wenx(i),i=netb,net) read(11)(wdfx(i),i=netb,net) read(11)(pcerr(i),i=netb,net) read(11)(((wkif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(11)(((wpif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(11)(((wvif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(11)(((wplx(i,m,n),i=netb,net),m=1,md),n=1,nd) read(11)(wakx(i),i=netb,net) read(11)(waky(i),i=netb,net) read(11)(therm(i),i=netb,net) read(11)(wxsp(i),i=netb,net) read(11)(((grtmx(i,m,n),i=netb,net),m=1,md),n=1,nd) ncnt=ncnt+ntotal read(11)(((uctim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(11)(((ustim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(11)(timf(i),i=1,ntotal) read(11)x netm=net-1 do 15 n=1,nd do 15 m=1,md grtmx(net,m,n)=grtmx(netm,m,n) 15 continue c nff=nff+1 if(nff.gt.nfile) goto 1000 netb=net+1 ncntb=ncnt+1 read(13)ld,ldb,nmin,nmax,nd,lin,ne,nfreq,ntotal,md read(13)(mmin(n),n=1,nd) read(13)(mmax(n),n=1,nd) read(13)tim,shr,qsf,epsn,etai,rmu1,rmu2,rkap,aex,rmu2pa,rkpara read(13)dt0,dt,x0,y0,z0,rkkmult net=net+ne read(13)(((ws(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((wc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((us(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((uc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((tpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((tparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((prs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((prc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((vs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((vc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((ts(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((tc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((qs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((qc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((qpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(((qparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(13)(timo(i),i=netb,net) read(13)(wpfx(i),i=netb,net) read(13)(gamx(i),i=netb,net) read(13)(wrhx(i),i=netb,net) read(13)(wrhy(i),i=netb,net) read(13)(wdux(i),i=netb,net) read(13)(wu2x(i),i=netb,net) read(13)(wp2x(i),i=netb,net) read(13)(wv2x(i),i=netb,net) read(13)(wlux(i),i=netb,net) read(13)(wdpx(i),i=netb,net) read(13)(wdvx(i),i=netb,net) read(13)(wenx(i),i=netb,net) read(13)(wdfx(i),i=netb,net) read(13)(pcerr(i),i=netb,net) read(13)(((wkif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(13)(((wpif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(13)(((wvif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(13)(((wplx(i,m,n),i=netb,net),m=1,md),n=1,nd) read(13)(wakx(i),i=netb,net) read(13)(waky(i),i=netb,net) read(13)(therm(i),i=netb,net) read(13)(wxsp(i),i=netb,net) read(13)(((grtmx(i,m,n),i=netb,net),m=1,md),n=1,nd) ncnt=ncnt+ntotal read(13)(((uctim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(13)(((ustim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(13)(timf(i),i=1,ntotal) read(13)x netm=net-1 do 16 n=1,nd do 16 m=1,md grtmx(net,m,n)=grtmx(netm,m,n) 16 continue c c nff=nff+1 if(nff.gt.nfile) goto 1000 netb=net+1 ncntb=ncnt+1 read(14)ld,ldb,nmin,nmax,nd,lin,ne,nfreq,ntotal,md read(14)(mmin(n),n=1,nd) read(14)(mmax(n),n=1,nd) read(14)tim,shr,qsf,epsn,etai,rmu1,rmu2,rkap,aex,rmu2pa,rkpara read(14)dt0,dt,x0,y0,z0,rkkmult net=net+ne read(14)(((ws(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((wc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((us(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((uc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((tpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((tparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((prs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((prc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((vs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((vc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((ts(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((tc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((qs(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((qc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((qpars(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(((qparc(l,m,n),l=1,ld),m=1,md),n=1,nd) read(14)(timo(i),i=netb,net) read(14)(wpfx(i),i=netb,net) read(14)(gamx(i),i=netb,net) read(14)(wrhx(i),i=netb,net) read(14)(wrhy(i),i=netb,net) read(14)(wdux(i),i=netb,net) read(14)(wu2x(i),i=netb,net) read(14)(wp2x(i),i=netb,net) read(14)(wv2x(i),i=netb,net) read(14)(wlux(i),i=netb,net) read(14)(wdpx(i),i=netb,net) read(14)(wdvx(i),i=netb,net) read(14)(wenx(i),i=netb,net) read(14)(wdfx(i),i=netb,net) read(14)(pcerr(i),i=netb,net) read(14)(((wkif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(14)(((wpif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(14)(((wvif(i,m,n),i=netb,net),m=1,md),n=1,nd) read(14)(((wplx(i,m,n),i=netb,net),m=1,md),n=1,nd) read(14)(wakx(i),i=netb,net) read(14)(waky(i),i=netb,net) read(14)(therm(i),i=netb,net) read(14)(wxsp(i),i=netb,net) read(14)(((grtmx(i,m,n),i=netb,net),m=1,md),n=1,nd) ncnt=ncnt+ntotal read(14)(((uctim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(14)(((ustim(i,m,n),i=ncntb,ncnt),m=1,md),n=1,nd) read(14)(timf(i),i=1,ntotal) read(14)x netm=net-1 1000 continue c c now write out to unit 7 c if(newg.ne.1) then ld1=ld ldb1=ldb nmin1=nmin nmax1=nmax nd1=nd do 10 i=1,nd mmax1(i)=mmax(i) mmin1(i)=mmin(i) md1=md 10 continue endif c write(7,1)ld1,ldb1,nmin1,nmax1,nd1,lin,ne write(7,1)(mmin1(n),n=1,nd1) write(7,1)(mmax1(n),n=1,nd1) write(7,2)tim,shr,qsf,epsn,etai,rmu1,rmu2,rkap,aex,rmu2pa,rkpara write(7,2)dt0,dt,x0,y0,z0 write(7,2)(((ws(l,m,n),l=1,ld),m=1,md1),n=1,nd1) write(7,2)(((wc(l,m,n),l=1,ld),m=1,md1),n=1,nd1) write(7,2)(((us(l,m,n),l=1,ld),m=1,md1),n=1,nd1) write(7,2)(((uc(l,m,n),l=1,ld),m=1,md1),n=1,nd1) write(7,2)(((prs(l,m,n),l=1,ld),m=1,md1),n=1,nd1) write(7,2)(((prc(l,m,n),l=1,ld),m=1,md1),n=1,nd1) write(7,2)(((vs(l,m,n),l=1,ld),m=1,md1),n=1,nd1) write(7,2)(((vc(l,m,n),l=1,ld),m=1,md1),n=1,nd1) write(7,2)(((tc(l,m,n),l=1,ld),m=1,md1),n=1,nd) write(7,2)(((ts(l,m,n),l=1,ld),m=1,md1),n=1,nd) write(7,2)(((qc(l,m,n),l=1,ld),m=1,md1),n=1,nd) write(7,2)(((qs(l,m,n),l=1,ld),m=1,md1),n=1,nd) write(7,2)(((qparc(l,m,n),l=1,ld),m=1,md1),n=1,nd) write(7,2)(((qpars(l,m,n),l=1,ld),m=1,md1),n=1,nd) write(7,2)(timo(i),i=1,ne) write(7,2)(wpfx(i),i=1,ne) write(7,2)(gamx(i),i=1,ne) write(7,2)(wrhx(i),i=1,ne) write(7,2)(wrhy(i),i=1,ne) write(7,2)(wdux(i),i=1,ne) write(7,2)(wu2x(i),i=1,ne) write(7,2)(wp2x(i),i=1,ne) write(7,2)(wv2x(i),i=1,ne) write(7,2)(wlux(i),i=1,ne) write(7,2)(wdpx(i),i=1,ne) write(7,2)(wdvx(i),i=1,ne) write(7,2)(wenx(i),i=1,ne) write(7,2)(wdfx(i),i=1,ne) write(7,2)(pcerr(i),i=1,ne) write(7,2)(((wkif(i,m,n),i=1,ne),m=1,md1),n=1,nd1) write(7,2)(((wpif(i,m,n),i=1,ne),m=1,md1),n=1,nd1) write(7,2)(((wvif(i,m,n),i=1,ne),m=1,md1),n=1,nd1) write(7,2)(((wplx(i,m,n),i=1,ne),m=1,md1),n=1,nd1) write(7,2)(wakx(i),i=1,ne) write(7,2)(waky(i),i=1,ne) write(7,2)(therm(i),i=1,ne) write(7,2)(wxsp(i),i=1,ne) write(7,2)(((grtmx(i,m,n),i=1,ne),m=1,md1),n=1,nd1) write(7,2)(((uctim(i,m,n),i=1,ntotal),m=1,md1),n=1,nd) write(7,2)(((ustim(i,m,n),i=1,ntotal),m=1,md1),n=1,nd) write(7,2)(timf(i),i=1,ntotal) write(7,2)avgflag 1 format(1x,3i10) 2 format(1x,3e22.14) stop end