source: codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/config/ppsrc/phys/su_watercycle.f90 @ 224

Last change on this file since 224 was 224, checked in by ymipsl, 10 years ago
File size: 5.6 KB
Line 
1      subroutine su_watercycle
2
3      use watercommon_h
4      implicit none
5!-----------------------------------------------------------------------
6! INCLUDE "comcstfi.h"
7
8      common/comcstfi/pi,rad,g,r,cpp,rcp,dtphys,daysec,mugaz,omeg
9      common/comcstfi/avocado!,molrad,visc
10     
11      real pi,rad,g,r,cpp,rcp,dtphys,daysec,mugaz,omeg
12      real avocado!,molrad,visc
13
14!
15! For Fortran 77/Fortran 90 compliance always use line continuation
16! symbols '&' in columns 73 and 6
17!
18! Group commons according to their type for minimal performance impact
19
20      COMMON/callkeys_l/callrad,corrk,calldifv,UseTurbDiff,calladj      &
21     &   , co2cond,callsoil                                             &
22     &   , season,diurnal,tlocked,rings_shadow,lwrite                   &
23     &   , callstats,calleofdump                                        &
24     &   , enertest                                                     &
25     &   , callgasvis,continuum,H2Ocont_simple,graybody                 &
26     &   , radfixed                                                     &
27     &   , meanOLR, specOLR                                             &
28     &   , kastprof                                                     &
29     &   , nosurf, oblate                                               &     
30     &   , newtonian, testradtimes                                      &
31     &   , check_cpp_match, force_cpp                                   &
32     &   , rayleigh                                                     &
33     &   , stelbbody                                                    &
34     &   , nearco2cond                                                  &
35     &   , tracer, mass_redistrib, varactive, varfixed                  &
36     &   , sedimentation,water,watercond,waterrain                      &
37     &   , aeroco2,aeroh2o,aeroh2so4,aeroback2lay                       &
38     &   , aerofixco2,aerofixh2o                                        &
39     &   , hydrology, sourceevol                                        &
40     &   , CLFvarying                                                   &
41     &   , strictboundcorrk                                             &                                       
42     &   , ok_slab_ocean                                                &
43     &   , ok_slab_sic                                                  &
44     &   , ok_slab_heat_transp                                         
45
46
47      COMMON/callkeys_i/iaervar,iddist,iradia,startype
48     
49      COMMON/callkeys_r/topdustref,Nmix_co2,dusttau,Fat1AU,stelTbb,     &
50     &                  Tstrat,tplanet,obs_tau_col_tropo,               &
51     &                  obs_tau_col_strato,pres_bottom_tropo,           &
52     &                  pres_top_tropo,pres_bottom_strato,              &
53     &                  pres_top_strato,size_tropo,size_strato,satval,  &
54     &                  CLFfixval,n2mixratio,co2supsat,pceil,albedosnow,&
55     &                  maxicethick,Tsaldiff,tau_relax,cloudlvl,        &
56     &                  icetstep,intheat,flatten,Rmean,J2,MassPlanet
57     
58      logical callrad,corrk,calldifv,UseTurbDiff                        &
59     &   , calladj,co2cond,callsoil                                     &
60     &   , season,diurnal,tlocked,rings_shadow,lwrite                   &
61     &   , callstats,calleofdump                                        &
62     &   , callgasvis,continuum,H2Ocont_simple,graybody                 &
63     &   , strictboundcorrk                                             
64
65      logical enertest
66      logical nonideal
67      logical meanOLR
68      logical specOLR
69      logical kastprof
70      logical newtonian
71      logical check_cpp_match
72      logical force_cpp
73      logical testradtimes
74      logical rayleigh
75      logical stelbbody
76      logical ozone
77      logical nearco2cond
78      logical tracer
79      logical mass_redistrib
80      logical varactive
81      logical varfixed
82      logical radfixed
83      logical sedimentation
84      logical water,watercond,waterrain
85      logical aeroco2,aeroh2o,aeroh2so4,aeroback2lay
86      logical aerofixco2,aerofixh2o
87      logical hydrology
88      logical sourceevol
89      logical CLFvarying
90      logical nosurf
91      logical oblate
92      logical ok_slab_ocean
93      logical ok_slab_sic
94      logical ok_slab_heat_transp
95
96      integer iddist
97      integer iaervar
98      integer iradia
99      integer startype
100
101      real topdustref
102      real Nmix_co2
103      real dusttau
104      real Fat1AU
105      real stelTbb
106      real Tstrat
107      real tplanet
108      real obs_tau_col_tropo
109      real obs_tau_col_strato
110      real pres_bottom_tropo
111      real pres_top_tropo
112      real pres_bottom_strato
113      real pres_top_strato
114      real size_tropo
115      real size_strato
116      real satval
117      real CLFfixval
118      real n2mixratio
119      real co2supsat
120      real pceil
121      real albedosnow
122      real maxicethick
123      real Tsaldiff
124      real tau_relax
125      real cloudlvl
126      real icetstep
127      real intheat
128      real flatten
129      real Rmean
130      real J2
131      real MassPlanet
132
133
134!==================================================================
135!
136!     Purpose
137!     -------
138!     Set up relevant constants and parameters for the water cycle, and water cloud properties
139!
140!     Authors
141!     -------
142!     Robin Wordsworth (2010)
143!     Jeremy Leconte (2012)
144!
145!==================================================================
146
147      epsi   = mH2O / mugaz
148      RCPD   = cpp 
149
150      !RV = 1000.*R/mH2O
151      RV = 1000.*8.314/mH2O ! caution! R is R/mugaz already!
152
153      RCPV   = 1.88e3 ! specific heat capacity of water vapor at 350K
154
155      RVTMP2 = RCPV/RCPD-1. ! not currently used...
156
157      end subroutine su_watercycle
Note: See TracBrowser for help on using the repository browser.