source: trunk/LMDZ.MARS/libf/phymars/tracer.h @ 1009

Last change on this file since 1009 was 635, checked in by emillour, 13 years ago

Mars GCM: Update of the chemistry package, including:

  • 93 reactions are accounted for (instead of 22); tracking 28 species (instead of 11)
  • computation of photoabsorption using raytracing
  • improved time stepping in the photochemistry
  • updated parameters (cross-sections); with this new version input files

are in 'EUV/param_v5' of "datafile" directory.

  • transition between lower and upper atmosphere chemistry set to 0.1 Pa (calchim.F90)
  • Lots of code clean-up: removed obsolete files column.F, param_v3.h, flujo.F, phdisrate.F, ch.F, interpfast.F, paramfoto.F, getch.F Converted chemtermos.F -> chemthermos.F90 and euvheat.F -> euvheat.F90. Added paramfoto_compact.F , param_v4.h and iono.h
  • Upadted surfacearea.F
  • Cleaned initracer.F and callkeys.h (removed use of obsolete "nqchem" and "oldnames" case when initializing tracers).
  • Minor correction in "callsedim": compute "rdust" and/or "rice" only when it makes sense.

FGG+FL+EM

File size: 3.8 KB
RevLine 
[38]1!-----------------------------------------------------------------------
2! INCLUDE 'tracer.h'
3
4      character*20  noms(nqmx)  ! name of the tracer
5      real mmol(nqmx)           ! mole mass of tracer (g/mol-1) 
6      real radius(nqmx)   ! dust and ice particle radius (m)
7      real rho_q(nqmx)    ! tracer densities (kg.m-3)
8      real alpha_lift(nqmx) ! saltation vertical flux/horiz flux ratio (m-1)
9      real alpha_devil(nqmx) ! lifting coeeficient by dust devil
10
11      real varian      ! Characteristic variance of log-normal distribution
12      real r3n_q     ! used to compute r0 from number and mass mixing ratio
13      real rho_dust     ! Mars dust density (kg.m-3)
14      real rho_ice     ! Water ice density (kg.m-3)
15      real nuice_ref   ! Effective variance of the water ice dist.
[358]16      real nuice_sed   ! Sedimentation effective variance of the water ice dist.
[38]17      real ref_r0        ! for computing reff=ref_r0*r0 (in log.n. distribution)
[420]18     
19      real ccn_factor  ! ratio of nuclei for water ice particles
[38]20
21      real dryness(ngridmx)!"Dryness coefficient" for grnd water ice sublimation
22     
23! tracer indexes: these are initialized in initracer and should be 0 if the
24!                 corresponding tracer does not exist
25      ! dust
26      integer :: igcm_dustbin(nqmx) ! for dustbin 'dust' tracers
27      ! dust, special doubleq case
28      integer :: igcm_dust_mass   ! dust mass mixing ratio
29                                  !   (for transported dust)
30      integer :: igcm_dust_number ! dust number mixing ratio
31                                  !   (transported dust)
[358]32      integer :: igcm_ccn_mass   ! CCN mass mixing ratio
33      integer :: igcm_ccn_number ! CCN number mixing ratio
[38]34      integer :: igcm_dust_submicron ! submicron dust mixing ratio
35                                     !   (transported dust)
36      ! water
37      integer :: igcm_h2o_vap ! water vapour
38      integer :: igcm_h2o_ice ! water ice
39      ! chemistry:
40      integer :: igcm_co2
41      integer :: igcm_co
42      integer :: igcm_o
43      integer :: igcm_o1d
44      integer :: igcm_o2
45      integer :: igcm_o3
46      integer :: igcm_h
47      integer :: igcm_h2
48      integer :: igcm_oh
49      integer :: igcm_ho2
50      integer :: igcm_h2o2
51      integer :: igcm_n2
52      integer :: igcm_ar
[171]53      integer :: igcm_n
54      integer :: igcm_no
55      integer :: igcm_no2
56      integer :: igcm_n2d
[324]57      integer :: igcm_ch4
[171]58      ! Ions
59      integer :: igcm_co2plus
60      integer :: igcm_oplus
61      integer :: igcm_o2plus
62      integer :: igcm_coplus
63      integer :: igcm_cplus
64      integer :: igcm_nplus
65      integer :: igcm_noplus
66      integer :: igcm_n2plus
67      integer :: igcm_hplus
[635]68      integer :: igcm_hco2plus
[171]69      integer :: igcm_elec
[38]70      ! other tracers
71      integer :: igcm_ar_n2 ! for simulations using co2 +neutral gas
72
73
74! NB: to keep commons aligned: 
75!     split them in groups (reals, integers and characters)
76      COMMON/tracer/radius,rho_q,alpha_lift,alpha_devil,mmol,           &
[358]77     & varian,r3n_q,rho_dust,rho_ice,nuice_ref,nuice_sed,               &
[420]78     & ref_r0,ccn_factor,dryness
[38]79      COMMON/tracer2/                                                   &
[358]80     & igcm_dustbin,igcm_dust_mass,igcm_dust_number,                    &
81     & igcm_ccn_mass,igcm_ccn_number,igcm_dust_submicron,               &
[38]82     & igcm_h2o_vap,igcm_h2o_ice,igcm_co2,igcm_co,igcm_o,igcm_o1d,      &
83     & igcm_o2,igcm_o3,igcm_h,igcm_h2,igcm_oh,igcm_ho2,igcm_h2o2,       &
[171]84     & igcm_n2,igcm_ar,igcm_n,igcm_no,igcm_no2,igcm_n2d,                &
[324]85     & igcm_ch4,                                                        &
[171]86     & igcm_co2plus,igcm_oplus,igcm_o2plus,igcm_coplus,igcm_cplus,      &
87     & igcm_nplus,igcm_noplus,igcm_n2plus,igcm_hplus,igcm_elec,         &
[635]88     & igcm_hco2plus,igcm_ar_n2!,nbqchem,niqchem
[38]89      COMMON/tracer3/noms
90!-----------------------------------------------------------------------
Note: See TracBrowser for help on using the repository browser.