|
Last change
on this file since 848 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:
876 bytes
|
| Line | |
|---|
| 1 | subroutine su_watercycle |
|---|
| 2 | |
|---|
| 3 | ! to use 'getin' |
|---|
| 4 | use ioipsl_getincom |
|---|
| 5 | use watercommon_h |
|---|
| 6 | implicit none |
|---|
| 7 | #include "comcstfi.h" |
|---|
| 8 | #include "callkeys.h" |
|---|
| 9 | |
|---|
| 10 | |
|---|
| 11 | !================================================================== |
|---|
| 12 | ! |
|---|
| 13 | ! Purpose |
|---|
| 14 | ! ------- |
|---|
| 15 | ! Set up relevant constants and parameters for the water cycle, and water cloud properties |
|---|
| 16 | ! |
|---|
| 17 | ! Authors |
|---|
| 18 | ! ------- |
|---|
| 19 | ! Robin Wordsworth (2010) |
|---|
| 20 | ! Jeremy Leconte (2012) |
|---|
| 21 | ! |
|---|
| 22 | !================================================================== |
|---|
| 23 | |
|---|
| 24 | epsi = mH2O / mugaz |
|---|
| 25 | RCPD = cpp |
|---|
| 26 | |
|---|
| 27 | !RV = 1000.*R/mH2O |
|---|
| 28 | RV = 1000.*8.314/mH2O ! caution! R is R/mugaz already! |
|---|
| 29 | |
|---|
| 30 | RCPV = 4. *RV ! could be more precise... |
|---|
| 31 | !RCPV = 5/2 *8.31*1000/mH2O ! BC modif here: doesn't appear to be correct |
|---|
| 32 | |
|---|
| 33 | RVTMP2 = RCPV/RCPD-1. ! not currently used... |
|---|
| 34 | |
|---|
| 35 | end subroutine su_watercycle |
|---|
Note: See
TracBrowser
for help on using the repository browser.
