| 1 | -------------------------------------------------------------------- |
|---|
| 2 | Changes in the LMD MARS GCM reference sources |
|---|
| 3 | (most recent version available in LMDZ.MARS.source.tar) |
|---|
| 4 | -------------------------------------------------------------------- |
|---|
| 5 | |
|---|
| 6 | 09/2008 dyn3d/physdem1.F |
|---|
| 7 | Correct possible bugs on some computers + english text & comments |
|---|
| 8 | |
|---|
| 9 | 09/2008 dyn3d/defrun_new.F |
|---|
| 10 | Corrected small bug (tapedef argument should not be altered) |
|---|
| 11 | |
|---|
| 12 | 09/2008 dyn3d/eofdump.F |
|---|
| 13 | Corrected small bug in the sampling time |
|---|
| 14 | |
|---|
| 15 | 09/2008 dyn3d/vanleer.F |
|---|
| 16 | removed inapropriate 'EXTERNAL' statement |
|---|
| 17 | |
|---|
| 18 | O9/2008 aeronomars/euvheat.F |
|---|
| 19 | removed print of 'ngrid','nlayer' (undefined and unused variables) |
|---|
| 20 | |
|---|
| 21 | O9/2008 phymars/mkstat.F90 |
|---|
| 22 | Cleaned up (removed unused variables) |
|---|
| 23 | |
|---|
| 24 | 09/2008 phymars/uniplot.F bibio/uniplot.F are removed because useless |
|---|
| 25 | |
|---|
| 26 | O9/2008 phymars/comcstfi.h |
|---|
| 27 | 09/2008 phymars/writediagfi.F |
|---|
| 28 | 08/2007 phymars/tabfi.F |
|---|
| 29 | Cosmetic: english text & comments |
|---|
| 30 | |
|---|
| 31 | 09/2008 dyn3d/iniwrite.F => phymars/iniwrite.F |
|---|
| 32 | - moved it from dyn3d to phymars |
|---|
| 33 | - cleanup & english |
|---|
| 34 | |
|---|
| 35 | 09/2008 phymars/physdem1.F |
|---|
| 36 | commented unused (wrong) 'external' qualifiers |
|---|
| 37 | |
|---|
| 38 | 09/2008 phymars/orosetup.F |
|---|
| 39 | Replaced logical 'xor' with 'neqv' (for standards compliance) |
|---|
| 40 | |
|---|
| 41 | 09/2008 phymars/physiq.F |
|---|
| 42 | - Change to ice albedo to fit new TES observation (Franck L.) |
|---|
| 43 | - cleanup + comments in english (Ehouarn) |
|---|
| 44 | |
|---|
| 45 | 09/2008 phymars/watercloud.F |
|---|
| 46 | cleaned version with reduced number of nuclei, and TES like output. |
|---|
| 47 | |
|---|
| 48 | 09/2008 phymars/growthrate.F |
|---|
| 49 | updated to improved version provided by F. Montmessin (Feb 2006...) |
|---|
| 50 | |
|---|
| 51 | 09/2008 phymars/vdifc.F |
|---|
| 52 | Change to remove an old quick fix to take into account condensation in |
|---|
| 53 | vertical diffusion |
|---|
| 54 | |
|---|
| 55 | 09/2008 |
|---|
| 56 | - phymars/mkstat.F90 |
|---|
| 57 | - phymars/def_var.F90 |
|---|
| 58 | Slight change to improve compilation (include => #include) |
|---|
| 59 | |
|---|
| 60 | O7/2008 MAJOR CHANGE to solve the "warm polar night" problem |
|---|
| 61 | - dyn3d/addfi.F |
|---|
| 62 | - dyn3d/gcm.F |
|---|
| 63 | Recompute mass after incrementation of ps in addfi.F (masse is added in the argument |
|---|
| 64 | list of addfi). |
|---|
| 65 | Also : the call of groupeun is now controlled by a flag. |
|---|
| 66 | |
|---|
| 67 | |
|---|
| 68 | 07/2008 phymars/physiq.F |
|---|
| 69 | - Change made (by Aymeric) to improve the models when iradia > 1 (iradia |
|---|
| 70 | is set is callphys.def). |
|---|
| 71 | (i.e. when one don't want to compute radiative transfer at every timestep): |
|---|
| 72 | The change has no impact when iradia=1 |
|---|
| 73 | => The net radiative surface flux "fluxrad" is computed outside the |
|---|
| 74 | test on iradia |
|---|
| 75 | => Add intermediate variable "fluxrad_sky" (rad. flux from sky absorbed |
|---|
| 76 | by surface) to be saved when rad. transfer is not computed. |
|---|
| 77 | - Cleaning of old commented lines about radiative transfer in clouds |
|---|
| 78 | |
|---|
| 79 | |
|---|
| 80 | 07/2008 Utility concatnc.F90 (not used by the gcm): improvement in order |
|---|
| 81 | to 1) concatenate 1D variable and 2) increase the number of input files |
|---|
| 82 | up to 1000 |
|---|
| 83 | |
|---|
| 84 | 12/2007 dyn3d/lect_start_archive.F |
|---|
| 85 | some old stuuf is commented to avoid bugs with some start_archive with |
|---|
| 86 | many tracers |
|---|
| 87 | |
|---|
| 88 | 12/2007 dyn3d/dump2d.F |
|---|
| 89 | Changed to improve crazy output to screen with newstart... |
|---|
| 90 | |
|---|
| 91 | 12/2007 bibio/mxva.F |
|---|
| 92 | Changed to avoid array bounds problems on some compilers |
|---|
| 93 | |
|---|
| 94 | 12/2007 |
|---|
| 95 | aeronomars/conduction.F |
|---|
| 96 | aeronomars/molvis.F |
|---|
| 97 | Changed to ensure variable initialisation on some compilers. |
|---|
| 98 | |
|---|
| 99 | 12/2007 phymars/convadj.F |
|---|
| 100 | Improvement regarding convection in the polar night |
|---|
| 101 | |
|---|
| 102 | 12/2007 |
|---|
| 103 | -phymars/callsedim.F |
|---|
| 104 | -phymars/initracer.F |
|---|
| 105 | -phymars/tracer.h |
|---|
| 106 | -phymars/callsedim2q.F |
|---|
| 107 | Changed to initialize tracer particle densities in initracer.F, and |
|---|
| 108 | transfer it with common in tracer.h |
|---|
| 109 | |
|---|
| 110 | 12/2007 script "makegcm" "create_make_gcm" |
|---|
| 111 | Improved to separate object compiled files in LIBO depending on the type |
|---|
| 112 | of machine you are (32bits or 64bits). Also allow to compute with LMDZ.2 |
|---|
| 113 | dynamics. |
|---|
| 114 | |
|---|
| 115 | 10/2007 Utility lslin.F90 (not used by the gcm) |
|---|
| 116 | changed evaluation of 'start_var' from hard-coded values to a computed value |
|---|
| 117 | |
|---|
| 118 | 09/2007 phymars/wstats.F90 and phymars/readtesassim.F90 |
|---|
| 119 | change " include" to "#include" to force recompilation by makegcm |
|---|
| 120 | |
|---|
| 121 | 08/2007 dyn3d/newstart.F90 |
|---|
| 122 | Translate french questions asked by program in english... |
|---|
| 123 | Improve options to remove topography... (for intercomparison) |
|---|
| 124 | |
|---|
| 125 | 08/2007 phymars/dustopacity.F |
|---|
| 126 | add possibility for scenario iaervar=5 (for paleoclimate...) |
|---|
| 127 | |
|---|
| 128 | 08/2007 phymars/tabfi.F |
|---|
| 129 | Translate french questions asked by program in english... |
|---|
| 130 | |
|---|
| 131 | ---------------------------------- |
|---|
| 132 | ---------------------------------- |
|---|
| 133 | ---------------------------------- |
|---|
| 134 | ---------------------------------- |
|---|
| 135 | |
|---|
| 136 | |
|---|
| 137 | 11/2006 phymars/convadj.F |
|---|
| 138 | Correction to remove a bug for compiler which did not "save" variable |
|---|
| 139 | by default (e.g. pgf on opteron, or g95) |
|---|
| 140 | |
|---|
| 141 | 11/2006 aeronomars/thermosphere.F |
|---|
| 142 | Correction to remove a bug if tracer=F and callthermos=T |
|---|
| 143 | |
|---|
| 144 | 10/2006 phymars/readtesassim.F90 |
|---|
| 145 | Correction to use dust scenario MY24 at high horizontal resolution |
|---|
| 146 | |
|---|
| 147 | 10/2006 dyn3d/dynredem.F and phymars/physdem1.F |
|---|
| 148 | Syntax changes to avoid problem when writing start files when compiling |
|---|
| 149 | on 64bits computers. |
|---|
| 150 | |
|---|
| 151 | 05/2006 : phymars/writediagfi.F |
|---|
| 152 | Add a test to be sure not to use writediagfi with 1D version of GCM |
|---|
| 153 | |
|---|
| 154 | 09/2005 : phymars/mucorr.F |
|---|
| 155 | Fix a weird bug that occured when running the model with "Diurnal=F" |
|---|
| 156 | on Linux, pgf90 |
|---|
| 157 | |
|---|
| 158 | ----------------------------------- |
|---|
| 159 | Update in june 2005 |
|---|
| 160 | ----------------------------------- |
|---|
| 161 | |
|---|
| 162 | 06/2005 : deftank/run0 |
|---|
| 163 | Improved so that the "run manager" run0 send you an E-mail when a run crash |
|---|
| 164 | (this can easily be commented) |
|---|
| 165 | |
|---|
| 166 | 06/2005 : deftank/run_mcd and deftank/callphys.def |
|---|
| 167 | Cosmetic changes to make thing clearer in english |
|---|
| 168 | |
|---|
| 169 | 06/2005 : makegcm |
|---|
| 170 | IMPORTANT FOR 1D GCM VERSION ON LINUX : modified to follow a change in |
|---|
| 171 | the Linux syntax of "sed" |
|---|
| 172 | |
|---|
| 173 | 06/2005 : aeronomars/calchim.F |
|---|
| 174 | NEW version provided by S. Lebonnois and F. Lefevre on May, 11, 2005 |
|---|
| 175 | ("Cosmetic changes")+syntax changes by FF (1 line > 72 characters) |
|---|
| 176 | |
|---|
| 177 | 06/2005 : phymars/convadj.F |
|---|
| 178 | Modified if one tracer is CO2, take into account the |
|---|
| 179 | Molecular mass variation (e.g. when CO2 condense) to trigger convection |
|---|
| 180 | |
|---|
| 181 | 06/2005 : phymars/newcondens.F |
|---|
| 182 | The CO2 deplation induced convection is now skipped, since it is properly |
|---|
| 183 | done in convadj.F (see above) |
|---|
| 184 | |
|---|
| 185 | 06/2005 : phymars/watercloud.F |
|---|
| 186 | IMPORTANT (especially for PALEOCLIMATE) : correction if a lot of subliming |
|---|
| 187 | CO2 fills the 1st layer. The number of input is modified (see below) |
|---|
| 188 | |
|---|
| 189 | 06/2005 : phymars/physiq.F |
|---|
| 190 | - Change in the input argument of watercloud.F : "pdpsrf" is added (see |
|---|
| 191 | above) |
|---|
| 192 | - If "caps", the cold trap at the south pole is only active for |
|---|
| 193 | obliquity lower than 27 deg (only useful for paleoclimate...) |
|---|
| 194 | |
|---|
| 195 | 06/2005 : dyn3d/groupe.F |
|---|
| 196 | Just cleaned to remove ambiguity of where "ngroup" is set (in groupeun, and |
|---|
| 197 | not in groupe). |
|---|
| 198 | |
|---|
| 199 | 06/2005 : phymars/inistats.F |
|---|
| 200 | remove declaration of "pseudoalt", already in comvert.h |
|---|
| 201 | (was a problem only for some compiler) |
|---|
| 202 | |
|---|
| 203 | 06/2005 : dyn3d/newstart.F |
|---|
| 204 | Cosmetic changes in some format for compatibility with gnu compiler |
|---|
| 205 | |
|---|
| 206 | 06/2005 : dyn3d/readhead_NC.F |
|---|
| 207 | Cosmetic changes in some format for compatibility with gnu compiler |
|---|
| 208 | |
|---|
| 209 | 06/2005 : dyn3d/sortvarc.F |
|---|
| 210 | Cosmetic changes in some format for compatibility with gnu compiler |
|---|
| 211 | |
|---|
| 212 | 06/2005 : phymars/testphys1d.F |
|---|
| 213 | Cleaning and changes in declaration for compatibility with gnu compiler |
|---|
| 214 | |
|---|
| 215 | 06/2005 : phymars/lwflux.F |
|---|
| 216 | Initialisation of "ksidb" for compiler that does not unitialized variable. |
|---|
| 217 | |
|---|
| 218 | 06/2005 : phymars/lwxd.F |
|---|
| 219 | Initialisation of "ksi_emis" for compiler that does not unitialized variable. |
|---|
| 220 | |
|---|
| 221 | 06/2005 : phymars/nirco2abs.F |
|---|
| 222 | Negligible change of one loop for compiler that does not like loop on real. |
|---|
| 223 | |
|---|
| 224 | 06/2005 : aeronomars/paramfoto.F |
|---|
| 225 | Negligible change of one loop for compiler that does not like loop on real. |
|---|
| 226 | |
|---|
| 227 | |
|---|
| 228 | -------------------------------------------------------------- |
|---|
| 229 | Change in December 2004 - January 2005 |
|---|
| 230 | -------------------------------------------------------------- |
|---|
| 231 | |
|---|
| 232 | 01/2005 : aeronomars/euvheat.F |
|---|
| 233 | IMPORTANT FOR THERMOSPHERE : correction to avoid EUV flux time variation |
|---|
| 234 | and keep it constant throughout the run |
|---|
| 235 | |
|---|
| 236 | 01/2005 : LMDZ.MARS/create_make_gcm |
|---|
| 237 | add phymars on line 187 to ensure recompilation of aeronomars subroutines |
|---|
| 238 | when any *.h file is updated in phymars |
|---|
| 239 | |
|---|
| 240 | 01/2005 : dyn3d/dynetat0.F |
|---|
| 241 | Add an option to reindex the tracer when using newstart if needed |
|---|
| 242 | |
|---|
| 243 | 01/2005 : phymars/phyetat0.F |
|---|
| 244 | Fix a bug when using newstart with more tracers than in old startfi |
|---|
| 245 | Add an option to reindex the surface tracer if needed |
|---|
| 246 | |
|---|
| 247 | 12/2004 : phymars/watercloud.F |
|---|
| 248 | IMPORTANT: Several Changes to fix a bug that made watercloud |
|---|
| 249 | to create negative mixing ratio in some conditions (past climate) |
|---|
| 250 | |
|---|
| 251 | 12/2004 : dyn3d/interp_horiz.F |
|---|
| 252 | IMPORTANT : add a test to fix a bug when using newstart with |
|---|
| 253 | a new grid size. |
|---|
| 254 | |
|---|
| 255 | 12/2004 : phymars/vdifc.F |
|---|
| 256 | FOR compilation on DEC machine : |
|---|
| 257 | Just change order of declaration of pqsurf(ngrid,nq |
|---|
| 258 | |
|---|
| 259 | 12/2004 : phymars/phyetat0.F |
|---|
| 260 | FOR compilation on DEC machine :Just change order of declaration |
|---|
| 261 | |
|---|
| 262 | 12/2004 : dyn3d/addfi.F |
|---|
| 263 | FOR compilation on DEC machine : qtestt = 1.0e-30 (before: qtestt = 1.0e-40) |
|---|
| 264 | |
|---|
| 265 | 12/2004 : LMDZ.MARS/makefile |
|---|
| 266 | Add new option to allow compilation on DEC machine |
|---|
| 267 | |
|---|
| 268 | -------------------------------------------------------------- |
|---|
| 269 | Change in August-september 2004 : before production of MCD V4.0 |
|---|
| 270 | -------------------------------------------------------------- |
|---|
| 271 | |
|---|
| 272 | 09/2004 : phymars/eofdump.F |
|---|
| 273 | Change to store density in EOF to build Mars Climate Database |
|---|
| 274 | |
|---|
| 275 | 09/2004 : phymars/physiq.F |
|---|
| 276 | add rho in the input of eofdump: |
|---|
| 277 | CALL eofdump(ngrid, nlayer, zu, zv, zt, rho, ps) |
|---|
| 278 | |
|---|
| 279 | 09/2004 : phymars/ineofdump.F |
|---|
| 280 | Change to store hybrid coordinates in EOF to build Mars Climate Database |
|---|
| 281 | |
|---|
| 282 | |
|---|
| 283 | 08/2004 : dyn3d/inidissip.F |
|---|
| 284 | Form change to simplify the choice of dissipation parameters |
|---|
| 285 | and use the same code for the thermosphere of only atmosphere model |
|---|
| 286 | |
|---|
| 287 | 08/2004 : dyn3d/newstart.F |
|---|
| 288 | Add an option "q=x" to initialise tracer manually |
|---|
| 289 | |
|---|
| 290 | 08/2004 : aeronomars/chimiedata.h |
|---|
| 291 | remove mmol,noms (now in tracer.h common) |
|---|
| 292 | |
|---|
| 293 | 08/2004 : In aeronomars: calchim.F, concentrations.F, euvheat.F, |
|---|
| 294 | inichim_newstart.F, moldiff.F, moldiffcoeff.F, perosat.F |
|---|
| 295 | Add : #include "tracer.h" (where tracer names are now) |
|---|
| 296 | |
|---|
| 297 | 08/2004 : phymars/inifis.F |
|---|
| 298 | Change some test on callphys.def options to allow any tracer transport... |
|---|
| 299 | |
|---|
| 300 | 08/2004 : phymars/initracer |
|---|
| 301 | Initialization of tracer name "noms(iq)" usable in many place |
|---|
| 302 | |
|---|
| 303 | 08/2004 : phymars/tracer.h |
|---|
| 304 | add noms,mmol to the common tracer |
|---|
| 305 | |
|---|
| 306 | 08/2004 : phymars/newcondens.F |
|---|
| 307 | IMPORTANT : major changes to allow proper treatment of non volatile |
|---|
| 308 | gaz enrichment when CO2 condense and depletion when CO2 sublime |
|---|
| 309 | |
|---|
| 310 | 08/2004 : phymars/writediagfi.F |
|---|
| 311 | remove the "abort" command so that the run does not stop due to a NETCDF bug |
|---|
| 312 | (NETCDF cannot write in file larger than 2098184 kb |
|---|
