LMDzt-INCA (Interaction of chemistry and aerosol)  
 
Science :: Documentation :: Projets :: Publications :: Talk :: People

PATH :: home  / documentation / Reprobus structure / liste routines


analytic_1401

subroutine analytic_1401(temp,ptot,hnm,ck,vsed,
      $      k1,k2,k3,k4,k5,k6,k7,k8,k9,
      $      parthno3,parthno3sl,parth2o,parthcl,
      $      parthbr,
      $      surfarealiq,surfareanat,surfareaice,
      $      ilat,ndrop1,xlatt)
USE MODULES:
CHEM_REP
      $      only: nlat=>nlatitude, ! 1
      $      nlon0, ! klon = 6818
      $      nlon, ! 6818*(nivbas0=36)
      $      niv=>naltitude, ! klev = 50
      $      nivbas0=>nivbas, ! 36
      $      nivhau0=>nivhau, ! 12
      $      nbcon=>nbqj1, ! 43
      $      ilinit, ! nlon0*(nivhau0-1)+1
      $      preshau, ! 10
      $      presbas, ! 150
      $      latana, ! 50
      $      deltat ! 5 K
 Reprobus
CALL SUBROUTINES:
density_1401 Reprobus
hetero_1401 Reprobus

chem3d_rlong_1401

chem3d_rlong_1401( ilat, daynum, dt, alt, trajlon, trajlat, trajpre, trajtem, irapp, qj1, hc, vsed3d, parth2osed, parthno3sed, surfarealiq, surfareanat, surfareaice, sza3d, o3t, ndrop2D, prodHOBr0, ghgco2v)
USE MODULES:
CHEM_REP
      $      only: nlat=>nlatitude, ! 1
      $      nlon=>nlon0, ! klon = 6818
      $      niv=>naltitude, ! klev = 50
      $      nivbas0=>nivbas, ! 36
      $      nivhau0=>nivhau, ! 12
      $      nbcon=>nbqj1, ! 43
      $      ncm=>nbnas, ! 15
      $      niter, ! 5
      $      nphot, ! 41
      $      nz, ! 81
      $      nsza, ! 27
      $      nozo, ! 7
      $      nbmino3up, nbmaxo3up,! 0-79
      $      presbas, ! 200
      $      nbrflux, ! 19 !JK
      $      nbjflux, ! 14 !JK
      $      rflux, ! Rate Flux
      $      jflux,
      $      parthno3_out, parth2o_out,
      $      ok_emission, ok_plume, iplume
 Reprobus
phys_local_var_mod, only : ql=>ql_seriLmdz

CALL SUBROUTINES:
total_1401 Reprobus
zenith2_1401 Reprobus
phot_1401 Reprobus
analytic_1401 Reprobus
plurckpa Reprobus

chemmain_rlong_1401

subroutine chemmain_rlong_1401( vmr, lon, lat, t_seri, pplay, pphi, pphis, presnivs, sh, debut)
USE MODULES:
dimphy Lmdz
infotrac Lmdz
phys_cal_mod, only : xjourmod=>xjour, year_cur, jH_1jan, mth_cur, day_cur, hour Lmdz
ioipsl
CHEM_REP, ONLY :
      $      nbnas, ! 15
      $      nbqj1, ! 43
      $      nbqinterp, ! 56
      $      vmr_s=>nas, ! esp.non.transp
      $      pdel, ! delta pression
      $      daynum, ! jour de l annee
      $      dt=>pdt_rep, ! pas.temps30 min
      CDC$      ql=>ql_seri, ! eau liquide mmr
      $      surf_aeroliq, ! densite.surf.aero
      $      surf_aeronat, ! densite.surf.aero
      $      surf_aeroice, ! densite.surf.aero
      $      ptrop, ! pression tropopause
      $      ttrop, ! temperature tropopause (DC)
      $      ztrop, ! altitude tropopause
      $      colo3, ! colonne ozone
      $      prodHOBr3, ! product HOBr
      $      xjour, ! jour
      $      ghgco2v, ! co2
      $       nphot, ! 40
      $      nz, ! 81
      $      nsza, ! 27
      $      nozo, ! 7
      $      nlat=>nlatitude, ! 1
      $      nivbas, ! 36
      $      nbmino3up, nbmaxo3up, ! 0 - 79
      $      sedpas, ! 24
      $      gravit, ! 9.8
      $      ilat, irapp,
      $      ok_debut_integ, ! clef.deb.long.run
      $      toto,
      $      ndrop2D_rep,
      $      iageaire, ! 44
      $      latageaire, ! 15
      $      preshau, ! 10
      $      presbas, ! 150
      $      latana, ! 50
      $      ok_aga, xjour_aga,
      $      totolinear,
      $      refyear,
      $      Concatene_Output,
      $      qjH2SO4,
      $      iplume, ok_plume, ok_emission		 Reprobus
