221 |
do L = 1,lm |
do L = 1,lm |
222 |
do i = 1,istrip |
do i = 1,istrip |
223 |
dtrad(i,L) = ( flx(i,L)- flx(i,L+1))*adelpl(i,L) |
dtrad(i,L) = ( flx(i,L)- flx(i,L+1))*adelpl(i,L) |
224 |
tmpstrip(i,L) = ( flx(i,L)- flx(i,L+1)) |
tmpstrip(i,L) = flx(i,L) |
225 |
dtdtg(i,L) = ( dfdts(i,L)- dfdts(i,L+1))*adelpl(i,L) |
dtdtg(i,L) = ( dfdts(i,L)- dfdts(i,L+1))*adelpl(i,L) |
226 |
dtradc(i,L) = (flxclr(i,L)-flxclr(i,L+1))*adelpl(i,L) |
dtradc(i,L) = (flxclr(i,L)-flxclr(i,L+1))*adelpl(i,L) |
227 |
enddo |
enddo |
815 |
c save c1,c2,c3,o1,o2,o3 |
c save c1,c2,c3,o1,o2,o3 |
816 |
c save h11,h12,h13,h21,h22,h23,h81,h82,h83 |
c save h11,h12,h13,h21,h22,h23,h81,h82,h83 |
817 |
|
|
818 |
c if (first) then |
if (first) then |
819 |
|
|
820 |
c-----tables co2 and h2o are only used with 'high' option |
c-----tables co2 and h2o are only used with 'high' option |
821 |
|
|
883 |
enddo |
enddo |
884 |
enddo |
enddo |
885 |
|
|
886 |
c first=.false. |
first=.false. |
887 |
|
|
888 |
c endif |
endif |
889 |
|
|
890 |
c-----set the pressure at the top of the model atmosphere |
c-----set the pressure at the top of the model atmosphere |
891 |
c to 1.0e-4 if it is zero |
c to 1.0e-4 if it is zero |
1493 |
dpe=0.2 |
dpe=0.2 |
1494 |
call tablup(k1,k2,m,n,np,nx,no,nt,sco3,scopre,scotem, |
call tablup(k1,k2,m,n,np,nx,no,nt,sco3,scopre,scotem, |
1495 |
* w1,p1,dwe,dpe,o1,o2,o3,trant) |
* w1,p1,dwe,dpe,o1,o2,o3,trant) |
1496 |
|
|
1497 |
endif |
endif |
1498 |
|
|
1499 |
c***** for trace gases ***** |
c***** for trace gases ***** |