Changeset 3196 for trunk/LMDZ.PLUTO

Timestamp:

Jan 31, 2024, 5:30:43 PM (10 months ago)

Author:

afalco

Message:

Pluto PCM:

Included N2 condensation from PLUTO.old & read haze aerosols.

AF

Location:

trunk/LMDZ.PLUTO

Files:

• changelog.txt (modified) (62 diffs)
• libf/phypluto/aeroptproperties.F90 (modified) (12 diffs)
• libf/phypluto/physiq_mod.F90 (modified) (2 diffs)

trunk/LMDZ.PLUTO/changelog.txt

@@ -39 +39 @@
     - adapted phyetat0.F for subterranean temperature & inertia
     - adapted physdem1.F to include new soil stuff
-    - physiq.F : added calls to writediagsoil 
+    - physiq.F : added calls to writediagsoil
     - soil.F : new routine (fixed vertical grid + variable thermal inertia)
     - added soil_settings.F (to read/initialize/interpolate soil properties)
@@ -82 +82 @@
     - modified vdifc.F and callsedim.F back
     - modified initracer.F (so that water names are h2o_vap & h2o_ice)
-    
+
 == 29/09/08 ==
     - modified aeronomars/init_chimie_B (cosmetics)
@@ -88 +88 @@
       "pause" error messages to 'stop' messages
     - modified phymars/watercloud.F to use tracers by name
-    - corrected aeronomars/molvis.F (undefined 'fac' and 'Akk' written to 
+    - corrected aeronomars/molvis.F (undefined 'fac' and 'Akk' written to
        output at first call)
-    
+
 == 30/09/08 ==
     - modified aeronomars/calchim.F to use tracers by name
@@ -171 +171 @@
      enforce (if set to .true.) that co2 mass mixing ratio remains bounded.
 == 07/11/08 ==
-    - corrected 'writediagfi' & 'writediagsoil' so that an error message is 
+    - corrected 'writediagfi' & 'writediagsoil' so that an error message is
       issued if called with a variable name which is too long.
 == 18/12/08 ==
@@ -230 +230 @@
         (removed dustopacity.F)
     - updated aeropacity.F with new tracer names
-    - changed the call to callradite.F in physiq.F, added the initialization 
+    - changed the call to callradite.F in physiq.F, added the initialization
         of reffrad and nueffrad (aerosol effective radius and variance)
     - removed all the lines relative to the old "activice" option, including
@@ -252 +252 @@
 
 == 06/07/09 ==
->>> Modified 'makegcm' and makegcm_g95' so that modules files are put 
+>>> Modified 'makegcm' and makegcm_g95' so that modules files are put
     with libraries (and not in current directory)
 
 ==
 == WORK FOR GENERIC MODEL STARTS HERE
-== 
-    
+==
+
 == 01/08/09 ==
 >>> Old Martian radiative transfer removed completely, new correlated-k
@@ -279 +279 @@
 >>> Many new options added to callphys.def.
 
-    
-    
-    
-    
+
+
+
+
 == 25/11/2009 ==
 >>>> List of changes and modifications
@@ -302 +302 @@
 -- adapted 'create_make_gcm' and 'makegcm' scripts so that we don't need
    to have and (unused) 'aeronomars' directory around.
- 
+
 -- shifted to reading traceur.def (dyn3d/iniadvtrac.F) in an Earth-LMDZ4-like
    fashion:
@@ -313 +313 @@
    such as rotation rate, molar mass of atmosphere, ...
    (NB: note however that we seem to lose dynamics/physics coherence ...
-    we should probably read these values from a planet.def file) 
+    we should probably read these values from a planet.def file)
 
 -- Changed condens_co2cloud.F : added a check for presence of CO2 gas and output
@@ -351 +351 @@
 >>> cosmetic change in callphys.def : separate "universal" parameters
     and planet-specific parameters (in a separate file, e.g. earlymars.def)
-    
+
 == 22/01/2010 ==
 ==> call this version LMDZ.GENERIC
@@ -388 +388 @@
    Two new functions were added: one to produce .nc files of the OLR
    by band, and the other to save the upgoing radiation by band in 1D
-   as a .txt file. 
+   as a .txt file.
 -- general structure --
    tabfi.F, newstart.F etc. cleaned up + it was checked that they pass
@@ -452 +452 @@
    Pressure is now distributed better in newstart.F. There was a bug before when the global mean of phi
    was not equal to zero and the surface pressure was changed with topography still present.
