chlml.f conmc.f dtast.f dvast.f fft99m.f fulldivt.f g2s0.f gdcrdn.f gdcrdp.f genqsat.f getbaln.f getlalo.f getpln.f glbsoi2.f 61,97c61,97 < dimension rpw(2*nlath+1,nlon+2),sigi(nsig+1) < dimension xhat((jcap+1)*(jcap+2),nsig,4) < dimension xhatp((jcap+1)*(jcap+2)) < dimension del2((jcap+1)*(jcap+2)) < dimension trigs(nlon*2),ifax(10) < dimension rlats(nlath*2),slat(nlath),clat(nlath) < dimension wgts(nlath*2),sigl(nsig),pwcon(nsig) < dimension isatv(nsig) < dimension rt(2*nlath+1,nlon+2,nsig) < dimension ru(2*nlath+1,nlon+2,nsig) < dimension rv(2*nlath+1,nlon+2,nsig) < dimension rq(2*nlath+1,nlon+2,nsig) < dimension rp(2*nlath+1,nlon+2) < dimension rus(2*nlath+1,nlon+2,nsig) < dimension rvs(2*nlath+1,nlon+2,nsig) < dimension rts(2*nlath+1,nlon+2,nsig) < dimension rvorts(2*nlath+1,nlon+2,nsig) < dimension rplons(2*nlath+1,nlon+2) < dimension rplats(2*nlath+1,nlon+2) < dimension rlsg(nsig),aibw(nsig,nsig),w(nsig),tbar(nsig) < dimension a3(nsig,nsig),dstlast((jcap+1)*(jcap+2),nsig) < dimension dstb((jcap+1)*(jcap+2),nsig) < dimension in((jcap+1)*(jcap+2)) < dimension pln((jcap+1)*(jcap+2),nlath) < dimension qln((jcap+1)*(jcap+2),nlath) < dimension rln((jcap+1)*(jcap+2),nlath) < dimension mlad(0:jcap,0:jcap) < dimension ml2lm((jcap+1)*(jcap+2)) < dimension factslm((jcap+1)*(jcap+2)) < dimension factvlm((jcap+1)*(jcap+2)) < dimension lmad(0:jcap,0:jcap) < dimension lm2ml((jcap+1)*(jcap+2)) < dimension factsml((jcap+1)*(jcap+2)) < dimension factvml((jcap+1)*(jcap+2)) < dimension qfile(17*nqdata),uvfile(18*nwdata),tfile(17*ntdata) < dimension sfile((4+(28+2)*30)*nsprof),pwfile(12*npwdat) < dimension psfile(11*npdata) --- > C-CRA dimension rpw(2*nlath+1,nlon+2),sigi(nsig+1) > C-CRA dimension xhat((jcap+1)*(jcap+2),nsig,4) > C-CRA dimension xhatp((jcap+1)*(jcap+2)) > C-CRA dimension del2((jcap+1)*(jcap+2)) > C-CRA dimension trigs(nlon*2),ifax(10) > C-CRA dimension rlats(nlath*2),slat(nlath),clat(nlath) > C-CRA dimension wgts(nlath*2),sigl(nsig),pwcon(nsig) > C-CRA dimension isatv(nsig) > C-CRA dimension rt(2*nlath+1,nlon+2,nsig) > C-CRA dimension ru(2*nlath+1,nlon+2,nsig) > C-CRA dimension rv(2*nlath+1,nlon+2,nsig) > C-CRA dimension rq(2*nlath+1,nlon+2,nsig) > C-CRA dimension rp(2*nlath+1,nlon+2) > C-CRA dimension rus(2*nlath+1,nlon+2,nsig) > C-CRA dimension rvs(2*nlath+1,nlon+2,nsig) > C-CRA dimension rts(2*nlath+1,nlon+2,nsig) > C-CRA dimension rvorts(2*nlath+1,nlon+2,nsig) > C-CRA dimension rplons(2*nlath+1,nlon+2) > C-CRA dimension rplats(2*nlath+1,nlon+2) > C-CRA dimension rlsg(nsig),aibw(nsig,nsig),w(nsig),tbar(nsig) > C-CRA dimension a3(nsig,nsig),dstlast((jcap+1)*(jcap+2),nsig) > C-CRA dimension dstb((jcap+1)*(jcap+2),nsig) > C-CRA dimension in((jcap+1)*(jcap+2)) > C-CRA dimension pln((jcap+1)*(jcap+2),nlath) > C-CRA dimension qln((jcap+1)*(jcap+2),nlath) > C-CRA dimension rln((jcap+1)*(jcap+2),nlath) > C-CRA dimension mlad(0:jcap,0:jcap) > C-CRA dimension ml2lm((jcap+1)*(jcap+2)) > C-CRA dimension factslm((jcap+1)*(jcap+2)) > C-CRA dimension factvlm((jcap+1)*(jcap+2)) > C-CRA dimension lmad(0:jcap,0:jcap) > C-CRA dimension lm2ml((jcap+1)*(jcap+2)) > C-CRA dimension factsml((jcap+1)*(jcap+2)) > C-CRA dimension factvml((jcap+1)*(jcap+2)) > C-CRA dimension qfile(17*nqdata),uvfile(18*nwdata),tfile(17*ntdata) > C-CRA dimension sfile((4+(28+2)*30)*nsprof),pwfile(12*npwdat) > C-CRA dimension psfile(11*npdata) 98a99,139 > dimension rpw(2*48+1,192+2),sigi(28+1) > dimension xhat((62+1)*(62+2),28,4) > dimension xhatp((62+1)*(62+2)) > dimension del2((62+1)*(62+2)) > dimension trigs(192*2),ifax(10) > dimension rlats(48*2),slat(48),clat(48) > dimension wgts(48*2),sigl(28),pwcon(28) > dimension isatv(28) > dimension rt(2*48+1,192+2,28) > dimension ru(2*48+1,192+2,28) > dimension rv(2*48+1,192+2,28) > dimension rq(2*48+1,192+2,28) > dimension rp(2*48+1,192+2) > dimension rus(2*48+1,192+2,28) > dimension rvs(2*48+1,192+2,28) > dimension rts(2*48+1,192+2,28) > dimension rvorts(2*48+1,192+2,28) > dimension rplons(2*48+1,192+2) > dimension rplats(2*48+1,192+2) > dimension rlsg(28),aibw(28,28) > dimension w(28),tbar(28) > dimension a3(28,28) > dimension dstlast((62+1)*(62+2),28) > dimension dstb((62+1)*(62+2),28) > dimension in((62+1)*(62+2)) > dimension pln((62+1)*(62+2),48) > dimension qln((62+1)*(62+2),48) > dimension rln((62+1)*(62+2),48) > dimension mlad(0:62,0:62) > dimension ml2lm((62+1)*(62+2)) > dimension factslm((62+1)*(62+2)) > dimension factvlm((62+1)*(62+2)) > dimension lmad(0:62,0:62) > dimension lm2ml((62+1)*(62+2)) > dimension factsml((62+1)*(62+2)) > dimension factvml((62+1)*(62+2)) > dimension qfile(17*15000),uvfile(18*112000) > dimension tfile(17*60000) > dimension sfile((4+(28+2)*30)*10000) > dimension pwfile(12*1) > dimension psfile(11*28000) 148c189,192 < factvlm=factslm --- > C-CRA factvlm=factslm > DO i=1,(jcap+1)*(jcap+2) > factvlm(i)=factslm(i) > ENDDO 160c204,207 < factvml=factsml --- > C-CRA factvml=factsml > DO i=1,(jcap+1)*(jcap+2) > factvml(i)=factsml(i) > ENDDO grad2s.f gtbdivt.f gtbhalf.f hoper.f hopers.f htoper.f inguess.f inguessv.f initps2.f initqpw2.f initsat6.f initt3.f initw2.f intps2.f intqpw2.f 33,36c33,36 < dimension pwcon(nsig) < dimension rt(2*nlath+1,nlon+2,nsig) < dimension qfile(*),pwfile(*) < dimension st(2*nlath+1,nlon+2,nsig) --- > C-CRA dimension pwcon(nsig) > C-CRA dimension rt(2*nlath+1,nlon+2,nsig) > C-CRA dimension qfile(*),pwfile(*) > C-CRA dimension st(2*nlath+1,nlon+2,nsig) 37a38,41 > dimension pwcon(28) > dimension rt(2*48+1,192+2,28) > dimension qfile(*),pwfile(*) > dimension st(2*48+1,192+2,28) 42,43c46,55 < st=rt < rt=0. --- > C-CRA st=rt > c dimension st(2*nlath+1,nlon+2,nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > st(ITMP,1,1)=rt(ITMP,1,1) > ENDDO > C-CRA rt=0. > c dimension rt(2*nlath+1,nlon+2,nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > rt(ITMP,1,1)=0. > ENDDO intrp2.f intrp3.f intt3.f intw2.f isrchfge.f isrchfle.f iw3jdn.f m1glat.f m1ipqr.f m1poly.f m1rcons2.f minv.f pcgsoi5.f 81,137c81,206 < dimension pwcon(nsig) < dimension trigs(nlon*2),ifax(10) < dimension isatv(nsig),as(4) < dimension rus(2*nlath+1,nlon+2,nsig) < dimension rvs(2*nlath+1,nlon+2,nsig) < dimension rts(2*nlath+1,nlon+2,nsig) < dimension rvorts(2*nlath+1,nlon+2,nsig) < dimension rplons(2*nlath+1,nlon+2) < dimension rplats(2*nlath+1,nlon+2) < dimension mlad(0:jcap,0:jcap) < dimension ml2lm((jcap+1)*(jcap+2)) < dimension factslm((jcap+1)*(jcap+2)) < dimension factvlm((jcap+1)*(jcap+2)) < dimension lmad(0:jcap,0:jcap) < dimension lm2ml((jcap+1)*(jcap+2)) < dimension factsml((jcap+1)*(jcap+2)) < dimension factvml((jcap+1)*(jcap+2)) < dimension qfile(17*nqdata),uvfile(18*nwdata) < dimension tfile(17*ntdata) < dimension sfile((4+(28+2)*30)*nsprof),pwfile(12*npwdat) < dimension psfile(11*npdata) < dimension sigl(nsig),sigi(nsig+1) < dimension vz(nsig,nsig),vd(nsig,nsig) < dimension vq(nsig,nsig),vh(nsig,nsig) < dimension agvz(0:jcap,nsig,nmdszh),wgvz(0:jcap,nmdszh) < dimension bvz(0:jcap,nsig,nmdszh) < dimension in((jcap+1)*(jcap+2)) < dimension c1(2),b1(2) < dimension a3(nsig,nsig) < dimension dsms(nsig),dbms(nsig),difms(nsig) < real rlats(2*nlath) < real del2((jcap+1)*(jcap+2)),wgts(nlath*2) < real rp((2*nlath+1)*(nlon+2)) < real rplon((2*nlath+1)*(nlon+2)) < real rplat((2*nlath+1)*(nlon+2)) < real rpw((2*nlath+1)*(nlon+2)) < real rq((2*nlath+1)*(nlon+2)*nsig) < real rt((2*nlath+1)*(nlon+2)*nsig) < real st((2*nlath+1)*(nlon+2)*nsig) < real ru((2*nlath+1)*(nlon+2)*nsig) < real rv((2*nlath+1)*(nlon+2)*nsig) < real rvort((2*nlath+1)*(nlon+2)*nsig) < real pln((jcap+1)*(jcap+2),nlath) < real qln((jcap+1)*(jcap+2),nlath) < real rln((jcap+1)*(jcap+2),nlath) < real cshat((jcap+1)*(jcap+2)) < real bhalf((jcap+1)*(jcap+2),nsig,4) < real bhalfp((jcap+1)*(jcap+2)) < real bdivt(0:jcap,nsig) < real xhat((jcap+1)*(jcap+2),nsig,4) < real xhatp((jcap+1)*(jcap+2)) < real phat((jcap+1)*(jcap+2),nsig,4) < real phatp((jcap+1)*(jcap+2)) < real fhat((jcap+1)*(jcap+2),nsig,4) < real fhatp((jcap+1)*(jcap+2)) < real ghat((jcap+1)*(jcap+2),nsig,4) < real ghatp((jcap+1)*(jcap+2)) --- > C-CRA dimension pwcon(nsig) > C-CRA dimension trigs(nlon*2),ifax(10) > C-CRA dimension isatv(nsig),as(4) > C-CRA dimension rus(2*nlath+1,nlon+2,nsig) > C-CRA dimension rvs(2*nlath+1,nlon+2,nsig) > C-CRA dimension rts(2*nlath+1,nlon+2,nsig) > C-CRA dimension rvorts(2*nlath+1,nlon+2,nsig) > C-CRA dimension rplons(2*nlath+1,nlon+2) > C-CRA dimension rplats(2*nlath+1,nlon+2) > C-CRA dimension mlad(0:jcap,0:jcap) > C-CRA dimension ml2lm((jcap+1)*(jcap+2)) > C-CRA dimension factslm((jcap+1)*(jcap+2)) > C-CRA dimension factvlm((jcap+1)*(jcap+2)) > C-CRA dimension lmad(0:jcap,0:jcap) > C-CRA dimension lm2ml((jcap+1)*(jcap+2)) > C-CRA dimension factsml((jcap+1)*(jcap+2)) > C-CRA dimension factvml((jcap+1)*(jcap+2)) > C-CRA dimension qfile(17*nqdata),uvfile(18*nwdata) > C-CRA dimension tfile(17*ntdata) > C-CRA dimension sfile((4+(28+2)*30)*nsprof),pwfile(12*npwdat) > C-CRA dimension psfile(11*npdata) > C-CRA dimension sigl(nsig),sigi(nsig+1) > C-CRA dimension vz(nsig,nsig),vd(nsig,nsig) > C-CRA dimension vq(nsig,nsig),vh(nsig,nsig) > C-CRA dimension agvz(0:jcap,nsig,nmdszh),wgvz(0:jcap,nmdszh) > C-CRA dimension bvz(0:jcap,nsig,nmdszh) > C-CRA dimension in((jcap+1)*(jcap+2)) > C-CRA dimension c1(2),b1(2) > C-CRA dimension a3(nsig,nsig) > C-CRA dimension dsms(nsig),dbms(nsig),difms(nsig) > C-CRA real rlats(2*nlath) > C-CRA real del2((jcap+1)*(jcap+2)),wgts(nlath*2) > C-CRA real rp((2*nlath+1)*(nlon+2)) > C-CRA real rplon((2*nlath+1)*(nlon+2)) > C-CRA real rplat((2*nlath+1)*(nlon+2)) > C-CRA real rpw((2*nlath+1)*(nlon+2)) > C-CRA real rq((2*nlath+1)*(nlon+2)*nsig) > C-CRA real rt((2*nlath+1)*(nlon+2)*nsig) > C-CRA real st((2*nlath+1)*(nlon+2)*nsig) > C-CRA real ru((2*nlath+1)*(nlon+2)*nsig) > C-CRA real rv((2*nlath+1)*(nlon+2)*nsig) > C-CRA real rvort((2*nlath+1)*(nlon+2)*nsig) > C-CRA real pln((jcap+1)*(jcap+2),nlath) > C-CRA real qln((jcap+1)*(jcap+2),nlath) > C-CRA real rln((jcap+1)*(jcap+2),nlath) > C-CRA real cshat((jcap+1)*(jcap+2)) > C-CRA real bhalf((jcap+1)*(jcap+2),nsig,4) > C-CRA real bhalfp((jcap+1)*(jcap+2)) > C-CRA real bdivt(0:jcap,nsig) > C-CRA real xhat((jcap+1)*(jcap+2),nsig,4) > C-CRA real xhatp((jcap+1)*(jcap+2)) > C-CRA real phat((jcap+1)*(jcap+2),nsig,4) > C-CRA real phatp((jcap+1)*(jcap+2)) > C-CRA real fhat((jcap+1)*(jcap+2),nsig,4) > C-CRA real fhatp((jcap+1)*(jcap+2)) > C-CRA real ghat((jcap+1)*(jcap+2),nsig,4) > C-CRA real ghatp((jcap+1)*(jcap+2)) > C-CRA real denom(2) > C-CRA real ps((jcap+1)*(jcap+2)) > C-CRA real dstlast((jcap+1)*(jcap+2),nsig) > C-CRA real dstb((jcap+1)*(jcap+2),nsig) > C-CRA real factor((jcap+1)*(jcap+2)) > C-CRA real factori((jcap+1)*(jcap+2)) > C-CRA real baln((jcap+1)*(jcap+2)) > C-CRA integer idateg(4) > C-CRA character*4 on85dt(8) > C-CRA character*4 on85(8) > > dimension pwcon(28) > dimension trigs(192*2),ifax(10) > dimension isatv(28),as(4) > dimension rus(2*48+1,192+2,28) > dimension rvs(2*48+1,192+2,28) > dimension rts(2*48+1,192+2,28) > dimension rvorts(2*48+1,192+2,28) > dimension rplons(2*48+1,192+2) > dimension rplats(2*48+1,192+2) > dimension mlad(0:62,0:62) > dimension ml2lm((62+1)*(62+2)) > dimension factslm((62+1)*(62+2)) > dimension factvlm((62+1)*(62+2)) > dimension lmad(0:62,0:62) > dimension lm2ml((62+1)*(62+2)) > dimension factsml((62+1)*(62+2)) > dimension factvml((62+1)*(62+2)) > dimension qfile(17*15000),uvfile(18*112000) > dimension tfile(17*60000) > dimension sfile((4+(28+2)*30)*10000) > dimension pwfile(12*1),psfile(11*28000) > dimension sigl(28),sigi(28+1) > dimension vz(28,28),vd(28,28) > dimension vq(28,28),vh(28,28) > dimension agvz(0:62,28,28) > dimension wgvz(0:62,28) > dimension bvz(0:62,28,28) > dimension in((62+1)*(62+2)) > dimension c1(2),_b1(2) > dimension a3(28,28) > dimension dsms(28),dbms(28),difms(28) > real rlats(2*48) > real del2((62+1)*(62+2)),wgts(48*2) > real rp((2*48+1)*(192+2)) > real rplon((2*48+1)*(192+2)) > real rplat((2*48+1)*(192+2)) > real rpw((2*48+1)*(192+2)) > real rq((2*48+1)*(192+2)*28) > real rt((2*48+1)*(192+2)*28) > real st((2*48+1)*(192+2)*28) > real ru((2*48+1)*(192+2)*28) > real rv((2*48+1)*(192+2)*28) > real rvort((2*48+1)*(192+2)*28) > real pln((62+1)*(62+2),48) > real qln((62+1)*(62+2),48) > real rln((62+1)*(62+2),48) > real cshat((62+1)*(62+2)) > real bhalf((62+1)*(62+2),28,4) > real bhalfp((62+1)*(62+2)) > real bdivt(0:62,28) > real xhat((62+1)*(62+2),28,4) > real xhatp((62+1)*(62+2)) > real phat((62+1)*(62+2),28,4) > real phatp((62+1)*(62+2)) > real fhat((62+1)*(62+2),28,4) > real fhatp((62+1)*(62+2)) > real ghat((62+1)*(62+2),28,4) > real ghatp((62+1)*(62+2)) 139,144c208,213 < real ps((jcap+1)*(jcap+2)) < real dstlast((jcap+1)*(jcap+2),nsig) < real dstb((jcap+1)*(jcap+2),nsig) < real factor((jcap+1)*(jcap+2)) < real factori((jcap+1)*(jcap+2)) < real baln((jcap+1)*(jcap+2)) --- > real ps((62+1)*(62+2)) > real dstlast((62+1)*(62+2),28) > real dstb((62+1)*(62+2),28) > real factor((62+1)*(62+2)) > real factori((62+1)*(62+2)) > real baln((62+1)*(62+2)) 148d216 < 192c260,264 < factor=factor/(16.*atan(1.)) --- > C-CRA factor=factor/(16.*atan(1.)) > c real factor((jcap+1)*(jcap+2)) > DO ITMP=1,(jcap+1)*(jcap+2) > factor(ITMP)=factor(ITMP)/(16.*atan(1.)) > ENDDO 195c267,271 < rp=0. --- > C-CRA rp=0. > c real rp((2*nlath+1)*(nlon+2)) > DO ITMP=1,(2*nlath+1)*(nlon+2) > rp(ITMP)=0. > ENDDO 198,199c274,283 < ru=0. < rv=0. --- > C-CRA ru=0. > c real ru((2*nlath+1)*(nlon+2)*nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > ru(ITMP)=0. > ENDDO > C-CRA rv=0. > c real rv((2*nlath+1)*(nlon+2)*nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > rv(ITMP)=0. > ENDDO 205c289,293 < rt=0. --- > C-CRA rt=0. > c real rt((2*nlath+1)*(nlon+2)*nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > rt(ITMP)=0. > ENDDO 209c297,301 < rq=0. --- > C-CRA rq=0. > c real rq((2*nlath+1)*(nlon+2)*nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > rq(ITMP)=0. > ENDDO 216c308,312 < rt=rt+st --- > C-CRA rt=rt+st > c real rt((2*nlath+1)*(nlon+2)*nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > rt(ITMP)=rt(ITMP)+st(ITMP) > ENDDO 231c327,333 < bdivt=ampdivt*bdivt --- > C-CRA bdivt=ampdivt*bdivt > c real bdivt(0:jcap,nsig) > DO JTMP=1,nsig > DO ITMP=0,jcap > bdivt(ITMP,JTMP)=ampdivt*bdivt(ITMP,JTMP) > ENDDO > ENDDO 233,235c335,349 < dbms=0. < dsms=0. < difms=0. --- > C-CRA dbms=0. > c dimension dsms(nsig),dbms(nsig),difms(nsig) > DO ITMP=1,nsig > dbms(ITMP)=0. > ENDDO > C-CRA dsms=0. > c dimension dsms(nsig),dbms(nsig),difms(nsig) > DO ITMP=1,nsig > dsms(ITMP)=0. > ENDDO > C-CRA difms=0. > c dimension dsms(nsig),dbms(nsig),difms(nsig) > DO ITMP=1,nsig > difms(ITMP)=0. > ENDDO 281,283c395,409 < rvort=0. < rplat=0. < rplon=0. --- > C-CRA rvort=0. > c real rvort((2*nlath+1)*(nlon+2)*nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > rvort(ITMP)=0. > ENDDO > C-CRA rplat=0. > c real rplat((2*nlath+1)*(nlon+2)) > DO ITMP=1,(2*nlath+1)*(nlon+2) > rplat(ITMP)=0. > ENDDO > C-CRA rplon=0. > c real rplon((2*nlath+1)*(nlon+2)) > DO ITMP=1,(2*nlath+1)*(nlon+2) > rplon(ITMP)=0. > ENDDO 316c442,446 < xhatp=0. --- > C-CRA xhatp=0. > c real xhatp((jcap+1)*(jcap+2)) > DO ITMP=1,(jcap+1)*(jcap+2) > xhatp(ITMP)=0. > ENDDO 321c451,455 < phatp=-ghatp --- > C-CRA phatp=-ghatp > c real phatp((jcap+1)*(jcap+2)) > DO ITMP=1,(jcap+1)*(jcap+2) > phatp(ITMP)=-ghatp(ITMP) > ENDDO 345,346c479,484 < end do < fhatp=phatp --- > end do > C-CRA fhatp=phatp > c real fhatp((jcap+1)*(jcap+2)) > DO ITMP=1,(jcap+1)*(jcap+2) > fhatp(ITMP)=phatp(ITMP) > ENDDO 412c550,554 < rt=rt+st --- > C-CRA rt=rt+st > c real rt((2*nlath+1)*(nlon+2)*nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > rt(ITMP)=rt(ITMP)+st(ITMP) > ENDDO 424,426c566,580 < rvort=0. < rplon=0. < rplat=0. --- > C-CRA rvort=0. > c real rvort((2*nlath+1)*(nlon+2)*nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > rvort(ITMP)=0. > ENDDO > C-CRA rplon=0. > c real rplon((2*nlath+1)*(nlon+2)) > DO ITMP=1,(2*nlath+1)*(nlon+2) > rplon(ITMP)=0. > ENDDO > C-CRA rplat=0. > c real rplat((2*nlath+1)*(nlon+2)) > DO ITMP=1,(2*nlath+1)*(nlon+2) > rplat(ITMP)=0. > ENDDO 504,505c658,667 < ghat=xhat < ghatp=xhatp --- > C-CRA ghat=xhat > c real ghat((jcap+1)*(jcap+2),nsig,4) > DO ITMP=1,((jcap+1)*(jcap+2))*nsig*4 > ghat(ITMP,1,1)=xhat(ITMP,1,1) > ENDDO > C-CRA ghatp=xhatp > c real ghatp((jcap+1)*(jcap+2)) > DO ITMP=1,(jcap+1)*(jcap+2) > ghatp(ITMP)=xhatp(ITMP) > ENDDO 509,510c671,680 < fhat=ghat < ps=ghatp --- > C-CRA fhat=ghat > c real fhat((jcap+1)*(jcap+2),nsig,4) > DO ITMP=1,((jcap+1)*(jcap+2))*nsig*4 > fhat(ITMP,1,1)=ghat(ITMP,1,1) > ENDDO > C-CRA ps=ghatp > c real ps((jcap+1)*(jcap+2)) > DO ITMP=1,(jcap+1)*(jcap+2) > ps(ITMP)=ghatp(ITMP) > ENDDO 531c701,705 < ps=factor+ps --- > C-CRA ps=factor+ps > c real ps((jcap+1)*(jcap+2)) > DO ITMP=1,(jcap+1)*(jcap+2) > ps(ITMP)=factor(ITMP)+ps(ITMP) > ENDDO prepp.f 41,51c41,51 < dimension drptps(mpsdat),dlons(mpsdat) < dimension dlats(mpsdat),dps(mpsdat),td(mpsdat) < dimension topogd(mpsdat),tg(nlat+1,nlon+2,nsig) < dimension psg(nlat+1,nlon+2) < dimension topogm(nlat+1,nlon+2) < dimension glats(nlat),glons(nlon),sigl(nsig) < dimension ttypes(mpsdat) < dimension rbtr(mpsdat) < dimension rdelz(mpsdat),rbpres(mpsdat),rbprs2(mpsdat) < dimension rbtps(mpsdat),rbtps2(mpsdat),rdp(mpsdat) < dimension sigll(nsig) --- > C-CRA dimension drptps(mpsdat),dlons(mpsdat) > C-CRA dimension dlats(mpsdat),dps(mpsdat),td(mpsdat) > C-CRA dimension topogd(mpsdat),tg(nlat+1,nlon+2,nsig) > C-CRA dimension psg(nlat+1,nlon+2) > C-CRA dimension topogm(nlat+1,nlon+2) > C-CRA dimension glats(nlat),glons(nlon),sigl(nsig) > C-CRA dimension ttypes(mpsdat) > C-CRA dimension rbtr(mpsdat) > C-CRA dimension rdelz(mpsdat),rbpres(mpsdat),rbprs2(mpsdat) > C-CRA dimension rbtps(mpsdat),rbtps2(mpsdat),rdp(mpsdat) > C-CRA dimension sigll(nsig) 52a53,70 > C dimension drptps(_mpsdat_),dlons(_mpsdat_) > dimension drptps(28000),dlons(28000) > C dimension dlats(_mpsdat_),dps(_mpsdat_),td(_mpsdat_) > dimension dlats(28000),dps(28000),td(28000) > C dimension topogd(_mpsdat_),tg(96+1,192+2,28) > dimension topogd(28000),tg(96+1,192+2,28) > dimension psg(96+1,192+2) > dimension topogm(96+1,192+2) > dimension glats(96),glons(192),sigl(28) > C dimension ttypes(_mpsdat_) > dimension ttypes(28000) > C dimension rbtr(_mpsdat_) > dimension rbtr(28000) > C dimension rdelz(_mpsdat_),rbpres(_mpsdat_),rbprs2(_mpsdat_) > dimension rdelz(28000),rbpres(28000),rbprs2(28000) > C dimension rbtps(_mpsdat_),rbtps2(_mpsdat_),rdp(_mpsdat_) > dimension rbtps(28000),rbtps2(28000),rdp(28000) > dimension sigll(28) 59c77,81 < sigll=log(sigl) --- > C-CRA sigll=log(sigl) > c dimension sigll(nsig) > DO ITMP=1,nsig > sigll(ITMP)=log(sigl(ITMP)) > ENDDO preppw.f prepq.f 41,54c41,54 < dimension drptps(mqdat),dlons(mqdat) < dimension dlats(mqdat),dpres(mqdat) < dimension psg(nlat+1,nlon+2) < dimension glats(nlat),glons(nlon),sigl(nsig) < dimension rtypes(mqdat) < dimension rspres(mqdat) < dimension rmaxerr(mqdat) < dimension temp(nlat+1,nlon+2,nsig) < dimension rbqs(mqdat) < dimension qscal(nlat+1,nlon+2,nsig) < dimension rbpres(mqdat) < dimension sigll(nsig) < dimension rlow(mqdat),rhgh(mqdat) < dimension rbq(mqdat) --- > C-CRA dimension drptps(mqdat),dlons(mqdat) > C-CRA dimension dlats(mqdat),dpres(mqdat) > C-CRA dimension psg(nlat+1,nlon+2) > C-CRA dimension glats(nlat),glons(nlon),sigl(nsig) > C-CRA dimension rtypes(mqdat) > C-CRA dimension rspres(mqdat) > C-CRA dimension rmaxerr(mqdat) > C-CRA dimension temp(nlat+1,nlon+2,nsig) > C-CRA dimension rbqs(mqdat) > C-CRA dimension qscal(nlat+1,nlon+2,nsig) > C-CRA dimension rbpres(mqdat) > C-CRA dimension sigll(nsig) > C-CRA dimension rlow(mqdat),rhgh(mqdat) > C-CRA dimension rbq(mqdat) 55a56,78 > C dimension drptps(_mqdat_),dlons(_mqdat_) > dimension drptps(15000),dlons(15000) > C dimension dlats(_mqdat_),dpres(_mqdat_) > dimension dlats(15000),dpres(15000) > dimension psg(96+1,192+2) > dimension glats(96),glons(192),sigl(28) > C dimension rtypes(_mqdat_) > dimension rtypes(15000) > C dimension rspres(_mqdat_) > dimension rspres(15000) > C dimension rmaxerr(_mqdat_) > dimension rmaxerr(15000) > dimension temp(96+1,192+2,28) > C dimension rbqs(_mqdat_) > dimension rbqs(15000) > dimension qscal(96+1,192+2,28) > C dimension rbpres(_mqdat_) > dimension rbpres(15000) > dimension sigll(28) > C dimension rlow(_mqdat_),rhgh(_mqdat_) > dimension rlow(15000),rhgh(15000) > C dimension rbq(_mqdat_) > dimension rbq(15000) 63c86,90 < sigll=log(sigl) --- > C-CRA sigll=log(sigl) > c dimension sigll(nsig) > DO ITMP=1,nsig > sigll(ITMP)=log(sigl(ITMP)) > ENDDO 78,79c105,114 < rspres=10.*exp(dpres) < dpres=dpres-rbpres --- > C-CRA rspres=10.*exp(dpres) > c dimension rspres(mqdat) > DO ITMP=1,mqdat > rspres(ITMP)=10.*exp(dpres(ITMP)) > ENDDO > C-CRA dpres=dpres-rbpres > c dimension dlats(mqdat),dpres(mqdat) > DO ITMP=1,mqdat > dpres(ITMP)=dpres(ITMP)-rbpres(ITMP) > ENDDO 97c132,136 < qscal=0. --- > C-CRA qscal=0. > c dimension qscal(nlat+1,nlon+2,nsig) > DO ITMP=1,(nlat+1)*(nlon+2)*nsig > qscal(ITMP,1,1)=0. > ENDDO preps.f 36,42c36,42 < dimension dlons(msdat) < dimension dlats(msdat),dpres(msdat) < dimension psg(nlat+1,nlon+2) < dimension glats(nlat),glons(nlon),sigl(nsig) < dimension rspres(msdat) < dimension rbpres(msdat) < dimension sigll(nsig) --- > C-CRA dimension dlons(msdat) > C-CRA dimension dlats(msdat),dpres(msdat) > C-CRA dimension psg(nlat+1,nlon+2) > C-CRA dimension glats(nlat),glons(nlon),sigl(nsig) > C-CRA dimension rspres(msdat) > C-CRA dimension rbpres(msdat) > C-CRA dimension sigll(nsig) 43a44,54 > C dimension dlons(_msdat_) > dimension dlons(120000) > C dimension dlats(_msdat_),dpres(_msdat_) > dimension dlats(120000),dpres(120000) > dimension psg(96+1,192+2) > dimension glats(96),glons(192),sigl(28) > C dimension rspres(_msdat_) > dimension rspres(120000) > C dimension rbpres(_msdat_) > dimension rbpres(120000) > dimension sigll(28) 50c61,65 < sigll=log(sigl) --- > C-CRA sigll=log(sigl) > c dimension sigll(nsig) > DO ITMP=1,nsig > sigll(ITMP)=log(sigl(ITMP)) > ENDDO 66,67c81,90 < rspres=10.*exp(dpres) < dpres=dpres-rbpres --- > C-CRA rspres=10.*exp(dpres) > c dimension rspres(msdat) > DO ITMP=1,msdat > rspres(ITMP)=10.*exp(dpres(ITMP)) > ENDDO > C-CRA dpres=dpres-rbpres > c dimension dlats(msdat),dpres(msdat) > DO ITMP=1,msdat > dpres(ITMP)=dpres(ITMP)-rbpres(ITMP) > ENDDO prept.f 38,47c38,47 < dimension drptps(mtdat),dlons(mtdat) < dimension dlats(mtdat),dpres(mtdat) < dimension psg(nlat+1,nlon+2) < dimension glats(nlat),glons(nlon),sigl(nsig) < dimension rtypes(mtdat) < dimension rspres(mtdat) < dimension rbpres(mtdat) < dimension rbt(mtdat) < dimension rlow(mtdat),rhgh(mtdat) < dimension sigll(nsig) --- > C-CRA dimension drptps(mtdat),dlons(mtdat) > C-CRA dimension dlats(mtdat),dpres(mtdat) > C-CRA dimension psg(nlat+1,nlon+2) > C-CRA dimension glats(nlat),glons(nlon),sigl(nsig) > C-CRA dimension rtypes(mtdat) > C-CRA dimension rspres(mtdat) > C-CRA dimension rbpres(mtdat) > C-CRA dimension rbt(mtdat) > C-CRA dimension rlow(mtdat),rhgh(mtdat) > C-CRA dimension sigll(nsig) 48a49,65 > C dimension drptps(_mtdat_),dlons(_mtdat_) > dimension drptps(60000),dlons(60000) > C dimension dlats(_mtdat_),dpres(_mtdat_) > dimension dlats(60000),dpres(60000) > dimension psg(96+1,192+2) > dimension glats(96),glons(192),sigl(28) > C dimension rtypes(_mtdat_) > dimension rtypes(60000) > C dimension rspres(_mtdat_) > dimension rspres(60000) > C dimension rbpres(_mtdat_) > dimension rbpres(60000) > C dimension rbt(_mtdat_) > dimension rbt(60000) > C dimension rlow(_mtdat_),rhgh(_mtdat_) > dimension rlow(60000),rhgh(60000) > dimension sigll(28) prepw.f 40,47c40,47 < dimension drptps(mwdat),dlons(mwdat),dlats(mwdat),dpres(mwdat) < dimension psg(nlat+1,nlon+2) < dimension glats(nlat),glons(nlon),sigl(nsig) < dimension rtypes(mwdat) < dimension rspres(mwdat),fact(mwdat) < dimension rbpres(mwdat) < dimension rlow(mwdat),rhgh(mwdat) < dimension sigll(nsig+1) --- > C-CRA dimension drptps(mwdat),dlons(mwdat),dlats(mwdat),dpres(mwdat) > C-CRA dimension psg(nlat+1,nlon+2) > C-CRA dimension glats(nlat),glons(nlon),sigl(nsig) > C-CRA dimension rtypes(mwdat) > C-CRA dimension rspres(mwdat),fact(mwdat) > C-CRA dimension rbpres(mwdat) > C-CRA dimension rlow(mwdat),rhgh(mwdat) > C-CRA dimension sigll(nsig+1) 48a49,63 > C dimension drptps(_mwdat_),dlons(_mwdat_) > dimension drptps(112000),dlons(112000) > C dimension dlats(_mwdat_),dpres(_mwdat_) > dimension dlats(112000),dpres(112000) > dimension psg(96+1,192+2) > dimension glats(96),glons(192),sigl(28) > C dimension rtypes(_mwdat_) > dimension rtypes(112000) > C dimension rspres(_mwdat_),fact(_mwdat_) > dimension rspres(112000),fact(112000) > C dimension rbpres(_mwdat_) > dimension rbpres(112000) > C dimension rlow(_mwdat_),rhgh(_mwdat_) > dimension rlow(112000),rhgh(112000) > dimension sigll(28+1) 74,78c89,98 < < rspres=10.*exp(dpres) < < dpres=dpres-rbpres < --- > C-CRA rspres=10.*exp(dpres) > c dimension rspres(mwdat),fact(mwdat) > DO ITMP=1,mwdat > rspres(ITMP)=10.*exp(dpres(ITMP)) > ENDDO > C-CRA dpres=dpres-rbpres > c dimension drptps(mwdat),dlons(mwdat),dlats(mwdat),dpres(mwdat) > DO ITMP=1,mwdat > dpres(ITMP)=dpres(ITMP)-rbpres(ITMP) > ENDDO prnon85.f qoper.f qtoper.f rdfact.f rdgesc.f rdprep7.f 54,62c54,62 < dimension tdata(ntdata,7),wdata(nwdata,8) < dimension sdata(nsdata,5) < dimension psdata(npdata,8),qdata(nqdata,7) < dimension pwdata(npwdat,6) < dimension ttype(ntdata),wtype(nwdata) < dimension stype(nsdata) < dimension pstype(npdata),qtype(nqdata) < dimension pwtype(npwdat),iqtflg(ntdata) < dimension qmaxerr(nqdata),pwmerr(npwdat) --- > C-CRA dimension tdata(ntdata,7),wdata(nwdata,8) > C-CRA dimension sdata(nsdata,5) > C-CRA dimension psdata(npdata,8),qdata(nqdata,7) > C-CRA dimension pwdata(npwdat,6) > C-CRA dimension ttype(ntdata),wtype(nwdata) > C-CRA dimension stype(nsdata) > C-CRA dimension pstype(npdata),qtype(nqdata) > C-CRA dimension pwtype(npwdat),iqtflg(ntdata) > C-CRA dimension qmaxerr(nqdata),pwmerr(npwdat) 63a64,72 > dimension tdata(60000,7),wdata(112000,8) > dimension sdata(120000,5) > dimension psdata(28000,8),qdata(15000,7) > dimension pwdata(1,6) > dimension ttype(60000),wtype(112000) > dimension stype(120000) > dimension pstype(28000),qtype(15000) > dimension pwtype(1),iqtflg(60000) > dimension qmaxerr(15000),pwmerr(1) rdtest9.f 321c321,348 < --- > IF(NTDATA.GT.60000) THEN > PRINT *,'NTDATA too large', NTDATA, 60000 > CALL ABORT > ENDIF > IF(NSDATA.GT.120000) THEN > PRINT *,'NSDATA too large', NSDATA, 120000 > CALL ABORT > ENDIF > IF(NWDATA.GT.112000) THEN > PRINT *,'NWDATA too large', NWDATA,112000 > CALL ABORT > ENDIF > IF(NPDATA.GT.28000) THEN > PRINT *,'NPDATA too large',NPDATA,28000 > CALL ABORT > ENDIF > IF(NQDATA.GT.15000) THEN > PRINT *,'NQDATA too large',NQDATA,15000 > CALL ABORT > ENDIF > IF(NPWDAT.GT.1) THEN > PRINT *,'NPWDAT too large',NPWDAT,1 > CALL ABORT > ENDIF > IF(NSPROF.GT.10000) THEN > PRINT *,'NSPROF too large',NSPROF,10000 > CALL ABORT > ENDIF residw.f respsf.f 26,28c26,28 < dimension dps(mpsdat) < dimension ttypes(mpsdat) < dimension pp(mpsdat) --- > C-CRA dimension dps(mpsdat) > C-CRA dimension ttypes(mpsdat) > C-CRA dimension pp(mpsdat) 29a30,35 > C dimension dps(_mpsdat_) > C dimension ttypes(_mpsdat_) > C dimension pp(_mpsdat_) > dimension dps(28000) > dimension ttypes(28000) > dimension pp(28000) 34c40,44 < pp=1000. --- > C-CRA pp=1000. > c dimension pp(mpsdat) > DO ITMP=1,mpsdat > pp(ITMP)=1000. > ENDDO respw.f 29,33c29,33 < dimension drptpw(mpwdat) < dimension dpw(mpwdat) < dimension ttypes(mpwdat) < dimension rmaxerr(mpwdat) < dimension rspres(mpwdat) --- > C-CRA dimension drptpw(mpwdat) > C-CRA dimension dpw(mpwdat) > C-CRA dimension ttypes(mpwdat) > C-CRA dimension rmaxerr(mpwdat) > C-CRA dimension rspres(mpwdat) 34a35,44 > C dimension drptpw(_mpwdat_) > C dimension dpw(_mpwdat_) > C dimension ttypes(_mpwdat_) > C dimension rmaxerr(_mpwdat_) > C dimension rspres(_mpwdat_) > dimension drptpw(1) > dimension dpw(1) > dimension ttypes(1) > dimension rmaxerr(1) > dimension rspres(1) resq.f ressat.f restmp.f s2g0.f s2grad.f s2gvec.f s2mg2x.f satc2.f satcov.f satop46.f sdot.f setuprhs5.f 70,122c70,180 < dimension trigs(nlon*2),ifax(10),rlats(nlath*2) < dimension isatv(nsig) < dimension rt(2*nlath+1,nlon+2,nsig) < dimension ru(2*nlath+1,nlon+2,nsig) < dimension rv(2*nlath+1,nlon+2,nsig) < dimension rq(2*nlath+1,nlon+2,nsig) < dimension rp(2*nlath+1,nlon+2) < dimension in((jcap+1)*(jcap+2)) < dimension dstlast((jcap+1)*(jcap+2),nsig) < dimension dstb((jcap+1)*(jcap+2),nsig) < dimension pln((jcap+1)*(jcap+2),nlath) < dimension qln((jcap+1)*(jcap+2),nlath) < dimension rln((jcap+1)*(jcap+2),nlath) < dimension mlad(0:jcap,0:jcap) < dimension ml2lm((jcap+1)*(jcap+2)) < dimension factslm((jcap+1)*(jcap+2)) < dimension factvlm((jcap+1)*(jcap+2)) < dimension lmad(0:jcap,0:jcap) < dimension lm2ml((jcap+1)*(jcap+2)) < dimension factsml((jcap+1)*(jcap+2)) < dimension factvml((jcap+1)*(jcap+2)) < dimension rts(2*nlath+1,nlon+2,nsig) < dimension rus(2*nlath+1,nlon+2,nsig) < dimension rvs(2*nlath+1,nlon+2,nsig) < dimension rvorts(2*nlath+1,nlon+2,nsig) < dimension rplons(2*nlath+1,nlon+2) < dimension rplats(2*nlath+1,nlon+2) < dimension rlons(nlon) < dimension tdata(ntdata,8),wdata(nwdata,8) < dimension sdata(nsdata,6) < dimension psdata(npdata,8),qdata(nqdata,7) < dimension pwdata(npwdat,6) < dimension rpw(2*nlath+1,nlon+2) < dimension sigi(nsig+1) < dimension tempt(ntdata),temps(nsdata),tempq(nqdata) < dimension tempp(npdata),temppw(npwdat),tempw(nwdata) < dimension tpres(ntdata),qpres(nqdata),wpres(nwdata) < dimension spres(nsdata) < dimension tges(ntdata),qges(nqdata),vges(nwdata),sges(nsdata) < dimension pges(npdata),uges(nwdata),pwges(npwdat) < dimension rqtype(nqdata),rwtype(nwdata),qtges(ntdata) < dimension iqtflg(ntdata) < dimension ptype(npdata),pwtype(npwdat) < dimension qmaxerr(nqdata),pwmerr(npwdat) < dimension fact(nwdata) < dimension rbqs(nqdata) < dimension factor(2*nlath+1,nlon+2) < dimension qfile(17*nqdata),uvfile(18*nwdata),tfile(17*ntdata) < dimension sfile((4+(28+2)*30)*nsprof),pwfile(12*npwdat) < dimension psfile(11*npdata) < real del2((jcap+1)*(jcap+2)) < real pwcon(nsig) < real wgts(nlath*2),sigl(nsig) --- > C-CRA dimension trigs(nlon*2),ifax(10),rlats(nlath*2) > C-CRA dimension isatv(nsig) > C-CRA dimension rt(2*nlath+1,nlon+2,nsig) > C-CRA dimension ru(2*nlath+1,nlon+2,nsig) > C-CRA dimension rv(2*nlath+1,nlon+2,nsig) > C-CRA dimension rq(2*nlath+1,nlon+2,nsig) > C-CRA dimension rp(2*nlath+1,nlon+2) > C-CRA dimension in((jcap+1)*(jcap+2)) > C-CRA dimension dstlast((jcap+1)*(jcap+2),nsig) > C-CRA dimension dstb((jcap+1)*(jcap+2),nsig) > C-CRA dimension pln((jcap+1)*(jcap+2),nlath) > C-CRA dimension qln((jcap+1)*(jcap+2),nlath) > C-CRA dimension rln((jcap+1)*(jcap+2),nlath) > C-CRA dimension mlad(0:jcap,0:jcap) > C-CRA dimension ml2lm((jcap+1)*(jcap+2)) > C-CRA dimension factslm((jcap+1)*(jcap+2)) > C-CRA dimension factvlm((jcap+1)*(jcap+2)) > C-CRA dimension lmad(0:jcap,0:jcap) > C-CRA dimension lm2ml((jcap+1)*(jcap+2)) > C-CRA dimension factsml((jcap+1)*(jcap+2)) > C-CRA dimension factvml((jcap+1)*(jcap+2)) > C-CRA dimension rts(2*nlath+1,nlon+2,nsig) > C-CRA dimension rus(2*nlath+1,nlon+2,nsig) > C-CRA dimension rvs(2*nlath+1,nlon+2,nsig) > C-CRA dimension rvorts(2*nlath+1,nlon+2,nsig) > C-CRA dimension rplons(2*nlath+1,nlon+2) > C-CRA dimension rplats(2*nlath+1,nlon+2) > C-CRA dimension rlons(nlon) > C-CRA dimension tdata(ntdata,8),wdata(nwdata,8) > C-CRA dimension sdata(nsdata,6) > C-CRA dimension psdata(npdata,8),qdata(nqdata,7) > C-CRA dimension pwdata(npwdat,6) > C-CRA dimension rpw(2*nlath+1,nlon+2) > C-CRA dimension sigi(nsig+1) > C-CRA dimension tempt(ntdata),temps(nsdata),tempq(nqdata) > C-CRA dimension tempp(npdata),temppw(npwdat),tempw(nwdata) > C-CRA dimension tpres(ntdata),qpres(nqdata),wpres(nwdata) > C-CRA dimension spres(nsdata) > C-CRA dimension tges(ntdata),qges(nqdata),vges(nwdata),sges(nsdata) > C-CRA dimension pges(npdata),uges(nwdata),pwges(npwdat) > C-CRA dimension rqtype(nqdata),rwtype(nwdata),qtges(ntdata) > C-CRA dimension iqtflg(ntdata) > C-CRA dimension ptype(npdata),pwtype(npwdat) > C-CRA dimension qmaxerr(nqdata),pwmerr(npwdat) > C-CRA dimension fact(nwdata) > C-CRA dimension rbqs(nqdata) > C-CRA dimension factor(2*nlath+1,nlon+2) > C-CRA dimension qfile(17*nqdata),uvfile(18*nwdata),tfile(17*ntdata) > C-CRA dimension sfile((4+(28+2)*30)*nsprof),pwfile(12*npwdat) > C-CRA dimension psfile(11*npdata) > C-CRA real del2((jcap+1)*(jcap+2)) > C-CRA real pwcon(nsig) > C-CRA real wgts(nlath*2),sigl(nsig) > C-CRA integer idateg(4) > C-CRA integer idateg2(4),idate5(5) > > dimension trigs(192*2),ifax(10),rlats(48*2) > dimension isatv(28) > dimension rt(2*48+1,192+2,28) > dimension ru(2*48+1,192+2,28) > dimension rv(2*48+1,192+2,28) > dimension rq(2*48+1,192+2,28) > dimension rp(2*48+1,192+2) > dimension in((62+1)*(62+2)) > dimension dstlast((62+1)*(62+2),28) > dimension dstb((62+1)*(62+2),28) > dimension pln((62+1)*(62+2),48) > dimension qln((62+1)*(62+2),48) > dimension rln((62+1)*(62+2),48) > dimension mlad(0:62,0:62) > dimension ml2lm((62+1)*(62+2)) > dimension factslm((62+1)*(62+2)) > dimension factvlm((62+1)*(62+2)) > dimension lmad(0:62,0:62) > dimension lm2ml((62+1)*(62+2)) > dimension factsml((62+1)*(62+2)) > dimension factvml((62+1)*(62+2)) > dimension rts(2*48+1,192+2,28) > dimension rus(2*48+1,192+2,28) > dimension rvs(2*48+1,192+2,28) > dimension rvorts(2*48+1,192+2,28) > dimension rplons(2*48+1,192+2) > dimension rplats(2*48+1,192+2) > dimension rlons(192) > dimension tdata(60000,8),wdata(112000,8) > dimension sdata(120000,6) > dimension psdata(28000,8),qdata(15000,7) > dimension pwdata(1,6) > dimension rpw(2*48+1,192+2) > dimension sigi(28+1) > dimension tempt(60000),temps(120000),tempq(15000) > dimension tempp(28000),temppw(1),tempw(112000) > dimension tpres(60000),qpres(15000),wpres(112000) > dimension spres(120000) > dimension tges(60000),qges(15000) > dimension vges(112000),sges(120000) > dimension pges(28000),uges(112000),pwges(1) > dimension rqtype(15000),rwtype(112000),qtges(60000) > dimension iqtflg(60000) > dimension ptype(28000),pwtype(1) > dimension qmaxerr(15000),pwmerr(1) > dimension fact(112000) > dimension rbqs(15000) > dimension factor(2*48+1,192+2) > dimension qfile(17*15000),uvfile(18*112000) > dimension tfile(17*60000) > dimension sfile((4+(28+2)*30)*10000),pwfile(12*1) > dimension psfile(11*28000) > real del2((62+1)*(62+2)) > real pwcon(28) > real wgts(48*2),sigl(28) 125d182 < 131c188,192 < isatv=1 --- > C-CRA isatv=1 > c dimension isatv(nsig) > DO ITMP=1,nsig > isatv(ITMP)=1 > ENDDO 260,262c321,335 < rus=ru < rvs=rv < rts=rt --- > C-CRA rus=ru > c dimension rus(2*nlath+1,nlon+2,nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > rus(ITMP,1,1)=ru(ITMP,1,1) > ENDDO > C-CRA rvs=rv > c dimension rvs(2*nlath+1,nlon+2,nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > rvs(ITMP,1,1)=rv(ITMP,1,1) > ENDDO > C-CRA rts=rt > c dimension rts(2*nlath+1,nlon+2,nsig) > DO ITMP=1,(2*nlath+1)*(nlon+2)*nsig > rts(ITMP,1,1)=rt(ITMP,1,1) > ENDDO 338c411,415 < qtges=0. --- > C-CRA qtges=0. > c dimension rqtype(nqdata),rwtype(nwdata),qtges(ntdata) > DO ITMP=1,ntdata > qtges(ITMP)=0. > ENDDO 509c586,590 < dstb=dstlast --- > C-CRA dstb=dstlast > c dimension dstb((jcap+1)*(jcap+2),nsig) > DO ITMP=1,(jcap+1)*(jcap+2)*nsig > dstb(ITMP,1)=dstlast(ITMP,1) > END DO sgemv.f sprp2.f 33,35c33,35 < dimension pdata(npsdta,7),pges(npsdta),ptype(npsdta) < dimension icount(nlat,nlon) < dimension ncount(npsdta),jl(128,4) --- > C-CRA dimension pdata(npsdta,7),pges(npsdta),ptype(npsdta) > C-CRA dimension icount(nlat,nlon) > C-CRA dimension ncount(npsdta),jl(128,4) 36a37,44 > C dimension pdata(_npsdta_,7),pges(_npsdta_) > dimension pdata(28000,7),pges(28000) > C dimension ptype(_npsdta_) > dimension ptype(28000) > dimension icount(96,192) > C dimension ncount(_npsdta_),jl(128,4) > dimension ncount(28000),jl(128,4) > 44c52,56 < ncount=0 --- > C-CRA ncount=0 > c dimension ncount(npsdta),jl(128,4) > DO ITMP=1,npsdta > ncount(ITMP)=0 > ENDDO 69c81,85 < icount=0 --- > C-CRA icount=0 > c dimension icount(nlat,nlon) > DO ITMP=1,nlat*nlon > icount(ITMP,1)=0 > ENDDO sprqpw2.f 43,50c43,50 < dimension pwdata(npwdta,5) < dimension pwges(npwdta),pwtype(npwdta) < dimension qdata(nqdta,6) < dimension qges(nqdta),qtype(nqdta) < dimension numq(nsig) < dimension qplty(nsig) < dimension icount(nlat,nlon,nsig),ncount(nqdta) < dimension rbqs(nqdta) --- > C-CRA dimension pwdata(npwdta,5) > C-CRA dimension pwges(npwdta),pwtype(npwdta) > C-CRA dimension qdata(nqdta,6) > C-CRA dimension qges(nqdta),qtype(nqdta) > C-CRA dimension numq(nsig) > C-CRA dimension qplty(nsig) > C-CRA dimension icount(nlat,nlon,nsig),ncount(nqdta) > C-CRA dimension rbqs(nqdta) 51a52,66 > C dimension pwdata(_npwdta_,5) > dimension pwdata(1,5) > C dimension pwges(_npwdta_),pwtype(_npwdta_) > dimension pwges(1),pwtype(1) > C dimension qdata(_nqdta_,6) > dimension qdata(15000,6) > C dimension qges(_nqdta_),qtype(_nqdta_) > dimension qges(15000),qtype(15000) > dimension numq(28) > dimension qplty(28) > C dimension icount(96,192,28),ncount(_nqdta_) > dimension icount(96,192,28),ncount(15000) > C dimension rbqs(_nqdta_) > dimension rbqs(15000) > 162c177,181 < qplty=0. --- > C-CRA qplty=0. > c dimension qplty(nsig) > DO ITMP=1,nsig > qplty(ITMP)=0. > ENDDO 172c191,195 < ncount=0 --- > C-CRA ncount=0 > c dimension icount(nlat,nlon,nsig),ncount(nqdta) > DO ITMP=1,nqdta > ncount(ITMP)=0 > ENDDO 179c202,206 < icount=0 --- > C-CRA icount=0 > c dimension icount(nlat,nlon,nsig),ncount(nqdta) > DO ITMP=1,nlat*nlon*nsig > icount(ITMP,1,1)=0 > ENDDO sprs6.f 35,47c35,67 < dimension tdata(ntdta,6) < dimension rr(nsig),saterr(nsig,nsig,2) < dimension saterrin(nsigsat,nsigsat,2) < dimension scrat(3*nsig) < dimension iscrat(3*nsig) < dimension numtemps(nsig),rr2(200),vloc(200),alsig(nsig) < dimension indl(nsig) < dimension el(nsig*nsig) < dimension aerr(nsig*nsig) < dimension ist(2*nlath+1,nlon+2) < dimension obserror(nsig,2) < dimension sigl(nsig),siglin(nsigsat) < dimension iold(nsig) --- > C-CRA dimension tdata(ntdta,6) > C-CRA dimension rr(nsig),saterr(nsig,nsig,2) > C-CRA dimension saterrin(nsigsat,nsigsat,2) > C-CRA dimension scrat(3*nsig) > C-CRA dimension iscrat(3*nsig) > C-CRA dimension numtemps(nsig),rr2(200),vloc(200),alsig(nsig) > C-CRA dimension indl(nsig) > C-CRA dimension el(nsig*nsig) > C-CRA dimension aerr(nsig*nsig) > C-CRA dimension ist(2*nlath+1,nlon+2) > C-CRA dimension obserror(nsig,2) > C-CRA dimension sigl(nsig),siglin(nsigsat) > C-CRA dimension iold(nsig) > C-CRA dimension sfile(*) > C-CRA real rrs(2*nlath+1,nlon+2,nsig) > C-CRA real satf(2*nlath+1,nlon+2,nsig) > C-CRA logical odiag > > C WARNING!!! dimension tdata(_ntdta_,6) is wrong. > C It must use the dimension of the calling program. > dimension tdata(120000,6) > dimension rr(28),saterr(28,28,2) > dimension saterrin(28,28,2) > dimension scrat(3*28) > dimension iscrat(3*28) > dimension numtemps(28),rr2(200),vloc(200),alsig(28) > dimension indl(28) > dimension el(28*28) > dimension aerr(28*28) > dimension ist(2*48+1,192+2) > dimension obserror(28,2) > dimension sigl(28),siglin(28) > dimension iold(28) 49,50c69,70 < real rrs(2*nlath+1,nlon+2,nsig) < real satf(2*nlath+1,nlon+2,nsig) --- > real rrs(2*48+1,192+2,28) > real satf(2*48+1,192+2,28) 52,53d71 < < sprt2.f 31,34c31,34 < dimension tdata(ntdta,8),tges(ntdta) < dimension numtemps(nsig) < dimension tplty(nsig) < dimension icount(nlat,nlon,nsig),ncount(ntdta+2*nlat*nlon) --- > C-CRA dimension tdata(ntdta,8),tges(ntdta) > C-CRA dimension numtemps(nsig) > C-CRA dimension tplty(nsig) > C-CRA dimension icount(nlat,nlon,nsig),ncount(ntdta+2*nlat*nlon) 35a36,45 > C Warning !! dimension tdata(_ntdta_,8),tges(_ntdta_) > dimension tdata(60000,8),tges(60000) > c > dimension numtemps(28) > dimension tplty(28) > dimension icount(96,192,28) > c > C dimension ncount(_ntdta_+2*96*192) > dimension ncount(60000+2*96*192) > c spruv2.f 35,40c35,40 < dimension wdata(nwdta,7) < dimension uges(nwdta),vges(nwdta),wtype(nwdta) < dimension numw(nsig) < dimension fact(nwdta) < dimension uplty(nsig),vplty(nsig) < dimension icount(nlat,nlon,nsig),ncount(nwdta) --- > C-CRA dimension wdata(nwdta,7) > C-CRA dimension uges(nwdta),vges(nwdta),wtype(nwdta) > C-CRA dimension numw(nsig) > C-CRA dimension fact(nwdta) > C-CRA dimension uplty(nsig),vplty(nsig) > C-CRA dimension icount(nlat,nlon,nsig),ncount(nwdta) 41a42,52 > C dimension wdata(_nwdta_,7) > dimension wdata(112000,7) > C dimension uges(_nwdta_),vges(_nwdta_),wtype(_nwdta_) > dimension uges(112000),vges(112000),wtype(112000) > dimension numw(28) > C dimension fact(_nwdta_) > dimension fact(112000) > dimension uplty(28),vplty(28) > C dimension icount(96,192,28),ncount(_nwdta_) > dimension icount(96,192,28),ncount(112000) > 50c61,65 < numw=0 --- > C-CRA numw=0 > c dimension numw(nsig) > DO ITMP=1,nsig > numw(ITMP)=0 > ENDDO 61a77,81 > C-CRA uplty=0. > c dimension uplty(nsig),vplty(nsig) > DO ITMP=1,nsig > uplty(ITMP)=0. > ENDDO 62a83,87 > C-CRA vplty=0. > c dimension uplty(nsig),vplty(nsig) > DO ITMP=1,nsig > vplty(ITMP)=0. > ENDDO 70c95,99 < ncount=0 --- > C-CRA ncount=0 > c dimension icount(nlat,nlon,nsig),ncount(nwdta) > DO ITMP=1,nwdta > ncount(ITMP)=0 > ENDDO 74c103,107 < icount=0 --- > C-CRA icount=0 > c dimension icount(nlat,nlon,nsig),ncount(nwdta) > DO ITMP=1,nlat*nlon*nsig > icount(ITMP,1,1)=0 > ENDDO ssi2.f tg2s0.f tgrad2s.f 31,42c31,42 < dimension ds((jcap+1)*(jcap+2)) < dimension u(2*nlath+1,nlon+2) < dimension v(2*nlath+1,nlon+2) < dimension qln((jcap+1)*(jcap+2),nlath) < dimension rln((jcap+1)*(jcap+2),nlath) < dimension trigs(2*nlon),ifax(10) < dimension wgts(2*nlath) < dimension del2((jcap+1)*(jcap+2)) < dimension work(2*(2*nlath+1)*(nlon+2)) < dimension ue(2*(jcap+1)),uo(2*(jcap+1)) < dimension ve(2*(jcap+1)),vo(2*(jcap+1)) < dimension factor(2*jcap+2,nlath) --- > C-CRA dimension ds((jcap+1)*(jcap+2)) > C-CRA dimension u(2*nlath+1,nlon+2) > C-CRA dimension v(2*nlath+1,nlon+2) > C-CRA dimension qln((jcap+1)*(jcap+2),nlath) > C-CRA dimension rln((jcap+1)*(jcap+2),nlath) > C-CRA dimension trigs(2*nlon),ifax(10) > C-CRA dimension wgts(2*nlath) > C-CRA dimension del2((jcap+1)*(jcap+2)) > C-CRA dimension work(2*(2*nlath+1)*(nlon+2)) > C-CRA dimension ue(2*(jcap+1)),uo(2*(jcap+1)) > C-CRA dimension ve(2*(jcap+1)),vo(2*(jcap+1)) > C-CRA dimension factor(2*jcap+2,nlath) 43a44,55 > dimension ds((62+1)*(62+2)) > dimension u(2*48+1,192+2) > dimension v(2*48+1,192+2) > dimension qln((62+1)*(62+2),48) > dimension rln((62+1)*(62+2),48) > dimension trigs(2*192),ifax(10) > dimension wgts(2*48) > dimension del2((62+1)*(62+2)) > dimension work(2*(2*48+1)*(192+2)) > dimension ue(2*(62+1)),uo(2*(62+1)) > dimension ve(2*(62+1)),vo(2*(62+1)) > dimension factor(2*62+2,48) ts2g0.f ts2grad.f ts2gvec.f w3fs03.f w3fs11.f w3fs21.f wranlc.f 34,45c34,45 < dimension zc((jcap+1)*(jcap+2),nsig) < dimension dc((jcap+1)*(jcap+2),nsig) < dimension tc((jcap+1)*(jcap+2),nsig) < dimension qc((jcap+1)*(jcap+2),nsig) < dimension pc((jcap+1)*(jcap+2)) < dimension rc((jcap+1)*(jcap+2)) < dimension idateg(4),sigi(nsig+1),sigl(nsig) < dimension lm2ml((jcap+1)*(jcap+2)) < dimension factsml((jcap+1)*(jcap+2)) < dimension factvml((jcap+1)*(jcap+2)) < dimension dummy(201-nsig-1-nsig) < dimension z((jcap+1)*(jcap+2)) --- > C-CRA dimension zc((jcap+1)*(jcap+2),nsig) > C-CRA dimension dc((jcap+1)*(jcap+2),nsig) > C-CRA dimension tc((jcap+1)*(jcap+2),nsig) > C-CRA dimension qc((jcap+1)*(jcap+2),nsig) > C-CRA dimension pc((jcap+1)*(jcap+2)) > C-CRA dimension rc((jcap+1)*(jcap+2)) > C-CRA dimension idateg(4),sigi(nsig+1),sigl(nsig) > C-CRA dimension lm2ml((jcap+1)*(jcap+2)) > C-CRA dimension factsml((jcap+1)*(jcap+2)) > C-CRA dimension factvml((jcap+1)*(jcap+2)) > C-CRA dimension dummy(201-nsig-1-nsig) > C-CRA dimension z((jcap+1)*(jcap+2)) 46a47,58 > dimension zc((62+1)*(62+2),28) > dimension dc((62+1)*(62+2),28) > dimension tc((62+1)*(62+2),28) > dimension qc((62+1)*(62+2),28) > dimension pc((62+1)*(62+2)) > dimension rc((62+1)*(62+2)) > dimension idateg(4),sigi(28+1),sigl(28) > dimension lm2ml((62+1)*(62+2)) > dimension factsml((62+1)*(62+2)) > dimension factvml((62+1)*(62+2)) > dimension dummy(201-28-1-28) > dimension z((62+1)*(62+2)) 68c80,84 < dummy=0. --- > C-CRA dummy=0. > c dimension dummy(201-nsig-1-nsig) > DO ITMP=1,201-nsig-1-nsig > dummy(ITMP)=0. > ENDDO