Context Navigation

clesphys.h @ 3556

Last change on this file since 3556 was 3035, checked in by amartinez, 16 months ago

======= VENUS PCM ===============

COMMIT BY ANTOINE MARTINEZ

Implementation of vertical ambipolar diffusion in physics

=================================

NEW KEYWORD OF PHYSIQ.DEF

=================================

NEW VERSION OF "physiq.def"

deftank/physiq-96x96x90-chemistry-IONOSPHERE-IONDIFF.def

ok_iondiff: keyword supposed to activate ion ambipolar diffusion

nbapp_chem: replaces nbapp_chim, in order to characterize the Number of calls to the chemistry routines (per Venusian day)

================

phyvenus

================

Iondiff_red.F

Calculation of the Ion Ambipolar Diffusion for 13 ions on 14:

O2+, O+, C+, N+, CO2+,
NO+, CO+, H+, H2O+, H3O+,

HCO+, N2+, OH+

The ion temperature is fixed as the half of the electron temperature
calculated in ion_h.F for the stability of the program and because the
ion temperature is lower than electron temperature.

plasmaphys_venus_mod.F90

parameters of the ambipolar diffusion scheme:

parameter (Pdiffion = 7.0D-4) ! pressure in Pa below which ion diffusion is computed
parameter (naltvert = 168) ! number of level on the altimetric subgrid
parameter (nsubvert = 24000) ! ptimephysiq/nsubvert - minimum sub-timestep allowed
parameter ( nsubmin = 40) ! ptimephysiq/nsubmin - maximum sub-timestep allowed

physic_mod.F

nbapp_chem is not fixed anymore here but deftank/in physiq.def
Ambipolar diffusion activated if (ok_iondiff) is true

conf_phys.f90

add ok_iondiff as parameters to read from physiq.def file set to .false. by default
add nbapp_chem as parameters to characterize the Number of calls to the chemistry routines (per Venusian day) instead of be fixed in physic_mod.F

to read from physiq.def file set to 24000 by default

cleshphys.h

added ok_iondiff & nbapp_chem in COMMON/clesphys_l/
removed nbapp_chim

phytrac_chemistry.F

Added security in the calculation of the sza_local in order to avoid the rare case where the |range| is above 1

File size: 1.8 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 ok_iondiff
14	LOGICAL callnlte,callnirco2,callthermos
15	LOGICAL ok_cloud, ok_chem, reinit_trac, ok_sedim
16	LOGICAL ok_clmain, physideal, startphy_file
17	INTEGER nbapp_rad, nbapp_chem, iflag_con, iflag_ajs
18	INTEGER lev_histins, lev_histday, lev_histmth
19	INTEGER tr_scheme, cl_scheme
20	INTEGER nircorr, nltemodel, solvarmod
21	INTEGER nb_mode
22	INTEGER solarchoice
23	REAL ecriphy
24	REAL z0, lmixmin
25	REAL ksta, inertie
26	REAL euveff, fixed_euv_value
27
28	COMMON/clesphys_l/ cycle_diurne, soil_model, &
29	& ok_orodr, ok_orolf, ok_gw_nonoro, ok_kzmin, &
30	& tuneupperatm,callnlte,callnirco2,callthermos, &
31	& ok_cloud, ok_chem, reinit_trac, ok_sedim, &
32	& ok_clmain, physideal, startphy_file, ok_ionchem, ok_jonline, &
33	& ok_iondiff
34
35	COMMON/clesphys_i/ nbapp_rad, nbapp_chem, &
36	& iflag_con, iflag_ajs, &
37	& lev_histins, lev_histday, lev_histmth, tr_scheme, &
38	& cl_scheme, nircorr, nltemodel, solvarmod, nb_mode, &
39	& solarchoice
40
41	COMMON/clesphys_r/ ecriphy, z0, lmixmin, &
42	& ksta, inertie, euveff, fixed_euv_value
43

Note: See TracBrowser for help on using the repository browser.