- Timestamp:
- Nov 3, 2011, 2:28:57 PM (14 years ago)
- Location:
- trunk
- Files:
-
- 8 added
- 1 deleted
- 9 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_les.F
r330 r341 4 4 !! sensible heat flux in W/m2 5 5 6 hfx(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)6 sensheat(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig) 7 7 8 8 !! u star in similarity theory in m/s 9 ust (ig) = 0.49 ustar(ig) = 0.4 10 10 . * sqrt( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) ) 11 11 . / log( 1.E+0 + zzlay(ig,1)/z0_default ) … … 16 16 DO ig=1,ngrid 17 17 18 ! New SL parametrization, correct formulation for hfx:18 ! New SL parametrization, correct formulation for sensheat : 19 19 20 hfx(ig) = (pplay(ig,1)/(r*pt(ig,1)))*cpp20 sensheat(ig) = (pplay(ig,1)/(r*pt(ig,1)))*cpp 21 21 & *sqrt(pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) 22 22 & + (1.2*wmax_th(ig))**2) 23 23 & *zcdh(ig)*(tsurf(ig)-zh(ig,1)) 24 24 25 ! New SL parametrization, ust is more accurately computed in vdif_cd :26 ust (ig) = sqrt(zcdv(ig)*25 ! New SL parametrization, ustar is more accurately computed in vdif_cd : 26 ustar(ig) = sqrt(zcdv(ig)* 27 27 & (pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) + (0.7*wmax_th(ig))**2) 28 28 & ) … … 32 32 endif !of if callrichsl 33 33 34 ! write (*,*) 'PHYS HFX cp zdts', hfx(100), zflubid(100),34 ! write (*,*) 'PHYS HFX cp zdts', sensheat(100), zflubid(100), 35 35 ! . capcal(100), 36 36 ! . zdtsdif(100) 37 ! write (*,*) 'PHYS UST', ust (100)37 ! write (*,*) 'PHYS UST', ustar(100) 38 38 39 39 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! -
trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_var.F
r308 r341 10 10 REAL output_tab2d(ngridmx,n2d) 11 11 REAL output_tab3d(ngridmx,nlayer,n3d) 12 REAL hfx(ngridmx) !! pour LES avec isfflx!=013 REAL ust (ngridmx) !! pour LES avec isfflx!=012 REAL sensheat(ngridmx) !! pour LES avec isfflx!=0 13 REAL ustar(ngridmx) !! pour LES avec isfflx!=0 14 14 LOGICAL flag_LES !! pour LES avec isfflx!=0 15 15 REAL qsurfice(ngridmx) !! pour diagnostics -
trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM
r170 r341 995 995 state real CDA ij misc 1 - r "CDA" "DRAG COEFF AT LOWEST MODEL LVL" "" 996 996 state real UST ij misc 1 - r "UST" "U* IN SIMILARITY THEORY" "m s-1" 997 state real USTM ij misc 1 - rh "USTM" "U* IN SIMILARITY THEORY WITHOUT VCONV" "m s-1" #SAVEMARS2 ust 997 state real USTM ij misc 1 - rh "USTM" "U* IN SIMILARITY THEORY WITHOUT VCONV" "m s-1" #SAVEMARS2 ustar 998 998 i1 real HOL ij misc 1 - - "HOL" "PBL HEIGHT OVER MONIN-OBUKHOV LENGTH" "" 999 999 state real RMOL ij misc 1 - r "RMOL" "1./Monin Ob. Length" "" … … 1002 1002 state real CAPG ij misc 1 - r "CAPG" "HEAT CAPACITY FOR SOIL" "J K-1 m-3" 1003 1003 state real THC ij misc 1 - r "THC" "THERMAL INERTIA" "Cal cm-1 K-1 s-0.5" 1004 state real HFX ij misc 1 - rh "HFX" "UPWARD HEAT FLUX AT THE SURFACE" "W m-2" #SAVEMARS2 hfx1004 state real HFX ij misc 1 - rh "HFX" "UPWARD HEAT FLUX AT THE SURFACE" "W m-2" #SAVEMARS2 sensheat 1005 1005 state real QFX ij misc 1 - r "QFX" "UPWARD MOISTURE FLUX AT THE SURFACE" "kg m-2 s-1" 1006 1006 state real LH ij misc 1 - r "LH" "LATENT HEAT FLUX AT THE SURFACE" "W m-2" -
trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM
r324 r341 110 110 state real HFMAX_TH ij misc 1 - rd "HFMAX_TH" "MAXIMUM TH HEAT FLUX" "m.K/s" #SAVEMARS2 hfmax_th 111 111 state real WMAX_TH ij misc 1 - rd "WMAX_TH" "MAXIMUM TH VERTICAL VELOCITY" "m/s" #SAVEMARS2 wmax_th 112 state real Z0SET ij misc 1 - r hd"Z0SET" "SET SURFACE ROUGHNESS" "m" #SAVEMARS2 z0112 state real Z0SET ij misc 1 - rd "Z0SET" "SET SURFACE ROUGHNESS" "m" #SAVEMARS2 z0 113 113 114 114 #### … … 946 946 ##### LES mettre h pour les deux variables ci-dessous 947 947 state real UST ij misc 1 - r "UST" "U* IN SIMILARITY THEORY" "m s-1" 948 state real USTM ij misc 1 - rh "USTM" "U* IN SIMILARITY THEORY - LES MODE" "m s-1" #SAVEMARS2 ust 948 state real USTM ij misc 1 - rh "USTM" "U* IN SIMILARITY THEORY - LES MODE" "m s-1" #SAVEMARS2 ustar 949 949 i1 real HOL ij misc 1 - - "HOL" "PBL HEIGHT OVER MONIN-OBUKHOV LENGTH" "" 950 950 state real RMOL ij misc 1 - r "RMOL" "1./Monin Ob. Length" "" … … 953 953 state real CAPG ij misc 1 - r "CAPG" "HEAT CAPACITY FOR SOIL" "J K-1 m-3" 954 954 state real THC ij misc 1 - r "THC" "THERMAL INERTIA" "Cal cm-1 K-1 s-0.5" 955 state real HFX ij misc 1 - rh "HFX" "UPWARD HEAT FLUX AT THE SURFACE" "W m-2" #SAVEMARS2 hfx955 state real HFX ij misc 1 - rh "HFX" "UPWARD HEAT FLUX AT THE SURFACE" "W m-2" #SAVEMARS2 sensheat 956 956 state real QFX ij misc 1 - r "QFX" "UPWARD MOISTURE FLUX AT THE SURFACE" "kg m-2 s-1" 957 957 state real LH ij misc 1 - r "LH" "LATENT HEAT FLUX AT THE SURFACE" "W m-2" -
trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/physiq.F
r333 r341 178 178 REAL sl_ra, sl_di0 179 179 REAL sky 180 REAL hfx(ngridmx) !! pour LES avec isfflx!=0181 REAL ust (ngridmx) !! pour LES avec isfflx!=0180 REAL sensheat(ngridmx) !! pour LES avec isfflx!=0 181 REAL ustar(ngridmx) !! pour LES avec isfflx!=0 182 182 LOGICAL flag_LES !! pour LES avec isfflx!=0 183 183 REAL qsurflast(ngridmx) !! pour diagnostics … … 925 925 DO ig=1,ngrid 926 926 !! sensible heat flux in W/m2 927 hfx(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)927 sensheat(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig) 928 928 !! u star in similarity theory in m/s 929 ust (ig) = 0.4929 ustar(ig) = 0.4 930 930 . * sqrt( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) ) 931 931 . / log( 1.E+0 + zzlay(ig,1)/z0 ) … … 944 944 write (*,*) '************************************************' 945 945 946 ! DO ig=1,ngrid947 ! !! sensible heat flux in W/m2948 ! hfx(ig) = zflubid(ig)-capcal(ig)*zdtsdif(ig)949 ! !! u star in similarity theory in m/s950 ! ust(ig) = 0.4951 ! . * sqrt( pu(ig,1)*pu(ig,1) + pv(ig,1)*pv(ig,1) )952 ! . / log( 1.E+0 + zzlay(ig,1)/z0 )953 946 ! 954 947 DO l=1,nlayer … … 966 959 !write (*,*) 'GRD ',fluxgrd(igout) 967 960 !write (*,*) 'dTs/dt ',capcal(igout)*zdtsurf(igout) 968 !write (*,*) 'HFX ', hfx(igout)969 !write (*,*) 'UST ', ust(igout)970 961 ENDIF 971 962 !!! LES LES -
trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/makegcm_mpi
r335 r341 177 177 ## set opt_link=" -Mfree -L/usr/local/pgi/linux86/lib -lpgf90 -lpgftnrtl -lpghpf -lpghpf2 -L$NCDFLIB -lnetcdf -Wl,-Bstatic" 178 178 # set opt_link=" -Mlarge_arrays -mcmodel=medium -Mfree -L/usr/local/pgi/linux86/lib -lpgf90 -lpgftnrtl -lpghpf -lpghpf2 -L$NCDFLIB -lnetcdf" 179 set optim90=" -O2 -I. -Munroll -Mcache_align" 180 set optimtru90=" -O2 -I. -Munroll -Mcache_align" 179 set optim90="-O2 -I. -Munroll -Mcache_align" # -Mbounds" 180 set optimtru90="-O2 -I. -Munroll -Mcache_align" # -Mbounds" 181 #set optim90=""#-Mvect=noaltcode -Msmartalloc -Mprefetch=distance:8 -Mfprelaxed" #-fast -I." 182 #set optimtru90=""#-Mvect=noaltcode -Msmartalloc -Mprefetch=distance:8 -Mfprelaxed"#-fast -I." 181 183 #set optim90="-O3 -I. -tp=nehalem-64 -fastsse -Mscalarsse -Mvect=sse -Mipa=fast" 182 184 #set optimtru90="-O3 -I. -tp=nehalem-64 -fastsse -Mscalarsse -Mvect=sse -Mipa=fast" 183 set opt_link=" -L$NCDFLIB -lnetcdf -Bstatic"185 set opt_link=" -L$NCDFLIB -lnetcdf" # -Bstatic" 184 186 set mod_loc_dir="" 185 187 set mod_suffix=mod … … 362 364 set optim90="$optim90"" -G1 " 363 365 else if $LINUX then 364 set optim90="-g -Mbounds "366 set optim90="-g -Mbounds -Ktrap=fp" 365 367 # set optim90="-mcmodel=medium" 366 368 else -
trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F
r324 r341 412 412 ENDIF 413 413 414 !!!! for 'subgrid' temporal diagnostics415 test2 = MODULO(elaps,history_interval*100.)414 !!!!! for 'subgrid' temporal diagnostics 415 !test2 = MODULO(elaps,history_interval*100.) 416 416 417 417 !!!******!! … … 478 478 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 479 479 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 480 IF (test2.EQ.0) THEN481 print *, 'compute stats'482 print *, 'RESET'483 uave = uave*0.484 vave = vave*0.485 tave = tave*0.486 wave = wave*0.487 ustd = ustd*0.488 vstd = vstd*0.489 tstd = tstd*0.490 wstd = wstd*0.491 ENDIF492 uave = uave + u * dt / (float(history_interval)*100.)493 vave = vave + v * dt / (float(history_interval)*100.)494 tave = tave + th * dt / (float(history_interval)*100.)495 wave = wave + w * dt / (float(history_interval)*100.)496 ustd = ustd + u * u * dt / (float(history_interval)*100.)497 vstd = vstd + v * v * dt / (float(history_interval)*100.)498 tstd = tstd + th * th * dt / (float(history_interval)*100.)499 wstd = wstd + w * w * dt / (float(history_interval)*100.)480 !IF (test2.EQ.0) THEN 481 ! print *, 'compute stats' 482 ! print *, 'RESET' 483 ! uave = uave*0. 484 ! vave = vave*0. 485 ! tave = tave*0. 486 ! wave = wave*0. 487 ! ustd = ustd*0. 488 ! vstd = vstd*0. 489 ! tstd = tstd*0. 490 ! wstd = wstd*0. 491 !ENDIF 492 ! uave = uave + u * dt / (float(history_interval)*100.) 493 ! vave = vave + v * dt / (float(history_interval)*100.) 494 ! tave = tave + th * dt / (float(history_interval)*100.) 495 ! wave = wave + w * dt / (float(history_interval)*100.) 496 ! ustd = ustd + u * u * dt / (float(history_interval)*100.) 497 ! vstd = vstd + v * v * dt / (float(history_interval)*100.) 498 ! tstd = tstd + th * th * dt / (float(history_interval)*100.) 499 ! wstd = wstd + w * w * dt / (float(history_interval)*100.) 500 500 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 501 501 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! -
trunk/MESOSCALE/LMD_MM_MARS/makemeso
r335 r341 264 264 cp ../../SRC/WRFV2/Registry/Registry.EM Registry/ 265 265 fi 266 if [[ "${config}" == "les" ]] 267 then 268 cp ../../SRC/LES/WRFV2/Registry/Registry.EM Registry/ 269 fi 266 270 ### ensure the model will be recompiled from scratch with correct registry 267 271 cd Registry ; mv Registry Registry.bak ; Registry.bash ; cd .. … … 452 456 sed s+"-L../mars_lmd/libo -llmd"+"-L../mars_lmd/libo -llmd $NETCDF/lib/libnetcdf.a"+g configure.wrf > yeah ; mv -f yeah configure.wrf 453 457 sed s+"mpif90"+"${WHERE_MPI}/mpif90"+g configure.wrf | sed s+"mpicc"+"${WHERE_MPI}/mpicc"+g | sed s+"mpicc"+"mpicc -DMPI2_SUPPORT"+g > yeah 458 #sed s+"mpif90"+"${WHERE_MPI}/mpif90"+g configure.wrf | sed s+"mpicc"+"${WHERE_MPI}/mpicc"+g | sed s+"mpicc"+"mpicc -DMPI2_SUPPORT -w -O3"+g > yeah 454 459 #| sed s+"fastsse"+"fast"+g > yeah 455 #mv -f yeah configure.wrf #;;460 mv -f yeah configure.wrf ;; 456 461 # #sed s+"-fastsse"+"-O2 -Munroll -Mcache_align"+g configure.wrf > yeah 457 462 # ##sed s+"-fastsse"+"-O2 -fpic"+g configure.wrf > yeah ## marche pas, ILM problem … … 459 464 #sed s+"-fastsse"+"-O3 -tp=nehalem-64 -fastsse -Mscalarsse -Mvect=sse -Mipa=fast"+g configure.wrf > yeah 460 465 #sed s+"-fastsse"+"-O2 -Munroll -Mcache_align"+g configure.wrf > yeah 461 mv -f yeah configure.wrf ;; 466 #sed s+"-fastsse"+""+g configure.wrf > yeah 467 #sed s+"-fastsse"+"-fast"+g configure.wrf > yeah 468 #sed s+"-fastsse"+"-Mvect=noaltcode -Msmartalloc -Mprefetch=distance:8 -Mfprelaxed"+g configure.wrf > yeah 469 #mv -f yeah configure.wrf ;; 462 470 # MPICH, 64 bits, OK with periodic BC but no positive definite 463 471 mpi_64_test) cd arch ; ln -sf ../configure.defaults.fix64nest configure.defaults ; cd .. … … 523 531 then 524 532 echo 'DEBUG DEBUG DEBUG DEBUG' 525 sed s+"#-g"+"-g"+g configure.wrf > yeah 533 sed s+"#-g"+"-g -Ktrap=fp -Mbounds"+g configure.wrf > yeah 534 #sed s+"#-g"+"-g"+g configure.wrf > yeah 526 535 mv -f yeah configure.wrf 527 536 fi -
trunk/MESOSCALE_DEV/NOTES.txt
r335 r341 10 10 .... avec fastsse partout (mm MPI), crash ds la sortie 2 11 11 .... option de base SVN avec mpich sans fast, marche pas non plus crashe entre 9 et 10 12 .... iradia=1 ne change rien, toujours le probleme 13 .... la simulation tourne lorsqu'on fait callrad = F 14 .... test avec Mbounds ne renvoie aucuune erreur 15 .... test avec juste -Ktrap=fp et rien d'autre renvoie aucune erreur 16 .... semble OK avec iaervar=1, probleme dans meso_dustopacity ??? ou avec la lecture NETCDF ??? 17 .... crashe direct avec WHERE_MPI=/donnees/aslmd/MPI/mpich2-1.3.1_normal/bin 18 .... alors que ca passe nickel avec iaervar=1 19 .... crashe direct avec mes librairies netcdf3.6.1 compilees avec le dernier pgf2011 20 .... ne crashe pas direct avec mes librairies netcdf4.0.1 compilees avec le 21 dernier pgf2011 et les options utilisees dans le modele ajoutees aux options 22 dans la librairie compilee dans /distrib/local 23 mais crashe entre la sortie 9 et 10 comme les autres !!!! 24 .... marche bien aussi avec iaervar=3... probleme de lecture netcdf ??? 25 .... experience: on a iaervar=4 mais juste apres readtesassim 26 on regle tauref a 0.3 et ca passe. donc ce n est pas un bug structurel. 27 les valeurs lues ne sont probablement les bonnes soit par probleme 28 dans le fichier soit par probleme structurel dans readtesassim 29 .... pourtant en affichant les valeurs on ne voit pas de pb particulier ! 30 .... en changeant le nom hfx en sensheat et en enlevant z0 qui pose un pb 31 avec l'ancienne physique, ca ne marche toujours pas. 32 .... crash: semble stocker les variables qui sortent de la physique OK 33 mais le reste, par exemple tsurf, est NaN c'est bizarre 34 .... avec ndomainsz=ngridmx le modele va plus loin et crashe le second jour 35 a la sortie 2 36 .... mm comportement sur ulrich que sur penn 37 .... avec mcmodel crashe tout de suite 38 .... idem en invalidant les options d optimisation dans WRF et LMD phys [non en fait il faut enlever dans MPI] 39 .... test avec netcdf3: marche pas. mais ne faut-il pas enlever toutes les options? 40 .... avec aucune option d'optimisation et netcdf3: Nan avant la sortie 2 41 .... avec aucune option d'optimisation et netcdf4: va plus loin mais NaN entre 9 et 10 42 .... options d'optimisation en -O3 toujours le mm probleme 43 .... toujours le mm probleme mm avec ulimit -s unlimited 44 45 .... test qui donne des sortie 2 des NaN en recompilant avec -fast partout 46 avec mpirun -np 1, aucun souci tout va bien 47 avec mpirun -np 8, souci egalement des la sortie 2 48 ... visiblement un souci avec readtesassim ??? 49 .... MAIS NON CAR SOUCI AUSSI AVEC iaervar=1 avec 8 procs 50 .... ALORS QUE PAS DE SOUCI AVEC iaervar=1 avec 4 procs 51 export NETCDF=/donnees/aslmd/MODELES/MESOSCALE_DEV/NETCDF/pgfortran_64_netcdf-4.0.1_fast 52 export WHERE_MPI=/donnees/aslmd/MPI/pgfortran_mpich2-1.3.1/bin 53 .... corrections readtesassim ne semblent rien changer... 54 .... sorties frequentes permettent de voir que le probleme est localisee 55 mais rempli tres vite le domaine 56 avec dt=40s probleme apparait au bout de 700s 57 avec dt=10s probleme apparait au bout de 300s 58 avec dt=100s problemen apparait au bout de 1200s 59 ... visiblement le probleme apparait aux jointures des domaines ? 60 ... probleme sur le vitesse verticale calculee ??? 61 ... visiblement non puisque mm comportement avec non_hydrostatic ou W est normal 62 ... apparemment il s'agit vraiment d'une instabilite numerique 63 ... mettre les tendances R..BLEN a 0 ne change rien... 64 ... changer dimradmars n'arrange pas en fait lorsquon met des sorties frequentes 65 ... bizarre un des 4 processes wrf.exe passe en D quelques fois ???? 66 ... ne marche pas avec les options de compilation de WRF3 !!! 67 (mais domain met moins de temps a compiler) 68 ... toujours le mm probleme en acces local disque 69 70 TEST AVEC DEBUG 71 .... s'arrete au beau milieu d integrations sans sortir de message d'erreur 72 TEST AVEC LES POUR VOIR SI PB CORRIGE AVEC WRF3 73 .... rsl_lite ne compile pas avec l'option -fast 74 .... OK avec nofast_nemesis version compilee de mpich2 75 TEST avec le vieux mpich2... CRASH aussi entre la sortie 9 et 10 76 77 memes erreurs avec RSL_LITE de WRF3 78 alors qu il compile sans souci chez LES 79 .... un probleme d'options de compilations ???? 80 .... pendre direct la librairie compilee chez WRF3 ??? 81 LES: run OK 82 juste des NaN a la toute fin... 83 84 peut etre faut il regler dans WRFV2 les warnings relies a la compilation de rsl_lite 85 ------- il y a probablement des choses a corriger 86 ------- coupler avec gcc [-CC=gcc] comme dans LES ???? 87 .... mais lorsqu on utilise le vieux mpi compile avec pgf7 pas de warnings ! 88 89 90 ...le debugger voir une floating exception sur lwtt dans la boucle avec kdlon 91 ...avec les options debug le modele semble aller loin OK --> a verifier?? 92 ...les warnins a la compilation ont ils de l importance ? 93 ...le fait que netcdf4 ne soit pas supporte ??? 94 ...longue compil sur module_domain.... 95 96 ...des pistes ici 97 http://www.mmm.ucar.edu/wrf/users/wrfv3/known-prob-3.0.1.html 98 99 fonctionne avec le vieux mpi dans pgf2011 [et netcdf4] 100 mais les jobs ne sont pas a 100% sur les procs 101 probleme donc... c est tres lent 102 103 test basique avec WRFV2.2.1 et le cas em_quarter_ss et mpipgf 104 memes resultats avec un proc ou plusieurs 105 pas de crash 12 106 13 107
Note: See TracChangeset
for help on using the changeset viewer.