source: trunk/LMDZ.VENUS/libf/phyvenus/clesphys.h @ 2937

Last change on this file since 2937 was 2836, checked in by abierjon, 3 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: 1.7 KB
Line 
1!
2!  ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
3!                 veillez  n utiliser que des ! pour les commentaires
4!                 et  bien positionner les & des lignes de continuation
5!                 (les placer en colonne 6 et en colonne 73)
6!
7!
8!..include cles_phys.h
9!
10       LOGICAL cycle_diurne,soil_model
11       LOGICAL ok_orodr,ok_orolf,ok_gw_nonoro
12       LOGICAL ok_kzmin,tuneupperatm, ok_ionchem, ok_jonline
13       LOGICAL callnlte,callnirco2,callthermos
14       LOGICAL ok_cloud, ok_chem, reinit_trac, ok_sedim
15       LOGICAL ok_clmain, physideal, startphy_file
16       INTEGER nbapp_rad, nbapp_chim, iflag_con, iflag_ajs
17       INTEGER lev_histins, lev_histday, lev_histmth
18       INTEGER tr_scheme, cl_scheme
19       INTEGER nircorr, nltemodel, solvarmod
20       INTEGER nb_mode
21       INTEGER solarchoice
22       REAL    ecriphy
23       REAL    z0, lmixmin
24       REAL    ksta, inertie
25       REAL    euveff, fixed_euv_value
26
27       COMMON/clesphys_l/ cycle_diurne, soil_model,                     &
28     &     ok_orodr, ok_orolf, ok_gw_nonoro, ok_kzmin,                  &
29     &     tuneupperatm,callnlte,callnirco2,callthermos,                &
30     &     ok_cloud, ok_chem, reinit_trac, ok_sedim,                    &
31     &     ok_clmain, physideal, startphy_file, ok_ionchem, ok_jonline
32
33       COMMON/clesphys_i/ nbapp_rad, nbapp_chim,                        &
34     &     iflag_con, iflag_ajs,                                        &
35     &     lev_histins, lev_histday, lev_histmth, tr_scheme,            &
36     &     cl_scheme, nircorr, nltemodel, solvarmod, nb_mode,           &
37     &     solarchoice
38
39       COMMON/clesphys_r/ ecriphy, z0, lmixmin,                         &
40     &     ksta, inertie, euveff, fixed_euv_value
41
Note: See TracBrowser for help on using the repository browser.