/[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.6 by molod, Tue Jul 13 21:18:41 2004 UTC revision 1.11 by molod, Tue Jul 27 04:10:56 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 770  c     o3 (band 5), and h2o (bands 1, 2, Line 771  c     o3 (band 5), and h2o (bands 1, 2,
771        logical first        logical first
772        data first /.true./        data first /.true./
773    
774  #include "h2o.tran3"  #include "h2o-tran3.h"
775  #include "co2.tran3"  #include "co2-tran3.h"
776  #include "o3.tran3"  #include "o3-tran3.h"
777    
778        save c1,c2,c3,o1,o2,o3        save c1,c2,c3,o1,o2,o3
779        save h11,h12,h13,h21,h22,h23,h81,h82,h83        save h11,h12,h13,h21,h22,h23,h81,h82,h83
# Line 1222  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 1403  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 1438  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
# Line 1715  c*************************************** Line 1715  c***************************************
1715    
1716  c---- input parameters -----  c---- input parameters -----
1717    
1718        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)
1719    
1720  c---- output parameters -----  c---- output parameters -----
1721    
1722        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)
1723    
1724  c*********************************************************************  c*********************************************************************
1725          do j=1,n          do j=1,n
# Line 1770  c*************************************** Line 1770  c***************************************
1770    
1771  c---- input parameters ------  c---- input parameters ------
1772    
1773        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)
1774    
1775  c---- output parameters -----  c---- output parameters -----
1776    
1777        real h2oexp(m,n,np,6)        _RL h2oexp(m,n,np,6)
1778    
1779  c---- static data -----  c---- static data -----
1780    
1781        integer mw(9)        integer mw(9)
1782        real xkw(9),aw(9),bw(9),pm(9)        _RL xkw(9),aw(9),bw(9),pm(9)
1783    
1784  c---- temporary arrays -----  c---- temporary arrays -----
1785    
1786        real xh,xh1        _RL xh,xh1
1787    
1788  c**********************************************************************  c**********************************************************************
1789  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 1886  c***************************************
1886    
1887  c---- input parameters ------  c---- input parameters ------
1888    
1889        real dcont(m,n,np)        _RL dcont(m,n,np)
1890    
1891  c---- updated parameters -----  c---- updated parameters -----
1892    
1893        real conexp(m,n,np,3)        _RL conexp(m,n,np,3)
1894    
1895  c---- static data -----  c---- static data -----
1896    
1897        real xke(9)        _RL xke(9)
1898    
1899  c**********************************************************************  c**********************************************************************
1900    
# Line 1947  c*************************************** Line 1947  c***************************************
1947    
1948  c---- input parameters -----  c---- input parameters -----
1949    
1950        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)
1951    
1952  c---- output parameters -----  c---- output parameters -----
1953    
1954        real co2exp(m,n,np,6,2)        _RL co2exp(m,n,np,6,2)
1955    
1956  c---- temporary arrays -----  c---- temporary arrays -----
1957    
1958        real xc        _RL xc
1959    
1960  c**********************************************************************  c**********************************************************************
1961    
# Line 2049  c*************************************** Line 2049  c***************************************
2049    
2050  c---- input parameters -----  c---- input parameters -----
2051    
2052        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)
2053    
2054  c---- output parameters -----  c---- output parameters -----
2055    
2056        real n2oexp(m,n,np,4)        _RL n2oexp(m,n,np,4)
2057    
2058  c---- temporary arrays -----  c---- temporary arrays -----
2059    
2060        real xc,xc1,xc2        _RL xc,xc1,xc2
2061    
2062  c**********************************************************************  c**********************************************************************
2063    
# Line 2125  c*************************************** Line 2125  c***************************************
2125    
2126  c---- input parameters -----  c---- input parameters -----
2127    
2128        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)
2129    
2130  c---- output parameters -----  c---- output parameters -----
2131    
2132        real ch4exp(m,n,np,4)        _RL ch4exp(m,n,np,4)
2133    
2134  c---- temporary arrays -----  c---- temporary arrays -----
2135    
2136        real xc        _RL xc
2137    
2138  c**********************************************************************  c**********************************************************************
2139    
# Line 2197  c*************************************** Line 2197  c***************************************
2197    
2198  c---- input parameters -----  c---- input parameters -----
2199    
2200        real dcom(m,n,np),dt(m,n,np)        _RL dcom(m,n,np),dt(m,n,np)
2201    
2202  c---- output parameters -----  c---- output parameters -----
2203    
2204        real comexp(m,n,np,2)        _RL comexp(m,n,np,2)
2205    
2206  c---- temporary arrays -----  c---- temporary arrays -----
2207    
2208        real xc,xc1,xc2        _RL xc,xc1,xc2
2209    
2210  c**********************************************************************  c**********************************************************************
2211    
# Line 2274  c*************************************** Line 2274  c***************************************
2274    
2275  c---- input parameters -----  c---- input parameters -----
2276    
2277        real dcfc(m,n,np),dt(m,n,np)        _RL dcfc(m,n,np),dt(m,n,np)
2278    
2279  c---- output parameters -----  c---- output parameters -----
2280    
2281        real cfcexp(m,n,np)        _RL cfcexp(m,n,np)
2282    
2283  c---- static data -----  c---- static data -----
2284    
2285        real a1,b1,fk1,a2,b2,fk2        _RL a1,b1,fk1,a2,b2,fk2
2286    
2287  c---- temporary arrays -----  c---- temporary arrays -----
2288    
2289        real xf        _RL xf
2290    
2291  c**********************************************************************  c**********************************************************************
2292    
# Line 2336  c*************************************** Line 2336  c***************************************
2336    
2337  c---- input parameters -----  c---- input parameters -----
2338    
2339        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)
2340        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)
2341    
2342  c---- output parameters -----  c---- output parameters -----
2343    
2344        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)
2345       *    ,n2oexp(m,n,np,4)       *    ,n2oexp(m,n,np,4)
2346    
2347  c---- temporary arrays -----  c---- temporary arrays -----
2348    
2349        real xx,xx1,xx2,xx3        _RL xx,xx1,xx2,xx3
2350    
2351  c**********************************************************************  c**********************************************************************
2352    
# Line 2487  c*************************************** Line 2487  c***************************************
2487    
2488  c---- input parameters -----  c---- input parameters -----
2489    
2490        real w1,p1,dwe,dpe        _RL w1,p1,dwe,dpe
2491        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)
2492        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)
2493    
2494  c---- update parameter -----  c---- update parameter -----
2495    
2496        real tran(m,n)        _RL tran(m,n)
2497    
2498  c---- temporary variables -----  c---- temporary variables -----
2499    
2500        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
2501        integer iw,ip,nn        integer iw,ip,nn
2502    
2503  c**********************************************************************  c**********************************************************************
# Line 2595  c*************************************** Line 2595  c***************************************
2595    
2596  c---- input parameters ------  c---- input parameters ------
2597    
2598        real conexp(m,n,np,3),h2oexp(m,n,np,6)        _RL conexp(m,n,np,3),h2oexp(m,n,np,6)
2599        integer ne(9)        integer ne(9)
2600        real  fkw(6,9),gkw(6,3)        _RL  fkw(6,9),gkw(6,3)
2601    
2602  c---- updated parameters -----  c---- updated parameters -----
2603    
2604        real th2o(m,n,6),tcon(m,n,3),tran(m,n)        _RL th2o(m,n,6),tcon(m,n,3),tran(m,n)
2605    
2606  c---- temporary arrays -----  c---- temporary arrays -----
2607    
2608        real trnth2o        _RL trnth2o
2609    
2610  c-----tco2 are the six exp factors between levels k1 and k2  c-----tco2 are the six exp factors between levels k1 and k2
2611  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 2729  c***************************************
2729    
2730  c---- input parameters -----  c---- input parameters -----
2731    
2732        real co2exp(m,n,np,6,2)        _RL co2exp(m,n,np,6,2)
2733    
2734  c---- updated parameters -----  c---- updated parameters -----
2735    
2736        real tco2(m,n,6,2),tran(m,n)        _RL tco2(m,n,6,2),tran(m,n)
2737    
2738  c---- temporary arrays -----  c---- temporary arrays -----
2739    
2740        real xc        _RL xc
2741    
2742  c-----tco2 is the 6 exp factors between levels k1 and k2.  c-----tco2 is the 6 exp factors between levels k1 and k2.
2743  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 2818  c***************************************
2818    
2819  c---- input parameters -----  c---- input parameters -----
2820    
2821        real n2oexp(m,n,np,4)        _RL n2oexp(m,n,np,4)
2822    
2823  c---- updated parameters -----  c---- updated parameters -----
2824    
2825        real tn2o(m,n,4),tran(m,n)        _RL tn2o(m,n,4),tran(m,n)
2826    
2827  c---- temporary arrays -----  c---- temporary arrays -----
2828    
2829        real xc        _RL xc
2830    
2831  c-----tn2o is the 2 exp factors between levels k1 and k2.  c-----tn2o is the 2 exp factors between levels k1 and k2.
2832  c     xc is the total n2o transmittance  c     xc is the total n2o transmittance
# Line 2895  c*************************************** Line 2895  c***************************************
2895    
2896  c---- input parameters -----  c---- input parameters -----
2897    
2898        real ch4exp(m,n,np,4)        _RL ch4exp(m,n,np,4)
2899    
2900  c---- updated parameters -----  c---- updated parameters -----
2901    
2902        real tch4(m,n,4),tran(m,n)        _RL tch4(m,n,4),tran(m,n)
2903    
2904  c---- temporary arrays -----  c---- temporary arrays -----
2905    
2906        real xc        _RL xc
2907    
2908  c-----tch4 is the 2 exp factors between levels k1 and k2.  c-----tch4 is the 2 exp factors between levels k1 and k2.
2909  c     xc is the total ch4 transmittance  c     xc is the total ch4 transmittance
# Line 2969  c*************************************** Line 2969  c***************************************
2969    
2970  c---- input parameters -----  c---- input parameters -----
2971    
2972        real comexp(m,n,np,2)        _RL comexp(m,n,np,2)
2973    
2974  c---- updated parameters -----  c---- updated parameters -----
2975    
2976        real tcom(m,n,2),tran(m,n)        _RL tcom(m,n,2),tran(m,n)
2977    
2978  c---- temporary arrays -----  c---- temporary arrays -----
2979    
2980        real xc        _RL xc
2981    
2982  c-----tcom is the 2 exp factors between levels k1 and k2.  c-----tcom is the 2 exp factors between levels k1 and k2.
2983  c     xc is the total co2-minor transmittance  c     xc is the total co2-minor transmittance
# Line 3037  c*************************************** Line 3037  c***************************************
3037    
3038  c---- input parameters -----  c---- input parameters -----
3039    
3040        real cfcexp(m,n,np)        _RL cfcexp(m,n,np)
3041    
3042  c---- updated parameters -----  c---- updated parameters -----
3043    
3044        real tcfc(m,n),tran(m,n)        _RL tcfc(m,n),tran(m,n)
3045    
3046  c-----tcfc is the exp factors between levels k1 and k2.  c-----tcfc is the exp factors between levels k1 and k2.
3047    
# Line 3093  c*************************************** Line 3093  c***************************************
3093    
3094  c---- input parameters -----  c---- input parameters -----
3095    
3096        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)
3097       *    ,n2oexp(m,n,np,4)       *    ,n2oexp(m,n,np,4)
3098    
3099  c---- updated parameters -----  c---- updated parameters -----
3100    
3101        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)
3102       *    ,tran(m,n)       *    ,tran(m,n)
3103    
3104  c---- temporary arrays -----  c---- temporary arrays -----
3105    
3106        real xx        _RL xx
3107    
3108  c-----initialize tran  c-----initialize tran
3109    
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22