phys_local_var_mod, only : ql=>ql_seri Lmdz

CALL SUBROUTINES:
init_debutlongsim Reprobus
weight_fraction_h2so4 Reprobus
const_1401 Reprobus
FDTROPOPAUSE Reprobus
plumemis Reprobus
chem3d_rlong_1401 Reprobus
sedimentation_1401 Reprobus
conserv_1401 Reprobus
relax_rlong_1401 Reprobus
Concatene_Output Reprobus

Concatene_Output

SUBROUTINE Concatene_Output

conserv_1401

subroutine conserv_1401( qj1, hc )
c correction for mass conservation
USE MODULES:
mod_phys_lmdz_paraLmdz
CHEM_REP
      $      only: nlat=>nlatitude, ! 1
      $      nlon0, ! klon = 6818
      $      nlon, ! 6818*(nivbas0=36)
      $      niv=>naltitude, ! klev = 50
      $      nivbas0=>nivbas, ! 36
      $      nivhau0=>nivhau, ! 12
      $      nbcon=>nbqj1, ! 44
      $      ncm=>nbnas, ! 15
      $      ok_emission, ok_plume, iplume
 Reprobus

const_1401

subroutine const_1401
USE MODULES:
mod_phys_lmdz_paraLmdz
CHEM_REP
      $      only: nlat=>nlatitude, ! 1
      $      nlon=>nlon0, ! klon = 6818
      $      niv=>naltitude, ! klev = 50
      $      nivbas, ! 36
      $      nivhau, ! 12
      $      nphot, ! 41
      $      nz, ! 81
      $      nsza, ! 27
      $      nozo, ! 7
      $      nbmino3up, nbmaxo3up ! 0-79
 Reprobus

density_1401

SUBROUTINE DENSITY_1401(WS,WN,TEMP,DENS,ptot,xlatt)
c DENSITY OF TERNARY SOLUTION IN G/CM3
C WS ,WN ARE WT FRACTION,
C FITTED TO 0.05<WS+WN<0.70 WT FRACTION, BUT EXTRAPOLATES WELL
C 185 < T (K)
USE MODULES:
CHEM_REP, ONLY :
      $      nlat=>nlatitude, ! 1
      $      nlon0, ! klon = 6818
      $      nlon, ! 6818*(nivbas0=36)
      $      niv=>naltitude, ! klev = 50
      $      nivbas0=>nivbas, ! 36
      $      nivhau0=>nivhau, ! 12
      $      ILINIT, ! NLON0*(NIVHAU0-1)+1
      $      preshau, ! 10
      $      presbas, ! 150
      $      latana ! 50
 Reprobus

emiclfra

subroutine emiclfra(nlon,nlat,niv,p,frclx)

emirock

subroutine emirock (nlon,nlat,niv,t, p, alt, rocno2, roch2o, rocco, rochcl, ilatr,slat)

CALL SUBROUTINES:
READROCKET Reprobus

error

SUBROUTINE error(tab)

FDTROPOPAUSE

subroutine FDTROPOPAUSE (nx, ny, nz, pmid, temp,itrop)
USE MODULES:
mod_grid_phy_lmdzLmdz
mod_phys_lmdz_paraLmdz
CHEM_REP
      $      nlon=>nlon0, ! klon = 6818
      $      niv=>naltitude, ! klev = 50
 Reprobus

CALL SUBROUTINES:
FDTROPOPAUSE_GLO Reprobus

FDTROPOPAUSE_GLO

subroutine FDTROPOPAUSE_GLO (nx, ny, nz, pmid, temp,itrop)
! Purpose
-------
Calculates 2D fields of thermal tropopause pressures and tropopause level indices for given 3D fields of temperature and pressure.
Methods
-------
- Subroutine stattrop calculates the thermal tropopause pressure (Pa) for a 1D column (sounding) of pressures (Pa) and temperatures (K).