-   In initracer, we no longer automatically set h2o_vap=h2o_ice.   
+   In initracer, we no longer automatically set h2o_vap=h2o_ice.
 
 -- general --
@@ -460 +460 @@
    I cleaned up the method of displaying energy balance and temperature range in the model a bit.
    Pure H2 atmospheres are now permitted, although the code is _untested_. Rayleigh scattering,
-   specific heat capacities and molar masses are varied, and the H2-H2 collision-induced data 
+   specific heat capacities and molar masses are varied, and the H2-H2 collision-induced data
    of Grushka et al. is used for the radiative transfer.
 
@@ -472 +472 @@
    for use in aeropacity. Surface wetness (beta) is now calculated in a simple way, and
    outputed as a diagnostic, along with the atmospheric relative humidity.
-  
+
    mol/mol to kg/kg bug in callcorrk corrected.
 
 -- radiative transfer --
-   A major bug in gfluxv.F was corrected. It involved BSURF, the radiation reflected upwards from 
+   A major bug in gfluxv.F was corrected. It involved BSURF, the radiation reflected upwards from
    the ground.
    In aeropacity, we now set the cloud opacity in the top layer of the atmosphere to zero. This solves
@@ -501 +501 @@
 
 -- startup --
-   a bug involving tsoil interpolation in lect_start_archive.F was corrected.   
+   a bug involving tsoil interpolation in lect_start_archive.F was corrected.
 
 -- aerosols --
@@ -507 +507 @@
 
 -- water cycle --
-   Big fat bug in rain.F90 corrected. The routine internally updated the temperature 
+   Big fat bug in rain.F90 corrected. The routine internally updated the temperature
    with new tendencies before it was supposed to happen!
    Bug in oborealis where phi not properly updated was corrected.
@@ -517 +517 @@
 
 -- convection --
-   I found a bug in convadj.F that breaks tracer conservation when convection stops at one layer 
+   I found a bug in convadj.F that breaks tracer conservation when convection stops at one layer
    and immediately restarts in the next one.
 
 -- diagnostics --
-   Some bugs corrected in the computation of 3D averages for the H2O / energy conservation 
+   Some bugs corrected in the computation of 3D averages for the H2O / energy conservation
    diagnostics.
 
@@ -541 +541 @@
   to use "SOURCE" to identify code and be compatible with latest
   'create_make_gcm'.
- 
+
 == 10/11/2011 == EM
 - Upgrade: The location of the 'datagcm' directory can now be given in the
@@ -557 +557 @@
 - Allocation is done in sugas_corrk which is within firstcall loop in callcorrk.
 - Note that N_LGAUSS is still a parameter. It is not supposed to change much.
-- Compiles fine. Tested with debugging options through pgdbg. Runs fine. Exact same results in Early Mars test case. 
+- Compiles fine. Tested with debugging options through pgdbg. Runs fine. Exact same results in Early Mars test case.
 
 Modified files:
@@ -597 +597 @@
 
 == 27/02/2012 == AS
-- Temperature grid for Planck calculations can now be refined through the parameter NTfac in radinc_h. 
+- Temperature grid for Planck calculations can now be refined through the parameter NTfac in radinc_h.
   Default is NTfac = 1.0D-1, i.e. Delta T = 0.1 K
 
@@ -639 +639 @@
 - Cleaned continuum part in optc*
 - Added .def files for a typical 1d earth case in deftank (dry case for the moment)
-- Corrects a bug on potential temperature calculation in physic 
+- Corrects a bug on potential temperature calculation in physic
 
 == 27/03/2012 == JL
@@ -651 +651 @@
 - As much as possible, only the value of the stephan boltzmann constant defined in racommon_h (and the
 corresponding variable, sigma) should be used. Now done in physics, vdifc and turbdiff.
-  
+
 == 16/04/2012 == JL
 - Added consistency checks for calculations including water and global1d+diurnal.
@@ -657 +657 @@
 
 == 04/05/2012 == JL
-- Correction a huge bug in newstart: rcp and cpp can now be changed in start.nc files and are the same as in startfi.nc; 
+- Correction a huge bug in newstart: rcp and cpp can now be changed in start.nc files and are the same as in startfi.nc;
        Even when starting from start and startfi files.
    - rcp, cpp and mugaz can now be computed using gases.def in newstart
