/[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.10 by molod, Mon Jul 26 18:45:17 2004 UTC revision 1.14 by molod, Fri Jul 30 18:59:32 2004 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "FIZHI_OPTIONS.h"  #include "FIZHI_OPTIONS.h"
5        subroutine lwrio (nymd,nhms,bi,bj,istrip,npcs,low_level,mid_level,        subroutine lwrio (nymd,nhms,bi,bj,istrip,npcs,low_level,mid_level,
6         .                  im,jm,lm,
7       .                  pz,plz,plze,dpres,pkht,pkz,tz,qz,oz,       .                  pz,plz,plze,dpres,pkht,pkz,tz,qz,oz,
8       .                  co2,cfc11,cfc12,cfc22,methane,n2o,emissivity,       .                  co2,cfc11,cfc12,cfc22,methane,n2o,emissivity,
9       .                  tgz,radlwg,st4,dst4,       .                  tgz,radlwg,st4,dst4,
10       .                  dtradlw,dlwdtg,dtradlwc,lwgclr,       .                  dtradlw,dlwdtg,dtradlwc,lwgclr,
11       .                  im,jm,lm,ptop,       .                  ptop,nlwcld,cldlw,clwmo,nlwlz,lwlz,
      .                  nlwcld,cldlw,clwmo,nlwlz,lwlz,  
12       .                  lpnt,imstturb,qliqave,fccave,landtype)       .                  lpnt,imstturb,qliqave,fccave,landtype)
13    
14        implicit none        implicit none
# Line 1223  c-----initialize fluxes Line 1223  c-----initialize fluxes
1223         enddo         enddo
1224        enddo        enddo
1225    
   
1226        do 2000 k1=1,np        do 2000 k1=1,np
1227    
1228  c-----initialize fclr, th2o, tcon, tco2, and tranal  c-----initialize fclr, th2o, tcon, tco2, and tranal
# Line 1404  c-----compute water vapor transmittance Line 1403  c-----compute water vapor transmittance
1403            if (ib.eq.1) then            if (ib.eq.1) then
1404             call tablup(k1,k2,m,n,np,nx,nh,nt,sh2o,swpre,swtem,             call tablup(k1,k2,m,n,np,nx,nh,nt,sh2o,swpre,swtem,
1405       *                 w1,p1,dwe,dpe,h11,h12,h13,trant)       *                 w1,p1,dwe,dpe,h11,h12,h13,trant)
   
1406            endif            endif
1407            if (ib.eq.2) then            if (ib.eq.2) then
1408             call tablup(k1,k2,m,n,np,nx,nh,nt,sh2o,swpre,swtem,             call tablup(k1,k2,m,n,np,nx,nh,nt,sh2o,swpre,swtem,
# Line 1439  c-----compute co2 transmittance using ta Line 1437  c-----compute co2 transmittance using ta
1437            dpe=0.2            dpe=0.2
1438            call tablup(k1,k2,m,n,np,nx,nc,nt,sco3,scopre,scotem,            call tablup(k1,k2,m,n,np,nx,nc,nt,sco3,scopre,scotem,
1439       *                w1,p1,dwe,dpe,c1,c2,c3,trant)       *                w1,p1,dwe,dpe,c1,c2,c3,trant)
1440    
1441         else         else
1442    
1443  c-----compute co2 transmittance using k-distribution method  c-----compute co2 transmittance using k-distribution method
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22