- The program tropo_test also corrects for non-sense occuring at high latitudes > 60N (60S).
Theese problems occur at high latitudes on the winter hemisphere because of the cold stratosphere and the consequently weak troposphere-stratosphere transition of the thermal lapse rate -dT/dz which is used to determine the thermal tropopause.
The problem is much less severe on the winter hemisphere.

Here a fix is implemented applying a factor which scales the tropopause pressures at latitudes LAT > 60N (60S) to the mean tropopause pressure at 60N (60S), whenever the mean tropopause pressure at LAT is lower than that at 60N (60S).
This idea was formulated in the ECHAM4/ routine xttphwmo by D. Nodorp, Ch. Land, B. Steil and R. Hein.

Programmed by Dominik Brunner, ETHZ, Switzerland V.01, 10 Aug 2000
Modified by Didier Hauglustaine, IPSL, for LMDZ/INCA, Oct 2000.
USE MODULES:
mod_grid_phy_lmdz, ONLY : nlon0=>klon_gloLmdz
CHEM_REP
      $      niv=>naltitude, ! klev = 50
Reprobus

CALL SUBROUTINES:
gr_fi_dynLmdz
stattropReprobus
gr_dyn_fiLmdz

hetero_1401

subroutine hetero_1401(aw,temp,ptot,hnm,ck,ws,
      $      hhcl,hhocl,hhobr,hclono2,hhbr,
      $      rice,nice,rnats,rnatl,nnats,nnatl,
      $      aliq,rmean,condliq,lnat,lice,
      $      k1,k2,k3,k4,k5,k6,k7,k8,k9,wt,xlatt)
C ROUTINE TO CALCULATE UPTAKE COEFFICIENTS (GAMMA VALUES).
C GAMMA VALUES ARE INDICATED BY VARIABLES WITH PREFIX 'G', FOR EXAMPLE GHOCLHCL IS THE GAMMA VALUE OF HOCL DUE TO REACTION WITH HCL IN THE DROPLETS.
C FROM THE GAMMA VALUES, SECOND ORDER RATE CONSTANTS ARE CALCULATED.
C THESE HAVE THE PREFIX 'R' AND HAVE UNITS CM3 MOLECULE-1 S-1. FOR EXAMPLE, THE LOSS OF CLNO3 AND HCL DUE TO THE HETEROGENEOUS REACTION
C CLNO3+HCL -> CL2+HNO3 IS D(CLNO3)/DT (UNITS MOLECULE CM-3 S-1) =-RCLNO3HCL.[CLNO3].[HCL],
cWHERE [CLNO3] AND [HCL] ARE THE ****TOTAL**** AMOUNTS OF THESE SPECIES IN UNITS MOLECULE CM-3.
USE MODULES:
CHEM_REP
      $      only: nlat=>nlatitude, ! 1
      $      nlon0, ! klon = 6818
      $      nlon, ! 6818*(nivbas0=36)
      $      niv=>naltitude, ! klev = 50
      $      nivbas0=>nivbas, ! 36
      $      nivhau0=>nivhau, ! 12
      $      nbcon=>nbqj1, ! 43
      $      ILINIT, ! NLON0*(NIVHAU0-1)+1)
      $      preshau, ! 10
      $      presbas, ! 150
      $      latana ! 50

Init_chem_rep_phys

SUBROUTINE Init_chem_rep_phys(klon0,klev0)
USE MODULES:
ioipsl

Init_chem_rep_trac

SUBROUTINE Init_chem_rep_trac(nbtr0)

Init_chem_rep_xjour

SUBROUTINE Init_chem_rep_xjour(xjour0)

init_debutlongsim

SUBROUTINE init_debutlongsim(pmid,qj1,xjour,debut)
USE MODULES:
phys_cal_mod, only : xjourmod=>xjour, year_cur, jH_1jan
CHEM_REP
      $      only: nlat=>nlatitude, ! 1
      $      nlon=>nlon0, ! klon = 6818
      $      nivbas, ! 36
      $      nbcon=>nbqj1, ! 44
      $      ncm=>nbnas, ! 15
      $      refyear,
      $      nF2D, ! 32
      $      nTOM, ! 3
      $      ok_Rtime,
      $      ok_Rtime2D,
      $      niv200mb,
      $      iplume, ok_emission, ok_plume