@@ -688 +688 @@
 - Some cleanup in start2archive.F and ini_archive.F to get them to work for
   the generic model (removed some "Martian" specificities).
-  
+
 == 17/07/2012 == JL for LK
 - Generalization of aerosol scheme:
@@ -694 +694 @@
     not important anymore.
     - addition of a module with the id numbers for aerosols (aerosol_mod.F90).
-    - initialization of aerosols id numbers in iniaerosol.F90 
+    - initialization of aerosols id numbers in iniaerosol.F90
     - compile with -s x where x *must* be equal to the number of aerosols turned on in callphys.def (either by a
     flag or by dusttau>0 for dust).
@@ -702 +702 @@
 - If starting from an old start file, recreate start file with the q=0 option in newstart.e.
 - update callphys.def with aeroXXX and aerofixXXX options (only XXX=co2,h2o supported for
-now). Dust is activated by setting dusttau>0. See the early mars case in deftank. 
+now). Dust is activated by setting dusttau>0. See the early mars case in deftank.
 - To add other aerosols, see Laura Kerber.
 
 == 18/07/2012 == JL
 - New water cycle scheme:
-   - largescale now in F90. Robustness increased by i) including evap inside largescale ii) computing the 
+   - largescale now in F90. Robustness increased by i) including evap inside largescale ii) computing the
         condensed water amount iteratively
    - same improvements in moistadj.
@@ -715 +715 @@
        - Saturation curves are now generalized for arbitrary water amount (not just q<<1)
        - The old watersat should be removed soon.
-   - The effect of water vapor on total (surface) pressure can be taken into account by setting 
+   - The effect of water vapor on total (surface) pressure can be taken into account by setting
       mass_redistrib=.true. in callphys.def (routine mass_redistribution inspired from co2_condense in martian
       model but with a different scheme as many routines evaporate/condense water vapor).
@@ -723 +723 @@
    - Totalcloudfraction now set the total cloud fraction to the fraction of the
       optically thickest cloud and totalcloudfrac is thus called in aeropacity.
