source: trunk/LMDZ.VENUS/deftank/traceur-chemistry-IONOSPHERE.def @ 3471

Last change on this file since 3471 was 2836, checked in by abierjon, 2 years ago

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

  • Normalization with Mars

A.M

File size: 602 bytes
RevLine 
[2836]155
210 10 co2
310 10 co
410 10 h2
510 10 h2o
610 10 o1d
710 10 o
810 10 o2
910 10 o2dg
1010 10 o3
1110 10 h
1210 10 oh
1310 10 ho2
1410 10 h2o2
1510 10 cl
1610 10 clo
1710 10 cl2
1810 10 hcl
1910 10 hocl
2010 10 clco
2110 10 clco3
2210 10 cocl2
2310 10 s
2410 10 so
2510 10 so2
2610 10 so3
2710 10 s2o2
2810 10 ocs
2910 10 hso3
3010 10 h2so4
3110 10 s2
3210 10 clso2
3310 10 oscl
3410 10 n2
3510 10 he
3610 10 n
3710 10 n2d
3810 10 no
3910 10 no2
4010 10 co2plus
4110 10 coplus
4210 10 oplus
4310 10 o2plus
4410 10 hplus
4510 10 n2plus
4610 10 nplus
4710 10 noplus
4810 10 h2oplus
4910 10 ohplus
5010 10 cplus
5110 10 h3oplus
5210 10 hcoplus
5310 10 hco2plus
5410 10 elec
5510 10 h2oliq
5610 10 h2so4liq
Note: See TracBrowser for help on using the repository browser.