interp
ioipsl

CALL SUBROUTINES:
interp_F2D Reprobus

interp_0D

SUBROUTINE interp_0D(mode,spc,y,vals)
USE MODULES:
mod_phys_lmdz_para
netcdf

CALL SUBROUTINES:
YMD2y Reprobus

interp_AER

SUBROUTINE interp_AER(filem,file,filep,spc,year,pmid,vals,debut)
! Purpose: Interpolation from 2D model grid to LMDZ 3D grid
USE MODULES:
mod_grid_phy_lmdzLmdz
mod_phys_lmdz_paraLmdz
CHEM_REP, ONLY :
      $      nlon=>nlon0, ! klon = 6818
      $      niv=>naltitude, ! klev = 50
      $      valsform, valsfor, valsforp, nF2D
 Reprobus

CALL SUBROUTINES:
get_time Reprobus
read_M2D Reprobus
gr_fi_dynLmdz
gr_dyn_fiLmdz

interp_F2D

SUBROUTINE interp_F2D(filem,file,filep,spc,year,pmid,vals,debut)
! Purpose: Interpolation from 2D model grid to LMDZ 3D grid
USE MODULES:
mod_phys_lmdz_paraLmdz
mod_grid_phy_lmdzLmdz
CHEM_REP, ONLY :
      $      nlon=>nlon0, ! klon = 6818
      $      niv=>naltitude, ! klev = 50
      $      valsform, valsfor, valsforp, nF2D
 Reprobus

CALL SUBROUTINES:
get_time Reprobus
read_M2D Reprobus
gr_fi_dynLmdz
gr_dyn_fiLmdz

interp_TOM

SUBROUTINE interp_TOM(filem,file,filep,spc,year,pmid,vals,debut)
! Purpose: Vertical interpolation from 3D TOM grid to LMDZ 3D grid
USE MODULES:
mod_phys_lmdz_paraLmdz
mod_grid_phy_lmdzLmdz
CHEM_REP, ONLY :
      $      nlon=>nlon0, ! klon = 6818
      $      niv=>naltitude, ! klev = 50
      $      valsform, valsfor, valsforp, nF2D
      $      nTOM, ! nTOM=3
      $      nlevt=>nlevTOM, ! nlevTOM=17
      $      nlont=>nlonTOM, ! nlonTOM=97
      $      nlatt=>nlatTOM ! nlatTOM=73
 Reprobus

CALL SUBROUTINES:
get_time Reprobus
read_TOM Reprobus
gr_fi_dynLmdz
gr_dyn_fiLmdz

jno

subroutine jno(cc_no, hnm, nlon, nivbas, o3t, sza, pm, tjno)
! parametrisation de la photodissociation de no ! d'apres minschwaner and siskind, a new calculation of nitric oxide photolysis in the stratosphere, mesosphere, and lower thermosphere, j. geophys. res., 98, 20401-20412, 1993

minmaxqfi_chimie

SUBROUTINE minmaxqfi_chimie(iqui,zq,qmin,qmax,comment)
USE MODULES:
dimphy Lmdz

phot_1401

subroutine phot_1401( pm, ilat, alt, sza, o3t, tj, cc_no, hnm )
!calculates photodissociation frequencies (J coefficients) [sec-1]
USE MODULES:
CHEM_REP
      $      only: nlat=>nlatitude, ! 1
      $      nlon0, ! klon = 6818
      $      nlon, ! 6818*(nivbas0=36)
      $      niv=>naltitude, ! klev = 50
      $      nivbas0=>nivbas, ! 36
      $      nivhau0=>nivhau, ! 12
      $      nbcon=>nbqj1, ! 43
      $      ncm=>nbnas, ! 15
      $      nphot, ! 40
      $      nz, ! 81
      $      nsza, ! 27
      $      nozo, ! 7
      $      nbmino3up, nbmaxo3up ! 0-79
 Reprobus

CALL SUBROUTINES:
jno Reprobus

plumemis

subroutine plumemis (t, p, alt, qj1, xlat, xlon)
USE MODULES:
mod_phys_lmdz_para Lmdz
mod_grid_phy_lmdz Lmdz
CHEM_REP, ONLY :
      $      nbcon=>nbqj1,
      $      nlat=>nlatitude, !1
      $      nlon=>nlon0,
      $      niv=>naltitude,
      $      ok_plume
 Reprobus

