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

Rev	Line
[3]	1	!
[97]	2	! ATTENTION!!!!: ce fichier include est compatible format fixe/format libre
[1310]	3	! veillez n utiliser que des ! pour les commentaires
[97]	4	! et bien positionner les & des lignes de continuation
	5	! (les placer en colonne 6 et en colonne 73)
[3]	6	!
[97]	7	!
	8	!..include cles_phys.h
	9	!
	10	LOGICAL cycle_diurne,soil_model
[3]	11	LOGICAL ok_orodr,ok_orolf,ok_gw_nonoro
[2836]	12	LOGICAL ok_kzmin,tuneupperatm, ok_ionchem, ok_jonline
[3035]	13	LOGICAL ok_iondiff
[1310]	14	LOGICAL callnlte,callnirco2,callthermos
[1642]	15	LOGICAL ok_cloud, ok_chem, reinit_trac, ok_sedim
[2135]	16	LOGICAL ok_clmain, physideal, startphy_file
[3035]	17	INTEGER nbapp_rad, nbapp_chem, iflag_con, iflag_ajs
[1305]	18	INTEGER lev_histins, lev_histday, lev_histmth
[1661]	19	INTEGER tr_scheme, cl_scheme
[1310]	20	INTEGER nircorr, nltemodel, solvarmod
[2135]	21	INTEGER nb_mode
[2560]	22	INTEGER solarchoice
[1120]	23	REAL ecriphy
[97]	24	REAL z0, lmixmin
	25	REAL ksta, inertie
[2464]	26	REAL euveff, fixed_euv_value
[3]	27
[1305]	28	COMMON/clesphys_l/ cycle_diurne, soil_model, &
	29	& ok_orodr, ok_orolf, ok_gw_nonoro, ok_kzmin, &
[2580]	30	& tuneupperatm,callnlte,callnirco2,callthermos, &
[2135]	31	& ok_cloud, ok_chem, reinit_trac, ok_sedim, &
[3035]	32	& ok_clmain, physideal, startphy_file, ok_ionchem, ok_jonline, &
	33	& ok_iondiff
[3]	34
[3035]	35	COMMON/clesphys_i/ nbapp_rad, nbapp_chem, &
[2135]	36	& iflag_con, iflag_ajs, &
[1310]	37	& lev_histins, lev_histday, lev_histmth, tr_scheme, &
[2560]	38	& cl_scheme, nircorr, nltemodel, solvarmod, nb_mode, &
	39	& solarchoice
[1160]	40
[2487]	41	COMMON/clesphys_r/ ecriphy, z0, lmixmin, &
[2464]	42	& ksta, inertie, euveff, fixed_euv_value
[1160]	43

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