/[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.11 by molod, Fri Jul 16 20:08:08 2004 UTC revision 1.12 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 swrio (nymd,nhms,bi,bj,ndswr,myid,istrip,npcs,        subroutine swrio (nymd,nhms,bi,bj,ndswr,myid,istrip,npcs,
6       .        low_level,mid_level,       .        low_level,mid_level,
7       .        pz,plz,plze,dpres,pkht,pkz,tz,qz,oz,co2,       .        pz,plz,plze,dpres,pkht,pkz,tz,qz,oz,co2,
# Line 25  c --------------- Line 24  c ---------------
24        integer nymd,nhms,bi,bj,ndswr,myid,istrip,npcs        integer nymd,nhms,bi,bj,ndswr,myid,istrip,npcs
25        integer mid_level,low_level        integer mid_level,low_level
26        integer im,jm,lm                integer im,jm,lm        
27        real  ptop        _RL  ptop
28        real pz(im,jm),plz(im,jm,lm),plze(im,jm,lm+1),dpres(im,jm,lm)        _RL pz(im,jm),plz(im,jm,lm),plze(im,jm,lm+1),dpres(im,jm,lm)
29        real pkht(im,jm,lm+1),pkz(im,jm,lm)        _RL pkht(im,jm,lm+1),pkz(im,jm,lm)
30        real tz(im,jm,lm),qz(im,jm,lm)        _RL tz(im,jm,lm),qz(im,jm,lm)
31        real oz(im,jm,lm)        _RL oz(im,jm,lm)
32        real co2        _RL co2
33        real albvisdr(im,jm),albvisdf(im,jm),albnirdr(im,jm)        _RL albvisdr(im,jm),albvisdf(im,jm),albnirdr(im,jm)
34        real albnirdf(im,jm)        _RL albnirdf(im,jm)
35        real radswg(im,jm),swgclr(im,jm),fdifpar(im,jm),fdirpar(im,jm)        _RL radswg(im,jm),swgclr(im,jm),fdifpar(im,jm),fdirpar(im,jm)
36        real osr(im,jm),osrclr(im,jm),dtradsw(im,jm,lm),dtswclr(im,jm,lm)        _RL osr(im,jm),osrclr(im,jm),dtradsw(im,jm,lm),dtswclr(im,jm,lm)
37        integer nswcld,nswlz            integer nswcld,nswlz    
38        real cldsw(im,jm,lm),cswmo(im,jm,lm),swlz(im,jm,lm)          _RL cldsw(im,jm,lm),cswmo(im,jm,lm),swlz(im,jm,lm)  
39        logical lpnt                    logical lpnt            
40        integer imstturb                integer imstturb        
41        real qliqave(im,jm,lm),fccave(im,jm,lm)          _RL qliqave(im,jm,lm),fccave(im,jm,lm)  
42        integer landtype(im,jm)        integer landtype(im,jm)
43        real xlats(im,jm),xlons(im,jm)        _RL xlats(im,jm),xlons(im,jm)
44    
45  c Local Variables  c Local Variables
46  c ---------------  c ---------------
47        integer   i,j,L,nn,nsecf        integer   i,j,L,nn,nsecf
48        integer   ntmstp,nymd2,nhms2        integer   ntmstp,nymd2,nhms2
49        real      getcon,grav,cp,undef        _RL      getcon,grav,cp,undef
50        real      ra,alf,reffw,reffi,tminv        _RL      ra,alf,reffw,reffi,tminv
51    
52        parameter ( reffw = 10.0 )          parameter ( reffw = 10.0 )  
53        parameter ( reffi = 65.0 )          parameter ( reffi = 65.0 )  
54    
55        real tdry(im,jm,lm)        _RL tdry(im,jm,lm)
56        real TEMP1(im,jm)        _RL TEMP1(im,jm)
57        real TEMP2(im,jm)        _RL TEMP2(im,jm)
58        real zenith (im,jm)        _RL zenith (im,jm)
59        real cldtot (im,jm,lm)        _RL cldtot (im,jm,lm)
60        real cldmxo (im,jm,lm)        _RL cldmxo (im,jm,lm)
61        real totcld (im,jm)        _RL totcld (im,jm)
62        real cldlow (im,jm)        _RL cldlow (im,jm)
63        real cldmid (im,jm)        _RL cldmid (im,jm)
64        real cldhi  (im,jm)        _RL cldhi  (im,jm)
65        real taulow (im,jm)        _RL taulow (im,jm)
66        real taumid (im,jm)        _RL taumid (im,jm)
67        real tauhi  (im,jm)        _RL tauhi  (im,jm)
68        real tautype(im,jm,lm,3)        _RL tautype(im,jm,lm,3)
69        real tau(im,jm,lm)        _RL tau(im,jm,lm)
70        real albedo(im,jm)            _RL albedo(im,jm)    
71    
72        real PK(ISTRIP,lm)        _RL PK(ISTRIP,lm)
73        real qzl(ISTRIP,lm),CLRO(ISTRIP,lm)        _RL qzl(ISTRIP,lm),CLRO(ISTRIP,lm)
74        real TZL(ISTRIP,lm)        _RL TZL(ISTRIP,lm)
75        real OZL(ISTRIP,lm)        _RL OZL(ISTRIP,lm)
76        real PLE(ISTRIP,lm+1)        _RL PLE(ISTRIP,lm+1)
77        real COSZ(ISTRIP)        _RL COSZ(ISTRIP)
78        real dpstrip(ISTRIP,lm)        _RL dpstrip(ISTRIP,lm)
79    
80        real albuvdr(ISTRIP),albuvdf(ISTRIP)        _RL albuvdr(ISTRIP),albuvdf(ISTRIP)
81        real albirdr(ISTRIP),albirdf(ISTRIP)        _RL albirdr(ISTRIP),albirdf(ISTRIP)
82        real difpar (ISTRIP),dirpar (ISTRIP)        _RL difpar (ISTRIP),dirpar (ISTRIP)
83    
84        real fdirir(istrip),fdifir(istrip)        _RL fdirir(istrip),fdifir(istrip)
85        real fdiruv(istrip),fdifuv(istrip)        _RL fdiruv(istrip),fdifuv(istrip)
86    
87        real flux(istrip,lm+1)        _RL flux(istrip,lm+1)
88        real fluxclr(istrip,lm+1)        _RL fluxclr(istrip,lm+1)
89        real dtsw(istrip,lm)        _RL dtsw(istrip,lm)
90        real dtswc(istrip,lm)        _RL dtswc(istrip,lm)
91    
92        real taul   (istrip,lm)        _RL taul   (istrip,lm)
93        real reff   (istrip,lm,2)        _RL reff   (istrip,lm,2)
94        real tauc   (istrip,lm,2)        _RL tauc   (istrip,lm,2)
95        real taua   (istrip,lm)        _RL taua   (istrip,lm)
96        real tstrip (istrip)        _RL tstrip (istrip)
97    
98        logical first        logical first
99        data first /.true./        data first /.true./
# Line 586  C                  tau(im,jm,lm,2):  Sus Line 585  C                  tau(im,jm,lm,2):  Sus
585  C                  tau(im,jm,lm,3):  Raindrops  C                  tau(im,jm,lm,3):  Raindrops
586  C  C
587  C***********************************************************************  C***********************************************************************
 C*                  GODDARD LABORATORY FOR ATMOSPHERES                 *  
 C***********************************************************************  
588    
589        implicit none        implicit none
590    
591        integer  im,jm,lm,i,j,L        integer  im,jm,lm,i,j,L
592    
593        real  tl(im,jm,lm)        _RL  tl(im,jm,lm)
594        real  pl(im,jm,lm)        _RL  pl(im,jm,lm)
595        real ple(im,jm,lm+1)        _RL ple(im,jm,lm+1)
596        real  lz(im,jm,lm)        _RL  lz(im,jm,lm)
597        real  cf(im,jm,lm)        _RL  cf(im,jm,lm)
598        real cfm(im,jm,lm)        _RL cfm(im,jm,lm)
599        real tau(im,jm,lm,3)        _RL tau(im,jm,lm,3)
600        integer lwi(im,jm)        integer lwi(im,jm)
601    
602        real dp, alf, fracls, fraccu        _RL dp, alf, fracls, fraccu
603        real tauice, tauh2o, tauras        _RL tauice, tauh2o, tauras
604    
605  c Compute Cloud Optical Depths  c Compute Cloud Optical Depths
606  c ----------------------------  c ----------------------------
# Line 776  c-----Explicit Inline Directives Line 773  c-----Explicit Inline Directives
773  cfpp$ expand (expmn)  cfpp$ expand (expmn)
774  #endif  #endif
775  #endif  #endif
776        real expmn        _RL expmn
777    
778  c-----input parameters  c-----input parameters
779    
780        integer m,n,ndim,np,ict,icb        integer m,n,ndim,np,ict,icb
781        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)
782        real  taucld(m,ndim,np,2),reff(m,ndim,np,2)        _RL  taucld(m,ndim,np,2),reff(m,ndim,np,2)
783        real  fcld(m,ndim,np),taual(m,ndim,np)        _RL  fcld(m,ndim,np),taual(m,ndim,np)
784        real  rsirbm(m,ndim),rsirdf(m,ndim),        _RL  rsirbm(m,ndim),rsirdf(m,ndim),
785       *     rsuvbm(m,ndim),rsuvdf(m,ndim),cosz(m,ndim),co2       *     rsuvbm(m,ndim),rsuvdf(m,ndim),cosz(m,ndim),co2
786    
787  c-----output parameters  c-----output parameters
788    
789        real  flx(m,ndim,np+1),flc(m,ndim,np+1)        _RL  flx(m,ndim,np+1),flc(m,ndim,np+1)
790        real  fdirir(m,ndim),fdifir(m,ndim)        _RL  fdirir(m,ndim),fdifir(m,ndim)
791        real  fdirpar(m,ndim),fdifpar(m,ndim)        _RL  fdirpar(m,ndim),fdifpar(m,ndim)
792        real  fdiruv(m,ndim),fdifuv(m,ndim)        _RL  fdiruv(m,ndim),fdifuv(m,ndim)
793    
794  c-----temporary array  c-----temporary array
795    
796        integer i,j,k        integer i,j,k
797        real  cc(m,n,3),tauclb(m,n,np),tauclf(m,n,np)        _RL  cc(m,n,3),tauclb(m,n,np),tauclf(m,n,np)
798        real  dp(m,n,np),wh(m,n,np),oh(m,n,np),scal(m,n,np)        _RL  dp(m,n,np),wh(m,n,np),oh(m,n,np),scal(m,n,np)
799        real  swh(m,n,np+1),so2(m,n,np+1),df(m,n,np+1)        _RL  swh(m,n,np+1),so2(m,n,np+1),df(m,n,np+1)
800        real  sdf(m,n),sclr(m,n),csm(m,n),x        _RL  sdf(m,n),sclr(m,n),csm(m,n),x
801    
802  c-----------------------------------------------------------------  c-----------------------------------------------------------------
803    
# Line 1013  c*************************************** Line 1010  c***************************************
1010  c-----input parameters  c-----input parameters
1011    
1012        integer m,n,ndim,np,ict,icb        integer m,n,ndim,np,ict,icb
1013        real  cosz(m,ndim),fcld(m,ndim,np),taucld(m,ndim,np,2)        _RL  cosz(m,ndim),fcld(m,ndim,np),taucld(m,ndim,np,2)
1014    
1015  c-----output parameters  c-----output parameters
1016    
1017        real  cc(m,n,3),tauclb(m,n,np),tauclf(m,n,np)        _RL  cc(m,n,3),tauclb(m,n,np),tauclf(m,n,np)
1018    
1019  c-----temporary variables  c-----temporary variables
1020    
1021        integer i,j,k,im,it,ia,kk        integer i,j,k,im,it,ia,kk
1022        real   fm,ft,fa,xai,taux        _RL   fm,ft,fa,xai,taux
1023    
1024  c-----pre-computed table  c-----pre-computed table
1025    
1026        integer   nm,nt,na        integer   nm,nt,na
1027        parameter (nm=11,nt=9,na=11)        parameter (nm=11,nt=9,na=11)
1028        real   dm,dt,da,t1,caib(nm,nt,na),caif(nt,na)        _RL   dm,dt,da,t1,caib(nm,nt,na),caif(nt,na)
1029        parameter (dm=0.1,dt=0.30103,da=0.1,t1=-0.9031)        parameter (dm=0.1,dt=0.30103,da=0.1,t1=-0.9031)
1030    
1031  c-----include the pre-computed table for cai  c-----include the pre-computed table for cai
# Line 1238  cfpp$ expand (sagpol) Line 1235  cfpp$ expand (sagpol)
1235  cfpp$ expand (expmn)  cfpp$ expand (expmn)
1236  #endif  #endif
1237  #endif  #endif
1238        real expmn        _RL expmn
1239    
1240  c-----input parameters  c-----input parameters
1241    
1242        integer m,n,ndim,np,ict,icb        integer m,n,ndim,np,ict,icb
1243        real  taucld(m,ndim,np,2),reff(m,ndim,np,2),fcld(m,ndim,np)        _RL  taucld(m,ndim,np,2),reff(m,ndim,np,2),fcld(m,ndim,np)
1244        real  tauclb(m,n,np),tauclf(m,n,np),cc(m,n,3)        _RL  tauclb(m,n,np),tauclf(m,n,np),cc(m,n,3)
1245        real  rsirbm(m,ndim),rsirdf(m,ndim)        _RL  rsirbm(m,ndim),rsirdf(m,ndim)
1246        real  wh(m,n,np),taual(m,ndim,np),csm(m,n)        _RL  wh(m,n,np),taual(m,ndim,np),csm(m,n)
1247    
1248  c-----output (updated) parameters  c-----output (updated) parameters
1249    
1250        real  flx(m,ndim,np+1),flc(m,ndim,np+1)        _RL  flx(m,ndim,np+1),flc(m,ndim,np+1)
1251        real  fdirir(m,ndim),fdifir(m,ndim)        _RL  fdirir(m,ndim),fdifir(m,ndim)
1252    
1253  c-----static parameters  c-----static parameters
1254    
1255        integer nk,nband        integer nk,nband
1256        parameter (nk=10,nband=3)        parameter (nk=10,nband=3)
1257        real  xk(nk),hk(nband,nk),ssaal(nband),asyal(nband)        _RL  xk(nk),hk(nband,nk),ssaal(nband),asyal(nband)
1258        real  aia(nband,3),awa(nband,3),aig(nband,3),awg(nband,3)        _RL  aia(nband,3),awa(nband,3),aig(nband,3),awg(nband,3)
1259    
1260  c-----temporary array  c-----temporary array
1261    
1262        integer ib,ik,i,j,k        integer ib,ik,i,j,k
1263        real  ssacl(m,n,np),asycl(m,n,np)        _RL  ssacl(m,n,np),asycl(m,n,np)
1264        real  rr(m,n,np+1,2),tt(m,n,np+1,2),td(m,n,np+1,2),        _RL  rr(m,n,np+1,2),tt(m,n,np+1,2),td(m,n,np+1,2),
1265       *       rs(m,n,np+1,2),ts(m,n,np+1,2)       *       rs(m,n,np+1,2),ts(m,n,np+1,2)
1266        real  fall(m,n,np+1),fclr(m,n,np+1)        _RL  fall(m,n,np+1),fclr(m,n,np+1)
1267        real  fsdir(m,n),fsdif(m,n)        _RL  fsdir(m,n),fsdif(m,n)
1268    
1269        real  tauwv,tausto,ssatau,asysto,tauto,ssato,asyto        _RL  tauwv,tausto,ssatau,asysto,tauto,ssato,asyto
1270        real  taux,reff1,reff2,w1,w2,g1,g2        _RL  taux,reff1,reff2,w1,w2,g1,g2
1271        real  ssaclt(m,n),asyclt(m,n)        _RL  ssaclt(m,n),asyclt(m,n)
1272        real  rr1t(m,n),tt1t(m,n),td1t(m,n),rs1t(m,n),ts1t(m,n)        _RL  rr1t(m,n),tt1t(m,n),td1t(m,n),rs1t(m,n),ts1t(m,n)
1273        real  rr2t(m,n),tt2t(m,n),td2t(m,n),rs2t(m,n),ts2t(m,n)        _RL  rr2t(m,n),tt2t(m,n),td2t(m,n),rs2t(m,n),ts2t(m,n)
1274    
1275  c-----water vapor absorption coefficient for 10 k-intervals.  c-----water vapor absorption coefficient for 10 k-intervals.
1276  c     unit: cm^2/gm  c     unit: cm^2/gm
# Line 1650  cfpp$ expand (sagpol) Line 1647  cfpp$ expand (sagpol)
1647  c-----input parameters  c-----input parameters
1648    
1649        integer m,n,ndim,np,ict,icb        integer m,n,ndim,np,ict,icb
1650        real  taucld(m,ndim,np,2),reff(m,ndim,np,2),fcld(m,ndim,np)        _RL  taucld(m,ndim,np,2),reff(m,ndim,np,2),fcld(m,ndim,np)
1651        real  tauclb(m,n,np),tauclf(m,n,np),cc(m,n,3)        _RL  tauclb(m,n,np),tauclf(m,n,np),cc(m,n,3)
1652        real  oh(m,n,np),dp(m,n,np),taual(m,ndim,np)        _RL  oh(m,n,np),dp(m,n,np),taual(m,ndim,np)
1653        real  rsuvbm(m,ndim),rsuvdf(m,ndim),csm(m,n)        _RL  rsuvbm(m,ndim),rsuvdf(m,ndim),csm(m,n)
1654    
1655  c-----output (updated) parameter  c-----output (updated) parameter
1656    
1657        real  flx(m,ndim,np+1),flc(m,ndim,np+1)        _RL  flx(m,ndim,np+1),flc(m,ndim,np+1)
1658        real  fdirpar(m,ndim),fdifpar(m,ndim)        _RL  fdirpar(m,ndim),fdifpar(m,ndim)
1659        real  fdiruv(m,ndim),fdifuv(m,ndim)        _RL  fdiruv(m,ndim),fdifuv(m,ndim)
1660    
1661  c-----static parameters  c-----static parameters
1662    
1663        integer nband        integer nband
1664        parameter (nband=8)        parameter (nband=8)
1665        real  hk(nband),xk(nband),ry(nband)        _RL  hk(nband),xk(nband),ry(nband)
1666        real  asyal(nband),ssaal(nband),aig(3),awg(3)        _RL  asyal(nband),ssaal(nband),aig(3),awg(3)
1667    
1668  c-----temporary array  c-----temporary array
1669    
1670        integer i,j,k,ib        integer i,j,k,ib
1671        real  taurs,tauoz,tausto,ssatau,asysto,tauto,ssato,asyto        _RL  taurs,tauoz,tausto,ssatau,asysto,tauto,ssato,asyto
1672        real  taux,reff1,reff2,g1,g2,asycl(m,n,np)        _RL  taux,reff1,reff2,g1,g2,asycl(m,n,np)
1673        real  td(m,n,np+1,2),rr(m,n,np+1,2),tt(m,n,np+1,2),        _RL  td(m,n,np+1,2),rr(m,n,np+1,2),tt(m,n,np+1,2),
1674       *       rs(m,n,np+1,2),ts(m,n,np+1,2)       *       rs(m,n,np+1,2),ts(m,n,np+1,2)
1675        real  fall(m,n,np+1),fclr(m,n,np+1),fsdir(m,n),fsdif(m,n)        _RL  fall(m,n,np+1),fclr(m,n,np+1),fsdir(m,n),fsdif(m,n)
1676        real  asyclt(m,n)        _RL  asyclt(m,n)
1677        real  rr1t(m,n),tt1t(m,n),td1t(m,n),rs1t(m,n),ts1t(m,n)        _RL  rr1t(m,n),tt1t(m,n),td1t(m,n),rs1t(m,n),ts1t(m,n)
1678        real  rr2t(m,n),tt2t(m,n),td2t(m,n),rs2t(m,n),ts2t(m,n)        _RL  rr2t(m,n),tt2t(m,n),td2t(m,n),rs2t(m,n),ts2t(m,n)
1679    
1680  c-----hk is the fractional extra-terrestrial solar flux.  c-----hk is the fractional extra-terrestrial solar flux.
1681  c     the sum of hk is 0.47074.  c     the sum of hk is 0.47074.
# Line 1960  c-----Explicit Inline Directives Line 1957  c-----Explicit Inline Directives
1957  cfpp$ expand (expmn)  cfpp$ expand (expmn)
1958  #endif    #endif  
1959  #endif  #endif
1960        real expmn        _RL expmn
1961    
1962        real  zero,one,two,three,four,fourth,seven,tumin        _RL  zero,one,two,three,four,fourth,seven,tumin
1963        parameter (one=1., three=3.)        parameter (one=1., three=3.)
1964        parameter (seven=7., two=2.)        parameter (seven=7., two=2.)
1965        parameter (four=4., fourth=.25)        parameter (four=4., fourth=.25)
1966        parameter (zero=0., tumin=1.e-20)        parameter (zero=0., tumin=1.e-20)
1967    
1968  c-----input parameters  c-----input parameters
1969        real  tau,ssc,g0,csm        _RL  tau,ssc,g0,csm
1970    
1971  c-----output parameters  c-----output parameters
1972        real  rr,tt,td        _RL  rr,tt,td
1973    
1974  c-----temporary parameters  c-----temporary parameters
1975    
1976        real  zth,ff,xx,taup,sscp,gp,gm1,gm2,gm3,akk,alf1,alf2,        _RL  zth,ff,xx,taup,sscp,gp,gm1,gm2,gm3,akk,alf1,alf2,
1977       *     all,bll,st7,st8,cll,dll,fll,ell,st1,st2,st3,st4       *     all,bll,st7,st8,cll,dll,fll,ell,st1,st2,st3,st4
1978  c  c
1979                  zth = one / csm                  zth = one / csm
# Line 2084  c-----Explicit Inline Directives Line 2081  c-----Explicit Inline Directives
2081  cfpp$ expand (expmn)  cfpp$ expand (expmn)
2082  #endif    #endif  
2083  #endif  #endif
2084        real expmn        _RL expmn
2085    
2086        real  one,three,four        _RL  one,three,four
2087        parameter (one=1., three=3., four=4.)        parameter (one=1., three=3., four=4.)
2088    
2089  c-----output parameters:  c-----output parameters:
2090    
2091        real  tau,ssc,g0        _RL  tau,ssc,g0
2092    
2093  c-----output parameters:  c-----output parameters:
2094    
2095        real  rll,tll        _RL  rll,tll
2096    
2097  c-----temporary arrays  c-----temporary arrays
2098    
2099        real  xx,uuu,ttt,emt,up1,um1,st1        _RL  xx,uuu,ttt,emt,up1,um1,st1
2100    
2101               xx  = one-ssc*g0               xx  = one-ssc*g0
2102               uuu = sqrt( xx/(one-ssc))               uuu = sqrt( xx/(one-ssc))
# Line 2123  c*************************************** Line 2120  c***************************************
2120  c compute exponential for arguments in the range 0> fin > -10.  c compute exponential for arguments in the range 0> fin > -10.
2121  c*******************************************************************  c*******************************************************************
2122        implicit none        implicit none
2123        real  fin,expmn        _RL  fin,expmn
2124    
2125        real one,expmin,e1,e2,e3,e4        _RL one,expmin,e1,e2,e3,e4
2126        parameter (one=1.0, expmin=-10.0)        parameter (one=1.0, expmin=-10.0)
2127        parameter (e1=1.0,        e2=-2.507213e-1)        parameter (e1=1.0,        e2=-2.507213e-1)
2128        parameter (e3=2.92732e-2, e4=-3.827800e-3)        parameter (e3=2.92732e-2, e4=-3.827800e-3)
# Line 2178  c-----input parameters Line 2175  c-----input parameters
2175    
2176        integer m,n,np,ict,icb        integer m,n,np,ict,icb
2177    
2178        real  rr(m,n,np+1,2),tt(m,n,np+1,2),td(m,n,np+1,2)        _RL  rr(m,n,np+1,2),tt(m,n,np+1,2),td(m,n,np+1,2)
2179        real  rs(m,n,np+1,2),ts(m,n,np+1,2)        _RL  rs(m,n,np+1,2),ts(m,n,np+1,2)
2180        real  cc(m,n,3)        _RL  cc(m,n,3)
2181    
2182  c-----temporary array  c-----temporary array
2183    
2184        integer i,j,k,ih,im,is        integer i,j,k,ih,im,is
2185        real  rra(m,n,np+1,2,2),tta(m,n,np+1,2,2),tda(m,n,np+1,2,2)        _RL  rra(m,n,np+1,2,2),tta(m,n,np+1,2,2),tda(m,n,np+1,2,2)
2186        real  rsa(m,n,np+1,2,2),rxa(m,n,np+1,2,2)        _RL  rsa(m,n,np+1,2,2),rxa(m,n,np+1,2,2)
2187        real  ch(m,n),cm(m,n),ct(m,n),flxdn(m,n,np+1)        _RL  ch(m,n),cm(m,n),ct(m,n),flxdn(m,n,np+1)
2188        real  fdndir(m,n),fdndif(m,n),fupdif        _RL  fdndir(m,n),fdndif(m,n),fupdif
2189        real  denm,xx        _RL  denm,xx
2190    
2191  c-----output parameters  c-----output parameters
2192    
2193        real  fclr(m,n,np+1),fall(m,n,np+1)        _RL  fclr(m,n,np+1),fall(m,n,np+1)
2194        real  fsdir(m,n),fsdif(m,n)        _RL  fsdir(m,n),fsdif(m,n)
2195    
2196  c-----initialize all-sky flux (fall) and surface downward fluxes  c-----initialize all-sky flux (fall) and surface downward fluxes
2197    
# Line 2494  c     due to co2 absorption. Line 2491  c     due to co2 absorption.
2491  c-----input parameters  c-----input parameters
2492    
2493        integer m,n,np        integer m,n,np
2494        real  csm(m,n),swc(m,n,np+1),swh(m,n,np+1),cah(22,19)        _RL  csm(m,n),swc(m,n,np+1),swh(m,n,np+1),cah(22,19)
2495    
2496  c-----output (undated) parameter  c-----output (undated) parameter
2497    
2498        real  df(m,n,np+1)        _RL  df(m,n,np+1)
2499    
2500  c-----temporary array  c-----temporary array
2501    
2502        integer i,j,k,ic,iw        integer i,j,k,ic,iw
2503        real  xx,clog,wlog,dc,dw,x1,x2,y2        _RL  xx,clog,wlog,dc,dw,x1,x2,y2
2504    
2505  c********************************************************************  c********************************************************************
2506  c-----include co2 look-up table  c-----include co2 look-up table
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22