CALL SUBROUTINES:
gr_fi_dynLmdz
emirock Reprobus
emiclfra Reprobus
gr_dyn_fiLmdz

plurckpa

subroutine plurckpa(p,sza,reffclo3,dilutim,deactim)
c calculation of rocket plume chemical parameters
USE MODULES:
CHEM_REP, ONLY : nlon! 6818*(nivbas0=36) Reprobus

POSITIO

Subroutine POSITIO(XC,N,X,JX,IER)

RAD_INTERACTIF

SUBROUTINE RAD_INTERACTIF(POZON,iof)
USE MODULES:
dimphy Lmdz
CHEM_REP, ONLY : nas, ok_Rtime2D, ok_o3interact, RCH42D, RN2O2D, RCFC112D, RCFC122D, ok_debut, ndimozon Reprobus
phys_local_var_mod, ONLY : tr_seri Lmdz

read_M2D

SUBROUTINE read_M2D(file,spc,vals,lat,pre)
! Purpose: 2D fields reading (no longitude)
USE MODULES:
netcdf

READROCKET

Subroutine READROCKET(nlon,nlat,niv,alt,fuel,fuel1)
CALL SUBROUTINES:
POSITIO Reprobus

read_TOM

SUBROUTINE read_TOM(file,spc,vals,lon,lat,pre)
! Purpose: 3D model files reading
USE MODULES:
netcdf
CHEM_REP, ONLY :
      $      nsp=>nTOM, ! nTOM=3
      $      nlev=>nlevTOM, ! nlevTOM=17
      $      nlon=>nlonTOM, ! nlonTOM=97
      $      nlat=>nlatTOM ! nlatTOM=73
 Reprobus

relax_rlong_1401

