Changeset 77 for trunk/mesoscale
- Timestamp:
- Feb 24, 2011, 4:14:33 PM (14 years ago)
- Location:
- trunk/mesoscale
- Files:
-
- 6 edited
-
LMDZ.MARS.new/myGCM/callphys.def (modified) (1 diff)
-
LMD_MM_MARS/SRC/ARWpost/idl/profile.pro (modified) (1 diff)
-
LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM (modified) (2 diffs)
-
LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F (modified) (2 diffs)
-
NOTES.txt (modified) (3 diffs)
-
TESTS/newphys_tracers/namelist.input (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/mesoscale/LMDZ.MARS.new/myGCM/callphys.def
r54 r77 1 link DEFS_JB/callphys.def 1 link DEFS_JB/callphys.def.orig -
trunk/mesoscale/LMD_MM_MARS/SRC/ARWpost/idl/profile.pro
r11 r77 124 124 ytitle=ytitle,$ 125 125 title=title_user,$ 126 /ylog127 ;subtitle=mention126 ;/ylog 127 subtitle=mention 128 128 129 129 ;------------- -
trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM
r76 r77 92 92 state real MTOT ij misc 1 - rhd "MTOT" "TOTAL MASS WATER VAPOR" "kg m-2" #SAVEMARS2 mtot 93 93 state real ICETOT ij misc 1 - rhd "ICETOT" "TOTAL MASS WATER ICE" "kg m-2" #SAVEMARS2 icetot 94 state real RAVE ij misc 1 - r hd"RAVE" "MEAN ICE RADIUS" "m" #SAVEMARS2 rave95 state real RICE ikj misc 1 - r d"RICE" "ICE RADIUS" "m" #SAVEMARS3 rice94 state real RAVE ij misc 1 - rd "RAVE" "MEAN ICE RADIUS" "m" #SAVEMARS2 rave 95 state real RICE ikj misc 1 - rhd "RICE" "ICE RADIUS" "m" #SAVEMARS3 rice 96 96 state real HR_SW ikj misc 1 - rd "HR_SW" "HEATING RATE SW" "K/s" #SAVEMARS3 zdtsw 97 97 state real HR_LW ikj misc 1 - rd "HR_LW" "HEATING RATE LW" "K/s" #SAVEMARS3 zdtlw 98 98 state real HR_SH ikj misc 1 - rd "HR_SH" "HEATING RATE sens. heat" "K/s" #SAVEMARS3 zdtdif 99 state real QSURF ij misc 1 - rhd "QSURF" "TRACER AT SURFACE" "kg m-2" #SAVEMARS2 qsurflast 99 state real QSURFICE ij misc 1 - rhd "QSURFICE" "WATER ICE AT SURFACE" "kg m-2" #SAVEMARS2 qsurfice 100 state real RDUST ikj misc 1 - rhd "RDUST" "DUST RADIUS" "m" #SAVEMARS3 rdust 101 state real HR_NIR ikj misc 1 - rd "HR_NIR" "HEATING RATE nirco2" "K/s" #SAVEMARS3 zdtnirco2 102 state real HR_NLTE ikj misc 1 - rd "HR_NLTE" "HEATING RATE nlte" "K/s" #SAVEMARS3 zdtnlte 100 103 #### 101 104 #### … … 1237 1240 rconfig logical pd_moist namelist,dynamics max_domains .false. rh "pd_moist" "positive-definite RK3 transport switch" "" 1238 1241 rconfig logical pd_chem namelist,dynamics max_domains .false. rh "pd_chem" "positive-definite RK3 transport switch" "" 1239 #rconfig logical pd_scalar namelist,dynamics max_domains .true. rh "pd_scalar" "positive-definite RK3 transport switch" ""1240 rconfig logical pd_scalar namelist,dynamics max_domains .false. rh "pd_scalar" "positive-definite RK3 transport switch" ""1242 rconfig logical pd_scalar namelist,dynamics max_domains .true. rh "pd_scalar" "positive-definite RK3 transport switch" "" 1243 #rconfig logical pd_scalar namelist,dynamics max_domains .false. rh "pd_scalar" "positive-definite RK3 transport switch" "" 1241 1244 rconfig logical pd_tke namelist,dynamics max_domains .false. rh "pd_tke" "positive-definite RK3 transport switch" "" 1242 1245 rconfig logical top_radiation namelist,dynamics max_domains .false. rh "top_radiation" "" "" -
trunk/mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
r76 r77 832 832 !! ----- [H2O ice is tracer nqmx in qsurf in LMD physics] 833 833 CASE(2) 834 qsurf_val(1)= 10000. !! infinite amount of dust, basically834 qsurf_val(1)=0. !! not coupled with lifting for the moment [non remobilise] 835 835 #ifdef NEWPHYS 836 836 CASE(3) … … 839 839 qsurf_val(1)=0. 840 840 qsurf_val(2)=MARS_WICE(i,j) !! logique avec wtnom(2) = 'h2o_ice' defini ci-dessus 841 !!!! petit probleme... alors QSURF sorti dans wrfout n'est pas h2o_ice ! 842 qsurf_val(3)=10000. !! infinite amount of dust, basically 841 qsurf_val(3)=0. !! not coupled with lifting for the moment [non remobilise] 843 842 #else 844 843 CASE(3:) -
trunk/mesoscale/NOTES.txt
r73 r77 1 1 2 changer le renormalisation dans aeropacity ???? 3 on ne laisse pas aerosol comme le lifting veut qu'il soit ! 4 tenter des taux de soulevement pour que taudust_tmp soit les obs 5 en prescivant une dust bomb fixe d opacite, on aura au moins la structure verticale 6 7 tester traceurs radiativement actifs avec la nouvelle physique ????? 8 9 A FAIRE: PB LES sur iDATAPLEX (les points HFX nuls) (pas de soucis sur ciclad) 10 METTRE SUR LE svn LA BASE d'ETATS INITIAUX ???? 2 11 3 12 more than 4 procs w/ nest ??? y reflechir 4 ----------------------------------------------------------------------- 5 -- si possible comment determiner taille ? 6 nproc doit diviser e_we-1 (1er nest) 7 grid_ratio doit diviser e_we-1 +4 (1er nest) 8 soit e_we=ye+1 9 grid_ratio divise ye+4 et nproc divise ye 10 soit nproc=8, ye=8*i 11 ainsi il existe j tel que 8i + 4 = 3j ou encore 4*[2i+1] = 3j 12 verifie par exemple si 2i+1 est multiple de 3 13 il suffit donc de trouver un multiple impair de 3 et de deduire i 14 par exemple 2i+1=33 >>>> i=16 15 >>>> e_we = 129 pour le 1er nest (et ajouter 4 pour les suivants) 16 ------------------------------------------------------------------------ 13 ----------------------------------------------------------------------- 14 -- si possible comment determiner taille ? 15 nproc doit diviser e_we-1 (1er nest) 16 grid_ratio doit diviser e_we-1 +4 (1er nest) 17 soit e_we=ye+1 18 grid_ratio divise ye+4 et nproc divise ye 19 soit nproc=8, ye=8*i 20 ainsi il existe j tel que 8i + 4 = 3j ou encore 4*[2i+1] = 3j 21 verifie par exemple si 2i+1 est multiple de 3 22 il suffit donc de trouver un multiple impair de 3 et de deduire i 23 par exemple 2i+1=33 >>>> i=16 24 >>>> e_we = 129 pour le 1er nest (et ajouter 4 pour les suivants) 25 ------------------------------------------------------------------------ 26 27 ne pas utiliser le FASTCASE avec traceurs (instabilites en haut) 28 ces instabilites sont cependant reglees si on augmente radt a 10 par exemple 17 29 18 30 pour le cycle de l'eau c'est OK de regler caps=F dans le mesoscale 19 31 sauf si on commence a devoiler la calotte permanente nord 20 32 21 NE SERAIT-CE PAS MIEUX DE TOUT TRANSMETTRE AUX BORNES ???22 tous les traceurs, pas seulement vapor33 NE SERAIT-CE PAS MIEUX DE TOUT TRANSMETTRE AUX BORNES ??? 34 tous les traceurs, pas seulement vapor 23 35 24 36 25 - attention il faut les trois MARS sinon il s arrete sans message clair 26 - attention a ne pas lancer le modele s il est deja lance 37 - attention il faut les trois MARS sinon il s arrete sans message clair 38 - attention a ne pas lancer le modele s il est deja lance 39 - important que pd_scalar soit a T ... le mettre par defaut ???? 27 40 28 41 … … 32 45 - il faut tester le nest pour verifier les lignes trop longues 33 46 34 (ok) lier gr_fi_dyn qui est dans dyn3d35 (ok) regler le pb du nouveau readtesassim (ou alors le lier tout simplement ou36 l'appeler meso_readtesassim)37 (ok) regler le pb meso_dustlift (le lier dans makemeso comme point precedent)38 (car le souci c que dustlift est appele dans vdifc)47 (ok) lier gr_fi_dyn qui est dans dyn3d 48 (ok) regler le pb du nouveau readtesassim (ou alors le lier tout simplement ou 49 l'appeler meso_readtesassim) 50 (ok) regler le pb meso_dustlift (le lier dans makemeso comme point precedent) 51 (car le souci c que dustlift est appele dans vdifc) 39 52 40 RESTE a ADAPTER le LES a la NOUVELLE PHYSIQUE41 il y a normalement peu a faire42 reste a faire egalement le -DNEWPHYS pour le LES53 RESTE a ADAPTER le LES a la NOUVELLE PHYSIQUE 54 il y a normalement peu a faire 55 reste a faire egalement le -DNEWPHYS pour le LES 43 56 44 attention pb d'affichage des valeurs dans le fichier texte avec LES ???45 bien que les valeurs du fichier soient tout a fait raisonnables46 ... n'est-ce pas un effet de bord cache ????57 attention pb d'affichage des valeurs dans le fichier texte avec LES ??? 58 bien que les valeurs du fichier soient tout a fait raisonnables 59 ... n'est-ce pas un effet de bord cache ???? 47 60 48 61 49 apres fusion, le LES est reconnu par module_lmd_driver lorsque diff_opt=2 km_opt=262 apres fusion, le LES est reconnu par module_lmd_driver lorsque diff_opt=2 km_opt=2 50 63 51 64 52 -attention PB si on ne sort pas HFX et USTM (note dans le Registry)53 -il faut run.def nouvelle physique [c est meme ce qui est utilise par runmeso]54 - IL FAUT SE PENCHER SUR LE FAIT QU'ON INDIQUE q2val=0 dans lmd_driver ....65 -attention PB si on ne sort pas HFX et USTM (note dans le Registry) 66 -il faut run.def nouvelle physique [c est meme ce qui est utilise par runmeso] 67 - IL FAUT SE PENCHER SUR LE FAIT QU'ON INDIQUE q2val=0 dans lmd_driver .... 55 68 69 ----------------------- 56 70 ATTENTION NOUVELLE PHYSIQUE 57 71 Oui, c'est quelque chose qu'il faut probablement changer partout … … 81 95 >> set optim90=" -fast" 82 96 >> set optimtru90=" -fast -c -Mfree " 97 ------------------------------ 98 99 100 - attention a cp et R, normaliser une bonne fois pour toutes 101 - il manque sur le SVN les cas idealises 102 - il manque sur le SVN les scripts MPI 103 - il faut recompiler les librairies NETCDF 104 - mettre la nouvelle physique 105 - mettre les DEF du meso-echelle 106 107 - modele ok sur auric 108 - modele pas ok sur ciclad avec pgf2010, erreur inedite un seul module manquant 109 - modele LES OK sur ciclad 110 - modele LES ok sur auric 111 112 24/01/2011 113 tests g95 en 64bits natif sur systeme Linux 114 -- modifications de makemeso, tests 115 -- tout est OK sauf les libraires NETCDF, probleme d'underscore 116 -- OK avec libraires maison compilees avec g95 standard sur flores [et tourne OK] 83 117 84 118 85 119 86 - attention a cp et R, normaliser une bonne fois pour toutes 87 - il manque sur le SVN les cas idealises 88 - il manque sur le SVN les scripts MPI 89 - il faut recompiler les librairies NETCDF 90 - mettre la nouvelle physique 91 - mettre les DEF du meso-echelle 92 93 - modele ok sur auric 94 - modele pas ok sur ciclad avec pgf2010, erreur inedite un seul module manquant 95 - modele LES OK sur ciclad 96 - modele LES ok sur auric 97 98 24/01/2011 99 tests g95 en 64bits natif sur systeme Linux 100 -- modifications de makemeso, tests 101 -- tout est OK sauf les libraires NETCDF, probleme d'underscore 102 -- OK avec libraires maison compilees avec g95 standard sur flores [et tourne OK] 103 104 105 106 mpi_64_pgf7_ncdf4_mpi1.2.txt 107 - probleme lors de la compilation de solve_em : LINUX runs out of memory [huchard] 108 - IL FAUT COMPILER SUR auric 109 110 nougaro est lent a la compilation, utiliser surtout auric 120 mpi_64_pgf7_ncdf4_mpi1.2.txt 121 - probleme lors de la compilation de solve_em : LINUX runs out of memory [huchard] 122 - IL FAUT COMPILER SUR auric 123 nougaro est lent a la compilation, utiliser surtout auric 111 124 112 125 113 126 114 127 115 128 ______________________________________________________ 116 129 117 130 -
trunk/mesoscale/TESTS/newphys_tracers/namelist.input
r76 r77 8 8 end_day = 03, !! (p1) End Martian Day 9 9 end_hour = 01, !! (p1) End Martian Hour (at longitude 0) 10 history_interval = 10, !! Frequency of outputs (37 --> 3700s = 1 Martian hour)10 history_interval = 37, !! Frequency of outputs (37 --> 3700s = 1 Martian hour) 11 11 frames_per_outfile = 24, !! Size of time dimension in files 12 12 restart = .false. !! (*) Output restart files ?
Note: See TracChangeset
for help on using the changeset viewer.
