------------------------------------- Doing diff of addw.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//addw.f 2,3c2 < include 'params.h' < parameter (imx=ima , jmx=jma, nmx=24) --- > parameter (imx=150 , jmx=150, nmx=24) ------------------------------------- Doing diff of amatrix.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//amatrix.f ------------------------------------- Doing diff of asym.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//asym.f 2c2 < include 'params.h' --- > parameter (isize=25,ires=2,ngprd=6) 4,6c4,6 < parameter (imx=ima, jmx=jma) < parameter (imx1=imx*ngprd, jmx1=jmx*ngprd) < parameter (imf =npts, jmf =npts) --- > parameter (imx=150, jmx=150) > parameter (imx1=imx*6, jmx1=jmx*6) > parameter (imf=npts, jmf=npts) ------------------------------------- Doing diff of bound.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//bound.f 3,4c3 < include 'params.h' < parameter (imx=ima , jmx=jma) --- > parameter (imx=150 , jmx=150) 15,16c14,15 < x=ro(i)/fact*cos(theta)+xc + 1./float(ires) < y=ro(i)*sin(theta)+yc + 1./float(ires) --- > x=ro(i)/fact*cos(theta)+xc + .5 > y=ro(i)*sin(theta)+yc + .5 ------------------------------------- Doing diff of bound2.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//bound2.f 2,3c2 < include 'params.h' < parameter(imx=ima, jmx=jma,nmx=24) --- > parameter(imx=150, jmx=150,nmx=24) ------------------------------------- Doing diff of calcr.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//calcr.f 2,4c2,3 < include 'params.h' < parameter (nmx=24) < parameter (imx=ima , jmx=jma) --- > parameter ( nmx=24) > parameter (imx=150 , jmx=150) 13c12 < pi = 4.*atan(1.0) --- > pi = 4.*atan(1.0) ------------------------------------- Doing diff of calcra.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//calcra.f 3,4c3 < include 'params.h' < parameter (imx=ima , jmx=jma) --- > parameter (imx=150 , jmx=150) 23,24d21 < < c rayr: why a +1 in the 1/2 degree version and not in the 1 degree version? 27d23 < ------------------------------------- Doing diff of calct.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//calct.f 5d4 < include 'params.h' 7c6 < parameter (imx=ima , jmx=jma) --- > parameter (imx=150 , jmx=150, iimx=110) 83c82 < cccc if(irend.eq.idst(i).or.iravg.eq.idst(i))then <---- remove dist as lower endpoint.... --- > cccc if(irend.eq.idst(i).or.iravg.eq.idst(i))then <---- remove di 111c110 < c make sure the search does not start less then 200 km from the center: --- > c make sure the search does not start less then 200 km from the cen ------------------------------------- Doing diff of center.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//center.f 3,5c3,4 < include 'params.h' < parameter (imx=ima, jmx=jma) < parameter ( kmax=kma, lgi=20 ) --- > parameter (imx=150, jmx=150) > parameter ( kmax=42, lgi=20 ) ------------------------------------- Doing diff of check.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//check.f ------------------------------------- Doing diff of compute.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//compute.f 3,5c3,4 < include 'params.h' < parameter (imx=ima, jmx=ima, nmx=24) < parameter ( kmax=kma) --- > parameter (imx=150, jmx=150, nmx=24) > parameter ( kmax=42) 160,164c159 < cc if(irort.eq.1)then < c < c do not allow code to modify surface winds (2013 test) < c < if(irort.eq.2)then --- > if(irort.eq.1)then ------------------------------------- Doing diff of fdump.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//fdump.f ------------------------------------- Doing diff of findra.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//findra.f 5,6c5 < include 'params.h' < parameter(imx=ima,jmx=jma,nmx=24) --- > parameter(imx=150, jmx=150,nmx=24,iimx=110) ------------------------------------- Doing diff of h12.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//h12.f ------------------------------------- Doing diff of i1mach.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//i1mach.f ------------------------------------- Doing diff of input.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//input.f 34,35c34 < include 'params.h' < parameter (imx=ima, jmx=jma,nmx=24) --- > parameter (imx=150 , jmx=150, nmx=24) 50,51c49,50 < c call start ! jim tuccillo's resource tracker < c call w3tagb('gfdl_rmenvir',2001,0081,0073,'np22') --- > call start ! jim tuccillo's resource tracker > call w3tagb('gfdl_rmenvir',2001,0081,0073,'np22') 641c640 < irort = 1 --- > irort = 1 877,878c876,877 < c call summary ! jim tuccillo's resource tracker < c call w3tage('gfdl_rmenvir') --- > call summary ! jim tuccillo's resource tracker > call w3tage('gfdl_rmenvir') ------------------------------------- Doing diff of interad.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//interad.f 3,4c3 < include 'params.h' < parameter (npts=isize*ires*ngprd) --- > parameter (isize=25,ires=2,ngprd=6,npts=isize*ires*ngprd) 7,8c6,7 < parameter (inmax = ima, jnmax = jma) < parameter (imaxn = ima, jmaxn = jma) --- > parameter (inmax = 150, jnmax = 150) > parameter (imaxn = 150, jmaxn = 150) ------------------------------------- Doing diff of interp.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//interp.f ------------------------------------- Doing diff of isamax.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//isamax.f NOTE: File isamax.f does not exist in directory /nwprod/sorc/gfdl_rmenvir.fd/ ------------------------------------- Doing diff of j4save.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//j4save.f ------------------------------------- Doing diff of maxmag.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//maxmag.f ------------------------------------- Doing diff of maxth.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//maxth.f 2,3c2 < include 'params.h' < parameter(nmx=8,imx=ima,jmx=jma,lgth=30) --- > parameter(nmx=8,imx=150,jmx=150,lgth=30,iimx=110) 53c52 < if(yv .lt. 0.0) then --- > if(yv .LT. 0.0) then 66c65 < do 54 ir=2,lgth-1 --- > do 54 ir=2,lgth 88,92c87,88 < c ycn=float(mod(ipos,1000))-2. < c xcn=float(ipos/1000)-2. < c rayr: not totally sure that these next two lines are correct... < ycn=float(mod(ipos,1000))-1.*float(ires) < xcn=float(ipos/1000)-1.*float(ires) --- > ycn=float(mod(ipos,1000))-2. > xcn=float(ipos/1000)-2. 114c110 < print*,'new center: ',xcn,ycn,ipos --- > print*,'new center: ',ycn,xcn,ipos 139c135 < if(yv .lt. 0.0) then --- > if(yv .LT. 0.0) then ------------------------------------- Doing diff of maxw.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//maxw.f 3,4c3,4 < include 'params.h' < parameter (imx=ima,jmx=jma,kmax=kma,nmx=24) --- > parameter (imx=150, jmx=150, nmx=24,kmax=42) > c 129,130c129 < windmx = 100.*1.25*windx - tfact*tanmf < cccc windmx = 100.*1.20*windx - tfact*tanmf --- > windmx = 100.*1.20*windx - tfact*tanmf ------------------------------------- Doing diff of modules.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//modules.f 6a7 > ------------------------------------- Doing diff of phase.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//phase.f ------------------------------------- Doing diff of readdropw.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//readdropw.f 5c5 < c this routine will read in the dropwinds and determine c --- > c this routine will read in the dropwinds and determine c 11,12c11 < include 'params.h' < parameter (imx=ima,jmx=jma,kmax=kma,nmx=24) --- > parameter (imx=150 , jmx=150, kmax = 18,nmx=24) 43,51c42,44 < data q /.0020416, .0076468, .0138312, .0209377, .0291166, < a .0385169, .0492901, .0615968, .0756123, .0915121, < b .1094567, .1296058, .1521000, .1770445, .2045146, < c .2345257, .2670232, .3018821, .3389029, .3778005, < d .4182100, .4597012, .5018039, .5440080, .5857977, < e .6266828, .6662129, .7039977, .7397218, .7731550, < f .8041373, .8325886, .8585089, .8819433, .9029819, < g .9217599, .9384272, .9531540, .9661055, .9774466, < h .9873524, .9959831/ --- > data q /.0207469,.0739862,.1244004, .1745733, .2246687,.2747291, > * .3247711,.3748014,.4248250, .4974484, .5935378,.6881255, > * .7772229,.8563145,.9204018, .9604809, .9814907,.9949968/ 237c230 < c determine epsoln as a function of the distance of the point from the storm --- > c determine epsoln as a function of the distance of the point from the s 338,342c331,335 < cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx < c x < c also compute the distance of each dropwind from the dopolar storm center x < c x < cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx --- > cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx > c > c also compute the distance of each dropwind from the dopolar storm cent > c > cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx ------------------------------------- Doing diff of readit.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//readit.f 3,5c3,4 < include 'params.h' < parameter (imx=ima, jmx=jma,nmx=24) < parameter ( kmax=kma,lgi=20) --- > parameter (imx=150 , jmx=150, nmx=24) > parameter ( kmax=42,lgi=20 ,iimx=110) 30c29 < dimension up(imx,jmx),vp(imx,jmx),q(kmax) --- > dimension up(imx,jmx),vp(imx,jmx),q(kmax) 108c107 < if(pcur .lt. 993.) ifl = 3 --- > if(pcur .lt. 993.) ifl = 3 127,130d125 < c rayr: changed this section so that center is not called if < c we don't have any vortex to remove. < < if (irem .eq. 1) then 132,138c127 < else < xcg=xv < ycg=yv < res = (del(2,1) - del(1,1))/pi180 < print*,'resolution of outer grid (degrees) :',res < end if < print*,'xcg,ycg=',xcg,ycg --- > c 148,158c137,142 < < c rayr: eliminate maxth call if we are not removing global vortex. < if (irem .eq. 1) then < call maxth(up,vp,xcgnew,ycgnew,rmxlim,tanp,res1) < xold=xcgnew * 1.0 + xcorn < yold=ycgnew * 1.0 + ycorn < dist= rmxlim < else < rfavg=.1 < write(15,*) rfavg,pcur < end if --- > call maxth(up,vp,xcgnew,ycgnew,rmxlim,tanp,res) > xold=xcgnew * 1.0 + xcorn > yold=ycgnew * 1.0 + ycorn > cc test xold= xcg / pi180 > cc test yold = ycg / pi180 > dist= rmxlim 193c177 < ccc write(77,7700)((tanp(i,j)/100.,i=ixc-5,ixc+5),j=iyc+5,iyc-5,-1) --- > ccc write(77,7700)((tanp(i,j)/100.,i=ixc-5,ixc+5),j=iyc+5,iyc-5,-1 219a204 > c call calcra(r,rtan,iang) 227c212 < c print*,'r and rtan: ',r,rtan --- > cc print*,'r and rtan: ',r,rtan 250d234 < c 255c239 < c call calcra(rzr,rtan,iang) --- > c call calcra(rzr,rtan,iang) 282d265 < c 302c285 < else ! if irem=0 --- > else 307c290 < endif ! end if irem .ne. 0 --- > endif 309,311c292,293 < xcg1 = xcg/pi180 - 360. < ycg1=ycg/pi180 < write(6,455)xcg1,ycg1 --- > xcg1 = xcg - 360. > write(6,455)xcg1,ycg ------------------------------------- Doing diff of rodist.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//rodist.f ------------------------------------- Doing diff of sasum.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//sasum.f NOTE: File sasum.f does not exist in directory /nwprod/sorc/gfdl_rmenvir.fd/ ------------------------------------- Doing diff of saxpy.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//saxpy.f NOTE: File saxpy.f does not exist in directory /nwprod/sorc/gfdl_rmenvir.fd/ ------------------------------------- Doing diff of scopy.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//scopy.f NOTE: File scopy.f does not exist in directory /nwprod/sorc/gfdl_rmenvir.fd/ ------------------------------------- Doing diff of sdot.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//sdot.f NOTE: File sdot.f does not exist in directory /nwprod/sorc/gfdl_rmenvir.fd/ ------------------------------------- Doing diff of separ.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//separ.f 5,6c5 < include 'params.h' < parameter (imx=ima,jmx=jma) --- > 7a7 > parameter (imx=150 , jmx=150) ------------------------------------- Doing diff of snrm2.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//snrm2.f NOTE: File snrm2.f does not exist in directory /nwprod/sorc/gfdl_rmenvir.fd/ ------------------------------------- Doing diff of srotm.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//srotm.f ------------------------------------- Doing diff of srotmg.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//srotmg.f ------------------------------------- Doing diff of sscal.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//sscal.f NOTE: File sscal.f does not exist in directory /nwprod/sorc/gfdl_rmenvir.fd/ ------------------------------------- Doing diff of sswap.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//sswap.f NOTE: File sswap.f does not exist in directory /nwprod/sorc/gfdl_rmenvir.fd/ ------------------------------------- Doing diff of wght.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//wght.f 3,4c3 < include 'params.h' < parameter (npts=isize*ires*ngprd) --- > parameter (isize=25,ires=2,ngprd=6,npts=isize*ires*ngprd) 7,8c6,7 < parameter (inmax = ima, jnmax = jma) < parameter (imaxn = ima, jmaxn = jma) --- > parameter (inmax = 150, jnmax = 150) > parameter (imaxn = 150, jmaxn = 150) ------------------------------------- Doing diff of wnlit.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//wnlit.f ------------------------------------- Doing diff of wnlsm.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//wnlsm.f ------------------------------------- Doing diff of wnnls.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//wnnls.f ------------------------------------- Doing diff of xerabt.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//xerabt.f ------------------------------------- Doing diff of xerctl.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//xerctl.f ------------------------------------- Doing diff of xerprt.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//xerprt.f ------------------------------------- Doing diff of xerror.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//xerror.f ------------------------------------- Doing diff of xerrwv.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//xerrwv.f ------------------------------------- Doing diff of xersav.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//xersav.f ------------------------------------- Doing diff of xgetua.f <<< ------ >>> /nwprod/sorc/gfdl_rmenvir.fd//xgetua.f