/[MITgcm]/MITgcm/pkg/fizhi/fizhi_lwrad.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_lwrad.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.18 by molod, Thu Aug 26 14:37:33 2004 UTC revision 1.19 by molod, Sun Aug 29 19:46:19 2004 UTC
# Line 221  C ************************************** Line 221  C **************************************
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
# Line 815  c     o3 (band 5), and h2o (bands 1, 2, Line 815  c     o3 (band 5), and h2o (bands 1, 2,
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    
# Line 883  c-----always use table look-up for ozone Line 883  c-----always use table look-up for ozone
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
# Line 1493  c-----All use table look-up to compute o Line 1493  c-----All use table look-up to compute o
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 *****
Legend:
Removed from v.1.18  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.22