1 | **** commit_v5 **** |
---|
2 | |
---|
3 | ** ADAPTATION GCM POUR CP(T) |
---|
4 | **-------------------------- |
---|
5 | |
---|
6 | * Ajout d'un fichier cpdet.F contenant toutes les routines nécessaires à |
---|
7 | l'utilisation de Cp(T) |
---|
8 | -> routines ini_cpdet, t2tpot, tpot2t |
---|
9 | -> fonction cpdet(t) |
---|
10 | |
---|
11 | * Du coup, suppression de tous les cpdet.* de phytitan et phyvenus |
---|
12 | |
---|
13 | * Modif de gcm.F : initialisation ini_cpdet ajoutée |
---|
14 | |
---|
15 | * Modif de comconst.h : ajout de nu_venus et t0_venus |
---|
16 | |
---|
17 | * Modif de caldyn0.F : |
---|
18 | - dimensions de pk (ip1jmp1 à la place de iip1,jjp1) |
---|
19 | - ajout de temp et tsurpk |
---|
20 | - calcul de temp (tpot2t) et tsurpk |
---|
21 | - passage de tsurpk à la place de teta dans sortvarc0 |
---|
22 | |
---|
23 | * Modif de caldyn.F : |
---|
24 | - passage de tsurpk en argument |
---|
25 | - dimensions de pk (ip1jmp1 à la place de iip1,jjp1) |
---|
26 | - ajout de temp |
---|
27 | - passage de tsurpk à la place de teta dans dudv2 |
---|
28 | - passage de tsurpk à la place de teta dans sortvarc |
---|
29 | |
---|
30 | * Modif de calfis.F : |
---|
31 | - ajout de zteta et zpk |
---|
32 | - calculs de zteta et zpk |
---|
33 | - calcul de ztfi via tpot2t |
---|
34 | - mise en commentaires de pcvgu/v/t/q car ne servent pas... |
---|
35 | - (ajout commentaire sur traceurs intensifs) |
---|
36 | - modification du calcul de pdhfi |
---|
37 | |
---|
38 | * Modif de diagedyn.F : |
---|
39 | (même si c'est une routine pure Terre, pour l'instant) |
---|
40 | - calcul de la temperature |
---|
41 | - remplacement de RCPD par cpdet(T) |
---|
42 | |
---|
43 | * Modif de leapfrog.F : |
---|
44 | - ajout de temp et tsurpk |
---|
45 | - calcul de temp (tpot2t) et tsurpk (3 endroits) |
---|
46 | - passage de tsurpk à la place de teta dans geopot (3 endroits) |
---|
47 | - passage de tsurpk en argument de caldyn |
---|
48 | - calcul de dtetaecdt |
---|
49 | |
---|
50 | * Modif de vlspltqs.F |
---|
51 | - modif taille de tempe |
---|
52 | - calcul de tempe (via tpot2t) |
---|
53 | |
---|
54 | |
---|
55 | ** APPEL DE LA PHYSIQUE |
---|
56 | **---------------------- |
---|
57 | |
---|
58 | * Modif de calfis.F : |
---|
59 | - on appelle physiq comme pour Venus et Titan quand on n'a pas la |
---|
60 | clé CPP_EARTH |
---|
61 | - j'ai enlevé le mot clé planet_type, ça fait double emploi, non ? |
---|
62 | - j'ai laissé nqtot dans les arguments d'appel à physiq... A reverifier. |
---|
63 | |
---|
64 | **** commit_v6 **** |
---|
65 | |
---|
66 | * Modif de defrun.F : commentaire indiquant qu'il ne sert plus donc pas a jour ! |
---|
67 | |
---|
68 | * Modif de dynredem.F : correction faute orthographe dans commentaire |
---|
69 | |
---|
70 | ** IOIPSL, real8 vs real4 |
---|
71 | **----------------------- |
---|
72 | |
---|
73 | J'élimine ce truc bancal. Il faudra compiler IOIPSL en double précision. |
---|
74 | Du coup, pas de modifs spécifiques dans dyn3d. Par contre, il faudra éliminer |
---|
75 | ça dans les physiques ! |
---|
76 | |
---|
77 | |
---|
78 | ** Nombre de traceurs |
---|
79 | **------------------- |
---|
80 | |
---|
81 | * Modif de infotrac.F90 : |
---|
82 | - definition de type_trac |
---|
83 | - definition de nqtrue depend de planet_type |
---|
84 | - lecture de traceur.def depend de planet_type |
---|
85 | |
---|
86 | * Modif de addfi.F : |
---|
87 | - separation des traceurs 1 et 2 pour qtestt/qtestw valable seulement |
---|
88 | pour planet_type=earth |
---|
89 | |
---|
90 | * Modif de caladvtrac.F : |
---|
91 | - deplacement du test sur planet_type pour virer tout ce qui concerne |
---|
92 | dq quand on n'est pas sur Terre |
---|
93 | |
---|
94 | * Modif de calfis.F : |
---|
95 | - ajout d'une remarque concernant pdq... |
---|
96 | |
---|
97 | |
---|
98 | ** Avec ou sans advection de traceurs |
---|
99 | **----------------------------------- |
---|
100 | |
---|
101 | Reglé grace à iadv lu dans traceur.def: si 0, pas advecté. |
---|
102 | Sans traceurs, on a donc un seul traceur dummy, non advecté, mais |
---|
103 | écrit dans start.nc et trimbalé partout... |
---|
104 | |
---|
105 | |
---|
106 | ** Homogeneisation des .def |
---|
107 | **------------------------- |
---|
108 | |
---|
109 | * Modif des deftanks/[titan/venus]/[gcm/physiq/run].def : |
---|
110 | - ajout planet_type |
---|
111 | - passage des cles physiques vers physiq.def |
---|
112 | - homogeneisation des parametres |
---|
113 | |
---|
114 | * Modif de conf_gcm.F : |
---|
115 | - ok_strato mis à true par defaut |
---|
116 | |
---|
117 | * Modif de clesph0.h : |
---|
118 | - ajout d'un commentaire: il faut virer ce fichier, ne sert plus a rien. |
---|
119 | |
---|
120 | |
---|
121 | ** Homogeneisation de bilan_dyn |
---|
122 | **----------------------------- |
---|
123 | |
---|
124 | * Modif de leapfrog.def : |
---|
125 | - premier appel: il y a first, deplacement de first=.false., elimination prem |
---|
126 | - ajout d'initialisations |
---|
127 | - ajout des tendances de la couche superieure (d*top) |
---|
128 | - elimination itau_w (ne sert pas) |
---|
129 | - couche superieure: ajout addfi séparé |
---|
130 | - passage des d*dis en */s |
---|
131 | - appel bilan_dyn |
---|
132 | |
---|
133 | * Modif de bilan_dyn.def : |
---|
134 | - modification appel (pour sortie des tendances en moment cinetique) |
---|
135 | - ajout tendances du* en m/s2 |
---|
136 | - ajout variables pour transport vertical |
---|
137 | - ajout variables pour tendances en moment cinetique |
---|
138 | - Champs advectés: base=6 (ovap est un traceur...), |
---|
139 | en commentaire pour l'instant, les traceurs. |
---|
140 | - introduction des calculs pour transport vertical |
---|
141 | - introduction des calculs pour tendances moment cinétique |
---|
142 | - présentation des déclarations |
---|
143 | - facteur fact_geovenus=-1 pour Vénus |
---|
144 | - petit ajustement possible pour ncum (problèmes de précision) |
---|
145 | - élminination des sorties grads |
---|
146 | - modif noms temp et ecin |
---|
147 | - modif sur incrementation itau et time |
---|
148 | - definitions et sorties: |
---|
149 | limité pour l'instant à temp et ang, sans les psiQ |
---|
150 | sans la somme en lat du transport vertical, à revoir (1D seulement?). |
---|
151 | - modification calcul ang |
---|
152 | - modification tpot2t |
---|
153 | - correction sur fonction de courant psiQ |
---|
154 | |
---|
155 | |
---|
156 | ** Choix de spécificités Terre dans gcm.F |
---|
157 | **--------------------------------------- |
---|
158 | ( certains appels ont été mis sous flag CPP_EARTH ou planet_type=earth alors ) |
---|
159 | ( qu'il sont aussi valables pour Venus et Titan... ) |
---|
160 | ( ils etaient sous flag IOIPSL avant... ) |
---|
161 | |
---|
162 | * Modif de gcm.F : |
---|
163 | - appel de dynetat0 modulé selon planet_type (mars est un cas à part) |
---|
164 | - appel de dynredem0 modulé selon planet_type (mars est un cas à part) |
---|
165 | - appel de iniphysiq: a mettre sous cle CPP_PHYS pour cas sans physique |
---|
166 | |
---|
167 | * Modif de leapfrog.F : |
---|
168 | - appel de dynredem1 modulé selon planet_type (mars est un cas à part) |
---|
169 | |
---|
170 | * Modif de calfis.F : |
---|
171 | - appel de physiq: a mettre sous cle CPP_PHYS pour cas sans physique |
---|
172 | - conservation de la possibilité de nsplit_phys |
---|
173 | - ajout commentaire sur unité des tendances |
---|
174 | |
---|
175 | * Modif de conf_gcm.F : |
---|
176 | - commentaires pour nsplit_phys |
---|
177 | |
---|
178 | |
---|
179 | **** commit_v7 **** |
---|
180 | |
---|
181 | * Mise à niveau (de la dynamique) par rapport à la version 1447 de LMDZ5 |
---|
182 | * ---------------------------------------------------------------------- |
---|
183 | |
---|
184 | * Pour l'instant, on n'a pas touché à dyn3dpar, donc pour l'instant on ne |
---|
185 | fait que recopier les fichiers de la version 1447 |
---|
186 | |
---|
187 | * Mise à jour de routines (dans dyn3d) qui n'ont changé que côté LMDZ5: |
---|
188 | - ener.h |
---|
189 | - grid_noro.F |
---|
190 | - iniacademic.F |
---|
191 | - integrd.F |
---|
192 | - limit_netcdf.F90 |
---|
193 | |
---|
194 | * Mise à jour des routines qui ont changé des deux côtés: |
---|
195 | - addfi.F |
---|
196 | * ajout d'un "USE control_mod, ONLY : planet_type" |
---|
197 | * mise en forme du if (planet_type=="earth") à l'identique de la revision 1447 (pour être peinard vis-à-vis de futures diff). |
---|
198 | - caladvtrac.F |
---|
199 | * ajout d'un "USE control_mod, ONLY : iapp_tracvl,planet_type" |
---|
200 | * mise en conformité de la taille du tableau dq(ip1jmp1,llm,nqtot) et non pas dq(ip1jmp1,llm,2) (qui oblige à avoir au moins 2 traceurs...) |
---|
201 | * mise entre if (planet_type=="earth") des calculs après advtrac (à nouveau, il faut avoir au moins 2 traceurs pour que ça ait un sens; et encore, ce morceau ne sert peut-être au final à rien). |
---|
202 | * A nouveau ces modifs étaient déjà en partie déjà faites; j'ai juste remis en forme comme dans la révision 1447 pour rester homogène avec le code terrestre. |
---|
203 | - gcm.F |
---|
204 | * Mise d'un #ifdef CPP_PHYS autoure de Init_Phys_lmdz et InitComgeomphy , au lieu du (bancal) #ifdef CPP_EARTH |
---|
205 | * De même il faut bel et bien un #ifdef CPP_PHYS autour de iniphysiq & co. |
---|
206 | - infotrac.F90 |
---|
207 | * Ajout des modifs de forme et correction sur l'utilisation de trim() et taille des tableaux str1(23),str2(23) |
---|
208 | * On garde un type_trac='plnt' pour les planètes (Attention, il faudra vérifier que ça ne casse rien par ailleurs et surtout n'est-ce finalement pas redoundant avec le planet_type?). |
---|
209 | - leapfrog.F |
---|
210 | * Ajout des modifs/améliorations pour le cas académique sans physique |
---|
211 | mais en utilisant dutop(),dvtop(), etc. pour le sponge layer |
---|
212 | |
---|
213 | **** commit_v8 **** |
---|
214 | |
---|
215 | ** ADAPTATION DYNAMIQUE PARALLELE AU CP(T) |
---|
216 | ------------------------------------------ |
---|
217 | |
---|
218 | Remarques: |
---|
219 | ---------- |
---|
220 | - Modifs non testées par une compilation/exécution du code; finalement s'organiser un mode de fonctionnement ou on peut compiler le code (au moins sans physique) est assez impératif. |
---|
221 | - On pourrait peut-être faire de cpdet.F90 un module... |
---|
222 | |
---|
223 | Retour sur leapfrog.F: (version séquentielle, dans dyn3d): |
---|
224 | ---------------------- |
---|
225 | - On revient en arrière sur l'utilisation de dutop(),dvtop(),dtetatop() au lieu des dufi(),dvfi() et dtetafi() pour les appels à top_bound. Dans la version terrestre, différente du sponge Martien (et Vénusien), ces tendances sont liées. |
---|
226 | => A 'hamoniser' plus tard. |
---|
227 | |
---|
228 | Modifications (dans dyn3dpar): |
---|
229 | ------------------------------ |
---|
230 | |
---|
231 | * infotrac.F90 : recopié à l'identique depuis dyn3d |
---|
232 | |
---|
233 | * comconst.h : recopié à l'identique depuis dyn3d |
---|
234 | |
---|
235 | * caldyn0.F : recopié à l'identique depuis dyn3d |
---|
236 | |
---|
237 | * cpdet.F : |
---|
238 | - création des routines t2tpot_p et tpot2t_p , versions // de t2tpot et tpot2t |
---|
239 | |
---|
240 | * gcm.F : |
---|
241 | - on encadre les appels à Init_Phys_lmdz et Init_interface_dyn_phys par des #ifdef CPP_PHYS (et non plus ifdef CPP_EARTH). |
---|
242 | - ajout de l'appel à ini_cpdet pour Cp(T) |
---|
243 | - appel de dynetat0 modulé selon planet_type (Mars est le cas pathologique) |
---|
244 | - appel à iniphysiq encadré par un #ifdef CPP_PHYS |
---|
245 | - appel de dynredem0_p modulé selon planet_type (Mars est le cas pathologique) |
---|
246 | |
---|
247 | * leapfrog_p.F : |
---|
248 | - ajout variables temp et tsurpk |
---|
249 | - ajout des calculs des hallos pour temp() et tsurkp() (NOTE POUR PLUS TARD: vérifier que c'est bien nécessaire...) |
---|
250 | - utilisation tpot2t_p et calcul de tsurkp() avant appel à geopot_p (sans oublier de modifier argument de geopot_p, teta => tsurpk). (3 fois) |
---|
251 | - ajout argument tsurpk dans l'appel à caldyn_p. |
---|
252 | - ajout des passages en (m/s)/s des dudis , dvdis et dtetadis. |
---|
253 | - calcul de dtetaecdt en passant par la température. |
---|
254 | - changement des arguments de l'appel à bilan_dyn_p. |
---|
255 | - appel de dynredem1_p modulé selon planet_type (Mars est le cas pathologique). |
---|
256 | |
---|
257 | * calfis_p.F : |
---|
258 | - utilisation de la clef CPP_PHYS (et non plus CPP_EARTH). |
---|
259 | - ajout des variables locales zteta et zpk. |
---|
260 | - ajout calcul de zteta et zpk. |
---|
261 | - calcul de ztfi via tpot2t_p (NOTE POUR PLUS TARD: vérifier que la taille "physique" du champ ne pose pas pb par rapport au découpage //). |
---|
262 | - appel de la routine physiq selon planet_type. |
---|
263 | - modification calcul de pdhfi. |
---|
264 | |
---|
265 | * caldyn_p.F : |
---|
266 | - ajout de tsurpk dans les arguments. |
---|
267 | - passage de tsurkp à la place de teta dans dudv2_p. |
---|
268 | - passage de tsurkp à la place de teta dans sortvarc. |
---|
269 | |
---|
270 | * vlspltqs_p.F : |
---|
271 | - modification de la taille de tempe (tempe(ip1jmp1,llm) au lieu de tempe(ip1jmp1). |
---|
272 | - calcul de tempe via tpot2t_p (NOTE POUR PLUS TARD: probablement à revoir! Etrangement vlspltqs_p ne contient aucune instruction OpenMP !?! Il y a là un truc qui m'échappe...) |
---|
273 | |
---|
274 | * bilan_dyn_p.F : A FAIRE... (mis de côté pour l'instant car pas sûr que l'ancien bilan_dyn_p fonctionne en mode mixte mpi/omp, à vérifier d'abord, avant de tenter des modifs...) |
---|
275 | |
---|
276 | **** commit_v11 **** |
---|
277 | |
---|
278 | 24/01/2011 ajout du modele meso-echelle martien, ancienne physique -- A. Spiga |
---|
279 | |
---|
280 | **** commit_v17 **** |
---|
281 | |
---|
282 | 24/01/2011 ajout du module LES, ancienne physique -- A. Spiga |
---|
283 | |
---|
284 | **** commit_v37 **** |
---|
285 | 27/01/2011 Remise en route chantier compilation -- Ehouarn |
---|
286 | - Modifs et corrections pour pouvoir compiler le gcm (en séquentiel, avec |
---|
287 | makelmdz_fcm pour l'instant): |
---|
288 | * ajout de fichiers 'arch' pour linux-64 (pour Bellonzi, avec ioipsl et en r8) |
---|
289 | * modification de makelmdz_fcm, ajout de la clé CPP_PHYS si on compile avec |
---|
290 | une physique |
---|
291 | * correction de quelques typos/bugs révélés à la compilation: |
---|
292 | - infotrac.F90 : supression des appels à 'lnblnk' (remplacés par len_trim) |
---|
293 | - bilan_dyn.F : déclaration des variables znom3,znom3l,zunites3, planet_type |
---|
294 | - cpdet.F : "use control_mod, ONLY: planet_type" mis aux bons endroits |
---|
295 | (idem sur cpdet.F dans dyn3dpar) |
---|
296 | - leapfrog.F : declaration de ztetaec(), dtec, cpdet , itau_w, duspg() |
---|
297 | - diagedyn.F : correction typo; attention dans diagedyn.F il y a du |
---|
298 | include "../phylmd/YOMCST.h" |
---|
299 | Ca va poser problème qd on change de physique.... |
---|
300 | |
---|
301 | - Avec ces modifs, la compilation marche sans physque, avec ou sans ioipsl |
---|
302 | et avec la physique terrestre phylmd. |
---|
303 | |
---|
304 | ******************************************** |
---|
305 | **** commit_v51 (et quelques autres...) **** |
---|
306 | ******************************************** |
---|
307 | |
---|
308 | 29/01/2011 --- A. SPIGA |
---|
309 | |
---|
310 | Nouveautes LMD_MM_MARS - LMD_LES_MARS - LMDZ.MARS [dossier mars] |
---|
311 | |
---|
312 | Ce commit fait le bilan des chantiers d'uniformisation des diverses sources du mesochelle |
---|
313 | [tous les commits par "spiga" depuis le debut sont relatifs a ces chantiers] |
---|
314 | |
---|
315 | Chantiers acheves avec compilation et execution testees: |
---|
316 | |
---|
317 | 1. le meme coeur dynamique peut maintenant etre compile avec ancienne ou nouvelle physique; option -p sur makemeso. |
---|
318 | |
---|
319 | 2. la meme interface module_lmd_driver est partagee par LMD_MM_MARS et LMD_LES_MARS |
---|
320 | |
---|
321 | 3. le meme libf est partage entre modele mesoechelle et gcm |
---|
322 | les quelques routines qui different sont notees meso_ et integrees aux sources du gcm |
---|
323 | |
---|
324 | =====> De 4 version du modele + 1 version de la nouvelle physique, on passe a une seule version du modele |
---|
325 | meso-echelle dont la physique est automatiquement indexee sur les modifications effectuees dans le GCM |
---|
326 | [quelques routines mesoechelle mises a part] |
---|
327 | |
---|
328 | * Reste a tester LMD_LES_MARS avec la nouvelle physique. Modifications mineures attendues suite au point 2. |
---|
329 | * Reste a tester le nesting qui devrait induire des modifications dues a des lignes trop longues... ce sera transparent du point de vue GCM. |
---|
330 | |
---|
331 | ******************** |
---|
332 | **** commit_v52 **** |
---|
333 | ******************** |
---|
334 | |
---|
335 | ********** chantier principal du commit |
---|
336 | ********** --- version LMDZ5 qui fonctionne pour tests geantes |
---|
337 | ********** --- prochaine etape, tests sur GNOME |
---|
338 | |
---|
339 | M libf/dyn3dpar/comconst.h |
---|
340 | M libf/dyn3dpar/conf_planete.F90 |
---|
341 | ajout du flux de chaleur intrinseque: ihf |
---|
342 | [par defaut il est nul] |
---|
343 | |
---|
344 | M libf/dyn3dpar/gcm.F |
---|
345 | changements cosmetiques |
---|
346 | [pour diff plus efficace avec version non par] |
---|
347 | |
---|
348 | M libf/dyn3dpar/iniacademic.F |
---|
349 | possibilites de variations latitudinales |
---|
350 | de temperature plus originales |
---|
351 | [seulement pour planet_type.eq."giant"] |
---|
352 | |
---|
353 | M libf/dyn3dpar/leapfrog_p.F |
---|
354 | 1. ajout d'une tendance causee par le flux de chaleur intrinseque |
---|
355 | (seulement prise en compte si planet_type.eq."giant") |
---|
356 | 2. correction bugs problematiques a la compilation et au run |
---|
357 | --> probleme dans les boucles (l'indice etait llm et non l) |
---|
358 | --> ajout de SAVE pour les variables paralleles |
---|
359 | --> correction des declarations de variables manquantes |
---|
360 | |
---|
361 | M libf/dyn3dpar/calfis_p.F |
---|
362 | correction d'une deuxieme parenthese manquante sur ALLOCATE(zteta(klon,llm)) |
---|
363 | |
---|
364 | M libf/phylmd/regr_lat_time_climoz_m.F90 |
---|
365 | erreur a la compilation avec FCM... il s'agit d'une routine terrestre |
---|
366 | il y a visiblement un probleme avec o3_in |
---|
367 | en attendant, les lignes sont commentees avec !AS |
---|
368 | |
---|
369 | A deftanks/giant 8 fichiers |
---|
370 | ajout de fichiers de configuration typiques pour les geantes gazeuses |
---|
371 | [experimental pour le moment... on est loin de jupiter] |
---|
372 | --> comparaisons entre un run ancien [avec LMDZ5-dev sur SVN ipsl sans cp var] |
---|
373 | et run avec version sur ce SVN planeto donne des resultats similaires |
---|
374 | |
---|
375 | ********** pratique |
---|
376 | |
---|
377 | A ioipsl |
---|
378 | A ioipsl/compile_ioipsl.bash |
---|
379 | A ioipsl/util 16 fichiers |
---|
380 | script et utilitaire pour compiler IOIPSL de facon independante |
---|
381 | il suffit d'executer ./compile_ioipsl.bash |
---|
382 | |
---|
383 | M arch/arch-AMD64_CICLAD.path |
---|
384 | si IOIPSL a ete compile avec la methode precedente, les bons |
---|
385 | PATH sont definis dans ce fichier [le NETCDF est aussi OK] |
---|
386 | |
---|
387 | M 000-README-svn |
---|
388 | mise a jour options "svn status" |
---|
389 | |
---|
390 | M mars/libf/phymars/meso_callkeys.h |
---|
391 | mise a jour mineure du fichier |
---|
392 | [ecri_phys etait defini mais pas dans la liste] |
---|
393 | |
---|
394 | |
---|
395 | ******************** |
---|
396 | **** commit_v54 **** |
---|
397 | ******************** |
---|
398 | |
---|
399 | |
---|
400 | -------------------------- |
---|
401 | LMD_LES_MARS: Tests with new physics in unified frame |
---|
402 | M 53 mesoscale/LMD_LES_MARS/modif_mars/module_model_constants.F |
---|
403 | M 53 mesoscale/LMD_LES_MARS/modif_mars/module_first_rk_step_part1.F |
---|
404 | M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F |
---|
405 | Modifications made to cope with new soil scheme |
---|
406 | |
---|
407 | M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/meso_physiq.F |
---|
408 | Commented aeronomars in old physics so that the model can compile with lots of grid points |
---|
409 | |
---|
410 | -------------------------- |
---|
411 | LMD_MM_MARS: Implement tracers with the new physics. Test with JBM runs with 5 tracers |
---|
412 | including radiatively active dust and water ice |
---|
413 | TODO: propagate reading "traceur.def" and few initialisations from |
---|
414 | testphys1d.F into module_lmd_driver.F |
---|
415 | |
---|
416 | M 53 mars/libf/phymars/dimradmars.h |
---|
417 | M 53 mars/libf/phymars/callradite.F |
---|
418 | Change it back to both dust and ice radiatively active |
---|
419 | |
---|
420 | M 53 mesoscale/LMD_MM_MARS/makemeso |
---|
421 | Allow to "fresh start" with the inclusion of "clean -a" in the script |
---|
422 | |
---|
423 | M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM |
---|
424 | M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F |
---|
425 | M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F |
---|
426 | M 53 mesoscale/LMD_MM_MARS/SRC/WRFV2/main/real_em.F |
---|
427 | M 53 mesoscale/LMD_MM_MARS/SIMU/runmeso |
---|
428 | Added the option mars==11 with 5 tracers: CO2 H2Ovap H2Oice DUST DUSTN |
---|
429 | -- Note that the order matters for H2Ovap and H2Oice |
---|
430 | |
---|
431 | LMD_MM_MARS: Create a folder to run LMD GCM with exact same physics as mesoscale |
---|
432 | |
---|
433 | A 0 mesoscale/LMDZ.MARS/libf/phymars/physiq.F |
---|
434 | All files in mesoscale/LMDZ.MARS/libf are links except this one |
---|
435 | --> Because specific WRITEDIAGFI commands are set to output what is needed PREP_MARS |
---|
436 | --> Something must be done so that qCO2, qdust, qdustN can be propagated too !!! |
---|
437 | |
---|
438 | A 0 mesoscale/LMDZ.MARS/libf_gcm |
---|
439 | A 0 mesoscale/LMDZ.MARS/in_lmdz_mars_newphys |
---|
440 | Two important links that must not be broken |
---|
441 | |
---|
442 | A 0 mesoscale/LMDZ.MARS/myGCM |
---|
443 | This folder is for runs. Only links here -- except temp files which are not synchronized. |
---|
444 | |
---|
445 | M 53 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/compile |
---|
446 | A 0 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/run.def |
---|
447 | A 0 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/traceur.def |
---|
448 | M 53 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/callphys.def |
---|
449 | A 0 mesoscale/LMD_MM_MARS/SIMU/in_lmdz_mars_newphys/myGCM/callphys.def.csttau |
---|
450 | Those files are needed to run the GCM in order to prepare initial & boundary |
---|
451 | conditions for the mesoscale... but basically it is the exact same GCM |
---|
452 | --> compile is a convenient script |
---|
453 | |
---|
454 | A 0 mesoscale/LMDZ.MARS/makegcm |
---|
455 | A 0 mesoscale/LMDZ.MARS/makegcm_g95 |
---|
456 | A 0 mesoscale/LMDZ.MARS/create_make_gcm |
---|
457 | This is supposed to be merged one day with files in trunk/mars/ |
---|
458 | |
---|
459 | A 0 mesoscale/LMDZ.MARS/myGCM/DEFS_JB |
---|
460 | Necessary files to restart a GCM run corresponding to new water cycle in JB's thesis |
---|
461 | |
---|
462 | **** commit_v65 **** |
---|
463 | -- 16/02/2011 -- EM , Mise a niveau par rapport à la version terrestre |
---|
464 | (LMDZ5V2.0-dev rev. 1487) |
---|
465 | * Mise à jour des scripts (terrestres) 'makegcm' et 'create_make_gcm' |
---|
466 | * Ajout du script 'makelmdz' (version "améliorée", en Bash, de makegcm) |
---|
467 | * Mise à jour des routines dans phylmd (sauf regr_lat_time_climoz_m.F) |
---|
468 | |
---|
469 | * disvert (dans dyn3d et dyn3dpar): passage au Fortran 90 |
---|
470 | * parallel.F90 (dyn3dpar): correction bug |
---|
471 | * etat0_netcdf.F90 (dyn3d et dyn3dpar) : mise à jour mineure |
---|
472 | * ce0l.F90 (dyn3dpar) : correction bug |
---|
473 | * abort_gcm.F (dyn3dpar) : correction bug |
---|
474 | * ugeostr.F90 (dyn3d et dyn3dpar) : passage au Fortran 90 |
---|
475 | * fluxstokenc_p.F (dyn3dpar) : correction bug |
---|
476 | * iniacademic.F90 (dyn3d et dyn3dpar) : passage au Fortran 90 |
---|
477 | * friction_p.F (dyn3dpar) : correction bug |
---|
478 | * infotrac.F90 (dyn3d et dyn3dpar) : correction bug mineur sur lecture traceurs |
---|
479 | * caladvtrac.F (dyn3d) : modifications cosmétiques |
---|
480 | |
---|
481 | ******************** |
---|
482 | **** commit_v76 **** etape importante d'une serie de commit sur les traceurs dans le mesoscale |
---|
483 | ******************** |
---|
484 | |
---|
485 | LMD_MM_MARS: remise a plat du traitement des traceurs terminee |
---|
486 | [en particulier en vue de l'utilisation nouvelle physique] |
---|
487 | --> fonctionnement concluant sur un cas Tharsis hydro 61x61x61 |
---|
488 | --> ne pas utiliser le cas FASTCASE trop instable avec traceurs |
---|
489 | --> reste a tester l'effet d'especes radiativement actives |
---|
490 | |
---|
491 | options dans namelist.input : |
---|
492 | mars = 0 ---> pas de traceurs |
---|
493 | mars = 1 ---> cycle de l'eau : water vapour + ice |
---|
494 | mars = 2 ---> cycle poussieres : un dust bin |
---|
495 | mars = 3 ---> cycle poussieres : dustq + dustn [NOUVELLE PHYS seulement] |
---|
496 | mars = 11 ---> cycle de l'eau + poussieres [1+3] [NOUVELLE PHYS seulement] |
---|
497 | NB: pour les deux derniers, reste un petit travail mineur sur qsurf |
---|
498 | (voir dans module_lmd_driver.F) |
---|
499 | il faut definir conjointement le callphys.def associe et relancer real.exe |
---|
500 | (sinon on transporte juste dynamiquement des 'dummy' traceurs) |
---|
501 | il n'est necessaire de recompiler que si le nombre total de traceurs change |
---|
502 | |
---|
503 | Fichiers a mettre a jour si l'on ajoute une option "mars" dans le Registry |
---|
504 | -------------------------------------------------------------------------- |
---|
505 | Attention suivant les inputs GCM, il faut peut etre egalement changer |
---|
506 | 1. readmeteo.F90 dans PREP_MARS et 2. METGRID.TBL dans WPS/metgrid |
---|
507 | |
---|
508 | M 75 mesoscale/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM |
---|
509 | ---> definition du scenario et de l'ordre des traceurs dans SCALAR |
---|
510 | |
---|
511 | M 75 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F |
---|
512 | ---> definition des interpolations verticales des champs de traceurs venant du GCM |
---|
513 | |
---|
514 | M 75 mesoscale/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F |
---|
515 | ---> definition du traitement a reserver aux bornes |
---|
516 | ---> --- dans les 4 scenarios precites, on passe aux bornes les champs du GCM |
---|
517 | [y compris QH2O_ICE contrairement a ce qui etait par defaut precedemment] |
---|
518 | ---> --- les lignes relatives a un flux nul aux bornes sont laissees a un cas hypothetique mars>50 |
---|
519 | |
---|
520 | M 75 mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F |
---|
521 | ---> definition de l'ordre correct des traceurs pour le passage a la physique |
---|
522 | ---> recuperation des tendances de la physique pour passage a la dynamique |
---|
523 | ---> NB: c'est dans ce module que sont presents des STOP si mars = 4-10 ou mars > 11 |
---|
524 | .... il faut donc modifier si l'on ajoute des options |
---|
525 | |
---|
526 | M 75 mesoscale/LMD_MM_MARS/SRC/WRFV2/main/real_em.F |
---|
527 | ---> definition et calcul des champs a appliquer aux bornes |
---|
528 | ---> generalise desormais, il n'y a plus qu'a ajouter d'eventuelles nouvelles options 'mars' |
---|
529 | ---> .... des modifications sont necessaires si on passe plus de 4 traceurs aux bornes |
---|
530 | |
---|
531 | M 75 mesoscale/LMD_MM_MARS/SIMU/runmeso |
---|
532 | ---> definition du bon nombre de traceurs dans la compilation de la physique puis l'execution |
---|
533 | [l'option mars est lue par le script dans namelist.input] |
---|
534 | |
---|
535 | Fichiers tests |
---|
536 | -------------- |
---|
537 | A 0 mesoscale/TESTS/newphys_tracers/* |
---|
538 | ---> pour la nouvelle physique (ici seulement les fichiers def) |
---|
539 | ---> toutes les options precitees ont ete testes avec succes a l'execution [pas de crash] |
---|
540 | ---> ... plausibilite physique verifiee rapidement, PAS d'ANALYSE APPROFONDIE pour le moment |
---|
541 | |
---|
542 | A 0 mesoscale/TESTS/LMD_MM_MARS_TESTCASE_water.tar.gz |
---|
543 | ---> pour l'ancienne physique (introduit precedemment mais n'avait pas ete synchronise) |
---|
544 | |
---|
545 | M 75 mars/libf/phymars/meso_dustlift.F |
---|
546 | NB: correction mineure, de facon a recuperer alpha_lift de initracer |
---|
547 | |
---|
548 | |
---|
549 | ******************** |
---|
550 | **** commit_v97 **** |
---|
551 | ******************** |
---|
552 | |
---|
553 | Serie de modifs Sebastien Lebonnois |
---|
554 | |
---|
555 | ** Correction bug |
---|
556 | **--------------- |
---|
557 | |
---|
558 | * Modif de bilan_dyn.F : |
---|
559 | - bug à la ligne 915: jjm et non jjp1... |
---|
560 | |
---|
561 | ** Modifs de dyn3d |
---|
562 | **--------------- |
---|
563 | |
---|
564 | * Modif de leapfrog.F : |
---|
565 | - ajout test pk monotone |
---|
566 | |
---|
567 | ** Sortie des opacites plutot que des aerosols dans phytitan |
---|
568 | **---------------------------------------------------------- |
---|
569 | |
---|
570 | * Modif de physiq.F : |
---|
571 | - ajout de variables et d'un common (pour variables venant de radtitan) |
---|
572 | |
---|
573 | * Modif de ini/write_histmth.h |
---|
574 | - sorties des opacites (thv* et khv*) |
---|
575 | - commentaire de la sortie des q* |
---|
576 | |
---|
577 | ** Reglage .def |
---|
578 | **--------------- |
---|
579 | |
---|
580 | * titan/physiq.def : |
---|
581 | - ajustement tx |
---|
582 | |
---|
583 | ** Sorties IOIPSL |
---|
584 | **----------------------- |
---|
585 | |
---|
586 | * Modif de phy[venus/titan]/physiq.F |
---|
587 | - suppression des REAL*4 lies aux sorties IOIPSL |
---|
588 | |
---|
589 | * Modif de phy[venus/titan]/[ini/write]_hist* |
---|
590 | - suppression des variables intermediaires en REAL*4 lies aux sorties IOIPSL |
---|
591 | |
---|
592 | ** nqmx remplace par nqtot |
---|
593 | **----------------------- |
---|
594 | |
---|
595 | * Modif de phytitan/ |
---|
596 | calchim.F,optci.F,optcv.F,pg3.F,physiq.F,phytrac.F, |
---|
597 | radlwsw.F,radtitan.F,rdf.F |
---|
598 | |
---|
599 | * Modif de phyvenus/ |
---|
600 | phytrac.F,testphys1d.F |
---|
601 | |
---|
602 | ** nettoyage |
---|
603 | **----------------------- |
---|
604 | |
---|
605 | * Suppression de phytitan/ |
---|
606 | calchim.dbleprec,calchim.spleprec |
---|
607 | |
---|
608 | ** calendriers |
---|
609 | **----------------------- |
---|
610 | |
---|
611 | * Modif de dyn3d[par]/gcm.F: |
---|
612 | - calendriers: il faut une option dans ioconf_calendar pour definir |
---|
613 | les calendriers titan et venus |
---|
614 | - pour definition jD_ref et jH_ref: avertissement pour Venus et Titan |
---|
615 | - avertissement pour separation en annees |
---|
616 | |
---|
617 | * Modif de deftank/[titan/venus]/run.def |
---|
618 | - ajout de calend => calend=titan ou venus |
---|
619 | - ajout de anneeref et raz_date |
---|
620 | |
---|
621 | ** Possibilite de tourner moins d'une journee |
---|
622 | **----------------------- |
---|
623 | |
---|
624 | * Modif de dyn3d[par]/leapfrog[_p].F |
---|
625 | - ajout de la possibilite de tourner moins d'1 jour |
---|
626 | |
---|
627 | * Modif de dyn3d[par]/control_mod.F90 |
---|
628 | - ajout de less1day et fractday |
---|
629 | |
---|
630 | * Modif de dyn3d[par]/conf_gcm.F |
---|
631 | - definition/lecture de less1day et fractday |
---|
632 | |
---|
633 | * Modif de deftank/venus/run.def |
---|
634 | - ajout de less1day et fractday |
---|
635 | |
---|
636 | ** suppression advtrac.h et control.h dans les phy* |
---|
637 | **----------------------- |
---|
638 | |
---|
639 | * Modif de : phy[venus/titan]/physiq.F |
---|
640 | - remplacement de #include "advtrac.h" par USE infotrac |
---|
641 | - remplacement de #include "control.h" par USE control_mod |
---|
642 | |
---|
643 | * Modif de : phy[venus/titan]/phytrac.F |
---|
644 | - remplacement de #include "advtrac.h" par USE infotrac |
---|
645 | - remplacement de #include "control.h" par USE control_mod |
---|
646 | |
---|
647 | * Modif de : phy[venus/titan]/phyredem.F |
---|
648 | - remplacement de #include "control.h" par USE control_mod |
---|
649 | |
---|
650 | ** ajustements lecture des .def |
---|
651 | **----------------------- |
---|
652 | |
---|
653 | * Modif de phy[titan/venus]/conf_phys.F90 |
---|
654 | - ajout de tous les parametres transferes de gcm.def dans physiq.def |
---|
655 | - utilisation de ecritphy: |
---|
656 | pour Venus, frequence de sortie dans mth et ins |
---|
657 | pour Titan, frequence de sortie dans ins |
---|
658 | - suppression de zzz |
---|
659 | |
---|
660 | * Modif de phy[titan/venus]/clesphys.h |
---|
661 | - mise a jour des declarations necessaires |
---|
662 | |
---|
663 | ** Elimination de clesphy0 |
---|
664 | **----------------------- |
---|
665 | |
---|
666 | * Modif de : phy[titan/venus]/phyetat0.F |
---|
667 | - elimination de clesphy0, longcles |
---|
668 | |
---|
669 | * Modif de : phy[titan/venus]/physiq.F |
---|
670 | - elimination de clesphy0, longcles |
---|
671 | |
---|
672 | * Modif de dyn3d[par]/ |
---|
673 | calfis.F, ce0l.F90, conf_gcm.F, gcm.F, leapfrog.F |
---|
674 | - elimination de clesphy0, longcles |
---|
675 | |
---|
676 | ** Ajustement format .h |
---|
677 | **----------------------- |
---|
678 | |
---|
679 | * Modif de phy[titan/venus]/clesphys.h (clesphys.inc elimine) |
---|
680 | |
---|
681 | * Modif de phy[titan/venus]/YOMCST.h (YOMCST.inc elimine) |
---|
682 | |
---|
683 | * Elimination de tous les YOETHF dans phy[titan/venus] et nettoyage de suphec.F |
---|
684 | |
---|
685 | * Elimination de tous les FCTTRE dans phy[titan/venus] |
---|
686 | |
---|
687 | * Modif de : phytitan/comorbit.h (comorbit.inc elimine) |
---|
688 | |
---|
689 | ** Controle de la remise a zero des varc |
---|
690 | **----------------------- |
---|
691 | |
---|
692 | * Modif de dyn3d[par]/sortvarc.F |
---|
693 | - ajout firstcal |
---|
694 | - ajout flag pour controle (resetvarc) |
---|
695 | - correction pour avoir le moment cinetique total absolu |
---|
696 | |
---|
697 | * Modif de dyn3d[par]/conf_gcm.F |
---|
698 | - ajout flag pour controle (resetvarc) |
---|
699 | |
---|
700 | * Modif de dyn3d[par]/control_mod.F90 |
---|
701 | - ajout de resetvarc |
---|
702 | |
---|
703 | * Modif de deftank/venus/run.def |
---|
704 | - ajout de resetvarc |
---|
705 | |
---|
706 | * Modif de dyn3d[par]/sortvarc0.F |
---|
707 | - correction pour avoir le moment cinetique total absolu |
---|
708 | |
---|
709 | |
---|
710 | ********************* |
---|
711 | **** commit_v101 **** |
---|
712 | ********************* |
---|
713 | |
---|
714 | Serie de modifs Sebastien Lebonnois |
---|
715 | |
---|
716 | Compilation de phyvenus, ça passe ! |
---|
717 | Modifs propagees a phytitan (sans compilation... a cause des routines C) |
---|
718 | dans la v102. |
---|
719 | |
---|
720 | Reste des chantiers a finir: |
---|
721 | - discretisation verticale |
---|
722 | - sponge layer |
---|
723 | |
---|
724 | ********************* |
---|
725 | **** commit_v104 **** |
---|
726 | ********************* |
---|
727 | |
---|
728 | SLebonnois: modification de makelmdz et create_make_gcm pour pouvoir |
---|
729 | compiler la chimie titan. Pas de raison que ca gene les autres. |
---|
730 | Dans cette version, les compilations de Venus et Titan fonctionnent. |
---|
731 | |
---|
732 | Phytitan: modifications pour pouvoir compiler correctement. |
---|
733 | Il ne manque plus que physiq.F a faire. |
---|
734 | --- |
---|
735 | Fait: |
---|
736 | Avec la v105, la compilation Titan se fait sans problemes. |
---|
737 | |
---|
738 | **----------- |
---|
739 | (commit 106: ajustements dans phytitan) |
---|
740 | **----------- |
---|
741 | - Sorties des opacites (*_hist*.h) |
---|
742 | - g(z) pris en compte (physiq.F, zzlay) |
---|
743 | - correction topo (physiq.F, zzlev(i,1)) |
---|
744 | - correction lell.F et lell_light.F (et aussi radtitan.F du coup) |
---|
745 | - remise a jour des sorties |
---|
746 | - suphec.F : RDAY=1.37889e6 et non RSIDAY |
---|
747 | **----------- |
---|
748 | |
---|
749 | ********************* |
---|
750 | **** commit_v108 **** |
---|
751 | ********************* |
---|
752 | |
---|
753 | Modifs Sebastien Lebonnois dans dyn3d[par] pour sponge layer et dissip horizontale. |
---|
754 | |
---|
755 | ** Correction bug |
---|
756 | **--------------- |
---|
757 | |
---|
758 | * Modif de dyn3d[par]/calfis[_p].F |
---|
759 | => ztfi doublement incremente... |
---|
760 | |
---|
761 | ** Sponge layer |
---|
762 | **------------- |
---|
763 | |
---|
764 | * Modif de dyn3d[par]/leapfrog[_p].F |
---|
765 | Utilisation de d*top dans appel a top_bound.F |
---|
766 | Initialisation des d* (dis, fi, top) |
---|
767 | |
---|
768 | * Modif de dyn3d[par]/top_bound[_p].F |
---|
769 | Implementation de mode_top_bound |
---|
770 | Pas d'incrementation |
---|
771 | |
---|
772 | * Modif de dyn3d[par]/conf_gcm.F |
---|
773 | Parametres de top_bound |
---|
774 | |
---|
775 | * Modif de dyn3d[par]/comconst.h |
---|
776 | Parametres de top_bound |
---|
777 | |
---|
778 | * Modif de deftank/[titan/venus]/gcm.def |
---|
779 | Ajout ok_strato |
---|
780 | Ajout parametres de top_bound |
---|
781 | |
---|
782 | * Ecriture doc |
---|
783 | Ajout top_bound.tex (et .pdf) dans la documentation. |
---|
784 | |
---|
785 | ** Dissipation horizontale |
---|
786 | **------------------------ |
---|
787 | |
---|
788 | * Modif de dyn3d[par]/inidissip.F |
---|
789 | passage a dissip_fac_mid automatique |
---|
790 | passage a dissip_fac_up sur modele martien si ok_strato |
---|
791 | |
---|
792 | * Modif de dyn3d[par]/conf_gcm.F |
---|
793 | Parametres dissip_* |
---|
794 | |
---|
795 | * Modif de dyn3d[par]/comconst.h |
---|
796 | Parametres dissip_* |
---|
797 | |
---|
798 | * Modif de deftank/[titan/venus]/gcm.def |
---|
799 | Ajout parametres dissip_* !!! REGLAGES A FAIRE !!! |
---|
800 | |
---|
801 | * Ecriture doc |
---|
802 | Ajout dissip_horiz.tex (et .pdf) dans la documentation. |
---|
803 | |
---|
804 | Reste la discretisation verticale a voir. |
---|
805 | |
---|