Index: trunk/LMDZ.PLUTO/deftank/README.md
===================================================================
--- trunk/LMDZ.PLUTO/deftank/README.md	(revision 3772)
+++ trunk/LMDZ.PLUTO/deftank/README.md	(revision 3845)
@@ -1,7 +1,7 @@
 In this folder, you will find examples .def files for the Pluto PCM.
 
-In case you use start*nc files from the old Pluto model, you can use `restart.old_pluto_to_new.def` to convert them (new model is in AU while old model in million of km).
+In case you use start*nc files from the old Pluto model, you can use `restart.old_pluto_to_new.def` (in newstart/) to convert them (new model is in AU while old model in million of km).
 
-In the 1D/ folder, callphys.def for various examples including hazes with profiles from Zhang et al 2017, a haze-generation routine (hazecloud) with a fixed particle radius, and without hazes.
+In the run_1D/ folder, callphys.def for various examples including hazes with profiles from Zhang et al 2017, a haze-generation routine (hazecloud) with a fixed particle radius, and without hazes.
 
 In the xios/ folder, you will find the necessary .xml files for xios.
@@ -9,5 +9,5 @@
 In the dynamico/ folder, you will find the necessary .xml and .def files for dynamico.
 
-In the 3D/ folder, TBD
+In the run_3D/ folder, TBD
 
 - newstart:
