VENUS GCM:
INCLUDING THE IONOSPHERE CODE IN VENUS GCM
ATTENTION: INCREASED MEMORY DEMAND
NEEDS AT LEAST 135 GB of allocated memory
==============================================================
===> LIST OF APPENDED FILES AND INTERNAL ADDITIONS <==========
==============================================================
NEW MODEL SPECIES
- deftank/traceur-chemistry-IONOSPHERE.def
Neutrals: 4 Species: N, NO, NO2, N(2D)
Ions: 15 species:
CO2+, CO+, O+, O2+, H+,
N2+, H2O+, OH+, C+, HCO+,
H3O+, HCO2+, N+, NO+, elec
FROM 36 to 55 chemical species
NEW KEYWORD OF PHYSIQ.DEF
- deftank/physiq-96x96x90-chemistry-IONOSPHERE.def
- ok_ionchem: keyword supposed to activate ion chemistry.
(be careful that n, no, no2, n2d and the ion species are in the deftank/traceur-chemistry-IONOSPHERE.def)
- ok_jonline: keyword supposed to activate the online photochemistry
(be careful that n, no, no2 and n2d are in the traceur-chemistry-IONOSPHERE.def)
==============================================================
===> LIST OF MODIFIED PROGRAMS <=========================
==============================================================
nonoro_gwd_ran_mod.F90
- Change EPFLUXMAX value from 5.E-3 to 1.E-3
photochemistry_venus.F90 (krates; photolysis_online; indices;)
- Import of keywords: ok_ionchem, tuneupperatm & ok_jonline
- addition of ion species in order
- Forcing electroneutrality
- Update of the reactions a001 and a002 with taking into account the other species
- Change of formula for a002 with the formulation of Baulch et al., 1976 (confirmed by Smith and Robertson, 2008)
photolysis_mod.F90
- Modification of the rdsolarflux subroutine to include interpolation with Atlas1 and Atlas3 in connection with E10.7
concentration2.F90
- Added chemical species n2d, no, no2 and n in the conductivity calculation.
The ions have been excluded because their sum is 105 times less dense than the neutrals and
their thermal conductivity is unknown
iono_h.F90
- Addition of the phdisrate routine
- replace the electronic temperature of Mars by that of
origin = 1: Theis et al. 1980 (Venus) with bilinear interpolation altitude/cos(SZA)
origin = 2: Theis et al. 1984 (Venus) with the formula of the electronic temperature at high solar activity
- addition of an ion temperature model based on VIRA
cleshphys.h
- added ok_ionchem & ok_jonline in COMMON/clesphys_l/
conf_phys.f90
- add ok_ionchem & ok_jonline as parameters to read from physiq.def file set to .false. by default
chemparam_mod.F90
- Add species in M_tr and corresponding i_X. Set all i_X to zero before reading traceur-chemistry-IONOSPHERE.def
- Added Type_tr table to differentiate species: 1 == neutral, 2 == ION, 3 == liquid, 10 == others
euvheat.F90; hrtherm.F; jthermcalc_e107.F; param_read_e107.F
A.M
|