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

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

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

revision 1.24 by ce107, Thu Jun 16 16:46:12 2005 UTC revision 1.25 by ce107, Fri Jun 17 16:51:24 2005 UTC
# Line 169  c -------------------------------------- Line 169  c --------------------------------------
169          do L =1,lm          do L =1,lm
170          do j =1,jm          do j =1,jm
171          do i =1,im          do i =1,im
172           cldtot(i,j,L)=min(1.0 _d 0,max(cldsw(i,j,L),fccave(i,j,L)/imstturb))             cldtot(i,j,L)=min(1.0 _d 0,max(cldsw(i,j,L),fccave(i,j,L)/
173           cldmxo(i,j,L)=min(1.0 _d 0,cswmo(i,j,L))       $          imstturb))
174               cldmxo(i,j,L)=min(1.0 _d 0,cswmo(i,j,L))
175             swlz(i,j,L)=swlz(i,j,L)+qliqave(i,j,L)/imstturb             swlz(i,j,L)=swlz(i,j,L)+qliqave(i,j,L)/imstturb
176          enddo          enddo
177          enddo          enddo
# Line 586  c ---------------------------- Line 587  c ----------------------------
587        do L=1,lm        do L=1,lm
588        do j=1,jm        do j=1,jm
589        do i=1,im        do i=1,im
590                   alf   =  min( max(( tl(i,j,L)-233.15)/20.,0. _d 0) ,1. _d 0)           alf   =  min( max(( tl(i,j,L)-233.15)/20.,0. _d 0), 1. _d 0)
591                    dp   =  ple(i,j,L+1)-ple(i,j,L)           dp   =  ple(i,j,L+1)-ple(i,j,L)
592           tau(i,j,L,1)  = 0.0           tau(i,j,L,1)  = 0.0
593           tau(i,j,L,2)  = 0.0           tau(i,j,L,2)  = 0.0
594           tau(i,j,L,3)  = 0.0           tau(i,j,L,3)  = 0.0
# Line 606  c -------------------------------------- Line 607  c --------------------------------------
607    
608  c Large-Scale Ice  c Large-Scale Ice
609  c ---------------  c ---------------
610           tauice = max( 0.0002 _d 0, 0.002*min(500*lz(i,j,L)*1000,1.0 _d 0) )           tauice = max( 0.0002 _d 0, 0.002*min(500*lz(i,j,L)*1000,
611         $        1.0 _d 0) )
612           tau(i,j,L,1) = fracls*(1-alf)*tauice*dp           tau(i,j,L,1) = fracls*(1-alf)*tauice*dp
613    
614  c Large-Scale Water  c Large-Scale Water
615  c -----------------  c -----------------
616  C Over Land  C Over Land
617           if( lwi(i,j).le.10 ) then           if( lwi(i,j).le.10 ) then
618            tauh2o = max( 0.0020 _d 0, 0.200*min(200*lz(i,j,L)*1000,1.0 _d 0) )                tauh2o = max( 0.0020 _d 0, 0.200*min(200*lz(i,j,L)*1000,
619         $           1.0 _d 0) )  
620            tau(i,j,L,3) = fracls*alf*tauh2o*dp            tau(i,j,L,3) = fracls*alf*tauh2o*dp
621           else           else
622  C Non-Precipitation Clouds Over Ocean  C Non-Precipitation Clouds Over Ocean
# Line 622  C Non-Precipitation Clouds Over Ocean Line 625  C Non-Precipitation Clouds Over Ocean
625             tau(i,j,L,2) = fracls*alf*tauh2o*dp             tau(i,j,L,2) = fracls*alf*tauh2o*dp
626            else            else
627  C Over Ocean  C Over Ocean
628             tauh2o = max( 0.0020 _d 0, 0.120*min( 20*lz(i,j,L)*1000,1.0 _d 0) )                 tauh2o = max( 0.0020 _d 0, 0.120*min( 20*lz(i,j,L)*1000,
629         $            1.0 _d 0) )
630             tau(i,j,L,3) = fracls*alf*tauh2o*dp             tau(i,j,L,3) = fracls*alf*tauh2o*dp
631            endif            endif
632           endif           endif
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22