Index: trunk/LMDZ.PLUTO/deftank/newstart/restart.old_pluto_to_new.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/newstart/restart.old_pluto_to_new.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/newstart/restart.old_pluto_to_new.def	(revision 3845)
@@ -0,0 +1,8 @@
+1
+n
+periastr
+29.658309835823086
+apoastr
+49.30504669275349
+
+
Index: trunk/LMDZ.PLUTO/deftank/run_1d/README
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_1d/README	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_1d/README	(revision 3845)
@@ -0,0 +1,2 @@
+
+ref_debatz_2025_07 : 1D_Gas+Haze_rp-5nm_rm-20nm_rho-800kgm-3
Index: trunk/LMDZ.PLUTO/deftank/run_1d/callphys.def.hazecloud_10
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_1d/callphys.def.hazecloud_10	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_1d/callphys.def.hazecloud_10	(revision 3845)
@@ -0,0 +1,187 @@
+## Datagcm path
+datadir=datadir
+## Orbit / general options
+## ~~~~~~~~~~~~~~~~~~~~~~~
+# Run with or without tracer transport ?
+tracer    = .true.
+# Diurnal cycle ?  if diurnal=false, diurnally averaged solar heating
+diurnal   = .false.
+# Seasonal cycle ? if season=false, Ls stays constant, to value set in "start"
+season    = .false.
+noseason_day = 1550
+# Fixed zenith angle in 1D
+global1d = .true.
+# Run 1D or 3D version (false) / fast version volatile transport (true)
+fast    = .false.
+# Write some more output on the screen ?
+lwrite    = .false.
+# Save statistics in file "stats.nc" ?
+callstats = .false.
+# Test energy conservation of model physics ?
+enertest  = .false.
+# run with triton orbit ?
+triton = .false.
+# conservation n2
+conservn2= .false.
+# check if cpp and mugaz from start.nc are consistent with values computed by comp_cpp_mugaz with gases.def
+check_cpp_match = .true.
+check_physics_inputs = .true.
+check_physics_outputs = .true.
+check_ps_min = 0.01
+
+## Radiative transfer options
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+# call radiative transfer?
+callrad    = .true.
+# the rad. transfer is computed every "iradia" physical timestep
+iradia     = 1
+# call multilayer correlated-k radiative transfer ?
+corrk      = .true.
+oldplutocorrk=.true.
+ch4fix = .true.
+vmrch4fix = 0.5
+nlte=.true.
+strobel=.false.
+cooling=.true.
+pref=.05
+alpha_top = 1e-11
+deltap = 0.7
+# folder in which correlated-k data is stored ?
+corrkdir   = pluton_ch4var_COfix
+# strictboundcorrk
+strictboundcorrk = .false.
+tplanckmin = 20
+dtplanck=0.1
+# call visible gaseous absorption in radiative transfer ?
+callgasvis = .true.
+# call continuum in radiative transfer ?
+continuum = .false.
+# Include Rayleigh scattering in the visible ?
+rayleigh   = .false.
+# Test physics timescale in 1D ?
+testradtimes = .false.
+# Output spectral OLR in 1D/3D?
+specOLR    = .true.
+# Output global radiative balance in file 'rad_bal.out' - slow for 1D!!
+meanOLR    = .false.
+# Variable gas species: Radiatively active ?
+varactive  = .true.
+# Variable gas species: Fixed vertical distribution ?
+#   (not to be used in time integration mode)
+varfixed   = .false.
+# Use fixed vertical profile, 1 step, no iteration ?
+kastprof     = .false.
+
+# Star type (inherited from generic model)
+startype = 1
+Fat1AU = 1366.0
+
+
+## Other physics options
+## ~~~~~~~~~~~~~~~~~~~~~
+# call turbulent vertical diffusion ?
+calldifv = .true.
+kmixmin=0.001
+vertdiff=.true.
+# use turbdiff instead of vdifc ? (if oldplutovdifc = .false.)
+# UseTurbDiff = .true.
+oldplutovdifc=.true.
+# call convective adjustment ?
+calladj  = .true.
+# call thermal conduction in the soil ?
+callsoil = .true.
+lay1_soil = 2e-4
+# geothermal flux at the bottom layer
+fluxgeo = 0.
+# call thermal conduction
+callconduct=.false.
+# call molecular viscosity 3D only
+callmolvis=.false.
+# call molecular diffusion
+callmoldiff=.false.
+# flux qui arrive au-dessus de latmosphere
+phitop=0.
+
+
+## Tracer and aerosol options
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+# atm mass update due to tracer evaporation/condensation?
+mass_redistrib = .false.
+# Gravitational sedimentation of tracers ?
+sedimentation = .true.
+oldsedimpluto= .true.
+naerkind = 1
+## Condensation
+## ~~~~~~~~~~~
+generic_condensation = .false.
+# call N2 condensation ?
+n2cond       = .false.
+# N2 condensation subtimestep ?
+nbcond = 20
+# Set initial temperature profile to 1 K above N2 condensation everywhere?
+nearn2cond   = .false.
+# Compute methane cycle ?
+methane = .true.
+condmetsurf = .false.
+# Compute CO cycle ?
+carbox = .true.
+condcosurf=.false.
+# Compute methane clouds ?
+metcloud = .false.
+# Compute CO clouds ?
+monoxcloud = .false.
+# Haze production ?
+haze = .true.
+hazeconservch4    = .true.
+# number mixing ratio of haze particles for ch4 clouds
+Nmix_ch4      = 100000.
+# number mixing ratio of haze particles for co clouds
+Nmix_co      = 100000.
+# Radiatively active haze ?
+aerohaze=.true.
+hazeprop_file=optprop_rannou_fractal_r010nm_N1_1e4_d2.dat
+# fractal particle ?
+fractal=.true.
+nb_monomer=10
+# haze particule mean radius
+rad_haze = 10e-9
+# Fixed haze radius
+haze_radproffix = .false.
+hazerad_file = rad_haze.txt
+# Fixed haze distribution? # t==  haze_prof (hazemmr_file) else  hazecloud()
+haze_proffix = .false.
+# hazemmr_file = hazemmr.txt
+hazedens_file = dens_haze.txt
+
+## Surface properties
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+mode_n2=0     ! 0=fixed, 1=change w thickness, 2=special for SP
+deltab=0.    ! if mod_n2=1, difference of albedo depending of thickness
+metlateq=35  ! latitude separating dark/bright methane
+tholateq=25  ! latitude separating eq/pole tholins
+mode_ch4=1
+specalb=.false.  ! Specific treatment for tholins (diff values)
+
+alb_n2b=0.7
+alb_n2a=0.7
+alb_ch4=0.5
+alb_ch4_eq=0.5
+alb_ch4_s=0.5
+alb_co=0.4
+alb_tho=0.2
+alb_tho_eq=0.2
+
+emis_n2b=0.8
+emis_n2a=0.8
+emis_ch4=0.8
+emis_co=0.8
+emis_tho=0.8
+emis_tho_eq=0.8
+
+# Change TI with time
+changeti=.false.
+changetid=.false.
+ITCH4=800.
+ITN2=800
+ITH2O=2000.
+
Index: trunk/LMDZ.PLUTO/deftank/run_1d/callphys.def.hazecloud_lora
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_1d/callphys.def.hazecloud_lora	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_1d/callphys.def.hazecloud_lora	(revision 3845)
@@ -0,0 +1,188 @@
+## Datagcm path
+datadir=datadir
+## Orbit / general options
+## ~~~~~~~~~~~~~~~~~~~~~~~
+# Run with or without tracer transport ?
+tracer    = .true.
+# Diurnal cycle ?  if diurnal=false, diurnally averaged solar heating
+diurnal   = .false.
+# Seasonal cycle ? if season=false, Ls stays constant, to value set in "start"
+season    = .false.
+noseason_day = 1550
+# Fixed zenith angle in 1D
+global1d = .true.
+# Run 1D or 3D version (false) / fast version volatile transport (true)
+fast    = .false.
+# Write some more output on the screen ?
+lwrite    = .false.
+# Save statistics in file "stats.nc" ?
+callstats = .false.
+# Test energy conservation of model physics ?
+enertest  = .false.
+# run with triton orbit ?
+triton = .false.
+# conservation n2
+conservn2= .false.
+# check if cpp and mugaz from start.nc are consistent with values computed by comp_cpp_mugaz with gases.def
+check_cpp_match = .true.
+check_physics_inputs = .true.
+check_physics_outputs = .true.
+check_ps_min = 0.01
+
+## Radiative transfer options
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+# call radiative transfer?
+callrad    = .true.
+# the rad. transfer is computed every "iradia" physical timestep
+iradia     = 1
+# call multilayer correlated-k radiative transfer ?
+corrk      = .true.
+oldplutocorrk=.true.
+ch4fix = .true.
+vmrch4fix = 0.5
+nlte=.true.
+strobel=.false.
+cooling=.true.
+pref=.05
+alpha_top = 1e-11
+deltap = 0.7
+# folder in which correlated-k data is stored ?
+corrkdir   = pluton_ch4var_COfix
+# strictboundcorrk
+strictboundcorrk = .false.
+tplanckmin = 20
+dtplanck=0.1
+# call visible gaseous absorption in radiative transfer ?
+callgasvis = .true.
+# call continuum in radiative transfer ?
+continuum = .false.
+# Include Rayleigh scattering in the visible ?
+rayleigh   = .false.
+# Test physics timescale in 1D ?
+testradtimes = .false.
+# Output spectral OLR in 1D/3D?
+specOLR    = .true.
+# Output global radiative balance in file 'rad_bal.out' - slow for 1D!!
+meanOLR    = .false.
+# Variable gas species: Radiatively active ?
+varactive  = .true.
+# Variable gas species: Fixed vertical distribution ?
+#   (not to be used in time integration mode)
+varfixed   = .false.
+# Use fixed vertical profile, 1 step, no iteration ?
+kastprof     = .false.
+
+# Star type (inherited from generic model)
+startype = 1
+Fat1AU = 1366.0
+
+
+## Other physics options
+## ~~~~~~~~~~~~~~~~~~~~~
+# call turbulent vertical diffusion ?
+calldifv = .true.
+kmixmin=0.001
+vertdiff=.true.
+# use turbdiff instead of vdifc ? (if oldplutovdifc = .false.)
+# UseTurbDiff = .true.
+oldplutovdifc=.true.
+# call convective adjustment ?
+calladj  = .true.
+# call thermal conduction in the soil ?
+callsoil = .true.
+lay1_soil = 2e-4
+# geothermal flux at the bottom layer
+fluxgeo = 0.
+# call thermal conduction
+callconduct=.false.
+# call molecular viscosity 3D only
+callmolvis=.false.
+# call molecular diffusion
+callmoldiff=.false.
+# flux qui arrive au-dessus de latmosphere
+phitop=0.
+
+
+## Tracer and aerosol options
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+# atm mass update due to tracer evaporation/condensation?
+mass_redistrib = .false.
+# Gravitational sedimentation of tracers ?
+sedimentation = .true.
+oldsedimpluto= .true.
+naerkind = 1
+## Condensation
+## ~~~~~~~~~~~
+generic_condensation = .false.
+# call N2 condensation ?
+n2cond       = .false.
+# N2 condensation subtimestep ?
+nbcond = 20
+# Set initial temperature profile to 1 K above N2 condensation everywhere?
+nearn2cond   = .false.
+# Compute methane cycle ?
+methane = .true.
+condmetsurf = .false.
+# Compute CO cycle ?
+carbox = .true.
+condcosurf=.false.
+# Compute methane clouds ?
+metcloud = .false.
+# Compute CO clouds ?
+monoxcloud = .false.
+# Haze production ?
+haze = .true.
+hazeconservch4    = .true.
+# number mixing ratio of haze particles for ch4 clouds
+Nmix_ch4      = 100000.
+# number mixing ratio of haze particles for co clouds
+Nmix_co      = 100000.
+# Radiatively active haze ?
+aerohaze=.true.
+# hazeprop_file=optprop_rannou_fractal_r010nm_N1_1e4_d2.dat
+hazeprop_file=optprop_lora_fractal_PH_r010nm_N1_5e5_d2.dat
+# fractal particle ?
+fractal=.true.
+nb_monomer=20
+# haze particule mean radius
+rad_haze = 20e-9
+# Fixed haze radius
+haze_radproffix = .false.
+hazerad_file = rad_haze.txt
+# Fixed haze distribution? # t==  haze_prof (hazemmr_file) else  hazecloud()
+haze_proffix = .false.
+# hazemmr_file = hazemmr.txt
+hazedens_file = dens_haze.txt
+
+## Surface properties
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+mode_n2=0     ! 0=fixed, 1=change w thickness, 2=special for SP
+deltab=0.    ! if mod_n2=1, difference of albedo depending of thickness
+metlateq=35  ! latitude separating dark/bright methane
+tholateq=25  ! latitude separating eq/pole tholins
+mode_ch4=1
+specalb=.false.  ! Specific treatment for tholins (diff values)
+
+alb_n2b=0.7
+alb_n2a=0.7
+alb_ch4=0.5
+alb_ch4_eq=0.5
+alb_ch4_s=0.5
+alb_co=0.4
+alb_tho=0.2
+alb_tho_eq=0.2
+
+emis_n2b=0.8
+emis_n2a=0.8
+emis_ch4=0.8
+emis_co=0.8
+emis_tho=0.8
+emis_tho_eq=0.8
+
+# Change TI with time
+changeti=.false.
+changetid=.false.
+ITCH4=800.
+ITN2=800
+ITH2O=2000.
+
Index: trunk/LMDZ.PLUTO/deftank/run_1d/callphys.def.nohaze
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_1d/callphys.def.nohaze	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_1d/callphys.def.nohaze	(revision 3845)
@@ -0,0 +1,187 @@
+## Datagcm path
+datadir=datadir
+## Orbit / general options
+## ~~~~~~~~~~~~~~~~~~~~~~~
+# Run with or without tracer transport ?
+tracer    = .true.
+# Diurnal cycle ?  if diurnal=false, diurnally averaged solar heating
+diurnal   = .false.
+# Seasonal cycle ? if season=false, Ls stays constant, to value set in "start"
+season    = .false.
+noseason_day = 1550
+# Fixed zenith angle in 1D
+global1d = .true.
+# Run 1D or 3D version (false) / fast version volatile transport (true)
+fast    = .false.
+# Write some more output on the screen ?
+lwrite    = .false.
+# Save statistics in file "stats.nc" ?
+callstats = .false.
+# Test energy conservation of model physics ?
+enertest  = .false.
+# run with triton orbit ?
+triton = .false.
+# conservation n2
+conservn2= .false.
+# check if cpp and mugaz from start.nc are consistent with values computed by comp_cpp_mugaz with gases.def
+check_cpp_match = .true.
+check_physics_inputs = .true.
+check_physics_outputs = .true.
+check_ps_min = 0.01
+
+## Radiative transfer options
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+# call radiative transfer?
+callrad    = .true.
+# the rad. transfer is computed every "iradia" physical timestep
+iradia     = 1
+# call multilayer correlated-k radiative transfer ?
+corrk      = .true.
+oldplutocorrk=.true.
+ch4fix = .true.
+vmrch4fix = 0.5
+nlte=.true.
+strobel=.false.
+cooling=.true.
+pref=.05
+alpha_top = 1e-11
+deltap = 0.7
+# folder in which correlated-k data is stored ?
+corrkdir   = pluton_ch4var_COfix
+# strictboundcorrk
+strictboundcorrk = .false.
+tplanckmin = 30
+dtplanck=0.1
+# call visible gaseous absorption in radiative transfer ?
+callgasvis = .true.
+# call continuum in radiative transfer ?
+continuum = .false.
+# Include Rayleigh scattering in the visible ?
+rayleigh   = .false.
+# Test physics timescale in 1D ?
+testradtimes = .false.
+# Output spectral OLR in 1D/3D?
+specOLR    = .true.
+# Output global radiative balance in file 'rad_bal.out' - slow for 1D!!
+meanOLR    = .false.
+# Variable gas species: Radiatively active ?
+varactive  = .true.
+# Variable gas species: Fixed vertical distribution ?
+#   (not to be used in time integration mode)
+varfixed   = .false.
+# Use fixed vertical profile, 1 step, no iteration ?
+kastprof     = .false.
+
+# Star type (inherited from generic model)
+startype = 1
+Fat1AU = 1366.0
+
+
+## Other physics options
+## ~~~~~~~~~~~~~~~~~~~~~
+# call turbulent vertical diffusion ?
+calldifv = .true.
+kmixmin=0.001
+vertdiff=.true.
+# use turbdiff instead of vdifc ? (if oldplutovdifc = .false.)
+UseTurbDiff = .true.
+oldplutovdifc=.true.
+# call convective adjustment ?
+calladj  = .true.
+# call thermal conduction in the soil ?
+callsoil = .true.
+lay1_soil = 2e-4
+# geothermal flux at the bottom layer
+fluxgeo = 0.
+# call thermal conduction
+callconduct=.false.
+# call molecular viscosity 3D only
+callmolvis=.false.
+# call molecular diffusion
+callmoldiff=.false.
+# flux qui arrive au-dessus de latmosphere
+phitop=0.
+
+
+## Tracer and aerosol options
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+# atm mass update due to tracer evaporation/condensation?
+mass_redistrib = .false.
+# Gravitational sedimentation of tracers ?
+sedimentation = .true.
+oldsedimpluto= .true.
+naerkind = 1
+## Condensation
+## ~~~~~~~~~~~
+generic_condensation = .false.
+# call N2 condensation ?
+n2cond       = .false.
+# N2 condensation subtimestep ?
+nbcond = 20
+# Set initial temperature profile to 1 K above N2 condensation everywhere?
+nearn2cond   = .false.
+# Compute methane cycle ?
+methane = .true.
+condmetsurf = .false.
+# Compute CO cycle ?
+carbox = .true.
+condcosurf=.false.
+# Compute methane clouds ?
+metcloud = .false.
+# Compute CO clouds ?
+monoxcloud = .false.
+# Haze production ?
+haze = .true.
+hazeconservch4    = .true.
+# haze particule mean radius
+rad_haze = 10e-9
+# number mixing ratio of haze particles for ch4 clouds
+Nmix_ch4      = 100000.
+# number mixing ratio of haze particles for co clouds
+Nmix_co      = 100000.
+# Radiatively active haze ?
+aerohaze=.false.
+hazeprop_file=optprop_rannou_r2-200nm_nu003.dat
+# fractal particle ?
+fractal=.false.
+nb_monomer=10
+# Fixed haze radius
+haze_radproffix = .false.
+hazerad_file = rad_haze.txt
+# Fixed haze distribution? # t==  haze_prof (hazemmr_file) else  hazecloud()
+haze_proffix = .false.
+hazemmr_file = hazemmr.txt
+hazedens_file = dens_haze.txt
+
+## Surface properties
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+mode_n2=0     ! 0=fixed, 1=change w thickness, 2=special for SP
+deltab=0.    ! if mod_n2=1, difference of albedo depending of thickness
+metlateq=35  ! latitude separating dark/bright methane
+tholateq=25  ! latitude separating eq/pole tholins
+mode_ch4=1
+specalb=.false.  ! Specific treatment for tholins (diff values)
+
+alb_n2b=0.7
+alb_n2a=0.7
+alb_ch4=0.5
+alb_ch4_eq=0.5
+alb_ch4_s=0.5
+alb_co=0.4
+alb_tho=0.2
+alb_tho_eq=0.2
+
+emis_n2b=0.8
+emis_n2a=0.8
+emis_ch4=0.8
+emis_co=0.8
+emis_tho=0.8
+emis_tho_eq=0.8
+
+# Change TI with time
+changeti=.false.
+changetid=.false.
+ITCH4=800.
+ITN2=800
+ITH2O=2000.
+
Index: trunk/LMDZ.PLUTO/deftank/run_1d/callphys.def.zhang
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_1d/callphys.def.zhang	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_1d/callphys.def.zhang	(revision 3845)
@@ -0,0 +1,187 @@
+## Datagcm path
+datadir=datadir
+## Orbit / general options
+## ~~~~~~~~~~~~~~~~~~~~~~~
+# Run with or without tracer transport ?
+tracer    = .true.
+# Diurnal cycle ?  if diurnal=false, diurnally averaged solar heating
+diurnal   = .false.
+# Seasonal cycle ? if season=false, Ls stays constant, to value set in "start"
+season    = .false.
+noseason_day = 1550
+# Fixed zenith angle in 1D
+global1d = .true.
+# Run 1D or 3D version (false) / fast version volatile transport (true)
+fast    = .false.
+# Write some more output on the screen ?
+lwrite    = .false.
+# Save statistics in file "stats.nc" ?
+callstats = .false.
+# Test energy conservation of model physics ?
+enertest  = .false.
+# run with triton orbit ?
+triton = .false.
+# conservation n2
+conservn2= .false.
+# check if cpp and mugaz from start.nc are consistent with values computed by comp_cpp_mugaz with gases.def
+check_cpp_match = .true.
+check_physics_inputs = .true.
+check_physics_outputs = .true.
+check_ps_min = 0.01
+
+## Radiative transfer options
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+# call radiative transfer?
+callrad    = .true.
+# the rad. transfer is computed every "iradia" physical timestep
+iradia     = 1
+# call multilayer correlated-k radiative transfer ?
+corrk      = .true.
+oldplutocorrk=.true.
+ch4fix = .true.
+vmrch4fix = 0.5
+nlte=.true.
+strobel=.false.
+cooling=.true.
+pref=.05
+alpha_top = 1e-11
+deltap = 0.7
+# folder in which correlated-k data is stored ?
+corrkdir   = pluton_ch4var_COfix
+# strictboundcorrk
+strictboundcorrk = .false.
+tplanckmin = 30
+dtplanck=0.1
+# call visible gaseous absorption in radiative transfer ?
+callgasvis = .true.
+# call continuum in radiative transfer ?
+continuum = .false.
+# Include Rayleigh scattering in the visible ?
+rayleigh   = .false.
+# Test physics timescale in 1D ?
+testradtimes = .false.
+# Output spectral OLR in 1D/3D?
+specOLR    = .true.
+# Output global radiative balance in file 'rad_bal.out' - slow for 1D!!
+meanOLR    = .false.
+# Variable gas species: Radiatively active ?
+varactive  = .true.
+# Variable gas species: Fixed vertical distribution ?
+#   (not to be used in time integration mode)
+varfixed   = .false.
+# Use fixed vertical profile, 1 step, no iteration ?
+kastprof     = .false.
+
+# Star type (inherited from generic model)
+startype = 1
+Fat1AU = 1366.0
+
+
+## Other physics options
+## ~~~~~~~~~~~~~~~~~~~~~
+# call turbulent vertical diffusion ?
+calldifv = .true.
+kmixmin=0.001
+vertdiff=.true.
+# use turbdiff instead of vdifc ? (if oldplutovdifc = .false.)
+UseTurbDiff = .true.
+oldplutovdifc=.true.
+# call convective adjustment ?
+calladj  = .true.
+# call thermal conduction in the soil ?
+callsoil = .true.
+lay1_soil = 2e-4
+# geothermal flux at the bottom layer
+fluxgeo = 0.
+# call thermal conduction
+callconduct=.false.
+# call molecular viscosity 3D only
+callmolvis=.false.
+# call molecular diffusion
+callmoldiff=.false.
+# flux qui arrive au-dessus de latmosphere
+phitop=0.
+
+
+## Tracer and aerosol options
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+# atm mass update due to tracer evaporation/condensation?
+mass_redistrib = .false.
+# Gravitational sedimentation of tracers ?
+sedimentation = .true.
+oldsedimpluto= .true.
+naerkind = 1
+## Condensation
+## ~~~~~~~~~~~
+generic_condensation = .false.
+# call N2 condensation ?
+n2cond       = .false.
+# N2 condensation subtimestep ?
+nbcond = 20
+# Set initial temperature profile to 1 K above N2 condensation everywhere?
+nearn2cond   = .false.
+# Compute methane cycle ?
+methane = .true.
+condmetsurf = .false.
+# Compute CO cycle ?
+carbox = .true.
+condcosurf=.false.
+# Compute methane clouds ?
+metcloud = .false.
+# Compute CO clouds ?
+monoxcloud = .false.
+# Haze production ?
+haze = .true.
+hazeconservch4    = .true.
+# haze particule mean radius
+rad_haze = 10e-9
+# number mixing ratio of haze particles for ch4 clouds
+Nmix_ch4      = 100000.
+# number mixing ratio of haze particles for co clouds
+Nmix_co      = 100000.
+# Radiatively active haze ?
+aerohaze=.true.
+hazeprop_file=optprop_rannou_fractal_r010nm_N1_1e4_d2.dat
+# fractal particle ?
+fractal=.true.
+nb_monomer=10
+# Fixed haze radius
+haze_radproffix = .true.
+hazerad_file = rad_haze.txt
+# Fixed haze distribution? # t==  haze_prof (hazemmr_file) else  hazecloud()
+haze_proffix = .true.
+# hazemmr_file = hazemmr.txt
+hazedens_file = dens_haze.txt
+
+## Surface properties
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+mode_n2=0     ! 0=fixed, 1=change w thickness, 2=special for SP
+deltab=0.    ! if mod_n2=1, difference of albedo depending of thickness
+metlateq=35  ! latitude separating dark/bright methane
+tholateq=25  ! latitude separating eq/pole tholins
+mode_ch4=1
+specalb=.false.  ! Specific treatment for tholins (diff values)
+
+alb_n2b=0.7
+alb_n2a=0.7
+alb_ch4=0.5
+alb_ch4_eq=0.5
+alb_ch4_s=0.5
+alb_co=0.4
+alb_tho=0.2
+alb_tho_eq=0.2
+
+emis_n2b=0.8
+emis_n2a=0.8
+emis_ch4=0.8
+emis_co=0.8
+emis_tho=0.8
+emis_tho_eq=0.8
+
+# Change TI with time
+changeti=.false.
+changetid=.false.
+ITCH4=800.
+ITN2=800
+ITH2O=2000.
+
Index: trunk/LMDZ.PLUTO/deftank/run_1d/rcm1d.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_1d/rcm1d.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_1d/rcm1d.def	(revision 3845)
@@ -0,0 +1,100 @@
+#
+#-----------------------------------------------------------------------
+# Run parameters for the 1D 'rcm1d.e' model
+#-----------------------------------------------------------------------
+
+#### Time integration parameters
+#
+# Initial date (in martian sols ; =0 at Ls=0)
+day0=1500
+# Initial local time (in hours, between 0 and 24)
+time=0
+# Number of time steps per sol
+day_step=10
+# Number of sols to run
+ndt = 2000
+ecritphy=20
+
+cpp_mugaz_mode = 2
+
+#### Physical parameters
+#
+# Day length, in s
+daysec = 551854.0625
+# length of a year, in days
+year_day = 14178.3037109375
+# Radius of the planet, in m
+rad = 1187000
+# Surface pressure (Pa)
+psurf=1.1
+# Gravity (ms^-2)
+g=0.617188990116119
+# Molar mass of atmosphere (g)
+mugaz=28
+# Specific heat capacity of atmosphere (J kg-1 K-1) # ignored if cpp_mugaz_mode = 1
+cpp=1040.0
+# latitude (in degrees)
+latitude=0.0
+# orbital distance at perihelion (AU)
+periastr=29.65831
+# orbital distance at aphelion (AU)
+apoastr=49.3050467
+# date of periastron, in days
+peri_day = 87.3600006103516
+# obliquity (degrees)
+obliquit=119.591
+
+# Albedo of bare ground
+albedo=0.5
+# Emissivity of bare ground
+emis=0.85
+# Soil thermal inertia (SI)
+inertia=20
+# zonal eastward component of the geostrophic wind (m/s)
+u=0
+# meridional northward component of the geostrophic wind (m/s)
+v=0.
+# Initial N2 ice on the surface (kg.m-2)
+n2ice=0
+# hybrid vertical coordinate ? (.true. for hybrid and .false. for sigma levels)
+hybrid=.false.
+
+# CH4 VMR (%)
+ch4ref=0.5
+# CO VMR (%)
+coref=0.046
+# haze ref (kg/kg)
+hazeref=0.000000385
+# prec haze ref (kg/kg)
+prechazeref=0.0000004
+
+# read haze from profhaze
+lecthaze = 1
+
+###### Initial atmospheric temperature profile
+#
+# Type of initial temperature profile
+#         ichoice=1   Constant Temperature:  T=tref
+#         ichoice=2   Savidjari profile (as Seiff but with dT/dz=cte)
+#         ichoice=3   Lindner (polar profile)
+#         ichoice=4   inversion
+#         ichoice=5   Seiff  (standard profile, based on Viking entry)
+#         ichoice=6   constant T  +  gaussian perturbation (levels)
+#         ichoice=7   constant T  + gaussian perturbation (km)
+#         ichoice=8   Read in an ascii file "profile"
+# ichoice=1
+ichoice=8
+# Reference temperature tref (K)
+tref=38
+# Add a perturbation to profile if isin=1
+isin=0
+# peak of gaussian perturbation (for ichoice=6 or 7)
+pic=26.522
+# width of the gaussian perturbation (for ichoice=6 or 7)
+largeur=10
+# height of the gaussian perturbation (for ichoice=6 or 7)
+hauteur=30.
+
+
+# some definitions for the physics, in file 'callphys.def'
+INCLUDEDEF=callphys.def
Index: trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/callphys.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/callphys.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/callphys.def	(revision 3845)
@@ -0,0 +1,273 @@
+## Orbit / general options
+## ~~~~~~~~~~~~~~~~~~~~~~~
+# Run with or without tracer transport ?
+tracer       = .true.
+# Diurnal cycle ?  if diurnal=false, diurnally averaged solar heating
+diurnal      = .false.
+# Seasonal cycle ? if season=false, Ls stays constant, to value set in "start"
+season       = .false.
+noseason_day = 1550
+# Fixed zenith angle in 1D
+global1d     = .true.
+# Tidally resonant orbit ? must have diurnal=false, correct rotation rate in newstart
+tlocked      = .false.
+# Tidal resonance ratio ? ratio T_orbit to T_rotation
+nres         = 1
+# Run 1D or 3D version (false) / fast version volatile transport (true)
+fast         = .false.
+
+# Write some more output on the screen ?
+lwrite    = .false.
+# Save statistics in file "stats.nc" ?
+callstats = .false.
+# Test energy conservation of model physics ?
+enertest  = .false.
+
+# Check the validity of input fields coming from the dynamics ?
+check_physics_inputs  = .false.
+# Check the validity of updated fields at the end of the physics step ?
+check_physics_outputs = .false.
+# Check minimum reasonable surface pressure (Pa) ?
+check_ps_min          = 1e-8
+
+# Conservation tracers
+conservn2  = .true.
+conservch4 = .true.
+
+# run with triton orbit ?
+triton = .false.
+
+## Directory where external input files are
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+datadir = datadir
+
+## Radiative transfer options
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+# call radiative transfer?
+callrad          = .true.
+# the rad. transfer is computed every "iradia" physical timestep
+iradia           = 1
+# call multilayer correlated-k radiative transfer ?
+corrk            = .true.
+oldplutocorrk    = .false.
+# Compute radiative methane?
+ch4fix           = .true.
+vmrch4fix        = 0.5
+vmrch4_proffix   = .false.
+nlte             = .true.
+strobel          = .false.
+cooling          = .true.
+pref             = 5.0e-2
+alpha_top        = 1.3e-11
+deltap           = 0.55
+# folder in which correlated-k data is stored ?
+corrkdir         = pluton_ch4var_COfix
+# strictboundcorrk / strictboundcia
+strictboundcorrk = .false.
+strictboundcia   = .false.
+tplanckmin       = 5
+# call visible gaseous absorption in radiative transfer ?
+callgasvis       = .true.
+# call continuum in radiative transfer ?
+continuum        = .true.
+# Include Rayleigh scattering in the visible ?
+rayleigh         = .true.
+# Test physics timescale in 1D ?
+testradtimes     = .false.
+# Output spectral OLR in 1D/3D?
+specOLR          = .false.
+# Output global radiative balance in file 'rad_bal.out' - slow for 1D!!
+meanOLR          = .false.
+# Variable gas species: Radiatively active ?
+varactive        = .true.
+# Variable gas species: Fixed vertical distribution ?
+# (not to be used in time integration mode)
+varfixed         = .false.
+# Use fixed vertical profile, 1 step, no iteration ?
+kastprof         = .false.
+# Output diagnostic optical thickness attenuation exp(-tau)?
+diagdtau         = .true.
+
+## Star type (startype = 1 <=> Sol)
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+startype = 1
+Fat1AU = 1366.0
+
+## Tracer options
+## ~~~~~~~~~~~~~~
+# atm mass update due to tracer evaporation/condensation?
+mass_redistrib = .false.
+
+### Volatile cycles
+### ~~~~~~~~~~~~~~~
+# call N2 condensation ?
+n2cond      = .false.
+# N2 condensation subtimestep ?
+nbsub       = 20
+# Compute methane cycle ?
+methane     = .true.
+condmetsurf = .false.
+# Compute CO cycle ?
+carbox      = .true.
+condcosurf  = .false.
+# Compute methane clouds ?
+metcloud    = .false.
+# Compute CO clouds ?
+monoxcloud  = .false.
+
+## Microphysics options
+## ~~~~~~~~~~~~~~~~~~~~
+# call microphysics moment model?
+callmufi    = .true.
+# Path to microphysical config file.
+config_mufi = datadir/microphysics/config.cfg
+
+### Radiatively active aerosols
+### ~~~~~~~~~~~~~~~~~~~~~~~~~~~
+# Number of radiatively active aerosols (if callmufi = true, must be = 2)?
+naerkind      = 2
+# Radiatively active haze aerosol?
+optichaze     = .true.
+aersprop_file = optprop_rannou_r2-200nm_nu003.dat
+aerfprop_file = optprop_rannou_fractal_r020nm_N1_1e4_d2.dat
+
+### Aerosol options (if callmufi = .true.)
+### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+# Equivalent radius production (m).
+haze_rc_prod        = 5.0e-9
+# Monomer radius (m).
+haze_rm             = 20e-9
+# Aerosol's fractal dimension.
+haze_df             = 2.0
+# Aerosol density (kg.m-3).
+haze_rho            = 800.
+# Radius of air molecule (m) ? (Pluto : N2)
+air_rad             = 1.75e-10
+# Use production from CH4 photolysis (if not, use production rate)?
+call_haze_prod_pCH4 = .true.
+hazeconservch4      = .true.
+# (call_haze_prod_pCH4 == false): 
+# Pressure level of aerosol production (Pa).
+haze_p_prod         = 1.0e-2
+# Aerosol production rate (kg.m-2.s-1).
+haze_tx_prod        = 9.8e-14
+
+### Old aerosol options (if callmufi = .false.)
+### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+# haze production?
+haze            = .false.
+# Gravitational sedimentation of tracers ?
+sedimentation   = .false.
+# Fixed haze distribution? (false == hazecloud, true == haze_prof)
+haze_proffix    = .false.
+hazemmr_file    = hazemmr.txt
+hazedens_file   = dens_haze.txt
+# Haze particule mean radius
+rad_haze        = 10e-9
+haze_radproffix = .false.
+hazerad_file    = rad_haze.txt
+fractal         = .true.
+nb_monomer      = 10
+hazeprop_file   = optprop_rannou_fractal_r010nm_N1_1e4_d2.dat
+# Number mixing ratio of haze particles for ch4 clouds
+Nmix_ch4        = 100000.
+# Number mixing ratio of haze particles for co clouds
+Nmix_co         = 100000.
+
+## Other physics options
+## ~~~~~~~~~~~~~~~~~~~~~
+# call turbulent vertical diffusion ?
+calldifv      = .true.
+oldplutovdifc = .false.
+# use turbdiff instead of vdifc ? (if oldplutovdifc = .false.)
+UseTurbDiff   = .true.
+kmixmin       = 0.001
+vertdiff      =.true.
+# call convective adjustment ?
+calladj       = .true.
+# call thermal conduction in the soil ?
+callsoil      = .true.
+lay1_soil     = 2e-4
+# call thermal conduction
+callconduct   = .true.
+# call molecular viscosity 
+callmolvis    = .false.
+# call molecular diffusion
+callmoldiff   = .false.
+# flux qui arrive au-dessus de latmosphere
+phitop        = 0.
+# geothermal flux at the bottom layer
+fluxgeo       = 0.
+
+## Surface properties
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+# options to change albedo with time
+#!  1) ALBEDOS and EMISSIVITY
+#!       A. N2
+#!                 CASE (0) ! fixed albedo
+#!                 CASE (1) ! Albedo decreases with thickness
+#!                 CASE (2) ! Special Sputnik differences of albedo
+#!                 CASE (3) ! Albedo increases (delta neg) or decreases (delta pos) with sublimationi rates
+#!                 CASE (4) ! Albedo Difference in N/S (e.g. used for Triton)
+#!                 CASE (5) ! Special Sputnik differences of albedo in small (1 pixel) patches (e.g. simulating dark patches / plumes)
+#!            --> EMISSIVITY N2: based on the alpha/beta transition
+#!       B. CO
+#!       C. CH4
+#!                 CASE (0) ! 2 albedos, one for the tropics, one for the poles
+#!                 CASE (1) ! 3 albedos, one for the tropics, 2 for the poles (north and south)
+#!                 CASE (2) ! 2 albedos + albedo feedback
+#!                          SELECT CASE (feedback_met)
+#!                            CASE (0) ! Default (linear from alb_ch4_eq)
+#!                            CASE (1) ! Hyperbolic tangent old
+#!                            CASE (2) ! hyperbolic tangent old
+#!                            CASE (3) ! hyperbolic tangent equation with parameters
+#!                 CASE (3) ! Eq, poles N, pole S + depending on Ls
+#!       D. Tholins
+#!                 CASE (0) ! Default, 2 albedos, one for the tropics, one for the poles
+#!                 CASE (1) ! Special mode one region with a different albedo
+#!       E. Tholins read from file
+#                 ! specalb
+
+# N2
+mode_n2     = 0  
+thres_n2ice = 0.001
+deltab      = 0. ! if mod_n2 = 1, difference of albedo depending of thickness
+
+# CH4
+mode_ch4     = 1
+thres_ch4ice = 0.01
+metlateq     = 35 ! latitude separating dark/bright methane
+
+# CO
+thres_coice = 1.
+
+# Tholins
+mode_tholins = 0
+tholateq     = 25      ! latitude separating eq/pole tholins
+specalb      = .false. ! Specific treatment for tholins (diff values)
+
+alb_n2b    = 0.7
+alb_n2a    = 0.7
+alb_ch4    = 0.5
+alb_ch4_eq = 0.5
+alb_ch4_s  = 0.5
+alb_co     = 0.4
+alb_tho    = 0.2
+alb_tho_eq = 0.2
+
+emis_n2b    = 0.8
+emis_n2a    = 0.8
+emis_ch4    = 0.8
+emis_co     = 0.8
+emis_tho    = 0.8
+emis_tho_eq = 0.8
+
+# Change TI with time
+changeti=.false.
+ITCH4  = 800.
+ITN2   = 800
+ITH2O  = 2000.
+changetid=.false.
+ITCH4d = 20.
+ITN2d  = 20.
+ITH2Od = 20.
Index: trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/diagfi.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/diagfi.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/diagfi.def	(revision 3845)
@@ -0,0 +1,52 @@
+tsurf
+ps
+ALB
+ISR
+ASR
+OLR
+GND
+DYN
+emis
+zdtsw
+zdtlw
+dtrad
+zdtdyn
+tau_col
+u
+v
+temperature
+p
+n2
+n2_col
+n2_surf
+n2_iceflux
+ch4_gas
+ch4_gas_col
+ch4_ice
+ch4_ice_col
+ch4_ice_surf
+ch4_iceflux
+co_gas
+co_gas_col
+co_ice
+co_ice_col
+co_ice_surf
+co_iceflux
+haze
+haze_surf
+haze_col
+prec_haze
+prec_haze_col
+vmr_ch4
+vmr_co
+m0as
+m3as
+m0af
+m3af
+rcs
+rcf
+dtauv01
+dtauv10
+dtauv23
+dtaui01
+dtaui17
Index: trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/gases.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/gases.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/gases.def	(revision 3845)
@@ -0,0 +1,11 @@
+# gases
+2
+CO_
+CH4
+1
+-1
+
+# First line is number of gases
+# Followed by gas names (always 3 characters)
+# and then mixing ratios.
+# mixing ratio -1 means the gas is variable.
Index: trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/rcm1d.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/rcm1d.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/rcm1d.def	(revision 3845)
@@ -0,0 +1,99 @@
+#-----------------------------------------------------------------------
+# Run parameters for the 1D 'rcm1d.e' model
+#-----------------------------------------------------------------------
+
+#### Time integration parameters
+#
+# Initial date (in pluto sols ; =1551 --> 2015)
+day0=1551
+# Initial local time (in hours, between 0 and 24)
+time=0
+# Number of time steps per sol
+day_step=10
+# Number of sols to run
+ndt = 5000
+ecritphy=20
+
+cpp_mugaz_mode = 2
+
+#### Physical parameters
+#
+# Day length, in s
+daysec = 551854.0625
+# length of a year, in days
+year_day = 14178.3037109375
+# Radius of the planet, in m
+rad = 1187000
+# Surface pressure (Pa)
+psurf=1.1
+# Gravity (ms^-2)
+g=0.617188990116119
+# Molar mass of atmosphere (g)
+mugaz=28
+# Specific heat capacity of atmosphere (J kg-1 K-1) # ignored if cpp_mugaz_mode = 1
+cpp=1040.0
+# latitude (in degrees)
+latitude=0.0
+# orbital distance at perihelion (AU)
+periastr=29.65831
+# orbital distance at aphelion (AU)
+apoastr=49.3050467
+# date of periastron, in days
+peri_day = 87.3600006103516
+# obliquity (degrees)
+obliquit=119.591
+
+# Albedo of bare ground
+albedo=0.5
+# Emissivity of bare ground
+emis=0.85
+# Soil thermal inertia (SI)
+inertia=20
+# zonal eastward component of the geostrophic wind (m/s)
+u=0
+# meridional northward component of the geostrophic wind (m/s)
+v=0.
+# Initial N2 ice on the surface (kg.m-2)
+n2ice=0
+# hybrid vertical coordinate ? (.true. for hybrid and .false. for sigma levels)
+hybrid=.false.
+
+# CH4 VMR (%)
+ch4ref=0.5
+# CO VMR (%)
+coref=0.046
+# haze ref (kg/kg)
+hazeref=0.000000385
+# prec haze ref (kg/kg)
+prechazeref=0.0000004
+
+# read haze from profhaze.out
+lecthaze = 1
+
+###### Initial atmospheric temperature profile
+#
+# Type of initial temperature profile
+#         ichoice=1   Constant Temperature:  T=tref
+#         ichoice=2   Savidjari profile (as Seiff but with dT/dz=cte)
+#         ichoice=3   Lindner (polar profile)
+#         ichoice=4   inversion
+#         ichoice=5   Seiff  (standard profile, based on Viking entry)
+#         ichoice=6   constant T  +  gaussian perturbation (levels)
+#         ichoice=7   constant T  + gaussian perturbation (km)
+#         ichoice=8   Read in an ascii file "profile"
+# ichoice=1
+ichoice=8
+# Reference temperature tref (K)
+tref=38
+# Add a perturbation to profile if isin=1
+isin=0
+# peak of gaussian perturbation (for ichoice=6 or 7)
+pic=26.522
+# width of the gaussian perturbation (for ichoice=6 or 7)
+largeur=10
+# height of the gaussian perturbation (for ichoice=6 or 7)
+hauteur=30.
+
+
+# some definitions for the physics, in file 'callphys.def'
+INCLUDEDEF=callphys.def
Index: trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/traceur.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/traceur.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/traceur.def	(revision 3845)
@@ -0,0 +1,12 @@
+11
+mu_m0as
+mu_m3as
+mu_m0af
+mu_m3af
+n2
+co_ice
+co_gas
+ch4_ice
+ch4_gas
+prec_haze
+haze
Index: trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/z2sig.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/z2sig.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_1d/ref_debatz_2025_07/z2sig.def	(revision 3845)
@@ -0,0 +1,48 @@
+18.00000     H: atmospheric scale height (km) (used as a reference only)
+0.010       Typical pseudo-altitude (m) for 1st layer (z=H*log(sigma))
+0.015       ,, ,,  ,, ,, ,, ,,  ,, ,, ,,  2nd layer, etc...
+0.025
+0.040
+0.060
+0.090
+0.1200
+0.1500
+0.1800
+0.2200
+0.2800
+0.3400
+0.5000
+0.70000
+0.90000
+1.20000
+1.50000
+1.80000
+2.10000
+2.50000
+3.00000
+3.60000
+5.00000
+6.50000
+8.00000
+9.50000
+11.0000
+13.0000
+16.0000
+19.0000
+22.0000
+25.0000
+29.0000
+32.0000
+36.0000
+40.0000
+45.0000
+50.0000
+55.0000
+60.0000
+66.0000
+72.0000
+80.0000
+95.0000
+115.000
+140.000
+180.000
Index: trunk/LMDZ.PLUTO/deftank/run_3d/README
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_3d/README	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_3d/README	(revision 3845)
@@ -0,0 +1,10 @@
+
+reference_debatz : 
+- with radiative haze
+- using XIOS
+- timestep 4800 (larger than the usual 3200)
+- 47 vertical levels
+- corrk : pluton_ch4var_COfix_17x27
+
+
+
Index: trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/callphys.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/callphys.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/callphys.def	(revision 3845)
@@ -0,0 +1,273 @@
+## Orbit / general options
+## ~~~~~~~~~~~~~~~~~~~~~~~
+# Run with or without tracer transport ?
+tracer       = .true.
+# Diurnal cycle ?  if diurnal=false, diurnally averaged solar heating
+diurnal      = .true.
+# Seasonal cycle ? if season=false, Ls stays constant, to value set in "start"
+season       = .true.
+noseason_day = 0
+# Tidally resonant orbit ? must have diurnal=false, correct rotation rate in newstart
+tlocked      = .false.
+# Tidal resonance ratio ? ratio T_orbit to T_rotation
+nres         = 1
+
+# Frequency (in physical steps) of diagfi file writing ?
+diagfi_output_rate = 32
+# Save statistics in file "stats.nc" ?
+callstats = .false.
+# Test energy conservation of model physics ?
+enertest  = .false.
+
+# Check the validity of input fields coming from the dynamics ?
+check_physics_inputs  = .false.
+# Check the validity of updated fields at the end of the physics step ?
+check_physics_outputs = .false.
+# Check minimum reasonable surface pressure (Pa) ?
+check_ps_min          = 1e-8
+
+# Conservation tracers diagnostics
+conservn2  = .false.
+conservch4 = .false.
+
+# run with triton orbit ?
+triton = .false.
+
+## Directory where external input files are
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+datadir = /scratchu/bdebatz/Pluto/datadir
+
+## Radiative transfer options
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+# call radiative transfer?
+callrad          = .true.
+# the rad. transfer is computed every "iradia" physical timestep
+iradia           = 5
+# call multilayer correlated-k radiative transfer ?
+corrk            = .true.
+oldplutocorrk    = .false.
+# Compute radiative methane?
+ch4fix           = .true.
+vmrch4fix        = 0.5
+vmrch4_proffix   = .false.
+nlte             = .false.
+strobel          = .false.
+cooling          = .false.
+pref             = 5.0e-2
+alpha_top        = 1e-11
+deltap           = 0.7
+# folder in which correlated-k data is stored ?
+corrkdir         = pluton_ch4var_COfix_17x27
+# strictboundcorrk
+strictboundcorrk = .false.
+strictboundcia   = .false.
+tplanckmin       = 5
+# call visible gaseous absorption in radiative transfer ?
+callgasvis       = .true.
+# call continuum in radiative transfer ?
+continuum        = .true.
+# Include Rayleigh scattering in the visible ?
+rayleigh         = .true.
+# Test physics timescale in 1D ?
+testradtimes     = .false.
+# Output spectral OLR in 1D/3D?
+specOLR          = .false.
+# Output global radiative balance in file 'rad_bal.out' - slow for 1D!!
+meanOLR          = .false.
+# Variable gas species: Radiatively active ?
+varactive        = .true.
+# Variable gas species: Fixed vertical distribution ?
+#   (not to be used in time integration mode)
+varfixed         = .false.
+# Use fixed vertical profile, 1 step, no iteration ?
+kastprof         = .false.
+
+## Star type (startype = 1 <=> Sol)
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+startype = 1
+Fat1AU = 1366.0
+
+
+## Other physics options
+## ~~~~~~~~~~~~~~~~~~~~~
+# call turbulent vertical diffusion ?
+calldifv      = .true.
+oldplutovdifc = .true.
+# use turbdiff instead of vdifc ? (if oldplutovdifc = .false.)
+UseTurbDiff   = .false.
+# call convective adjustment ?
+calladj       = .true.
+# call thermal conduction in the soil ?
+callsoil      = .true.
+lay1_soil     = 2e-4
+# call thermal conduction
+callconduct   = .true.
+# call molecular viscosity
+callmolvis    = .true.
+# call molecular diffusion
+callmoldiff   = .false.
+# flux qui arrive au-dessus de latmosphere
+phitop        = 0.
+# geothermal flux at the bottom layer
+fluxgeo       = 0.
+# calldrag_noro_mod
+calllott      = .false.
+
+
+## Tracer options
+## ~~~~~~~~~~~~~~
+# atm mass update due to tracer evaporation/condensation?
+mass_redistrib = .false.
+
+### Volatile cycles
+### ~~~~~~~~~~~~~~~
+# call N2 condensation ?
+n2cond      = .true.
+# N2 condensation subtimestep ?
+nbsub       = 20
+# No formation of N2 frost ?
+no_n2frost  = .true.
+# Compute methane cycle ?
+methane     = .true.
+condmetsurf = .true.
+# Compute CO cycle ?
+carbox      = .true.
+condcosurf  = .false.
+# Compute methane clouds ?
+metcloud    = .false.
+# Compute CO clouds ?
+monoxcloud  = .false.
+
+## Microphysics options
+## ~~~~~~~~~~~~~~~~~~~~
+# call microphysics moment model?
+callmufi    = .true.
+# Path to microphysical config file.
+config_mufi = /scratchu/bdebatz/Pluto/datadir/microphysics/config.cfg
+
+### Radiatively active aerosols
+### ~~~~~~~~~~~~~~~~~~~~~~~~~~~
+# Number of radiatively active aerosols (if callmufi = true, must be = 2)?
+naerkind      = 2
+# Radiatively active haze aerosol?
+optichaze     = .true.
+##aersprop_file = optprop_rannou_r2-200nm_nu003.dat
+aersprop_file = optprop_lora_pluto_r4-1000nm_nu003.dat
+aerfprop_file = optprop_rannou_fractal_r020nm_N1_1e4_d2.dat
+##aerfprop_file = optprop_lora_fractal_PH_r010nm_N1_5e5_d2.dat
+
+### Aerosol options (if callmufi = .true.)
+### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+# Equivalent radius production (m).
+haze_rc_prod        = 8.0e-9
+# Monomer radius (m).
+haze_rm             = 20e-9
+# Aerosol's fractal dimension.
+haze_df             = 2.0
+# Aerosol density (kg.m-3).
+haze_rho            = 800.
+# Radius of air molecule (m) ? (Pluto : N2)
+air_rad             = 1.75e-10
+# Use production from CH4 photolysis (if not, use production rate)?
+call_haze_prod_pCH4 = .true.
+hazeconservch4      = .true.
+# (call_haze_prod_pCH4 == false): 
+# Pressure level of aerosol production (Pa).
+haze_p_prod         = 1.0e-2
+# Aerosol production rate (kg.m-2.s-1).
+haze_tx_prod        = 9.8e-14
+
+### Old aerosol options (if callmufi = .false.)
+### ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
+# haze production?
+haze            = .false.
+# Gravitational sedimentation of tracers ?
+sedimentation   = .false.
+# Fixed haze distribution? (false == hazecloud, true == haze_prof)
+haze_proffix    = .false.
+hazemmr_file    = hazemmr.txt
+hazedens_file   = dens_haze.txt
+# Haze particule mean radius
+rad_haze        = 10e-9
+haze_radproffix = .false.
+hazerad_file    = rad_haze.txt
+fractal         = .false.
+nb_monomer      = 10
+hazeprop_file   = optprop_rannou_fractal_r010nm_N1_1e4_d2.dat
+# Number mixing ratio of haze particles for ch4 clouds
+Nmix_ch4        = 100000.
+# Number mixing ratio of haze particles for co clouds
+Nmix_co         = 100000.
+
+
+## Surface properties
+## ~~~~~~~~~~~~~~~~~~~~~~~~~~
+# options to change albedo with time
+#!  1) ALBEDOS and EMISSIVITY
+#!       A. N2
+#!                 CASE (0) ! fixed albedo
+#!                 CASE (1) ! Albedo decreases with thickness
+#!                 CASE (2) ! Special Sputnik differences of albedo
+#!                 CASE (3) ! Albedo increases (delta neg) or decreases (delta pos) with sublimationi rates
+#!                 CASE (4) ! Albedo Difference in N/S (e.g. used for Triton)
+#!                 CASE (5) ! Special Sputnik differences of albedo in small (1 pixel) patches (e.g. simulating dark patches / plumes)
+#!            --> EMISSIVITY N2: based on the alpha/beta transition
+#!       B. CO
+#!       C. CH4
+#!                 CASE (0) ! 2 albedos, one for the tropics, one for the poles
+#!                 CASE (1) ! 3 albedos, one for the tropics, 2 for the poles (north and south)
+#!                 CASE (2) ! 2 albedos + albedo feedback
+#!                          SELECT CASE (feedback_met)
+#!                            CASE (0) ! Default (linear from alb_ch4_eq)
+#!                            CASE (1) ! Hyperbolic tangent old
+#!                            CASE (2) ! hyperbolic tangent old
+#!                            CASE (3) ! hyperbolic tangent equation with parameters
+#!                 CASE (3) ! Eq, poles N, pole S + depending on Ls
+#!       D. Tholins
+#!                 CASE (0) ! Default, 2 albedos, one for the tropics, one for the poles
+#!                 CASE (1) ! Special mode one region with a different albedo
+#!       E. Tholins read from file
+#                 ! specalb
+
+# N2
+mode_n2     = 0
+thres_n2ice = 0.001
+deltab      = 0. ! if mod_n2 = 1, difference of albedo depending of thickness
+
+# CH4
+mode_ch4     = 1
+thres_ch4ice = 0.01
+metlateq     = 35 ! latitude separating dark/bright methane
+
+# CO
+thres_coice = 1.
+
+# Tholins
+mode_tholins = 0
+tholateq     = 25      ! latitude separating eq/pole tholins
+specalb      = .false. ! Specific treatment for tholins (diff values)
+
+alb_n2b    = 0.74
+alb_n2a    = 0.74
+alb_ch4    = 0.7
+alb_ch4_eq = 0.55
+alb_ch4_s  = 0.6
+alb_co     = 0.4
+alb_tho    = 0.1
+alb_tho_eq = 0.1
+emis_n2b    = 0.8
+emis_n2a    = 0.8
+emis_ch4    = 0.8
+emis_co     = 0.8
+emis_tho    = 1.
+emis_tho_eq = 1.
+
+# Change TI with time
+changeti=.false.
+ITCH4  = 800.
+ITN2   = 800.
+ITH2O  = 800.
+changetid=.false.
+ITCH4d = 20.
+ITN2d  = 20.
+ITH2Od = 20.
Index: trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/context_pcm_physics.xml
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/context_pcm_physics.xml	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/context_pcm_physics.xml	(revision 3845)
@@ -0,0 +1,157 @@
+<context id="LMDZ">
+
+    <!-- Definition of a calendar is done in code via wxios_set_cal -->
+
+    <!-- Definition of domains and groups of domains -->
+    <domain_definition>
+        <domain id="dom_glo" data_dim="2" />
+
+        <domain id="dom_32_24" type="rectilinear" ni_glo="32" nj_glo="24"   >
+          <generate_rectilinear_domain/>
+          <interpolate_domain order="1"/>
+        </domain>
+	
+	<domain id="dom_64_48" type="rectilinear" ni_glo="64" nj_glo="48"   >
+          <generate_rectilinear_domain/>
+          <interpolate_domain order="1"/>
+        </domain>
+
+        <domain id="dom_96_95" type="rectilinear"  ni_glo="96" nj_glo="95"  >
+          <generate_rectilinear_domain/>
+          <interpolate_domain order="1"/>
+        </domain>
+
+        <domain id="dom_144_142" type="rectilinear" ni_glo="144" nj_glo="142"   >
+          <generate_rectilinear_domain/>
+          <interpolate_domain order="1"/>
+        </domain>
+
+        <domain id="dom_512_360" type="rectilinear" ni_glo="512" nj_glo="360"   >
+          <generate_rectilinear_domain/>
+          <interpolate_domain order="1"/>
+        </domain>
+
+        <domain id="dom_720_360" type="rectilinear" ni_glo="720" nj_glo="360" >
+          <generate_rectilinear_domain/>
+          <interpolate_domain order="1"/>
+        </domain>
+
+        <domain id="dom_128_96" type="rectilinear" ni_glo="128" nj_glo="96" >
+          <generate_rectilinear_domain/>
+          <interpolate_domain order="1"/>
+        </domain>
+
+        <domain id="dom_256_192" type="rectilinear" ni_glo="256" nj_glo="192" >
+          <generate_rectilinear_domain/>
+          <interpolate_domain order="1"/>
+        </domain>
+
+        <domain id="dom_out" domain_ref="dom_32_24"/>
+    </domain_definition>
+    <!-- <domain_definition src="./domain_def.xml"/> -->
+
+    <!-- Definition of vertical axes
+                           NB: these must be defined and set in code via
+         xios_output_mod -->
+    <axis_definition>
+        <axis id="presnivs"
+              standard_name="Pseudo-pressure of model vertical levels"
+              unit="Pa">
+        </axis>
+        <axis id="altitude"
+              standard_name="Pseudo-altitude of model vertical levels"
+              unit="km">
+        </axis>
+        <axis id="soil_layers"
+              standard_name="soil layer depth at intermediate level"
+              unit="m">
+        </axis>
+        <axis id="interlayer"
+              standard_name="Levels between layers (altitude+1)"
+              unit="km">
+        </axis>
+        <axis id="controle_axe"
+              standard_name="Controle"
+              n_glo="100">
+        </axis>
+        <!--Added for diagspecIR-->
+        <axis id="IR_Wavenumber"
+               standard_name="Band mid frequency in the infrared"
+               unit="cm-1">
+        </axis>
+        <!--Added for diagspecVI-->
+        <axis id="VI_Wavenumber"
+              standard_name="Band mid frequency in the visible"
+              unit=
+            "cm-1">
+        </axis>
+
+    </axis_definition>
+
+    <scalar_definition>
+           <scalar id="scalar" />
+    </scalar_definition>
+
+    <grid_definition>
+        <!-- toggle axis id below to change output vertical axis -->
+        <grid id="grid_2D">
+            <domain domain_ref="dom_glo" />
+        </grid>
+        <grid id="grid_3D">
+            <domain domain_ref="dom_glo" />
+            <axis axis_ref="altitude" />
+        </grid>
+        <grid id="grid_soil">
+            <domain domain_ref="dom_glo" />
+            <axis axis_ref="soil_layers" />
+        </grid>       
+        <!--Added for diagspecIR-->
+        <grid id="grid_IR_3D">
+            <domain domain_ref="dom_glo" />
+            <axis axis_ref='IR_Wavenumber'/>
+        </grid>
+        <!-- Added for diagspecVI-->
+        <grid id="grid_VI_3D">
+          <domain domain_ref="dom_glo" />
+          <axis axis_ref="VI_Wavenumber"/>
+        </grid>
+        <!-- grid for interpolated lon-lat outputs -->
+        <!-- to output scalars -->
+        <grid id="scalarpoint" >
+            <scalar scalar_ref="scalar" />
+        </grid>
+        <grid id="altitude_grid" >
+            <axis axis_ref="altitude" />
+        </grid>
+        <!-- output grids -->
+        <grid id="grid_3D_out">
+            <domain domain_ref="dom_out" />
+            <axis axis_ref="altitude" />
+        </grid>
+        <grid id="grid_2D_out">
+            <domain domain_ref="dom_out" />
+        </grid>
+        <!--added for diagspecIR-->
+        <grid id="grid_IR_3D_out">
+            <domain domain_ref="dom_out"/>
+            <axis axis_ref="IR_Wavenumber"/>
+        </grid>
+        <grid id="grid_IR">
+            <axis axis_ref="IR_Wavenumber"/>
+        </grid>
+        <!--Added for diagspecVI-->
+        <grid id="grid_VI_3D_out">
+            <domain domain_ref="dom_out"/>
+            <axis axis_ref="VI_Wavenumber"/>
+        </axis>
+        <grid id="grid_VI">
+          <axis axis_ref="VI_Wavenumber"/>
+        </grid>
+    </grid_definition>
+
+    <!-- Definition of model variables -->
+    <field_definition src="./field_def_physics.xml" />
+
+    <!-- Definition of output files -->
+    <file_definition src="./file_def_physics.xml" />
+</context>
Index: trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/field_def_physics.xml
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/field_def_physics.xml	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/field_def_physics.xml	(revision 3845)
@@ -0,0 +1,513 @@
+<!-- ============================================================== -->
+<!-- field_def_physics.xml                                          -->
+<!-- Definition of all existing variables that can be outputed
+     from Generic PCM physics                                       -->
+<!-- ============================================================== -->
+
+    <field_definition prec="4"
+                      freq_op="1ts"
+                      enabled="true">
+
+        <!-- 0D variables -->
+        <field_group id="fields_0D" grid_ref="scalarpoint">
+            <field id="Ls"
+                   long_name="Solar longitude"
+                   unit="degrees" />
+            <field id="RA"
+                   long_name="Right Ascension"
+                   unit="degrees" />
+            <field id="Declin"
+                   long_name="Solar declination"
+                   unit="degrees" />
+            <field id="dist_star"
+                   long_name="Distance to star"
+                   unit="AU" />
+            <field id="globave"
+                   long_name="Surface pressure"
+                   unit="Pa" />
+        </field_group>
+
+        <field_group id="fields_controle" grid_ref="controle_grid">
+            <field id="controle"
+                   long_name="controle" />
+        </field_group>
+
+        <field_group id="fields_altitude" grid_ref="altitude_grid">
+            <field id="aps"
+                   long_name="aps" />
+            <field id="bps"
+                   long_name="bps" />
+        </field_group>
+
+        <field_group id="fields_interlayer" grid_ref="interlayer_grid">
+            <field id="ap"
+                   long_name="ap" />
+            <field id="bp"
+                   long_name="bp" />
+        </field_group>
+
+        <!-- 1D variables (spectral)-->
+        <field_group id="IR_fields" grid_ref="grid_IR">
+              <field id="IR_Bandwidth"
+                     long_name="Bandwidth in the infrared"
+                     unit="cm-1"/>
+        </field_group>
+        <field_group id="VI_field" grid_ref="grid_VI">
+              <field id="VI_Bandwidth"
+                     long_name="Bandwith in the visible"
+                     unit="cm-1"/>
+       </field_group>
+        <!-- 2D variables -->
+        <field_group id="fields_2D" grid_ref="grid_2D">
+            <field id="area"
+                   long_name="Mesh area"
+                   unit="m2" />
+            <field id="phisfi"
+                   long_name="Surface geopotential"
+                   unit="m2.s-2" />
+            <field id="emis"
+                   long_name="Surface emissivity"
+                   unit="" />
+            <field id="ALB"
+                   long_name="Surface Albedo"
+                    unit="" />
+
+            <field id="ps"
+                   long_name="Surface Pressure"
+                   unit="Pa" />
+            <field id="tsurf"
+                   long_name="Surface Temperature"
+                   unit="K" />
+            <field id="zdtsdif"
+                   long_name="tendancy Ts diff"
+                   unit="K.s-1" />
+
+            <field id="capcal"
+                   long_name="Surface Heat Capacity"
+                   unit="W.s.m-2.K-1" />
+
+            <!-- Tracers -->
+            <field id="mu_m0as_surf"
+                   long_name="moment of order 0 at surface (sph. aer.)"
+                   unit="m-2" />
+            <field id="mu_m3as_surf"
+                   long_name="moment of order 3 at surface (sph. aer.)"
+                   unit="m3.m-2" />
+            <field id="mu_m0af_surf"
+                   long_name="moment of order 0 at surface (fra. aer.)"
+                   unit="m-2" />
+            <field id="mu_m3af_surf"
+                   long_name="moment of order 3 at surface (fra. aer.)"
+                   unit="m3.m-2" />
+            <field id="n2_surf"
+                   long_name="N2 ice thickness"
+                   unit="kg/m2" />
+            <field id="ch4_gas_surf"
+                   long_name="CH4 vapor at surface"
+                   unit="kg/m2" />
+            <field id="ch4_ice_surf"
+                   long_name="CH4 ice at surface"
+                   unit="kg/m2" />
+            <field id="co_gas_surf"
+                   long_name="co vapor at surface"
+                   unit="kg/m2" />
+            <field id="co_ice_surf"
+                   long_name="co ice at surface"
+                   unit="kg/m2" />
+            <field id="haze_surf"
+                   long_name="haze at surface"
+                   unit="kg/m2" />
+            <field id="prec_haze_surf"
+                   long_name="precursor haze at surface"
+                   unit="kg/m2" />
+
+            <field id="mu_m0as_col"
+                   long_name="moment of order 0 column (sph. aer.)"
+                   unit="m-2" />
+            <field id="mu_m3as_col"
+                   long_name="moment of order 3 column (sph. aer.)"
+                   unit="m3.m-2" />
+            <field id="mu_m0af_col"
+                   long_name="moment of order 0 column (fra. aer.)"
+                   unit="m-2" />
+            <field id="mu_m3af_col"
+                   long_name="moment of order 3 column (fra. aer.)"
+                   unit="m3.m-2" />
+            <field id="aers_prec"
+                   long_name="Spherical aerosols precipitations"
+                   unit="kg.m-2.s-1" />
+            <field id="aerf_prec"
+                   long_name="Fractal aerosols precipitations"
+		    unit="kg.m-2.s-1" />
+            <field id="n2_col"
+                   long_name="N2 vapor column"
+                   unit="kg/m2" />
+            <field id="ch4_gas_col"
+                   long_name="CH4 vapor column"
+                   unit="kg/m2" />
+            <field id="ch4_ice_col"
+                   long_name="CH4 ice column"
+                   unit="kg/m2" />
+            <field id="co_gas_col"
+                   long_name="co vapor column"
+                   unit="kg/m2" />
+            <field id="co_ice_col"
+                   long_name="co ice column"
+                   unit="kg/m2" />
+            <field id="haze_col"
+                   long_name="haze column"
+                   unit="kg/m2" />
+            <field id="prec_haze_col"
+                   long_name="precursor haze column"
+                   unit="kg/m2" />
+
+            <field id="n2_iceflux"
+                   long_name="N2 flux"
+                   unit="kg/m2/s" />
+            <field id="ch4_iceflux"
+                   long_name="CH4 ice flux"
+                   unit="kg/m2/s" />
+            <field id="co_iceflux"
+                   long_name="co ice flux"
+                   unit="kg/m2/s" />
+
+            <field id="vmr_ch4"
+                   long_name="CH4 volume mixing ratio"
+                   unit="" />
+            <field id="vmr_co"
+                   long_name="CO volume mixing ratio"
+                   unit="" />
+
+            <field id="tau_col"
+                   long_name="Total aerosol optical depth"
+                   unit="opacity" />
+
+            <field id="zdtsurfc"
+                   long_name="tendancy Ts cond N2"
+                   unit="K.s-1" />
+            <field id="zdqsc_n2"
+                   long_name="tendancy tracer surf cond N2"
+                   unit="kg.kg-1.s-1" />
+            <field id="zdqsdif_n2"
+                   long_name="tendancy tracer surf diff"
+                   unit="kg.kg-1.s-1" />
+	    
+	    <field id="zdqsdif_ch4_ice"
+                   long_name="zdqsdif ch4"
+                   unit="" />
+            <field id="zdqssed_ch4"
+                   long_name="zdqssed ch4"
+                   unit="" />
+            <field id="zdqssed_haze"
+                   long_name="zdqssed haze"
+                   unit="" />
+
+            <!-- Radiation -->
+            <field id="ASR"
+                   long_name="Absorbed stellar radiation"
+                    unit="W/m2" />
+            <field id="ISR"
+                   long_name="incoming solar radiation"
+                   unit="W/m2" />
+            <field id="OLR"
+                   long_name="outgoing longwave radiation"
+                   unit="W/m2" />
+            <field id="GND"
+                   long_name="Heat flux from ground"
+                   unit="W/m2" />
+            <field id="DYN"
+                   long_name="Dynamical heat input"
+                   unit="W/m2" />
+            <field id="fluxsurf_lw"
+                   long_name="Longwave radiation at the surface"
+                   unit="W.m-2" />
+            <field id="sensibFlux"
+                   long_name="Sensible heat flux"
+                   unit="W/m2" />
+            <field id="lscaledE"
+                   long_name="heat from largescale"
+                   unit="W/m2" />
+            <field id="madjdE"
+                   long_name="heat from moistadj"
+                   unit="W/m2" />
+            <field id="genericconddE"
+                   long_name="heat from generic condensation"
+                   unit="W/m2" />
+
+        </field_group>
+
+        <!-- 3D variables -->
+        <field_group id="fields_3D" grid_ref="grid_3D">
+
+            <field id="zzlay"
+                   long_name="Midlayer altitude"
+                   unit="m" />
+            <field id="zzlev"
+                   long_name="Interlayer altitude"
+                   unit="m" />
+            <field id="pphi"
+                   long_name="Geopotential"
+                   unit="m2.s-2" />
+
+            <field id="temperature"
+                   long_name="Atmospheric temperature"
+                   unit="K" />
+            <field id="teta"
+                   long_name="Potential temperature"
+                   unit="K" />
+            <field id="u"
+                   long_name="Eastward Zonal Wind"
+                   unit="m/s" />
+            <field id="v"
+                   long_name="Northward Meridional Wind"
+                   unit="m/s" />
+            <field id="w"
+                   long_name="Downward Vertical Wind"
+                   unit="m/s" />
+            <field id="omega"
+                   long_name="Downward omega velocity"
+                   unit="Pa/s" />
+            <field id="p"
+                   long_name="Atmospheric pressure"
+                   unit="Pa" />
+            <field id="detr"
+                   long_name="Detrainment"
+                   unit="kg/m2/s" />
+            <field id="entr"
+                   long_name="Entrainment"
+                   unit="kg/m2/s" />
+            <field id="w_plm"
+                   long_name="Plume vertical velocity"
+                   unit="m/s" />
+            <field id="q2"
+                   long_name="turbulent kinetic energy"
+                   unit="J.kg^-1" />
+            <field id="dEzradsw"
+                   long_name="radiative heating"
+                   unit="w.m^-2" />
+            <field id="dEzradlw"
+                   long_name="radiative heating"
+                   unit="w.m^-2" />
+            <field id="qsatatm"
+                   long_name="atm qsat"
+                   unit=" " />
+            <field id="dt_generic_condensation"
+                   long_name="heating from generic condensation"
+                   unit="K s-1" />
+            <field id="zdudyn"
+                   long_name="Dyn. U"
+                   unit="m.s-2" />
+            <field id="zdumolvis"
+                   long_name="tendancy molvis"
+                   unit="m.s-2" />
+            <field id="zdvmolvis"
+                   long_name="tendancy molvis"
+                   unit="m.s-2" />
+            <field id="zdtsw"
+                   long_name="SW heating"
+                   unit="K.s-1" />
+            <field id="zdtlw"
+                   long_name="LW heating"
+                   unit="K.s-1" />
+            <field id="dtrad"
+                   long_name="radiative heating"
+                   unit="K.s-1" />
+            <field id="zdtdif"
+                   long_name="tendancy T diff"
+                   unit="K.s-1" />
+            <field id="zdtadj"
+                   long_name="tendancy T adj"
+                   unit="K.s-1" />
+            <field id="zdtconduc"
+                   long_name="tendancy conduc"
+                   unit="K.s-1" />
+            <field id="zdtdyn"
+                   long_name="Dyn. heating"
+                   unit="K.s-1" />
+            <field id="dtmoistadj"
+                   long_name="moist adj heating"
+                   unit="K.s-1" />
+            <field id="dtdryadj"
+                   long_name="dry adj heating"
+                   unit="K.s-1" />
+
+            <!--Tracers-->
+            <field id="mu_m0as"
+                   long_name="moment of order 0 for spherical aerosols"
+                   unit="kg-1" />
+            <field id="mu_m3as"
+                   long_name="moment of order 3 for spherical aerosols"
+                   unit="m3.kg-1" />
+            <field id="mu_m0af"
+                   long_name="moment of order 0 for fractal aerosols"
+                   unit="kg-1" />
+            <field id="mu_m3af"
+                   long_name="moment of order 0 for fractal aerosols"
+                   unit="m3.kg-1" />
+            <field id="n2"
+                   long_name="N2 mass mixing ratio"
+                   unit="kg/kg" />
+            <field id="ch4_gas"
+                   long_name="CH4 vapor mass mixing ratio"
+                   unit="kg/kg" />
+            <field id="ch4_ice"
+                   long_name="CH4 ice mass mixing ratio"
+                   unit="kg/kg" />
+            <field id="co_gas"
+                   long_name="co vapor mass mixing ratio"
+                   unit="kg/kg" />
+            <field id="co_ice"
+                   long_name="co ice mass mixing ratio"
+                   unit="kg/kg" />
+            <field id="haze"
+                   long_name="haze mass mixing ratio"
+                   unit="kg/kg" />
+            <field id="prec_haze"
+                   long_name="precursor haze mass mixing ratio"
+                   unit="kg/kg" />
+            <field id="zdtc"
+                   long_name="tendancy T cond N2"
+                   unit="K" />
+
+            <field id="zrho_ch4"
+                   long_name="zrho_ch4"
+                   unit="kg/m3" />
+            <field id="zrho_co"
+                   long_name="zrho_co"
+                   unit="kg/m3" />
+            <field id="zrho_haze"
+                   long_name="zrho_haze"
+                   unit="kg/m3" />
+            <field id="zdqcn2_ch4"
+                   long_name="zdq condn2 ch4"
+                   unit="" />
+            <field id="zdqdif_ch4"
+                   long_name="zdqdif_ch4"
+                   unit="" />
+            <field id="zdqadj_ch4"
+                   long_name="zdqadj_ch4"
+                   unit="" />
+            <field id="zdqsed_ch4"
+                   long_name="zdqsed_ch4"
+                   unit="" />
+            <field id="zdtch4cloud"
+                   long_name="ch4 cloud"
+                   unit="T s-1" />
+            <field id="zdqch4cloud"
+                   long_name="ch4 cloud"
+                   unit="T s-1" />
+            <field id="zdqrho_photprec"
+                   long_name="zdqrho_photprec"
+                   unit="kg/m3/s" />
+            <field id="zdqphot_prec"
+                   long_name="zdqphot_prec"
+                   unit="" />
+            <field id="zdqphot_ch4"
+                   long_name="zdqphot_ch4"
+                   unit="" />
+
+            <field id="zduc"
+                   long_name="tendancy U cond N2"
+                   unit="m.s-2" />
+            <field id="zdvc"
+                   long_name="tendancy V cond N2"
+                   unit="m.s-1" />
+            <field id="zdqc_n2"
+                   long_name="tendancy tracer cond N2"
+		   unit="kg.kg-1.s-1" />
+            <field id="zdqdif_n2"
+                   long_name="tendancy tracer diff"
+		   unit="kg.kg-1.s-1" />
+            <field id="zdqadj_n2"
+                   long_name="tendancy tracer adj"
+		   unit="kg.kg-1.s-1" />
+
+            <field id="zdqhaze_ch4"
+                   long_name="zdqhaze_ch4"
+                   unit="" />
+            <field id="zdqhaze_prec"
+                   long_name="zdqhaze_prec"
+                   unit="" />
+            <field id="zdqhaze_haze"
+                   long_name="zdqhaze_haze"
+                   unit="" />
+            <field id="zdqconv_prec"
+                   long_name="zdqconv_prec"
+                   unit="" />
+
+            <field id="m0as"
+                   long_name="Density number of spherical aerosols"
+                   unit="m-3" />
+            <field id="m3as"
+                   long_name="Volume of spherical aerosols"
+		    unit="m3.m-3" />
+            <field id="m0af"
+                   long_name="Density number of fractal aerosols"
+                   unit="m-3" />
+            <field id="m3af"
+                   long_name="Volume of fractal aerosols"
+		    unit="m3.m-3" />
+            <field id="aers_w"
+                   long_name="Spherical aerosol settling velocity"
+                   unit="m.s-1" />
+            <field id="aerf_w"
+                   long_name="Fractal aerosol settling velocity"
+		    unit="m.s-1" />
+            <field id="rcs"
+                   long_name="Characteristic radius of spherical aerosols"
+                   unit="m" />
+            <field id="rcf"
+                   long_name="Characteristic radius of fractal aerosols"
+		    unit="m" />
+
+            <field id="dtauv_01"
+                   long_name="Layer optical thickness attenuation in VI band"
+		    unit="" />
+            <field id="dtauv_23"
+                   long_name="Layer optical thickness attenuation in VI band"
+		    unit="" />
+            <field id="dtauv_24"
+                   long_name="Layer optical thickness attenuation in VI band"
+		    unit="" />
+            <field id="dtauv_27"
+                   long_name="Layer optical thickness attenuation in VI band"
+		    unit="" />
+            <field id="dtaui_01"
+                   long_name="Layer optical thickness attenuation in IR band"
+		    unit="" />
+            <field id="dtaui_17"
+                   long_name="Layer optical thickness attenuation in IR band"
+		    unit="" />
+
+        </field_group>
+
+        <field_group id="fields_soil" grid_ref="grid_soil">
+            <field id="tsoil"
+                   long_name="Soil Temperature"
+                   unit="K" />
+            <field id="therm_inertia"
+                   long_name="Soil Thermal Inertia"
+                   unit="S.I." />
+
+        </field_group>
+
+
+        <!--Spectral IR 3D variables-->
+        <field_group id="spectral_IR_3D" grid_ref="grid_IR_3D">
+              <field id="OLR3D"
+                     long_name="OLR(lon,lat,band)"
+                     unit="W/m^2/cm^-1" />
+        </field_group>
+
+        <!-- SPectral VI 3D variables-->
+        <field_group id="spectral_VI_3D" grid_ref="grid_VI_3D">
+              <field id="OSR3D"
+                     long_name="OSR(lon,lat,band"
+                     unit="W/m^2/cm^-1"/>
+              <field id="GSR3D"
+                     long_name="GSR(lon,lat,band)"
+                     unit="w/m^2/cm^-1"/>
+        </field_group>
+    </field_definition>
+
+
Index: trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/file_def_physics.xml
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/file_def_physics.xml	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/file_def_physics.xml	(revision 3845)
@@ -0,0 +1,126 @@
+<!-- =================================================================== -->
+<!-- file_def_physics.xml                                                -->
+<!-- Definition of output files for Pluto PCM physics                    -->
+<!-- =================================================================== -->
+
+  <file_definition>
+
+        <!-- ======================================================= -->
+        <!-- Instantaneous outputs:
+             1 ts = 1 physical time step.
+             If day_step=3200, iphys=10 and output_freq=20 ==> 16/Pd -->
+        <!-- ======================================================= -->
+	<file id="histins"
+	    name="Xhistins"
+     output_freq="60ts"
+       sync_freq="480ts"
+            type="one_file"
+    output_level="1"
+         enabled=".true.">
+
+
+	   <!-- VARS 0D -->
+       	   <field_group operation="instant"
+			 grid_ref="scalarpoint"
+                          freq_op="1ts">
+
+		  <field field_ref="Ls" level="1" />
+           </field_group>
+
+           <field_group operation="once"
+			 grid_ref="altitude_grid"
+                          freq_op="1ts">
+
+                  <field field_ref="aps" level="1"/>
+                  <field field_ref="bps" level="1"/>
+           </field_group>
+
+
+           <!-- VARS 2D -->
+	   <field_group operation="once"
+			 grid_ref="grid_2D"
+                          freq_op="1ts">
+
+                  <field field_ref="area" level="1" />  
+                  <field field_ref="phisfi" level="1" />
+	   </field_group>
+
+	   <field_group operation="instant"
+			 grid_ref="grid_2D"
+                          freq_op="1ts">
+
+                  <field field_ref="ps" level="1" />
+                  <field field_ref="tsurf" level="1" />
+                  <field field_ref="ISR" level="1" />
+                  <field field_ref="OLR" level="1" />
+                  <field field_ref="ASR" level="1" />
+		  
+		  <field field_ref="n2_surf" level="1" />
+                  <field field_ref="ch4_ice_surf" level="1" />
+                  <field field_ref="co_ice_surf" level="1" />
+                  <field field_ref="n2_iceflux" level="1" />
+                  <field field_ref="ch4_iceflux" level="1" />
+                  <field field_ref="ch4_ice_col" level="1" />
+                  <field field_ref="ch4_gas_col" level="1" />
+                  <field field_ref="co_iceflux" level="1" />
+                  <field field_ref="co_ice_col" level="1" />
+		  <field field_ref="co_gas_col" level="1" />
+                  <field field_ref="haze_surf" level="999" />
+                  <field field_ref="haze_col" level="999" />
+                  <field field_ref="prec_haze_col" level="999" />
+                  <field field_ref="vmr_ch4" level="1" />
+                  <field field_ref="vmr_co" level="1" />
+		  
+		  <field field_ref="tau_col" level="1" />
+		  
+                  <field field_ref="aers_prec" level="1" />
+		  <field field_ref="aerf_prec" level="1" />
+           </field_group>
+
+
+	   <!-- VARS 3D -->
+	   <field_group operation="instant"
+			 grid_ref="grid_3D"
+                          freq_op="1ts">
+
+                  <field field_ref="pphi" level="1" />
+                  <field field_ref="p" level="1" />
+                  <field field_ref="u" level="1" />
+                  <field field_ref="v" level="1" />
+                  <field field_ref="w" level="1" />
+                  <field field_ref="omega" level="999" />
+		  <field field_ref="temperature" level="1" />
+
+                  <field field_ref="zdtsw" level="1" />
+                  <field field_ref="zdtlw" level="1" />
+                  <field field_ref="zdtdyn" level="1" />
+		  
+		  <field field_ref="n2" level="1" />
+                  <field field_ref="ch4_gas" level="1" />
+                  <field field_ref="ch4_ice" level="1" />
+                  <field field_ref="co_gas" level="1" />
+                  <field field_ref="co_ice" level="1" />
+                  <field field_ref="haze" level="1" />
+                  <field field_ref="prec_haze" level="1" />
+                  <field field_ref="zdqphot_prec" level="1" />
+                  <field field_ref="zdqrho_photprec" level="1" />
+
+		  <field field_ref="m0as" level="1" />
+                  <field field_ref="m3as" level="1" />
+                  <field field_ref="m0af" level="1" />
+                  <field field_ref="m3af" level="1" />
+                  <field field_ref="aers_w" level="1" />
+		  <field field_ref="aerf_w" level="1" />
+                  <field field_ref="rcs" level="1" />
+		  <field field_ref="rcf" level="1" />
+		  
+		  <field field_ref="dtauv_01" level="1" />
+		  <field field_ref="dtauv_23" level="1" />
+		  <field field_ref="dtauv_24" level="1" />
+		  <field field_ref="dtauv_27" level="1" />
+		  <field field_ref="dtaui_01" level="1" />
+		  <field field_ref="dtaui_17" level="1" />
+           </field_group>
+
+	</file>
+  </file_definition>
Index: trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/gases.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/gases.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/gases.def	(revision 3845)
@@ -0,0 +1,11 @@
+# gases
+2
+CO_
+CH4
+1
+-1
+
+# First line is number of gases
+# Followed by gas names (always 3 characters)
+# and then mixing ratios.
+# mixing ratio -1 means the gas is variable.
Index: trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/iodef.xml
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/iodef.xml	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/iodef.xml	(revision 3845)
@@ -0,0 +1,27 @@
+<?xml version="1.0"?>
+<simulation>
+        <context  id="xios">
+            <variable_definition>
+                <!-- set print_file to true to have XIOS write its
+                     outputs in separate files -->
+                <variable id="print_file" type="bool"> true </variable>
+
+                <variable id="buffer_size_factor" type="double">4</variable>
+                <variable_group id="buffer">
+                buffer_size = 85000000
+                buffer_server_factor_size = 2
+                </variable_group>
+
+                <variable_group id="parameters" >
+                    <!-- set using_server to true to run with XIOS
+                         in detached server/client mode -->
+                    <variable id="using_server" type="bool">false</variable>
+                    <!-- info_level is the verbosity level of XIOS
+                         0 for no messages; 100 for lots -->
+                    <variable id="info_level" type="int">10</variable>
+                </variable_group>
+            </variable_definition>
+        </context>
+
+        <context id="LMDZ" src="./context_pcm_physics.xml"/>
+</simulation>
Index: trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/run.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/run.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/run.def	(revision 3845)
@@ -0,0 +1,99 @@
+#
+#-----------------------------------------------------------------------
+#Parametres de controle du run:
+#------------------------------
+
+ planet_type = generic
+
+# Nombre de jours d'integration
+ nday = 57
+
+# nombre de pas par jour (multiple de iperiod) ( ici pour dt = 1 min )
+ day_step = 4800
+
+ day_ini = 0
+
+# periode de la physique (en pas)
+  iphysiq=10
+
+# periode pour le pas Matsuno (en pas)
+ iperiod=5
+
+# periode de stockage fichier histmoy (en jour)
+ periodav=70.
+
+# periode de la dissipation (en pas)
+  idissip=5
+
+# choix de l'operateur de dissipation (star ou  non star )
+ lstardis=.true.
+
+# avec ou sans coordonnee hybrides
+ hybrid=.true.
+
+# nombre d'iterations de l'operateur de dissipation   gradiv
+nitergdiv=1
+
+# nombre d'iterations de l'operateur de dissipation  nxgradrot
+nitergrot=2
+
+# nombre d'iterations de l'operateur de dissipation  divgrad
+niterh=2
+
+# temps de dissipation des plus petites long.d ondes pour u,v (gradiv)
+tetagdiv=80000. # weak 50000 strong 5000
+
+# temps de dissipation des plus petites long.d ondes pour u,v(nxgradrot)
+tetagrot=200000. # weak 20000 strong 2000
+
+# temps de dissipation des plus petites long.d ondes pour  h ( divgrad)
+tetatemp=200000. # weak 20000 strong 2000
+
+# Values of dissip_factz are all defaults. Used ?
+dissip_factz =    2.00000000000000
+
+# Values of dissip_deltaz are all defaults. Used ?
+dissip_deltaz =    10.
+
+# Values of dissip_zref are all defaults.
+dissip_zref =    100.
+
+# Values of dissip_fac_mid are all defaults.
+dissip_fac_mid =    1
+
+# Values of dissip_fac_up are all defaults.
+dissip_fac_up =    8.
+
+# Values of dissip_hdelta are all defaults. Higher values = faster transition between mid and up
+dissip_hdelta =    10.
+
+# Values of dissip_pupstart are all defaults. Where is the transition (in Pa)
+dissip_pupstart =   1000.
+
+# coefficient pour gamdissip
+  coefdis=0.
+
+# choix du shema d'integration temporelle (Matsuno ou Matsuno-leapfrog)
+  purmats=.false.
+
+# avec ou sans physique
+   physic=.true.
+
+# choix d'une grille reguliere
+  grireg=.true.
+
+# Avec sponge layer
+  callsponge  = .true.
+
+# Sponge:  mode0(u=v=0), mode1(u=umoy,v=0), mode2(u=umoy,v=vmoy)
+  mode_sponge= 2
+
+# Sponge:  number of sponge layers
+  nsponge= 3
+
+# Sponge:  tetasponge (secondes)
+  tetasponge = 50000
+
+# some definitions for the physics, in file 'callphys.def'
+INCLUDEDEF=callphys.def
+
Index: trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/traceur.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/traceur.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/traceur.def	(revision 3845)
@@ -0,0 +1,12 @@
+11
+mu_m0as
+mu_m3as
+mu_m0af
+mu_m3af
+n2
+co_ice
+co_gas
+ch4_ice
+ch4_gas
+prec_haze
+haze
Index: trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/z2sig.def
===================================================================
--- trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/z2sig.def	(revision 3845)
+++ trunk/LMDZ.PLUTO/deftank/run_3d/reference_deBatz/z2sig.def	(revision 3845)
@@ -0,0 +1,48 @@
+18.00000     H: atmospheric scale height (km) (used as a reference only)
+0.010       Typical pseudo-altitude (m) for 1st layer (z=H*log(sigma))
+0.025
+0.040
+0.060
+0.090
+0.1200
+0.1500
+0.1800
+0.2200
+0.2800
+0.3400
+0.4100
+0.5000
+0.70000
+0.90000
+1.20000
+1.50000
+1.80000
+2.10000
+2.50000
+3.00000
+3.60000
+4.30000
+5.00000
+6.50000
+8.00000
+9.50000
+11.0000
+13.0000
+16.0000
+19.0000
+22.0000
+25.0000
+28.0000
+32.0000
+36.0000
+40.0000
+45.0000
+50.0000
+55.0000
+60.0000
+66.0000
+72.0000
+80.0000
+95.0000
+115.000
+140.000