source: trunk/LMDZ.MARS/libf/phymars/nltedefs.h @ 414

Last change on this file since 414 was 414, checked in by aslmd, 13 years ago

LMDZ.MARS : NEW NLTE MODEL FROM GRANADA AMIGOS

23/11/11 == FGG + MALV

New parameterization of the NLTE 15 micron cooling. The old parameterization is kept as an option, including or not variable atomic oxygen concentration. A new flag is introduced in callphys.def, nltemodel, to select which parameterization wants to be used (new one, old one with variable [O], or old one with fixed [O], see below). Includes many new subroutines and commons in phymars. Some existing routines are also modified:

-physiq.F. Call to the new subroutine NLTE_leedat in first call. Call to nltecool modified to allow for variable atomic oxygen. Depending on the value of nltemodel, the new subroutine NLTEdlvr09_TCOOL is called instead of nltecool.

-inifis.F. Reading of nltemodel is added.

-callkeys.h Declaration of nltemodel is added.

The following lines should be added to callphys.def (ideally after setting callnlte):

# NLTE 15um scheme to use.
# 0-> Old scheme, static oxygen
# 1-> Old scheme, dynamic oxygen
# 2-> New scheme
nltemodel = 2

A new directory, NLTEDAT, has to be included in datagcm.

Improvements into NLTE NIR heating parameterization to take into account variability with O/CO2 ratio and SZA. A new subroutine, NIR_leedat.F, and a new common, NIRdata.h, are included. A new flag, nircorr, is added in callphys.def, to include or not these corrections. The following files are modified:

-nirco2abs.F: nq and pq are added in the arguments. The corrections factors are interpolated to the GCM grid and included in the clculation. A new subroutine, interpnir, is added at the end of the file.

-physiq.F: Call to NIR_leedat added at first call. Modified call to nirco2abs

-inifis: Reading new flag nircorr.

-callkeys.h: Declaration of nircorr.

The following lines have to be added to callphys.def (ideally after callnirco2):

# NIR NLTE correction for variable SZA and O/CO2?
# matters only if callnirco2=T
# 0-> no correction
# 1-> include correction
nircorr=1

A new data file, NIRcorrection_feb2011.dat, has to be included in datagcm.

Small changes to the molecular diffusion scheme to fix the number of species considered, to avoid problems when compiling with more than 15 tracers (for example, when CH4 is included). Modified subroutines: aeronomars/moldiff.F and aeronomars/moldiffcoeff.F
File size: 2.0 KB
Line 
1! NLTE parameters:
2
3        integer nl              ! actual # alt in NLTE module
4        parameter ( nl=20 )
5        !real deltaz
6        !parameter ( deltaz=5.0 )   ! Good option for Mars
7                                  ! Asegurarse que el conjunto de valores
8                                  ! (zmin,nl,deltaz) no se salen de los
9                                  ! límites del grid (zxmin,nz,deltazx)
10        integer nl2             ! = nl-2, needed for matrix inversion (mmh2)
11        parameter ( nl2=nl-2 ) 
12
13        integer nzy
14        parameter ( nzy = (nl-1)*4 + 1 )  ! Fine grid for mztud.f
15        !real deltazy
16        !parameter ( deltazy = deltaz*4.0 )
17
18!  Other NLTE parameters:
19        integer         nisot           ! number of isotopes considered
20        integer         nb              ! number of bands included
21        integer         nbmax_jt        ! number of bands for solar heating
22        parameter ( nisot=4, nb=41, nbmax_jt=47 )
23
24        integer         nhist                   ! # of temps in histogr.
25        parameter       ( nhist = 36 )          ! (get it from histograms!)
26
27        integer         nmax_freq               ! maximum number of points
28        parameter       ( nmax_freq = 40000 )   ! for overlap-freq integration
29                                                ! If dv (tcr_15um.drv) is very
30                                                ! small this must be increased.
31
32        integer         nbox_max
33        parameter       ( nbox_max = 70 )       ! max.# boxes in histogram
34
35!  Parameter for the composite vectors and matrixes of the compact-LU method
36        integer         nl5
37        parameter       ( nl5 = nl * 5 )
38
39! NLTE parameters that were needed by solar10 and are still in use
40        integer         ngroup                  ! normally we use 0 & 10 only
41        parameter       ( ngroup=10 )
42        integer         nisos                   ! #isotps para los vectores jt
43        parameter       ( nisos = 9 )
44        integer         nfile,nfile1            ! reminiscencia de cuando se
45        parameter       ( nfile=5, nfile1=4 )   ! leian las diff.hr desde
46                                                ! ficheros externos
47! NLTE Parameters for JT calculations
48        integer         nsublayers              ! subdivision para el fine grid
49        parameter       ( nsublayers = 50 )     ! called jj in the fot_alljt it
50                                                ! is usually 10 for vert. paths
51                                                ! or if you wish, for smallSZA,
52                                                ! and 50 for oblique (SZA>90)
Note: See TracBrowser for help on using the repository browser.