- Timestamp:
- Feb 24, 2011, 4:14:33 PM (14 years ago)
- Location:
- trunk
- Files:
-
- 10 edited
-
mars/libf/phymars/callradite.F (modified) (1 diff)
-
mars/libf/phymars/dimradmars.h (modified) (1 diff)
-
mars/libf/phymars/meso_physiq.F (modified) (2 diffs)
-
mars/libf/phymars/vdifc.F (modified) (1 diff)
-
mesoscale/LMDZ.MARS.new/myGCM/callphys.def (modified) (1 diff)
-
mesoscale/LMD_MM_MARS/SRC/ARWpost/idl/profile.pro (modified) (1 diff)
-
mesoscale/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM (modified) (2 diffs)
-
mesoscale/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F (modified) (2 diffs)
-
mesoscale/NOTES.txt (modified) (3 diffs)
-
mesoscale/TESTS/newphys_tracers/namelist.input (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/mars/libf/phymars/callradite.F
r55 r77 263 263 c PLEASE MAKE SURE that you set up the right number of 264 264 c scatterers in dimradmars.h (naerkind); 265 name_iaer(1) = "dust_conrath" !! poussiere classique266 cname_iaer(1) = "dust_doubleq"265 c name_iaer(1) = "dust_conrath" !! poussiere classique 266 name_iaer(1) = "dust_doubleq" 267 267 cc name_iaer(2) = "dust_submicron" !! JB: experimental 268 cname_iaer(2) = "h2o_ice"268 name_iaer(2) = "h2o_ice" 269 269 c ---------------------------------------------------------- 270 270 -
trunk/mars/libf/phymars/dimradmars.h
r55 r77 27 27 ! (ex: naerkind=2 if you use one dust mode and active ice ...) 28 28 integer naerkind 29 parameter (naerkind=1)30 !parameter (naerkind=2)29 ! parameter (naerkind=1) 30 parameter (naerkind=2) 31 31 32 32 ! Reference wavelengths used to compute reference optical depth (m) -
trunk/mars/libf/phymars/meso_physiq.F
r73 r77 192 192 REAL ust(ngridmx) !! pour LES avec isfflx!=0 193 193 LOGICAL flag_LES !! pour LES avec isfflx!=0 194 REAL qsurf last(ngridmx) !! pour diagnostics194 REAL qsurfice(ngridmx) !! pour diagnostics 195 195 real alpha,lay1 ! coefficients for building layers 196 196 integer iloop … … 1577 1577 wco2ice(ig) = co2ice(ig) !! co2 ice 1578 1578 1579 !!! TEMP TEMP TEMP TEMP TEMP TEMP TEMP1580 1579 !!! specific to WRF WRF WRF 1581 1580 !!! just to output water ice on surface 1582 !!! [it might not be water ice on surface but OK]1583 1581 !!! uncomment the Registry entry 1584 qsurflast(ig) = qsurf(ig,nqmx)1582 IF (igcm_h2o_ice .ne. 0) qsurfice(ig) = qsurf(ig,igcm_h2o_ice) 1585 1583 1586 1584 enddo -
trunk/mars/libf/phymars/vdifc.F
r38 r77 474 474 else if (doubleq) then 475 475 do ig=1,ngrid 476 !!! soulevement constant 476 477 pdqsdif(ig,igcm_dust_mass) = 477 478 & -alpha_lift(igcm_dust_mass) -
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.
