dimension dt(20),de(20),rqh(20,20),rqe(20,20) open (unit=22,file='gra.res',status='unknown') do 1 i=1,20 dt(i)=-5.0+float(i-1) 1 de(i)=-2.+float(i-1) c dt(i)=-3. c de(i)=-0.1 ta=10. c ta(i)=-5.+float(i-1)/4. c ZH=(0.377*13./1000.+1)*(2.87E2*(TA(i)+273.15)) c RO=1000.*100./ZH ro=1.24 at=0.144 ae=0.159 vis=(9.E-4*ta+0.132)/10000. alfa=1./(ta+273.15) beta=0.61 rkt=2.1389E-5 rke=2.3741E-5 rm=0.61*1004.8*(ta+273.15)/2468000. c************ do 22 j=1,20 do 22 i=1,20 if(dt(i).lt.0.0001.and.dt(i).gt.-0.0001) dt(i)=dt(i)+0.01 if(de(j).lt.0.0001.and.de(j).gt.-0.0001) de(j)=de(j)+0.01 rbo=0.662*(dt(i))/(de(j)) boh=1+rm/rbo boe=1+rbo/rm C******************************************** C f l u x e s * C******************************************** dt4=dt(i)*dt(i)*dt(i)*dt(i) dth=dt4**0.33333333 bb=(alfa*rkt*rkt*9.81*boh/vis) bb13=abs(bb)**0.3333333 if(bb.lt.0.) bb13=bb13*(-1.) rqh(i,j)=ro*at*1004.8*(dth)*bb13 dq=de(j) dq1=dq/1000. dq4=dq1*dq1*dq1*dq1 dqh=dq4**0.33333333 be=(beta*rke*rke*9.81*boe/vis) be13=abs(be)**0.3333333 if(be.lt.0.) be13=be13*(-1.) rqe(i,j)=ro*ae*2468000.*(dqh)*be13 c*********** 22 continue write (22,*) 'rqe(i,j),i=1,20' do 2 j=1,20 2 write (22,3) (rqe(i,j),i=1,20) 3 format (20f10.2) stop end