source: trunk/LMDZ.GENERIC/libf/phystd/callkeys.h @ 728

Last change on this file since 728 was 728, checked in by jleconte, 12 years ago

18/07/2012 == JL

  • New water cycle scheme:
    • largescale now in F90. Robustness increased by i) including evap inside largescale ii) computing the

condensed water amount iteratively

  • same improvements in moistadj.
  • Water thermodynamical data and saturation curves centralized in module watercommn_h
    • The saturation curves used are now Tetens formula as they are analyticaly inversible (Ts(P)-> Ps(T)).

New saturation curve yields very good agreement with the former one.

  • Saturation curves are now generalized for arbitrary water amount (not just q<<1)
  • The old watersat should be removed soon.
  • The effect of water vapor on total (surface) pressure can be taken into account by setting

mass_redistrib=.true. in callphys.def (routine mass_redistribution inspired from co2_condense in martian
model but with a different scheme as many routines evaporate/condense water vapor).

  • New cloud and precipitation scheme (JL + BC):
    • The default recovery assumption for computing the total cloud fraction has been changed (total random gave too

large cloud fractions). See totalcloudfrac.F90 for details and to change this.

  • Totalcloudfraction now set the total cloud fraction to the fraction of the

optically thickest cloud and totalcloudfrac is thus called in aeropacity.

  • Only the total cloud fraction is used to compute optical depth in aeropacity (no more effective

optical depth with exponential formula).

  • 4 precipitation schemes are now available (see rain.F90 for details). The choice can be made using precip_scheme

in callphys.def. Usage of the more physically based model of Boucher et al 95 (precip_scheme=4) is recommended.
default behavior is set to the former "simple scheme" (precip_scheme=1).

  • See rain.f90 to determine the parameter to be defined in callphys.def as a function of the precipitation scheme used.
  • Physiq.F90 now written in a matricial (more F90) way.
  • Radii (H2O and CO2 cloud particles, aerosols, duts, ...) calculations now centralized in module radii_mod.F90

and work with the new aerosol scheme implemented by Laura K. Some inconsistency may remain in callsedim.


Implementation compiled with ifort and pgf90.
gcm.e runs in Earth and Early Mars case with CO2 and H2O cycle + dust.

File size: 3.5 KB
Line 
1!
2! For Fortran 77/Fortran 90 compliance always use line continuation
3! symbols '&' in columns 73 and 6
4!
5      COMMON/callkeys/callrad,corrk,calldifv,UseTurbDiff,calladj        &
6     &   , co2cond,callsoil                                             &
7     &   , season,diurnal,tlocked,iradia,lwrite                         &
8     &   , iaervar,iddist,topdustref,callstats,calleofdump              &
9     &   , enertest                                                     &
10     &   , callgasvis,Continuum,H2Ocont_simple,graybody                 &
11     &   , Nmix_co2, radfixed, dusttau                                  &
12     &   , meanOLR, specOLR                                             &
13     &   , kastprof, noradsurf, Tstrat                                  &
14     &   , newtonian, tau_relax, testradtimes                           &
15     &   , check_cpp_match, force_cpp                                   &
16     &   , rayleigh                                                     &
17     &   , stelbbody, stelTbb                                           &
18     &   , tplanet                                                      &
19     &   , startype, Fat1AU                                             &
20     &   , nearco2cond                                                  &
21     &   , tracer, mass_redistrib, varactive, varfixed, satval          &
22     &   , sedimentation,water,watercond,waterrain                      &
23     &   , aeroco2,aeroh2o,aeroh2so4                                    &
24     &   , aerofixco2,aerofixh2o                                        &
25     &   , hydrology, sourceevol, icetstep, albedosnow                  &
26     &   , maxicethick, Tsaldiff                                        &
27     &   , CLFfixval, CLFvarying                                        &
28     &   , n2mixratio                                                   &
29     &   , co2supsat                                                    &
30     &   , cloudlvl                                                     &
31     &   , pceil                                                   
32
33      logical callrad,corrk,calldifv,UseTurbDiff                        &
34     &   , calladj,co2cond,callsoil                                     &
35     &   , season,diurnal,tlocked,lwrite                                &
36     &   , callstats,calleofdump                                        &
37     &   , callgasvis,Continuum,H2Ocont_simple,graybody
38
39      logical enertest
40      logical nonideal
41      logical meanOLR
42      logical specOLR
43      logical kastprof
44      logical newtonian
45      logical check_cpp_match
46      logical force_cpp
47      logical testradtimes
48      logical rayleigh
49      logical stelbbody
50      logical ozone
51      logical nearco2cond
52      logical tracer
53      logical mass_redistrib
54      logical varactive
55      logical varfixed
56      logical radfixed
57      logical sedimentation
58      logical water,watercond,waterrain
59      logical aeroco2,aeroh2o,aeroh2so4
60      logical aerofixco2,aerofixh2o
61      logical hydrology
62      logical sourceevol
63      logical CLFvarying
64      logical noradsurf
65
66      integer iddist
67      integer iaervar
68      integer iradia
69      integer startype
70
71      real topdustref
72      real Nmix_co2
73      real dusttau
74      real Fat1AU
75      real stelTbb
76      real Tstrat
77      real tplanet
78      real satval
79      real CLFfixval
80      real n2mixratio
81      real co2supsat
82      real pceil
83      real albedosnow
84      real maxicethick
85      real Tsaldiff
86      real tau_relax
87      real cloudlvl
88      real icetstep
89
Note: See TracBrowser for help on using the repository browser.