/[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.7 by molod, Tue Jul 13 21:26:32 2004 UTC revision 1.10 by molod, Mon Jul 26 18:45:17 2004 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "FIZHI_OPTIONS.h"
 #include "PACKAGES_CONFIG.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       .                  pz,plz,plze,dpres,pkht,pkz,tz,qz,oz,       .                  pz,plz,plze,dpres,pkht,pkz,tz,qz,oz,
7       .                  co2,cfc11,cfc12,cfc22,methane,n2o,emissivity,       .                  co2,cfc11,cfc12,cfc22,methane,n2o,emissivity,
# Line 14  C $Name$ Line 13  C $Name$
13    
14        implicit none        implicit none
15  #ifdef ALLOW_DIAGNOSTICS  #ifdef ALLOW_DIAGNOSTICS
16    #include "SIZE.h"
17    #include "diagnostics_SIZE.h"
18  #include "diagnostics.h"  #include "diagnostics.h"
19  #endif  #endif
20    
# Line 22  c --------------- Line 23  c ---------------
23        integer nymd,nhms,istrip,npcs,bi,bj        integer nymd,nhms,istrip,npcs,bi,bj
24        integer mid_level,low_level        integer mid_level,low_level
25        integer im,jm,lm                integer im,jm,lm        
26        real  ptop                      _RL  ptop              
27        real pz(im,jm),plz(im,jm,lm),plze(im,jm,lm+1)        _RL pz(im,jm),plz(im,jm,lm),plze(im,jm,lm+1)
28        real dpres(im,jm,lm),pkht(im,jm,lm+1),pkz(im,jm,lm)        _RL dpres(im,jm,lm),pkht(im,jm,lm+1),pkz(im,jm,lm)
29        real tz(im,jm,lm),qz(im,jm,lm),oz(im,jm,lm)        _RL tz(im,jm,lm),qz(im,jm,lm),oz(im,jm,lm)
30        real co2,cfc11,cfc12,cfc22,methane(lm),n2o(lm)            _RL co2,cfc11,cfc12,cfc22,methane(lm),n2o(lm)    
31        real emissivity(im,jm,10)        _RL emissivity(im,jm,10)
32        real tgz(im,jm),radlwg(im,jm),st4(im,jm),dst4(im,jm)            _RL tgz(im,jm),radlwg(im,jm),st4(im,jm),dst4(im,jm)    
33        real dtradlw(im,jm,lm),dlwdtg (im,jm,lm)        _RL dtradlw(im,jm,lm),dlwdtg (im,jm,lm)
34        real dtradlwc(im,jm,lm),lwgclr(im,jm)            _RL dtradlwc(im,jm,lm),lwgclr(im,jm)    
35        integer nlwcld,nlwlz        integer nlwcld,nlwlz
36        real cldlw(im,jm,lm),clwmo(im,jm,lm),lwlz(im,jm,lm)        _RL cldlw(im,jm,lm),clwmo(im,jm,lm),lwlz(im,jm,lm)
37        logical lpnt        logical lpnt
38        integer imstturb        integer imstturb
39        real qliqave(im,jm,lm),fccave(im,jm,lm)        _RL qliqave(im,jm,lm),fccave(im,jm,lm)
40        integer landtype(im,jm)        integer landtype(im,jm)
41    
42  c Local Variables  c Local Variables
43  c ---------------  c ---------------
44        integer i,j,l,n,nn        integer i,j,l,n,nn
45    
46        real cldtot (im,jm,lm)        _RL cldtot (im,jm,lm)
47        real cldmxo (im,jm,lm)        _RL cldmxo (im,jm,lm)
48    
49        real pl(istrip,lm)        _RL pl(istrip,lm)
50        real pk(istrip,lm)        _RL pk(istrip,lm)
51        real pke(istrip,lm)        _RL pke(istrip,lm)
52        real ple(istrip,lm+1)        _RL ple(istrip,lm+1)
53    
54        real ADELPL(ISTRIP,lm)        _RL ADELPL(ISTRIP,lm)
55        real dtrad(istrip,lm),dtradc(istrip,lm)        _RL dtrad(istrip,lm),dtradc(istrip,lm)
56        real OZL(ISTRIP,lm),TZL(ISTRIP,lm)        _RL OZL(ISTRIP,lm),TZL(ISTRIP,lm)
57        real SHZL(ISTRIP,lm),CLRO(ISTRIP,lm)        _RL SHZL(ISTRIP,lm),CLRO(ISTRIP,lm)
58        real CLMO(ISTRIP,lm)        _RL CLMO(ISTRIP,lm)
59        real flx(ISTRIP,lm+1),flxclr(ISTRIP,lm+1)        _RL flx(ISTRIP,lm+1),flxclr(ISTRIP,lm+1)
60        real cldlz(istrip,lm)        _RL cldlz(istrip,lm)
61        real dfdts(istrip,lm+1),dtdtg(istrip,lm)        _RL dfdts(istrip,lm+1),dtdtg(istrip,lm)
62    
63        real emiss(istrip,10)        _RL emiss(istrip,10)
64        real taual(istrip,lm,10)        _RL taual(istrip,lm,10)
65        real ssaal(istrip,lm,10)        _RL ssaal(istrip,lm,10)
66        real asyal(istrip,lm,10)        _RL asyal(istrip,lm,10)
67        real cwc(istrip,lm,3)        _RL cwc(istrip,lm,3)
68        real reff(istrip,lm,3)        _RL reff(istrip,lm,3)
69        real tauc(istrip,lm,3)        _RL tauc(istrip,lm,3)
70    
71        real SGMT4(ISTRIP),TSURF(ISTRIP),dsgmt4(ISTRIP)        _RL SGMT4(ISTRIP),TSURF(ISTRIP),dsgmt4(ISTRIP)
72        integer lwi(istrip)        integer lwi(istrip)
73    
74        real getcon,secday,convrt,pcheck        _RL getcon,secday,convrt
75    
76        logical high,  trace, cldwater        logical high,  trace, cldwater
77        data high /.true./        data high /.true./
# Line 549  cfpp$ expand (b10kdis) Line 550  cfpp$ expand (b10kdis)
550  c---- input parameters ------  c---- input parameters ------
551    
552        integer m,n,ndim,np,ict,icb        integer m,n,ndim,np,ict,icb
553        real pl(m,ndim,np+1),ta(m,ndim,np),wa(m,ndim,np),oa(m,ndim,np),        _RL pl(m,ndim,np+1),ta(m,ndim,np),wa(m,ndim,np),oa(m,ndim,np),
554       *     ts(m,ndim)       *     ts(m,ndim)
555        real co2,n2o(np),ch4(np),cfc11,cfc12,cfc22,emiss(m,ndim,10)        _RL co2,n2o(np),ch4(np),cfc11,cfc12,cfc22,emiss(m,ndim,10)
556        real cwc(m,ndim,np,3),taucl(m,ndim,np,3),reff(m,ndim,np,3),        _RL cwc(m,ndim,np,3),taucl(m,ndim,np,3),reff(m,ndim,np,3),
557       *     fcld(m,ndim,np)       *     fcld(m,ndim,np)
558        real taual(m,ndim,np,10),ssaal(m,ndim,np,10),asyal(m,ndim,np,10)        _RL taual(m,ndim,np,10),ssaal(m,ndim,np,10),asyal(m,ndim,np,10)
559        logical cldwater,high,trace        logical cldwater,high,trace
560    
561  c---- output parameters ------  c---- output parameters ------
562    
563        real flx(m,ndim,np+1),flc(m,ndim,np+1),dfdts(m,ndim,np+1),        _RL flx(m,ndim,np+1),flc(m,ndim,np+1),dfdts(m,ndim,np+1),
564       *     st4(m,ndim)       *     st4(m,ndim)
565    
566  c---- static data -----  c---- static data -----
567    
568        real cb(5,10)        _RL cb(5,10)
569        real xkw(9),aw(9),bw(9),pm(9),fkw(6,9),gkw(6,3),xke(9)        _RL xkw(9),aw(9),bw(9),pm(9),fkw(6,9),gkw(6,3),xke(9)
570        real aib(3,10),awb(4,10),aiw(4,10),aww(4,10),aig(4,10),awg(4,10)        _RL aib(3,10),awb(4,10),aiw(4,10),aww(4,10),aig(4,10),awg(4,10)
571        integer ne(9),mw(9)        integer ne(9),mw(9)
572    
573  c---- temporary arrays -----  c---- temporary arrays -----
574    
575        real pa(m,n,np),dt(m,n,np)        _RL pa(m,n,np),dt(m,n,np)
576        real sh2o(m,n,np+1),swpre(m,n,np+1),swtem(m,n,np+1)        _RL sh2o(m,n,np+1),swpre(m,n,np+1),swtem(m,n,np+1)
577        real sco3(m,n,np+1),scopre(m,n,np+1),scotem(m,n,np+1)        _RL sco3(m,n,np+1),scopre(m,n,np+1),scotem(m,n,np+1)
578        real dh2o(m,n,np),dcont(m,n,np),dco2(m,n,np),do3(m,n,np)        _RL dh2o(m,n,np),dcont(m,n,np),dco2(m,n,np),do3(m,n,np)
579        real dn2o(m,n,np),dch4(m,n,np)        _RL dn2o(m,n,np),dch4(m,n,np)
580        real df11(m,n,np),df12(m,n,np),df22(m,n,np)        _RL df11(m,n,np),df12(m,n,np),df22(m,n,np)
581        real th2o(m,n,6),tcon(m,n,3),tco2(m,n,6,2)        _RL th2o(m,n,6),tcon(m,n,3),tco2(m,n,6,2)
582        real tn2o(m,n,4),tch4(m,n,4),tcom(m,n,2)        _RL tn2o(m,n,4),tch4(m,n,4),tcom(m,n,2)
583        real tf11(m,n),tf12(m,n),tf22(m,n)        _RL tf11(m,n),tf12(m,n),tf22(m,n)
584        real h2oexp(m,n,np,6),conexp(m,n,np,3),co2exp(m,n,np,6,2)        _RL h2oexp(m,n,np,6),conexp(m,n,np,3),co2exp(m,n,np,6,2)
585        real n2oexp(m,n,np,4),ch4exp(m,n,np,4),comexp(m,n,np,2)        _RL n2oexp(m,n,np,4),ch4exp(m,n,np,4),comexp(m,n,np,2)
586        real f11exp(m,n,np),  f12exp(m,n,np),  f22exp(m,n,np)        _RL f11exp(m,n,np),  f12exp(m,n,np),  f22exp(m,n,np)
587        real clr(m,n,0:np+1),fclr(m,n)        _RL clr(m,n,0:np+1),fclr(m,n)
588        real blayer(m,n,0:np+1),dlayer(m,n,np+1),dbs(m,n)        _RL blayer(m,n,0:np+1),dlayer(m,n,np+1),dbs(m,n)
589        real clrlw(m,n),clrmd(m,n),clrhi(m,n)        _RL clrlw(m,n),clrmd(m,n),clrhi(m,n)
590        real cwp(m,n,np,3)        _RL cwp(m,n,np,3)
591        real trant(m,n),tranal(m,n),transfc(m,n,np+1),trantcr(m,n,np+1)        _RL trant(m,n),tranal(m,n),transfc(m,n,np+1),trantcr(m,n,np+1)
592        real flxu(m,n,np+1),flxd(m,n,np+1),flcu(m,n,np+1),flcd(m,n,np+1)        _RL flxu(m,n,np+1),flxd(m,n,np+1),flcu(m,n,np+1),flcd(m,n,np+1)
593        real rflx(m,n,np+1),rflc(m,n,np+1)        _RL rflx(m,n,np+1),rflc(m,n,np+1)
594    
595        logical oznbnd,co2bnd,h2otbl,conbnd,n2obnd        logical oznbnd,co2bnd,h2otbl,conbnd,n2obnd
596        logical ch4bnd,combnd,f11bnd,f12bnd,f22bnd,b10bnd        logical ch4bnd,combnd,f11bnd,f12bnd,f22bnd,b10bnd
597    
598        real c1 (nx,nc,nt),c2 (nx,nc,nt),c3 (nx,nc,nt)        _RL c1 (nx,nc,nt),c2 (nx,nc,nt),c3 (nx,nc,nt)
599        real o1 (nx,no,nt),o2 (nx,no,nt),o3 (nx,no,nt)        _RL o1 (nx,no,nt),o2 (nx,no,nt),o3 (nx,no,nt)
600        real h11(nx,nh,nt),h12(nx,nh,nt),h13(nx,nh,nt)        _RL h11(nx,nh,nt),h12(nx,nh,nt),h13(nx,nh,nt)
601        real h21(nx,nh,nt),h22(nx,nh,nt),h23(nx,nh,nt)        _RL h21(nx,nh,nt),h22(nx,nh,nt),h23(nx,nh,nt)
602        real h81(nx,nh,nt),h82(nx,nh,nt),h83(nx,nh,nt)        _RL h81(nx,nh,nt),h82(nx,nh,nt),h83(nx,nh,nt)
603    
604        real dp,xx,p1,dwe,dpe,a1,b1,fk1,a2,b2,fk2        _RL dp,xx,p1,dwe,dpe,a1,b1,fk1,a2,b2,fk2
605        real w1,w2,w3,g1,g2,g3,ww,gg,ff,taux,reff1,reff2        _RL w1,w2,w3,g1,g2,g3,ww,gg,ff,taux,reff1,reff2
606    
607  c-----the following coefficients (equivalent to table 2 of  c-----the following coefficients (equivalent to table 2 of
608  c     chou and suarez, 1995) are for computing spectrally  c     chou and suarez, 1995) are for computing spectrally
# Line 1715  c*************************************** Line 1716  c***************************************
1716    
1717  c---- input parameters -----  c---- input parameters -----
1718    
1719        real pa(m,n,np),dt(m,n,np),sabs0(m,n,np)        _RL pa(m,n,np),dt(m,n,np),sabs0(m,n,np)
1720    
1721  c---- output parameters -----  c---- output parameters -----
1722    
1723        real sabs(m,n,np+1),spre(m,n,np+1),stem(m,n,np+1)        _RL sabs(m,n,np+1),spre(m,n,np+1),stem(m,n,np+1)
1724    
1725  c*********************************************************************  c*********************************************************************
1726          do j=1,n          do j=1,n
# Line 1770  c*************************************** Line 1771  c***************************************
1771    
1772  c---- input parameters ------  c---- input parameters ------
1773    
1774        real dh2o(m,n,np),pa(m,n,np),dt(m,n,np)        _RL dh2o(m,n,np),pa(m,n,np),dt(m,n,np)
1775    
1776  c---- output parameters -----  c---- output parameters -----
1777    
1778        real h2oexp(m,n,np,6)        _RL h2oexp(m,n,np,6)
1779    
1780  c---- static data -----  c---- static data -----
1781    
1782        integer mw(9)        integer mw(9)
1783        real xkw(9),aw(9),bw(9),pm(9)        _RL xkw(9),aw(9),bw(9),pm(9)
1784    
1785  c---- temporary arrays -----  c---- temporary arrays -----
1786    
1787        real xh,xh1        _RL xh,xh1
1788    
1789  c**********************************************************************  c**********************************************************************
1790  c    note that the 3 sub-bands in band 3 use the same set of xkw, aw,  c    note that the 3 sub-bands in band 3 use the same set of xkw, aw,
# Line 1886  c*************************************** Line 1887  c***************************************
1887    
1888  c---- input parameters ------  c---- input parameters ------
1889    
1890        real dcont(m,n,np)        _RL dcont(m,n,np)
1891    
1892  c---- updated parameters -----  c---- updated parameters -----
1893    
1894        real conexp(m,n,np,3)        _RL conexp(m,n,np,3)
1895    
1896  c---- static data -----  c---- static data -----
1897    
1898        real xke(9)        _RL xke(9)
1899    
1900  c**********************************************************************  c**********************************************************************
1901    
# Line 1947  c*************************************** Line 1948  c***************************************
1948    
1949  c---- input parameters -----  c---- input parameters -----
1950    
1951        real dco2(m,n,np),pa(m,n,np),dt(m,n,np)        _RL dco2(m,n,np),pa(m,n,np),dt(m,n,np)
1952    
1953  c---- output parameters -----  c---- output parameters -----
1954    
1955        real co2exp(m,n,np,6,2)        _RL co2exp(m,n,np,6,2)
1956    
1957  c---- temporary arrays -----  c---- temporary arrays -----
1958    
1959        real xc        _RL xc
1960    
1961  c**********************************************************************  c**********************************************************************
1962    
# Line 2049  c*************************************** Line 2050  c***************************************
2050    
2051  c---- input parameters -----  c---- input parameters -----
2052    
2053        real dn2o(m,n,np),pa(m,n,np),dt(m,n,np)        _RL dn2o(m,n,np),pa(m,n,np),dt(m,n,np)
2054    
2055  c---- output parameters -----  c---- output parameters -----
2056    
2057        real n2oexp(m,n,np,4)        _RL n2oexp(m,n,np,4)
2058    
2059  c---- temporary arrays -----  c---- temporary arrays -----
2060    
2061        real xc,xc1,xc2        _RL xc,xc1,xc2
2062    
2063  c**********************************************************************  c**********************************************************************
2064    
# Line 2125  c*************************************** Line 2126  c***************************************
2126    
2127  c---- input parameters -----  c---- input parameters -----
2128    
2129        real dch4(m,n,np),pa(m,n,np),dt(m,n,np)        _RL dch4(m,n,np),pa(m,n,np),dt(m,n,np)
2130    
2131  c---- output parameters -----  c---- output parameters -----
2132    
2133        real ch4exp(m,n,np,4)        _RL ch4exp(m,n,np,4)
2134    
2135  c---- temporary arrays -----  c---- temporary arrays -----
2136    
2137        real xc        _RL xc
2138    
2139  c**********************************************************************  c**********************************************************************
2140    
# Line 2197  c*************************************** Line 2198  c***************************************
2198    
2199  c---- input parameters -----  c---- input parameters -----
2200    
2201        real dcom(m,n,np),dt(m,n,np)        _RL dcom(m,n,np),dt(m,n,np)
2202    
2203  c---- output parameters -----  c---- output parameters -----
2204    
2205        real comexp(m,n,np,2)        _RL comexp(m,n,np,2)
2206    
2207  c---- temporary arrays -----  c---- temporary arrays -----
2208    
2209        real xc,xc1,xc2        _RL xc,xc1,xc2
2210    
2211  c**********************************************************************  c**********************************************************************
2212    
# Line 2274  c*************************************** Line 2275  c***************************************
2275    
2276  c---- input parameters -----  c---- input parameters -----
2277    
2278        real dcfc(m,n,np),dt(m,n,np)        _RL dcfc(m,n,np),dt(m,n,np)
2279    
2280  c---- output parameters -----  c---- output parameters -----
2281    
2282        real cfcexp(m,n,np)        _RL cfcexp(m,n,np)
2283    
2284  c---- static data -----  c---- static data -----
2285    
2286        real a1,b1,fk1,a2,b2,fk2        _RL a1,b1,fk1,a2,b2,fk2
2287    
2288  c---- temporary arrays -----  c---- temporary arrays -----
2289    
2290        real xf        _RL xf
2291    
2292  c**********************************************************************  c**********************************************************************
2293    
# Line 2336  c*************************************** Line 2337  c***************************************
2337    
2338  c---- input parameters -----  c---- input parameters -----
2339    
2340        real dh2o(m,n,np),dcont(m,n,np),dn2o(m,n,np)        _RL dh2o(m,n,np),dcont(m,n,np),dn2o(m,n,np)
2341        real dco2(m,n,np),pa(m,n,np),dt(m,n,np)        _RL dco2(m,n,np),pa(m,n,np),dt(m,n,np)
2342    
2343  c---- output parameters -----  c---- output parameters -----
2344    
2345        real h2oexp(m,n,np,6),conexp(m,n,np,3),co2exp(m,n,np,6,2)        _RL h2oexp(m,n,np,6),conexp(m,n,np,3),co2exp(m,n,np,6,2)
2346       *    ,n2oexp(m,n,np,4)       *    ,n2oexp(m,n,np,4)
2347    
2348  c---- temporary arrays -----  c---- temporary arrays -----
2349    
2350        real xx,xx1,xx2,xx3        _RL xx,xx1,xx2,xx3
2351    
2352  c**********************************************************************  c**********************************************************************
2353    
# Line 2487  c*************************************** Line 2488  c***************************************
2488    
2489  c---- input parameters -----  c---- input parameters -----
2490    
2491        real w1,p1,dwe,dpe        _RL w1,p1,dwe,dpe
2492        real sabs(m,n,np+1),spre(m,n,np+1),stem(m,n,np+1)        _RL sabs(m,n,np+1),spre(m,n,np+1),stem(m,n,np+1)
2493        real coef1(nx,nh,nt),coef2(nx,nh,nt),coef3(nx,nh,nt)        _RL coef1(nx,nh,nt),coef2(nx,nh,nt),coef3(nx,nh,nt)
2494    
2495  c---- update parameter -----  c---- update parameter -----
2496    
2497        real tran(m,n)        _RL tran(m,n)
2498    
2499  c---- temporary variables -----  c---- temporary variables -----
2500    
2501        real x1,x2,x3,we,pe,fw,fp,pa,pb,pc,ax,ba,bb,t1,ca,cb,t2        _RL x1,x2,x3,we,pe,fw,fp,pa,pb,pc,ax,ba,bb,t1,ca,cb,t2
2502        integer iw,ip,nn        integer iw,ip,nn
2503    
2504  c**********************************************************************  c**********************************************************************
# Line 2595  c*************************************** Line 2596  c***************************************
2596    
2597  c---- input parameters ------  c---- input parameters ------
2598    
2599        real conexp(m,n,np,3),h2oexp(m,n,np,6)        _RL conexp(m,n,np,3),h2oexp(m,n,np,6)
2600        integer ne(9)        integer ne(9)
2601        real  fkw(6,9),gkw(6,3)        _RL  fkw(6,9),gkw(6,3)
2602    
2603  c---- updated parameters -----  c---- updated parameters -----
2604    
2605        real th2o(m,n,6),tcon(m,n,3),tran(m,n)        _RL th2o(m,n,6),tcon(m,n,3),tran(m,n)
2606    
2607  c---- temporary arrays -----  c---- temporary arrays -----
2608    
2609        real trnth2o        _RL trnth2o
2610    
2611  c-----tco2 are the six exp factors between levels k1 and k2  c-----tco2 are the six exp factors between levels k1 and k2
2612  c     tran is the updated total transmittance between levels k1 and k2  c     tran is the updated total transmittance between levels k1 and k2
# Line 2729  c*************************************** Line 2730  c***************************************
2730    
2731  c---- input parameters -----  c---- input parameters -----
2732    
2733        real co2exp(m,n,np,6,2)        _RL co2exp(m,n,np,6,2)
2734    
2735  c---- updated parameters -----  c---- updated parameters -----
2736    
2737        real tco2(m,n,6,2),tran(m,n)        _RL tco2(m,n,6,2),tran(m,n)
2738    
2739  c---- temporary arrays -----  c---- temporary arrays -----
2740    
2741        real xc        _RL xc
2742    
2743  c-----tco2 is the 6 exp factors between levels k1 and k2.  c-----tco2 is the 6 exp factors between levels k1 and k2.
2744  c     xc is the total co2 transmittance given by eq. (53).  c     xc is the total co2 transmittance given by eq. (53).
# Line 2818  c*************************************** Line 2819  c***************************************
2819    
2820  c---- input parameters -----  c---- input parameters -----
2821    
2822        real n2oexp(m,n,np,4)        _RL n2oexp(m,n,np,4)
2823    
2824  c---- updated parameters -----  c---- updated parameters -----
2825    
2826        real tn2o(m,n,4),tran(m,n)        _RL tn2o(m,n,4),tran(m,n)
2827    
2828  c---- temporary arrays -----  c---- temporary arrays -----
2829    
2830        real xc        _RL xc
2831    
2832  c-----tn2o is the 2 exp factors between levels k1 and k2.  c-----tn2o is the 2 exp factors between levels k1 and k2.
2833  c     xc is the total n2o transmittance  c     xc is the total n2o transmittance
# Line 2895  c*************************************** Line 2896  c***************************************
2896    
2897  c---- input parameters -----  c---- input parameters -----
2898    
2899        real ch4exp(m,n,np,4)        _RL ch4exp(m,n,np,4)
2900    
2901  c---- updated parameters -----  c---- updated parameters -----
2902    
2903        real tch4(m,n,4),tran(m,n)        _RL tch4(m,n,4),tran(m,n)
2904    
2905  c---- temporary arrays -----  c---- temporary arrays -----
2906    
2907        real xc        _RL xc
2908    
2909  c-----tch4 is the 2 exp factors between levels k1 and k2.  c-----tch4 is the 2 exp factors between levels k1 and k2.
2910  c     xc is the total ch4 transmittance  c     xc is the total ch4 transmittance
# Line 2969  c*************************************** Line 2970  c***************************************
2970    
2971  c---- input parameters -----  c---- input parameters -----
2972    
2973        real comexp(m,n,np,2)        _RL comexp(m,n,np,2)
2974    
2975  c---- updated parameters -----  c---- updated parameters -----
2976    
2977        real tcom(m,n,2),tran(m,n)        _RL tcom(m,n,2),tran(m,n)
2978    
2979  c---- temporary arrays -----  c---- temporary arrays -----
2980    
2981        real xc        _RL xc
2982    
2983  c-----tcom is the 2 exp factors between levels k1 and k2.  c-----tcom is the 2 exp factors between levels k1 and k2.
2984  c     xc is the total co2-minor transmittance  c     xc is the total co2-minor transmittance
# Line 3037  c*************************************** Line 3038  c***************************************
3038    
3039  c---- input parameters -----  c---- input parameters -----
3040    
3041        real cfcexp(m,n,np)        _RL cfcexp(m,n,np)
3042    
3043  c---- updated parameters -----  c---- updated parameters -----
3044    
3045        real tcfc(m,n),tran(m,n)        _RL tcfc(m,n),tran(m,n)
3046    
3047  c-----tcfc is the exp factors between levels k1 and k2.  c-----tcfc is the exp factors between levels k1 and k2.
3048    
# Line 3093  c*************************************** Line 3094  c***************************************
3094    
3095  c---- input parameters -----  c---- input parameters -----
3096    
3097        real h2oexp(m,n,np,6),conexp(m,n,np,3),co2exp(m,n,np,6,2)        _RL h2oexp(m,n,np,6),conexp(m,n,np,3),co2exp(m,n,np,6,2)
3098       *    ,n2oexp(m,n,np,4)       *    ,n2oexp(m,n,np,4)
3099    
3100  c---- updated parameters -----  c---- updated parameters -----
3101    
3102        real th2o(m,n,6),tcon(m,n,3),tco2(m,n,6,2),tn2o(m,n,4)        _RL th2o(m,n,6),tcon(m,n,3),tco2(m,n,6,2),tn2o(m,n,4)
3103       *    ,tran(m,n)       *    ,tran(m,n)
3104    
3105  c---- temporary arrays -----  c---- temporary arrays -----
3106    
3107        real xx        _RL xx
3108    
3109  c-----initialize tran  c-----initialize tran
3110    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.22