/[MITgcm]/MITgcm/pkg/fizhi/fizhi_driver.F
ViewVC logotype

Diff of /MITgcm/pkg/fizhi/fizhi_driver.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.9 by molod, Thu Jul 8 20:24:08 2004 UTC revision 1.13 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"
5        subroutine fizhi_driver (myid,im,jm,lm,bi,bj,ptracer,ntracer,        subroutine fizhi_driver (myid,im,jm,lm,bi,bj,ptracer,ntracer,
6       . xlats,xlons,       . xlats,xlons,
7       . p,u,v,t,q,pl,ple,dpres,pkht,pkl,fracland,landtype,radswt,       . p,u,v,t,q,pl,ple,dpres,pkht,pkl,fracland,landtype,radswt,
8       . phis_var,tgz,sea_ice,       . phis_var,tgz,sea_ice,
9       . nchp,chlat,chlon,igrd,nchpland,chfr,ityp,       . nchp,chlat,chlon,igrd,nchptot,nchpland,chfr,ityp,
10       . tcanopy,tdeep,ecanopy,swetshal,swetroot,swetdeep,capac,snodep,       . tcanopy,tdeep,ecanopy,swetshal,swetroot,swetdeep,capac,snodep,
11       . ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke,       . ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke,
12       . albvisdr,albvisdf,albnirdr,albnirdf,emiss,alai,agrn,       . albvisdr,albvisdf,albnirdr,albnirdf,emiss,alai,agrn,
13       . chemq,chemo3,co2,cfc11,cfc12,cfc22,methane,n2o,       . chemq,chemo3,co2,cfc11,cfc12,cfc22,methane,n2o,
14       . lwdt,lwdtclr,swdt,swdtclr,turbu,turbv,turbt,turbq,       . lwdt,lwdtclr,swdt,swdtclr,turbu,turbv,turbt,turbq,
15       . moistu,moistv,moistt,moistq,       . moistu,moistv,moistt,moistq,
16       . radswg,swgclr,fdirpar,fdifpar,osr,osrclr,tg0,radlwg,       . radswg,swgclr,fdirpar,fdifpar,osr,osrclr,tg0,radlwg,lwgclr,
17       . st4,dst4,dlwdtg,rainlsp,raincon,snowfall,iras,nlwcld,       . st4,dst4,dlwdtg,rainlsp,raincon,snowfall,iras,nlwcld,
18       . cldtot_lw,clras_lw,cldlsp_lw,nlwlz,lwlz,nswcld,cldtot_sw,       . cldtot_lw,cldras_lw,cldlsp_lw,nlwlz,lwlz,nswcld,cldtot_sw,
19       . clras_sw,cldlsp_sw,nswlz,swlz,imstturbsw,imstturblw,qliqavesw,       . cldras_sw,cldlsp_sw,nswlz,swlz,imstturbsw,imstturblw,qliqavesw,
20       . qliqavelw,fccavesw,fccavelw,qq)       . qliqavelw,fccavesw,fccavelw,qq)
21  C***********************************************************************          C***********************************************************************        
22  C  Purpose                                                                        C  Purpose                                                                      
# Line 42  c ----------------- Line 42  c -----------------
42    
43  c Timers Common  c Timers Common
44  c -------------  c -------------
45  #include 'chronos.h'  #include "chronos.h"
46    
47  c Input Parameters  c Input Parameters
48  c ----------------  c ----------------
49        integer myid,im,jm,lm,bi,bj,ptracer,ntracer        integer myid,im,jm,lm,bi,bj,ptracer,ntracer
50        integer nchp,igrd(nchp),nchpland,ityp(nchp)        integer nchp,igrd(nchp),nchptot,nchpland,ityp(nchp)
51        integer iras,nlwcld,nlwlz,nswcld,nswlz        integer iras,nlwcld,nlwlz,nswcld,nswlz
52        integer imstturbsw,imstturblw        integer imstturbsw,imstturblw
53    
54        real xlats(im,jm), xlons(im,jm)        _RL xlats(im,jm), xlons(im,jm)
55        real p(im,jm)        _RL p(im,jm)
56        real u(im,jm,lm), v(im,jm,lm), t(im,jm,lm)        _RL u(im,jm,lm),v(im,jm,lm),t(im,jm,lm)
57        real q(im,jm,lm,ntracer)        _RL q(im,jm,lm,ntracer)
58        real pl(im,jm,lm),ple(im,jm,lm+1),dpres(im,jm,lm),pkht(im,jm,lm)        _RL pl(im,jm,lm),ple(im,jm,lm+1),dpres(im,jm,lm)
59        real pkl(im,jm,lm)        _RL pkht(im,jm,lm)
60        real radswt(im,jm), fracland(im,jm), landtype(im,jm)        _RL pkl(im,jm,lm)
61        real phis_var(im,jm), sea_ice(im,jm)        _RL radswt(im,jm), fracland(im,jm), landtype(im,jm)
62        real chlat(nchp),chlon(nchp),chfr(nchp)        _RL phis_var(im,jm), sea_ice(im,jm)
63        real tcanopy(nchp),tdeep(nchp),ecanopy(nchp),swetshal(nchp)        _RL chlat(nchp),chlon(nchp),chfr(nchp)
64        real swetroot(nchp),swetdeep(nchp),capac(nchp),snodep(nchp)        _RL tcanopy(nchp),tdeep(nchp),ecanopy(nchp),swetshal(nchp)
65        real ctmt(nchp),xxmt(nchp),yymt(nchp),zetamt(nchp)        _RL swetroot(nchp),swetdeep(nchp),capac(nchp),snodep(nchp)
66        real xlmt(nchp,lm),khmt(nchp,lm),tke(nchp,lm)        _RL ctmt(nchp),xxmt(nchp),yymt(nchp),zetamt(nchp)
67        real co2,cfc11,cfc12,cfc22,methane(lm),n2o(lm)        _RL xlmt(nchp,lm),khmt(nchp,lm),tke(nchp,lm)
68        real lwdt(im,jm,lm),lwdtclr(im,jm,lm)        _RL co2,cfc11,cfc12,cfc22,methane(lm),n2o(lm)
69        real swdt(im,jm,lm),swdtclr(im,jm,lm)        _RL lwdt(im,jm,lm),lwdtclr(im,jm,lm)
70        real turbu(im,jm,lm),turbv(im,jm,lm),turbt(im,jm,lm)        _RL swdt(im,jm,lm),swdtclr(im,jm,lm)
71        real turbq(im,jm,lm)        _RL turbu(im,jm,lm),turbv(im,jm,lm),turbt(im,jm,lm)
72        real moistu(im,jm,lm),moistv(im,jm,lm),moistt(im,jm,lm)        _RL turbq(im,jm,lm)
73        real moistq(im,jm,lm)        _RL moistu(im,jm,lm),moistv(im,jm,lm),moistt(im,jm,lm)
74        real chemo3(im,jm,lm),chemq(im,jm,lm)        _RL moistq(im,jm,lm)
75        real albvisdr(im,jm),albvisdf(im,jm),albnirdr(im,jm)        _RL chemo3(im,jm,lm),chemq(im,jm,lm)
76        real albnirdf(im,jm),emiss(im,jm)        _RL albvisdr(im,jm),albvisdf(im,jm),albnirdr(im,jm)
77        real alai(nchp),agrn(nchp)        _RL albnirdf(im,jm),emiss(im,jm)
78        real radswg(im,jm),swgclr(im,jm)        _RL alai(nchp),agrn(nchp)
79        real fdirpar(im,jm),fdifpar(im,jm),osr(im,jm),osrclr(im,jm)        _RL radswg(im,jm),swgclr(im,jm)
80        real tg0(im,jm),radlwg(im,jm),st4(im,jm),dst4(im,jm)        _RL fdirpar(im,jm),fdifpar(im,jm),osr(im,jm),osrclr(im,jm)
81        real dlwdtg(im,jm,lm)        _RL tg0(im,jm),radlwg(im,jm),lwgclr(im,jm),st4(im,jm)
82        real rainlsp(im,jm),raincon(im,jm),snowfall(im,jm)        _RL dst4(im,jm)
83        real cldtot_lw(im,jm,lm),clras_lw(im,jm,lm),cldlsp_lw(im,jm,lm)        _RL dlwdtg(im,jm,lm)
84        real lwlz(im,jm,lm)        _RL rainlsp(im,jm),raincon(im,jm),snowfall(im,jm)
85        real cldtot_sw(im,jm,lm),clras_sw(im,jm,lm),cldlsp_sw(im,jm,lm)        _RL cldtot_lw(im,jm,lm),cldras_lw(im,jm,lm)
86        real swlz(im,jm,lm)        _RL cldlsp_lw(im,jm,lm)
87        real qliqavesw(im,jm,lm),qliqavelw(im,jm,lm)        _RL lwlz(im,jm,lm)
88        real fccavesw(im,jm,lm),fccavelw(im,jm,lm)        _RL cldtot_sw(im,jm,lm),cldras_sw(im,jm,lm)
89        real qq(im,jm,lm)        _RL cldlsp_sw(im,jm,lm)
90        real tgz(im,jm)        _RL swlz(im,jm,lm)
91          _RL qliqavesw(im,jm,lm),qliqavelw(im,jm,lm)
92          _RL fccavesw(im,jm,lm),fccavelw(im,jm,lm)
93          _RL qq(im,jm,lm)
94          _RL tgz(im,jm)
95    
96  c Local Variables  c Local Variables
97  c ---------------  c ---------------
# Line 96  c --------------- Line 100  c ---------------
100        integer numpcheck        integer numpcheck
101        parameter (numpcheck = 5)        parameter (numpcheck = 5)
102        integer pchecklevs(numpcheck)        integer pchecklevs(numpcheck)
103        real pcheckpress(numpcheck)        _RL pcheckpress(numpcheck)
104        data pcheckpress/950.,750.,700.,400.,10./        data pcheckpress/950.,750.,700.,400.,10./
105    
106        integer low_level,mid_level,nltop,nsubmin,nsubmax,Lup        integer low_level,mid_level,nltop,nsubmin,nsubmax,Lup
107        integer ndmoist,ndturb,ndlw,ndsw        integer ndmoist,ndturb,ndlw,ndsw
108        integer istrip,npcs        integer istrip,npcs
109        integer i,j,L,n        integer i,j
110        integer ndum,ndpnt        integer ndum,ndpnt
111        real dum,akap,getcon        _RL akap,getcon
112        real ptop        _RL ptop
113        logical lpnt        logical lpnt
114    
115        istrip = min( im*jm,71 )        istrip = min( im*jm,71 )
# Line 182  C ========== Line 186  C ==========
186        IF ( alarm('turb') ) THEN        IF ( alarm('turb') ) THEN
187    
188        call turbio (im,jm,lm,1,nymd,nhms,bi,bj,ndturb,        call turbio (im,jm,lm,1,nymd,nhms,bi,bj,ndturb,
189       .   ptop,p,u,v,t,q,ntracer,ptracer,pl,ple,dpres,pkht,pkl,       .  ptop,p,u,v,t,q,ntracer,ptracer,pl,ple,dpres,pkht,pkl,
190       .   ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke,       .  ctmt,xxmt,yymt,zetamt,xlmt,khmt,tke,
191       .   tgz,fracland,landtype,       .  tgz,fracland,landtype,
192       .   tcanopy,ecanopy,tdeep,swetshal,swetroot,swetdeep,snodep,capac,       .  tcanopy,ecanopy,tdeep,swetshal,swetroot,swetdeep,snodep,capac,
193       .   nchp,nchpland,chfr,chlt,chlon,igrd,ityp,alai,agrn,sea_ice,lpnt,       .  nchp,nchptot,nchpland,chfr,chlat,chlon,igrd,ityp,
194       .   turbu,turbv,turbt,turbq,radlwg,st4,dst4,radswg,radswt,       .  alai,agrn,sea_ice,lpnt,
195       .   fdifpar,fdirpar,rainlsp,rainconv,snowfall,tg0,       .  turbu,turbv,turbt,turbq,radlwg,st4,dst4,radswg,radswt,
196       .   imstturblw,imstturbsw,qliqavelw,qliqavesw,fccavelw,fccavesw,qq,       .  fdifpar,fdirpar,rainlsp,raincon,snowfall,tg0,
197       .   myid,comm )       .  imstturblw,imstturbsw,qliqavelw,qliqavesw,fccavelw,fccavesw,qq,
198         .  myid)
199    
200  c Add Gravity-Wave Drag Tendency  c Add Gravity-Wave Drag Tendency
201  c ------------------------------  c ------------------------------
# Line 210  C =============== Line 215  C ===============
215       .    low_level,mid_level,nltop,nsubmin,nsubmax,Lup,       .    low_level,mid_level,nltop,nsubmin,nsubmax,Lup,
216       .    p,pl,ple,dpres,pkht,pkl,t,q,bi,bj,ntracer,ptracer,qq,       .    p,pl,ple,dpres,pkht,pkl,t,q,bi,bj,ntracer,ptracer,qq,
217       .    moistu,moistv,moistt,moistq,im,jm,lm,ptop,iras,       .    moistu,moistv,moistt,moistq,im,jm,lm,ptop,iras,
218       .    rainlsp,rainconv,snowfall,       .    rainlsp,raincon,snowfall,
219       .    nswcld,cldtot_sw,cldras_sw,cldlsp_sw,nswlz,swlz,       .    nswcld,cldtot_sw,cldras_sw,cldlsp_sw,nswlz,swlz,
220       .    nlwcld,cldtot_lw,cldras_lw,cldlsp_lw,nlwlz,lwlz,       .    nlwcld,cldtot_lw,cldras_lw,cldlsp_lw,nlwlz,lwlz,
221       .              .false.,myid)       .              .false.,myid)
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22