SUBROUTINE relax_rlong_1401(pmid,qj1,hc,xjour,debut)
! Some comments:
! * The list of qj1 (REPROBUS array) index of the species available in forcing file (read using read_forcing(..'F2D'..) is iF2D.
! * qj1(:,:,:,k) is updated if lrelax(k)=.true. So, if this species(k) was available in the forcing file, then qj1 is updated to forcing field.
! If not, it is set to zero (or the initial value of Qdyn if you changed it)
USE MODULES:
phys_cal_mod, only : xjourmod=>xjour, year_cur,jH_1janLmdz
CHEM_REP, ONLY :
      $      nlat=>nlatitude, ! 1
      $      nlon=>nlon0, ! klon = 6818
      $      nivbas, ! 36
      $      nbcon=>nbqj1, ! 44
      $      ncm=>nbnas, ! 15
      $      refyear,
      $      nF2D, ! 32
      $      nTOM, ! 3
      $      ok_Rtime,
      $      ok_Rtime2D,
      $      iplume, ok_emission, ok_plume
 Reprobus
interp Reprobus
ioipsl, ONLY: ioget_year_len
CALL SUBROUTINES:
interp_F2D Reprobus
interp_TOM Reprobus
interp_0D Reprobus

total_1401

subroutine total_1401( trajpre, irapp, ilat, hc, hnm, o3t )
! calcul de la colonne d ozone [cm-2]
USE MODULES:
CHEM_REP, ONLY :
      $      nlat=>nlatitude, ! 1
      $      nlon0, ! klon = 6818
      $      nlon, ! 6818*(nivbas0=36)
      $      niv=>naltitude, ! klev = 50
      $      nivbas0=>nivbas, ! 36
      $      nivhau0=>nivhau, ! 12
      $      nbcon=>nbqj1, ! 43
      $      ncm=>nbnas ! 15
 Reprobus

sedimentation_1401

subroutine sedimentation_1401( trajpre, trajtem, irapp, qj1, vsed3d, parth2osed, parthno3sed, dt,itrop,lat )
!calcule les quantites de h2o et hno3 sedimentees
USE MODULES:
CHEM_REP, ONLY :
      $      nlat=>nlatitude, ! 1
      $      nlon=>nlon0, ! klon = 6818
      $      niv=>naltitude, ! klev = 50
      $      nivbas, ! 36
      $      nivhau, ! 12
      $      nbcon=>nbqj1, ! 43
      $      ncm=>nbnas, ! 15
      $      lolani, ! (nlon0=6818)*nlat*(niv=50)
      $      loni, ! (nlon0=6818)*(nivbas=36)
      $      presbas, ! 150.
      $      latana ! 50.
 Reprobus

stattrop

subroutine stattrop(pfull, tfull, nlev, ptropd, itropd)
c programmed by Dominik Brunner V1.0 Aug 2000
c built upon routine stattrop by Peter van Velthoven, KNMI, The Netherlands and on the ECHAM4/ routine xttphwmo by Thomas Reichle, Christine Land, B. Steil and R. Hein, DLR
c purpose
c -------
c stattrop computes the pressure (Pa) at the thermal (static) tropopause (TP) for a 1D column (sounding) of pressures and temperatures following the definition of the height of the tropopause as postulated by WMO (1957).
c ATTENTION: In the current formulation of the program the first level (index 1) must be at the top of the atmosphere and the last level (index nlev) at the surface
c interface
c ---------
c call stattrop(pfull, tfull, nlev, ptropd, itropd)
c - Input
c nlev : number of vertical levels
c pfull: pressure in 1D column at nlev full levels (layers)
c tfull: temperature in 1D column at nlev full levels
c - Output
c ptropd: height of the tropopause in Pa
c itropd: index of layer containing the tropopause
c c methods
c -------
c - Lapse rate gamma = -dT/dz
c Using the hydrostatic approximation
c
c dz = -dp/(g*rho) = -dp/p * R/g * T = -dlnp * R/g * T
c
c we get -dT/dz = dT/T * g/R *1/dlnp = dlnT/dlnp
c
c - Variables are assumed to vary linearly in log(pressure)
c - The tropopause is the lowest level above 450 hPa fullfilling the WMO criterium. If ptropd is less than 85 hPa it is set to 85 hPa.
c If no tropopause is found ptropd is set to -999.
c
c references
c ----------
c - WMO (1992): International meteorological vocabulary, Genf, 784pp.:
c
c 1. The first tropopause is defined as the lowest level at which the lapse rate decreases to 2 deg C per kilometer or less, provided also the average lapse rate between this level and all higher levels within 2 kilometers does not exceed 2 deg C
c
c - Randel WJ, Wu F, Gaffen DJ, Interannual variability of the tropical tropopause derived from radiosonde data and NCEP reanalyses,
c JOURNAL OF GEOPHYSICAL RESEARCH, 105, 15509-15523, 2000.
c
c The following webpage clearifies the calculation of the tropopause in the NCEP reanalysis: http://dss.ucar.edu/pub/reanalysis/FAQ.html
c
c For comparison NCEP reanalysis tropopause pressures can be obtained
c from http://www.cdc.noaa.gov/cdc/reanalysis/reanalysis.shtml

test_leap

SUBROUTINE test_leap(y,calendar,ndays,months)
! Purpose: This function tests if 'y' is a leap year or not.

get_time

SUBROUTINE get_time(file,year)
USE MODULES:
netcdf

weight_fraction_h2so4

SUBROUTINE weight_fraction_h2so4(trajpre, trajtem, qj1, ndrop1, debut)
USE MODULES:
phys_cal_mod, ONLY : year_cur, xjour, jH_1jan Lmdz
netcdf, ONLY : ioget_year_len
CHEM_REP, ONLY : nlon, nbcon=>nbqj1, ilinit, nlon0, nlat=>nlatitude, niv=>naltitude, nivbas0=>nivbas Reprobus
interp Reprobus
CALL SUBROUTINES:
interp_AER Reprobus

YMD2y

SUBROUTINE YMD2y(YMD,calendar,year,hms)
! Purpose: Converts a date (YYYYMMDD format) + a hours hms (hhmmss, optional) into a year (double format). Works for all calendars implemented in test_leap.

zenith2_1401

subroutine zenith2_1401(trajlon, trajlat, irapp, ilat, day, sza )
! calcul de l'angle zenithal solaire
USE MODULES:
CHEM_REP, ONLY :
      $      nlat=>nlatitude, ! 1
      $      nlon0, ! klon = 6818
      $      nlon, ! 6818*(nivbas0=36)
      $      niv=>naltitude, ! klev = 50
      $      nivbas0=>nivbas, ! 36
      $      nivhau0=>nivhau, ! 12
      $      nbcon=>nbqj1, ! 43
      $      ncm=>nbnas ! 15
 Reprobus
 

 

 
           © Laboratoire des Sciences du Climat et de l'Environnement