-   - Only the total cloud fraction is used to compute optical depth in aeropacity (no more effective 
+   - Only the total cloud fraction is used to compute optical depth in aeropacity (no more effective
       optical depth with exponential formula).
    - 4 precipitation schemes are now available (see rain.F90 for details). The choice can be made using precip_scheme
@@ -730 +730 @@
    - See rain.f90 to determine the parameter to be defined in callphys.def as a function of the precipitation scheme used.
 - Physiq.F90 now written in a matricial (more F90) way.
-- Radii (H2O and CO2 cloud particles, aerosols, duts, ...) calculations now centralized in module radii_mod.F90 
+- Radii (H2O and CO2 cloud particles, aerosols, duts, ...) calculations now centralized in module radii_mod.F90
    and work with the new aerosol scheme implemented by Laura K. Some inconsistency may remain in callsedim.
 - Corrected gray gas mode. The use of kappa_VI or kappa_IR far a given band is now based on its wavenumber independently of
-   SW and LW calculations 
-   
+   SW and LW calculations
+
 Implementation compiled with ifort and pgf90.
-gcm.e runs in Earth and Early Mars case with CO2 and H2O cycle + dust. 
+gcm.e runs in Earth and Early Mars case with CO2 and H2O cycle + dust.
 
 == 19/07/2012 == JL
@@ -755 +755 @@
     -Can now handle any prograde resonance with nres=omega_rot/omega_orb.
     -Sun now goes westward for the standard 2:1 case, as expected.
-- In the gray case, the separation between kappa_IR and VI is now set by 
+- In the gray case, the separation between kappa_IR and VI is now set by
     wave number, independently of the usual IR/VISIBLE calculation separation.
     i.e. kappa_IR can be used in the calculation of the downward stellar flux
@@ -787 +787 @@
   --> debug mode : perfect match. bit by bit (diff command). checked with plots
   --> O1 mode : close match (checked with plots)
-  --> O2 mode : sometimes up to 0.5 K departure.... 
+  --> O2 mode : sometimes up to 0.5 K departure....
   BUT in this new version O2 and O1 are quite close while in previous version O1 and O2 differed by about, well, typically 0.5 K
   (pictures available on request)
@@ -802 +802 @@
     tracer.h >> tracer_h.F90
     comsaison.h >> comsaison_h.F90
-    comgeomfi.h >> comgeomfi_h.F90   
+    comgeomfi.h >> comgeomfi_h.F90
     comsoil.h >> comsoil_h.F90
     comdiurn.h >> comdiurn_h.F90
@@ -811 +811 @@
 - Suppressing any mention to advtrac.h which is a common in the dynamics and needs nqmx
   This was easily solved by adding an argument with tracer names, coming from the dynamics
-  This is probably not a definitive solution, 
+  This is probably not a definitive solution,
   ... but this allows for generic physics to work easily with either LMDZ.GENERIC or LMDZ dynamical cores
 - Removing consistency tests between nq and nqmx ; and ngrid and ngridmx. No use now!
@@ -819 +819 @@
 - Now written so that a slice of horizontal size 'ngrid' starting at grid point 'cursor' is read in startfi.nc
   'cursor' is defined in dimphys.h and initialized by inifis (or in newstart)
-  this is useful for parallel computations. 
+  this is useful for parallel computations.
   default behavior is the usual one : sequential runs, cursor is 1, size ngrid is the whole global domain
 
-A note on an additional change : 
+A note on an additional change :
 - nueffrad is now an argument to callcorrk as is the case for reffrad
   both are saved in physiq
   this is for consistency and lisibility (previously nueffrad was saved in callcorrk)
   ... but there is a call to a function which modifies nueffrad in physiq
-  ... previously this was not modifying nueffrad (although it was quite cumbersome to detect this) 
+  ... previously this was not modifying nueffrad (although it was quite cumbersome to detect this)
   ... to be conservative I kept this behaviour and highlighted it with an array nueffrad_dummy
   ... I added a comment because someone might want to change this
 
 == 20/09/2012 == AS
-A more robust way to count lines in setspi and setspv. 
+A more robust way to count lines in setspi and setspv.
 bandlen.txt file is no longer used. This was causing problems with MPI computations.
 
@@ -844 +844 @@
 
 == 20/12/2012 == EM
-- Fixed sedimentation issue: ensure in callsedim that the correct radii are 
+- Fixed sedimentation issue: ensure in callsedim that the correct radii are
   provided to newsedim and also that the updated temperature and tracer
   mixing ratios are used to compute sedimentation.
@@ -854 +854 @@
 
 == 10/01/2013 == AS
-- Fixed an allocating bug which arises from previous modifications 
+- Fixed an allocating bug which arises from previous modifications
   owing to the double use of callcorrk with CLF_varying
 - Fixed a small bug in a diagnostic in the end of calc_rayleigh. Some picky compilers complain.
@@ -870 +870 @@
 
 == 24/01/2013 == AS + JL
-A more robust way to refer to gas type. 
+A more robust way to refer to gas type.
 - Gas names with an arbitrary number of characters (<20) can be used
   This is good for C2H2, C2H6, H2SO4, C17H21NO4, etc...
   !!! Remember this must be compliant with Q.dat in corrk_data !!!
-- igas_... labels are assigned once for all in su_gases 
+- igas_... labels are assigned once for all in su_gases
   Then using igas_... everywhere instead of gnom (except for kcm stuff)
 - Users can still use e.g. H2_ but H2 also works
@@ -894 +894 @@
 - Corrected a bug: Continuum in inifis instead of continuum
     ... until now, most users (unbeknownst to them) were running with the continuum by default!
-- Cosmetic changes in optcv (mostly spaces and line breaks) 
+- Cosmetic changes in optcv (mostly spaces and line breaks)
     ... so that comparisons with optci are easy e.g. through vimdiff
-    
+
 == 11/02/2013 == JL
-- Updated moist convection scheme to handle situations with a large water vapor content 
+- Updated moist convection scheme to handle situations with a large water vapor content
 - Added a keyword to enable ocean runoff in callphys.def (activerunoff)
 
@@ -907 +907 @@
 - because N2 and H2O continua databases are smaller, improvement around 15% for
 an earth case.
- 
+
 == 25/02/2013 == EM
 - added "-fp-model precise" option when compiling with ifort (fixes unexplained
   mass loss of tracer in the dynamics)
-  
+
 == 13/03/2013 == AS + EM
 - corrected a bug in bilinearbig (found with debug mode)
@@ -924 +924 @@
 - optimization of optci and optcv routines. 15to 25% gain on these routines.
      around 10% on the whole code with 1 scatterer.
-     - No changes on output (byte to byte) 
-- corrected bug in gray case in callcorrk. 
+     - No changes on output (byte to byte)
+- corrected bug in gray case in callcorrk.
 - added profiling option in makegcm_ifort. See the file for details
 - changed capcal back to default in physics
 
 == 01/05/2013 == AS + JL
-- added possibility to remove surface 
+- added possibility to remove surface
   (particularly exchanges heat+momentum exchanges between surface and atmosphere)
 - added possibility to add an internal heat flux
@@ -956 +956 @@
   file name)
 - Some general cleanup:
