/[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.3 by molod, Thu Jun 24 19:57:02 2004 UTC revision 1.4 by molod, Wed Jul 7 19:33:48 2004 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5        subroutine lwrio (nymd,nhms,bi,bj,istrip,npcs,        subroutine lwrio (nymd,nhms,bi,bj,istrip,npcs,
6         .                  low_level,mid_level,
7       .                  pz,tz,qz,plz,plze,pkz,pkht,oz,co2,       .                  pz,tz,qz,plz,plze,pkz,pkht,oz,co2,
8       .                  cfc11,cfc12,cfc22,       .                  cfc11,cfc12,cfc22,
9       .                  methane,n2o,emissivity,       .                  methane,n2o,emissivity,
# Line 20  C $Name$ Line 21  C $Name$
21  c Input Variables  c Input Variables
22  c ---------------  c ---------------
23        integer nymd,nhms,istrip,npcs,bi,bj        integer nymd,nhms,istrip,npcs,bi,bj
24          integer mid_level,low_level
25    
26        integer im,jm,lm                integer im,jm,lm        
27        real  ptop                      real  ptop              
# Line 60  c --------------- Line 62  c ---------------
62  c Local Variables  c Local Variables
63  c ---------------  c ---------------
64        integer i,j,l,n,nn        integer i,j,l,n,nn
       integer mid_level,low_level  
65    
66        real          PLZ(im,jm,lm)        real          PLZ(im,jm,lm)
67        real          PKZ(im,jm,lm)        real          PKZ(im,jm,lm)
# Line 113  C ************************************** Line 114  C **************************************
114        SECDAY = GETCON('SDAY')        SECDAY = GETCON('SDAY')
115        CONVRT = GETCON('GRAVITY') / ( 100.0 * GETCON('CP') )        CONVRT = GETCON('GRAVITY') / ( 100.0 * GETCON('CP') )
116    
 c Determine Level Indices for Low-Mid-High Cloud Regions  
 c ------------------------------------------------------  
       low_level = lm  
       mid_level = lm  
       do L = lm-1,1,-1  
       pcheck = plz(  
   
       if (pcheck.gt.700.0) low_level = L  
       if (pcheck.gt.400.0) mid_level = L  
       enddo  
   
117  c Adjust cloud fractions and cloud liquid water due to moist turbulence  c Adjust cloud fractions and cloud liquid water due to moist turbulence
118  c ---------------------------------------------------------------------  c ---------------------------------------------------------------------
119        if(imstturb.ne.0) then        if(imstturb.ne.0) then
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22