-  - in bibio: removed unused lmdstd.h readstd.F writestd.F mywrite.F 
+  - in bibio: removed unused lmdstd.h readstd.F writestd.F mywrite.F
               readcoord.F scatter.F gather.F ini36.F from36.F to36.F
               lnblnk.F (F90 len_trim() should be used instead)
@@ -1033 +1033 @@
 
 == 10/07/2014 == EM
-Some cleanup to simplify dynamics/physics interactions by getting rid 
+Some cleanup to simplify dynamics/physics interactions by getting rid
 of dimphys.h (i.e. the nlayermx parameter) and minimizing use of
 dimension.h in the physics.
@@ -1066 +1066 @@
   interactions, e.g. routines gr_dyn_fi or gr_fi_dyn and calfis
 - Moreover the dynlonlat_phylonlat contains directory "phystd".
-  This subdirectories should only contain specific interfaces (e.g. 
+  This subdirectories should only contain specific interfaces (e.g.
   iniphysiq) or main programs (e.g. newstart).
 * phystd/dyn1d: this subdirectory contains the 1D models.
@@ -1076 +1076 @@
 
 == 16/09/2015 == EM
-- reorganizing the "datadir" structure: aerosol properties should now be in 
+- reorganizing the "datadir" structure: aerosol properties should now be in
   subdirectory 'aerosol_properties' of datadir, and surface.nc files should
   be in subdirectory 'surface_data'. These subdirectory names are stored in
@@ -1155 +1155 @@
 == 05/04/2016 == EM
 - Got rid of references to "dimensions.h" from physics packages:
-  use nbp_lon (=iim), nbp_lat (=jjp1) and nbp_lev 
+  use nbp_lon (=iim), nbp_lat (=jjp1) and nbp_lev
   from module mod_grid_phy_lmdz (in phy_common) instead.
 - Removed module "comhdiff_mod.F90", as it is only used by module
   surf_heat_transp_mod.F90, moved module variables there.
 - Added in "surf_heat_transp_mod" local versions of some arrays
-  and routines (from dyn3d) required to compute gradient, divergence, etc. 
+  and routines (from dyn3d) required to compute gradient, divergence, etc.
   on the global dynamics grid.
   As before, the slab ocean only works in serial.
@@ -1178 +1178 @@
 
 == 08/04/2016 == EM
-- Made nsoilmx be no longer a "parameter" and thus added the possibility to define the 
+- Made nsoilmx be no longer a "parameter" and thus added the possibility to define the
   number of subsurface layers nsoilmx, along with first layer thickness "lay1_soil" and
   companion coefficient "alpha_soil", in callphys.def at run time.
@@ -1191 +1191 @@
 
 == 18/04/2016 == EM
-- fix for 1D in writediagfi to enable writing at 
+- fix for 1D in writediagfi to enable writing at
   "ecritphy" rate.
 - move iniprint.h to "misc"
 - Some code cleanup in anticipation of future updates:
  - changed variable names in comgeomphy.F90: give them more
-   explicit names: rlond => longitude , 
+   explicit names: rlond => longitude ,
    rlatd => latitude, airephy => cell_area,
    cuphy => dx , cvphy => dy
@@ -1234 +1234 @@
 == 12/07/2016 == EM
 - move initialization of dimphy from inigeomphy to iniphysiq (initializations
-related to routines in phy_common or dynphy_lonlat can be done in 
+related to routines in phy_common or dynphy_lonlat can be done in
 inigeomphy, but any initialization for modules/routines in a physics
 package (directory phystd) must be done in the related phystd/iniphysiq
@@ -1260 +1260 @@
 == 15/09/2016 == TB
 - bug fix in tpindex.F : when computing weights for abundance of species one
-  should not return the maximum index (leads to out-of-bounds index use in 
+  should not return the maximum index (leads to out-of-bounds index use in
   the caller).
 
@@ -1283 +1283 @@
 
 == 21/02/2017 == EM
-Added possibility to run without a startfi.nc file (mainly usefull for 
+Added possibility to run without a startfi.nc file (mainly usefull for
 tests with coupling with dynamico dynamical core):
 - added flag "startphy_file" flag (.false. if doing an "academic" start
@@ -1304 +1304 @@
 - xios_output_mod.F90 : update initialize_xios_output initialization
                         of the horizontal domain
-In dynphy_lonlat : 
+In dynphy_lonlat :
 - inigeomphy_mod.F90 : add ind_cell_glo computation and  transfer
                        to init_geometry
@@ -1376 +1376 @@
 Start a series of commits to change the upper boundary conditions in the radiative transfer to solve some issues with the last two layers.
 It seems to be good to have aerosols in the first "radiative layer" of the gcm in the IR but visible does not handle very well diffusion in first layer.
-Tauaero and tauray are set to 0 (a small value for rayleigh because the code crashes otherwise) in the 4 first semilayers in optcv, but not optci. 
-This solves random variations of the sw heating at the model top. 
+Tauaero and tauray are set to 0 (a small value for rayleigh because the code crashes otherwise) in the 4 first semilayers in optcv, but not optci.
+This solves random variations of the sw heating at the model top.
 
 == 28/08/2018 == JL
-We now shift the radiative model top from p=0 to the middle of the last physical layer. This is done by changing pmid and plevrad in callcorrk and some corrections need to be done in gfluxv. 
+We now shift the radiative model top from p=0 to the middle of the last physical layer. This is done by changing pmid and plevrad in callcorrk and some corrections need to be done in gfluxv.
 This seems to get rid of the aratic temperature behavior in the last two layers of the model (especially on the night side on synchronous planets).
 Additional speedup corrections have been made in gfluxi that change nothing to the result.
@@ -1479 +1479 @@
 
 == 29/04/2019 == JVO
-+ Add diagnostics of optical thickness, if 'diagdtau' key is activated, it 
++ Add diagnostics of optical thickness, if 'diagdtau' key is activated, it
  outputs dtaui/v(altitude) in diagfi.nc for every narrowband (could be done with one var
  but would require to be able to have writediag in 5D)
@@ -1508 +1508 @@
 - the thermal plume model is able to manage several plumes in the same column and work without the convective adjustment.
 
-== 24/02/2020 (r2243) == JVO 
+== 24/02/2020 (r2243) == JVO
 - Get rid of the old 'sparadrap' in iniorbit.F assuming Mkm instead of AU if periastre or apostre gt 25, (otherwise Neptune is a pb)
   Everybody should be using AU, that's all.
@@ -1537 +1537 @@
  It can be called using 'aeronlay=.true.' in callphys.def, and set the number of layers (up to 4) with 'nlayaero'.
  Then, the following parameters are read as arrays of size nlayaero in callphys.def (separated by blank space)
- 
+
 *aeronlay_tauref  (Optical depth of aerosol layer at ref wavelenght)
 *aeronlay_lamref  (Ref wavelenght (m))
@@ -1561 +1561 @@
    + Following the last point some comment is added in suaer_corrk about the 'not-really-dummy'ness of IR lamref..
 
-== 27/04/2020 (r2299) == DB 
-Add non-orographic gravity waves drag parameterization (Flott scheme adpated from the Earth GCM) 
+== 27/04/2020 (r2299) == DB
+Add non-orographic gravity waves drag parameterization (Flott scheme adpated from the Earth GCM)
  It can be called using 'calllott_nonoro=.true.' in callphys.def, and set the maximum value of the Eliassen-Plam flux 'epflux_max'.
  Cumulated output fields are du_nonoro, dv_nonoro (winds tendency due to GW drag), east_gwstress and west_gwstress (stress profile in esatward and westward direction due to GW drag)
@@ -1570 +1570 @@
 When using the "season=.false." option then the orbital parameters are frozen and set by day_ini which is stored in startfi.nc.
 But at the end of the run "day_ini" written in the restartfi.nc file is updated to match current calendar evolution (to match dynamics), so in chained simulations orbital parameters in fact change with each new simulation...
-Add a parameter to lock a no seasonal cycle simulation during restart. 
-One sets the initial day using 'season=.false.' and set 'noseason_day'. This commit fixed the tickets #42 on planeto tack 
+Add a parameter to lock a no seasonal cycle simulation during restart.
+One sets the initial day using 'season=.false.' and set 'noseason_day'. This commit fixed the tickets #42 on planeto tack
 BE CARREFUL, there is a mistake in the previous commit message: it is actually SEASON=.FALSE.
 
 == 05/05/2020 (r2308) == JVO
-Fix a nasty copy-paste bug from r2297 in n-layer aerosol scheme 
+Fix a nasty copy-paste bug from r2297 in n-layer aerosol scheme
 
 == 05/06/2020 == EM
@@ -1589 +1589 @@
 
 == 09/07/2020 == DB
-Add kstar parameter to control kmin value (set value in the parameterization). kmin=1/lambda_max, to ensure the "subgrid scale" characteristic, we have to constrain the maximum GW's wavelength (lambda_max) by the size of the mesh. So, kmin=max(kstar, kmin prescribed value) 
+Add kstar parameter to control kmin value (set value in the parameterization). kmin=1/lambda_max, to ensure the "subgrid scale" characteristic, we have to constrain the maximum GW's wavelength (lambda_max) by the size of the mesh. So, kmin=max(kstar, kmin prescribed value)
 
 == 21/10/2020 == YJ
@@ -1630 +1630 @@
 
 == 16/03/2021 == AB
-Commit the last changes in the thermal plume model which was waiting for one year. 
+Commit the last changes in the thermal plume model which was waiting for one year.
 
 == 18/03/2021 == YJ
@@ -1655 +1655 @@
 - Update of the water continuum (H2O-H2O and H2O-AIR): MT_CKD v3.3
 - New file in phystd/: interpolateH2O_self_foreign.F90
-- New continua (H2O-AIR_continuum_MT_CKD3.3.cia and H2O-H2O_continuum_MT_CKD3.3.cia) 
+- New continua (H2O-AIR_continuum_MT_CKD3.3.cia and H2O-H2O_continuum_MT_CKD3.3.cia)
 can be downloaded from  https://www.lmd.jussieu.fr/~lmdz/planets/LMDZ.GENERIC/datagcm/continuum_data/
 - Now the water continnum is defined between 0.1 and 20000cm-1
-- We also removed the old (obsolete) parameterizations (CKD and PPC) from the code 
+- We also removed the old (obsolete) parameterizations (CKD and PPC) from the code
 
 ==16/06/2021 == MT
@@ -1713 +1713 @@
 
 == 22/02/2022 == GM
-Minor bug fix for aerave_new.F when input wavelenght in data file are in 
+Minor bug fix for aerave_new.F when input wavelenght in data file are in
 descending order.
 
@@ -1832 +1832 @@
 Add output for atmospheric escape fluxes at the top of the model,
 available for any tracers, this is calculated in deposition routine.
-Corection: add hydrology flag to calculate surface albedo used for the 
+Corection: add hydrology flag to calculate surface albedo used for the
 radiative transfer in photolysis online calculation.
 This was needed to be consistent with the way to calculate the surface albedo.
 
 == 25/01/2024 == AF
-Added LMDZ.PLUTO, a copy of the generic model cleaned of some 
+Added LMDZ.PLUTO, a copy of the generic model cleaned of some
 unnecessary modules
+
+== 31/01/2024 == AF
+Included N2 condensation from PLUTO.old & read haze aerosols.

trunk/LMDZ.PLUTO/libf/phypluto/aeroptproperties.F90

@@ -38 +38 @@
 !     ==============================================================
 
-!     Local variables 
+!     Local variables
 !     ---------------
 
@@ -198 +198 @@
 !      REAL :: omegaREFvis3d(ngrid,nlayer,naerkind)
 !      REAL :: omegaREFir3d(ngrid,nlayer,naerkind)
+
+      REAL :: minrad    ! minimal radius in table .dat radiustab (outside 0)
+      REAL :: maxrad
 
 !     0. Allocate local saved arrays at firstcall
@@ -246 +249 @@
         ALLOCATE(qsqrefIRb(L_NSPECTI,ngau,naerkind))
         ALLOCATE(qrefIRb(ngau,naerkind))
-        
+
         ALLOCATE(omegIRa(L_NSPECTI,ngau,naerkind))
         ALLOCATE(omegrefIRa(ngau,naerkind))
@@ -253 +256 @@
         ALLOCATE(gIRa(L_NSPECTI,ngau,naerkind))
         ALLOCATE(gIRb(L_NSPECTI,ngau,naerkind))
-        
+
         first_allocate=.false.
       ENDIF ! of IF (first_allocate)
@@ -317 +320 @@
            enddo
         endif
+
+        ! check that radiustab (table.dat) cover the reffrad used for haze : iaer=1
+        minrad=min(MINVAL(radiustab(1,1,1:nsize(1,1))),MINVAL(radiustab(1,2,1:nsize(1,2))))
+        maxrad=min(MAXVAL(radiustab(1,1,1:nsize(1,1))),MAXVAL(radiustab(1,2,1:nsize(1,2))))
+        IF ((MINVAL(reffrad).LE.minrad).OR.(MAXVAL(reffrad).GE.maxrad)) then
+           WRITE(*,*) 'Warning: particle size in grid box #'
+           WRITE(*,*) ig,' is too large to be used by the '
+           WRITE(*,*) 'radiative transfer; please extend the '
+           WRITE(*,*) 'interpolation grid to larger grain sizes.'
+           WRITE(*,*) 'radiustab=',minrad,'-',maxrad
+           WRITE(*,*) 'reffrad=',MINVAL(reffrad),'-',MAXVAL(reffrad)
+           stop
+        ENDIF
 
         firstcall = .false.
@@ -412 +428 @@
             qsqrefVISb(m,gausind,iaer)=                         &
                     (1-kint)*QVISsQREF(m,iaer,radius_id) +      &
-                    kint*QVISsQREF(m,iaer,radius_id+1)    
+                    kint*QVISsQREF(m,iaer,radius_id+1)
             omegVISb(m,gausind,iaer)=                           &
                     (1-kint)*omegaVIS(m,iaer,radius_id) +       &
@@ -500 +516 @@
                     EXP(-dista(j,1,iaer,idomain,gausind) *            &
                     dista(j,1,iaer,idomain,gausind) *                 &
-                    0.5e0/sizedistk2)/(radiusm-drad)                  
+                    0.5e0/sizedistk2)/(radiusm-drad)
                   dista(j,1,iaer,idomain,gausind) =                   &
                     dista(j,1,iaer,idomain,gausind) /                 &
@@ -540 +556 @@
 !                     Convolution:
                       qextVISgrid(j,1,m,iaer) =              &
-                        qextVISgrid(j,1,m,iaer) +            & 
+                        qextVISgrid(j,1,m,iaer) +            &
                         weightgaus(gausind) *                &
                         (                                    &
@@ -609 +625 @@
                       (                                      &
                       omegrefVISb(gausind,iaer) *            &
-                      qrefVISb(gausind,iaer) *               & 
+                      qrefVISb(gausind,iaer) *               &
                       pi*radGAUSb(gausind,iaer,idomain) *    &
                       radGAUSb(gausind,iaer,idomain) *       &
@@ -622 +638 @@
 
                   qrefVISgrid(j,1,iaer)=qrefVISgrid(j,1,iaer) /          &
-                                normd(j,1,iaer,idomain)       
+                                normd(j,1,iaer,idomain)
                   qscatrefVISgrid(j,1,iaer)=qscatrefVISgrid(j,1,iaer) /  &
                                 normd(j,1,iaer,idomain)
@@ -736 +752 @@
                       )
                   ENDDO
- 
+
                   qrefIRgrid(j,1,iaer)=qrefIRgrid(j,1,iaer) /          &
                                 normd(j,1,iaer,idomain)
@@ -789 +805 @@
             omegaIR3d(ig,lg,m,iaer) =                              &
                         k1*omegIRgrid(grid_i,1,m,iaer) +           &
-                        k2*omegIRgrid(grid_i+1,1,m,iaer) 
-            gIR3d(ig,lg,m,iaer) =                                  & 
+                        k2*omegIRgrid(grid_i+1,1,m,iaer)
+            gIR3d(ig,lg,m,iaer) =                                  &
                         k1*gIRgrid(grid_i,1,m,iaer) +              &
                         k2*gIRgrid(grid_i+1,1,m,iaer)

trunk/LMDZ.PLUTO/libf/phypluto/physiq_mod.F90

@@ -1474 +1474 @@
          enddo
 
-         ! Generalised for arbitrary aerosols now. By LK
+         if (aerohaze) then
+            ! Generalised for arbitrary aerosols now. By LK
          reffcol(1:ngrid,1:naerkind)=0.0
             ! call n2_reffrad(ngrid,nlayer,nq,zq,reffrad(1,1,iaero_haze))
@@ -1496 +1497 @@
             ENDDO
          ENDDO
+         endif ! end of aerohaze
       endif ! end of 'tracer'