Changeset 34 for trunk/mesoscale/LMD_LES_MARS
- Timestamp:
- Jan 26, 2011, 3:37:26 AM (14 years ago)
- Location:
- trunk/mesoscale/LMD_LES_MARS
- Files:
-
- 1 added
- 2 deleted
- 7 edited
-
LMD_LES_MARS_install (modified) (2 diffs)
-
modif_mars/DEF/B_START8_REAL.tar.gz_is_a_good_quick_try (added)
-
modif_mars/Registry.EM (modified) (11 diffs)
-
modif_mars/home (deleted)
-
modif_mars/makeles (modified) (2 diffs)
-
modif_mars/makeles_ciclad (modified) (1 diff)
-
modif_mars/module_big_step_utilities_em.F (modified) (5 diffs)
-
modif_mars/module_first_rk_step_part1.F (modified) (4 diffs)
-
modif_mars/module_lmd_driver.F (deleted)
-
modif_mars/module_model_constants.F (modified) (8 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/mesoscale/LMD_LES_MARS/LMD_LES_MARS_install
r17 r34 22 22 #cd LMD_LES_MARS 23 23 #ln -sf ../LMD_LES_MARS_patches modif_mars 24 25 cd modif_mars 26 ln -sf ../../LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F 27 cd .. 24 28 25 29 cd dyn_em … … 56 60 57 61 62 ### create a fake directory 63 mkdir -p mars_lmd/libf/phymars 64 -
trunk/mesoscale/LMD_LES_MARS/modif_mars/Registry.EM
r17 r34 71 71 ################################################################################ 72 72 73 73 74 ####MARS MARS MARS MARS MARS 74 75 ####MARS MARS MARS MARS MARS … … 115 116 ####MARS MARS MARS MARS MARS 116 117 ####MARS MARS MARS MARS MARS 118 119 120 117 121 118 122 … … 185 189 endif 186 190 187 188 ###MARS MARS MARS MARS MARS 191 ####MARS MARS MARS MARS MARS 189 192 ####MARS MARS MARS MARS MARS 190 193 #### … … 197 200 state real slpx_gc ij dyn_em 1 - i1 "SLPX" "slope x direction" "rad" 198 201 state real slpy_gc ij dyn_em 1 - i1 "SLPY" "slope y direction" "rad" 199 state real hv_gc igj dyn_em 1 Z i1 "HV" "water vapor" "kg/kg" 202 state real hv_gc igj dyn_em 1 Z i1 "HV" "water vapor" "kg/kg" 200 203 state real hi_gc igj dyn_em 1 Z i1 "HI" "water ice" "kg/kg" 201 204 state real tsoil_gc igj dyn_em 1 Z i1 "TSOIL" "soil temperatures" "K" 205 state real isoil_gc igj dyn_em 1 Z i1 "ISOIL" "soil thermal inertia" "tiu" 206 state real dsoil_gc igj dyn_em 1 Z i1 "DSOIL" "soil depths" "m" 202 207 #### 203 208 #### variables dans wrf … … 210 215 state real mars_cice ij misc 1 - i012rd "MARS_CICE" "co2 ice" "kg" 211 216 state real mars_tsoil ilj misc 1 - i012rd "MARS_TSOIL" "soil temperatures" "K" 217 state real mars_isoil ilj misc 1 - i012rd "MARS_ISOIL" "soil thermal inertia" "tiu" 218 state real mars_dsoil ilj misc 1 - i012rd "MARS_DSOIL" "soil depths" "m" 212 219 state real mars_gw ilj misc 1 - i012rd "MARS_GW" "gwparam" " " 213 220 #### … … 218 225 219 226 220 221 227 #----------------------------------------------------------------------------------------------------------------------------------------------------------------- 222 228 # 223 229 # Variables for Eulerian mass coordinate dynamics 224 230 # 225 226 231 227 232 # Velocities … … 260 265 state real ph ikjb dyn_em 2 Z \ 261 266 irusdf=(bdy_interp:dt) "ph" "perturbation geopotential" "m2 s-2" 262 state real phb ikj dyn_em 1 Z irdus "phb" "base-state geopotential" "m2 s-2"267 state real phb ikj dyn_em 1 Z irdus "phb" "base-state geopotential" "m2 s-2" 263 268 state real phb_fine ikj dyn_em 1 Z - "phb_fine" "for nesting, temp holding interpolated coarse grid phb" "m2 s-2" 264 269 state real ph0 ikj dyn_em 1 Z r "ph0" "initial geopotential" "m2 s-2" … … 282 287 i1 real t_2save ikj dyn_em 1 - 283 288 state real t_save ikj dyn_em 1 - "t_save" 284 # 289 # 285 290 286 291 … … 288 293 state real mu ijb dyn_em 2 - \ 289 294 irusdf=(bdy_interp:dt) "mu" "perturbation dry air mass in column" "Pa" 290 state real mub ij dyn_em 1 - irdus "mub" "base state dry air mass in column" "Pa"295 state real mub ij dyn_em 1 - irdus "mub" "base state dry air mass in column" "Pa" 291 296 state real mub_fine ij dyn_em 1 - - "mub_fine" "nest temp, holds interpolated coarse grid mub" "Pa" 292 297 state real mub_save ij dyn_em 1 - - "mub_save" "nest temp, holds orig fine grid mub" "Pa" … … 315 320 316 321 # Pressure and Density 317 state real p ikj dyn_em 1 - r "p" "perturbation pressure" "Pa"322 state real p ikj dyn_em 1 - r "p" "perturbation pressure" "Pa" 318 323 state real al ikj dyn_em 1 - r "al" "inverse perturbation density" "m3 kg-1" 319 324 state real alt ikj dyn_em 1 - r "alt" "inverse density" "m3 kg-1" … … 323 328 state real rdz ikj dyn_em 1 Z - " " " " " " 324 329 state real rdzw ikj dyn_em 1 Z - " " " " " " 325 state real pb ikj dyn_em 1 - rdus "pb" "BASE STATE PRESSURE " "Pa"330 state real pb ikj dyn_em 1 - rdus "pb" "BASE STATE PRESSURE " "Pa" 326 331 ## 327 332 ## MARS: TO GET LIGHTER HISTORY FILES 328 333 ## 329 334 state real ptot ikj dyn_em 1 - irh "PTOT" "total pressure" "Pa" 335 ####### if you choose h here, you will probably want not to output u,v,w,t,etc... 336 state real uave ikj dyn_em 1 - ir "UAVE" "average zonal wind" "m/s" 337 state real ustd ikj dyn_em 1 - ir "USTD" "std dev zonal wind" "m/s" 338 state real vave ikj dyn_em 1 - ir "VAVE" "average meridional wind" "m/s" 339 state real vstd ikj dyn_em 1 - ir "VSTD" "std dev meridional wind" "m/s" 340 state real wave ikj dyn_em 1 - ir "WAVE" "average vertical wind" "m/s" 341 state real wstd ikj dyn_em 1 - ir "WSTD" "std dev vertical wind" "m/s" 342 state real tave ikj dyn_em 1 - ir "TAVE" "average pert. potential temperature" "K" 343 state real tstd ikj dyn_em 1 - ir "TSTD" "std dev pert. potential temperature" "K" 344 ####### 345 ## ajouter tk, ptot, phtot ??? les deux derniers sont assez smooth... 330 346 331 347 # -
trunk/mesoscale/LMD_LES_MARS/modif_mars/makeles
r24 r34 8 8 # A. Spiga 01/2009 9 9 #------------------------------ 10 11 12 echo "Use -n if you don't want to recompile physics" 10 13 11 14 … … 50 53 echo $libphys_folder 51 54 cp -rf ${libphys_folder}/liblmd.a ${whereliblmd}/ 52 cp -rf ${phymars_folder}/dimphys.h ${here}/modif_mars/ 53 cp -rf ${folder}/call_meso_inifis.inc ${here}/modif_mars/ 54 cp -rf ${folder}/call_meso_physiq.inc ${here}/modif_mars/ 55 cp -rf ${folder}/inc/module_lmd_driver_output*.inc ${here}/modif_mars/ 56 cp -rf ${phymars_folder}/wrf_output_2d.h ${here}/modif_mars/ 57 cp -rf ${phymars_folder}/wrf_output_3d.h ${here}/modif_mars/ 58 55 #cp -rf ${phymars_folder}/dimphys.h ${here}/modif_mars/ 56 cp -rf ${phymars_folder}/dimphys.h ${here}/mars_lmd/libf/phymars/ 57 #cp -rf ${folder}/call_meso_inifis.inc ${here}/modif_mars/ 58 #cp -rf ${folder}/call_meso_physiq.inc ${here}/modif_mars/ 59 cp -rf ${folder}/call_meso_inifis.inc ${here}/ 60 cp -rf ${folder}/call_meso_physiq.inc ${here}/ 61 #cp -rf ${folder}/inc/module_lmd_driver_output*.inc ${here}/modif_mars/ 62 cp -rf ${folder}/inc/module_lmd_driver_output*.inc ${here}/inc/ 63 #cp -rf ${phymars_folder}/wrf_output_2d.h ${here}/modif_mars/ 64 #cp -rf ${phymars_folder}/wrf_output_3d.h ${here}/modif_mars/ 65 cp -rf ${phymars_folder}/wrf_output_2d.h ${here}/mars_lmd/libf/phymars/ 66 cp -rf ${phymars_folder}/wrf_output_3d.h ${here}/mars_lmd/libf/phymars/ 59 67 fi 60 68 ############################################################## -
trunk/mesoscale/LMD_LES_MARS/modif_mars/makeles_ciclad
r17 r34 50 50 echo $libphys_folder 51 51 cp -rf ${libphys_folder}/liblmd.a ${whereliblmd}/ 52 cp -rf ${phymars_folder}/dimphys.h ${here}/modif_mars/ 53 cp -rf ${folder}/call_meso_inifis.inc ${here}/modif_mars/ 54 cp -rf ${folder}/call_meso_physiq.inc ${here}/modif_mars/ 55 cp -rf ${folder}/inc/module_lmd_driver_output*.inc ${here}/modif_mars/ 56 cp -rf ${phymars_folder}/wrf_output_2d.h ${here}/modif_mars/ 57 cp -rf ${phymars_folder}/wrf_output_3d.h ${here}/modif_mars/ 58 52 #cp -rf ${phymars_folder}/dimphys.h ${here}/modif_mars/ 53 cp -rf ${phymars_folder}/dimphys.h ${here}/mars_lmd/libf/phymars/ 54 #cp -rf ${folder}/call_meso_inifis.inc ${here}/modif_mars/ 55 #cp -rf ${folder}/call_meso_physiq.inc ${here}/modif_mars/ 56 cp -rf ${folder}/call_meso_inifis.inc ${here}/ 57 cp -rf ${folder}/call_meso_physiq.inc ${here}/ 58 #cp -rf ${folder}/inc/module_lmd_driver_output*.inc ${here}/modif_mars/ 59 cp -rf ${folder}/inc/module_lmd_driver_output*.inc ${here}/inc/ 60 #cp -rf ${phymars_folder}/wrf_output_2d.h ${here}/modif_mars/ 61 #cp -rf ${phymars_folder}/wrf_output_3d.h ${here}/modif_mars/ 62 cp -rf ${phymars_folder}/wrf_output_2d.h ${here}/mars_lmd/libf/phymars/ 63 cp -rf ${phymars_folder}/wrf_output_3d.h ${here}/mars_lmd/libf/phymars/ 59 64 fi 60 65 ############################################################## -
trunk/mesoscale/LMD_LES_MARS/modif_mars/module_big_step_utilities_em.F
r17 r34 4347 4347 4348 4348 SUBROUTINE phy_prep ( config_flags, & ! input 4349 mu, muu, muv, u, v, p, pb, alt, ph,& ! input4349 mu, muu, muv, u, v, w, p, pb, alt, ph, & ! input 4350 4350 phb, t, tsk, moist, n_moist, & ! input 4351 4351 mu_3d, rho, th_phy, p_phy , pi_phy , & ! output 4352 u_phy, v_phy, p8w, t_phy, t8w,& ! output4352 u_phy, v_phy, w_phy, p8w, t_phy, t8w, & ! output 4353 4353 z, z_at_w, dz8w, & ! output 4354 4354 fzm, fzp, & ! params … … 4383 4383 INTENT( OUT) :: u_phy, & 4384 4384 v_phy, & 4385 w_phy, & 4385 4386 pi_phy, & 4386 4387 p_phy, & … … 4400 4401 u, & 4401 4402 v, & 4403 w, & 4402 4404 alt, & 4403 4405 ph, & … … 4512 4514 do k = k_start, k_end 4513 4515 do i = i_start, i_end 4514 z(i,k,j) = 0.5*(z_at_w(i,k,j) +z_at_w(i,k+1,j) ) 4516 z(i,k,j) = 0.5*(z_at_w(i,k,j) + z_at_w(i,k+1,j) ) 4517 !!!! MARS MARS ajout aymeric (ainsi que les arguments de cette routine) 4518 w_phy(i,k,j) = 0.5*(w(i,k,j) + w(i,k+1,j) ) 4515 4519 enddo 4516 4520 enddo … … 4640 4644 IF ( (config_flags%bl_pbl_physics .gt. 0) & 4641 4645 .OR. (config_flags%modif_wrf) ) THEN 4642 4646 !****MARS 4643 4647 DO J=j_start,j_end 4644 4648 DO K=k_start,k_end -
trunk/mesoscale/LMD_LES_MARS/modif_mars/module_first_rk_step_part1.F
r17 r34 78 78 REAL ,DIMENSION(ims:ime,kms:kme,jms:jme), INTENT(INOUT) :: u_phy 79 79 REAL ,DIMENSION(ims:ime,kms:kme,jms:jme), INTENT(INOUT) :: v_phy 80 REAL ,DIMENSION(ims:ime,kms:kme,jms:jme) :: w_phy !!MARS aymeric var loc 80 81 REAL ,DIMENSION(ims:ime,kms:kme,jms:jme), INTENT(INOUT) :: dz8w 81 82 REAL ,DIMENSION(ims:ime,kms:kme,jms:jme), INTENT(INOUT) :: p8w … … 147 148 CALL phy_prep ( config_flags, & 148 149 grid%mut, grid%muu, grid%muv, grid%u_2, & 149 grid%v_2, grid% p, grid%pb, grid%alt, &150 grid%v_2, grid%w_2, grid%p, grid%pb, grid%alt, & 150 151 grid%ph_2, grid%phb, grid%t_2, grid%tsk, moist, num_moist, & 151 152 mu_3d, rho, & 152 th_phy, p_phy, pi_phy, u_phy, v_phy, &153 th_phy, p_phy, pi_phy, u_phy, v_phy, w_phy, & 153 154 p8w, t_phy, t8w, grid%z, z_at_w, & 154 155 dz8w, grid%fnm, grid%fnp, & … … 594 595 ! ... (cf phy_prep in module_big_step_utilities_em) 595 596 & ,HT=grid%ht & 596 & ,U=u_phy ,V=v_phy 597 & ,U=u_phy ,V=v_phy, W=w_phy & 597 598 & ,TH=th_phy ,T=t_phy & 598 599 & ,P=p_phy ,EXNER=pi_phy ,RHO=rho & … … 621 622 & ,CST_TI=config_flags%init_TI & 622 623 & ,ISFFLX=config_flags%isfflx & 624 & ,DIFF_OPT=config_flags%diff_opt & 625 & ,KM_OPT=config_flags%km_opt & 626 ! additional diagnostics 627 ! & ,UAVE=grid%em_uave, USTD=grid%em_ustd & 628 ! & ,VAVE=grid%em_vave, VSTD=grid%em_vstd & 629 ! & ,WAVE=grid%em_wave, WSTD=grid%em_wstd & 630 ! & ,TAVE=grid%em_tave, TSTD=grid%em_tstd & 631 & ,UAVE=grid%uave, USTD=grid%ustd & 632 & ,VAVE=grid%vave, VSTD=grid%vstd & 633 & ,WAVE=grid%wave, WSTD=grid%wstd & 634 & ,TAVE=grid%tave, TSTD=grid%tstd & 635 & ,HISTORY_INTERVAL=model_config_rec%history_interval(1) & 623 636 !------------------! 624 637 ! OUTPUT VARIABLES ! 625 638 !------------------! 626 #include " ../modif_mars/module_lmd_driver_output4.inc"639 #include "module_lmd_driver_output4.inc" 627 640 & ,SLPX=grid%slpx,SLPY=grid%slpy) 628 641 ENDIF -
trunk/mesoscale/LMD_LES_MARS/modif_mars/module_model_constants.F
r17 r34 1 1 !WRF:MODEL_LAYER:CONSTANTS 2 2 ! 3 !PLANET MARS 3 4 4 5 MODULE module_model_constants … … 18 19 19 20 #if ( NMM_CORE == 1 ) 20 REAL , PARAMETER :: r_d = 192. 21 REAL , PARAMETER :: cp = 844.6 21 REAL , PARAMETER :: r_d = 192. ! gas constant m2 s-2 K-1 22 REAL , PARAMETER :: cp = 844.6 ! r= 8.314511E+0 *1000.E+0/mugaz 22 23 #else 23 24 REAL , PARAMETER :: r_d = 192. … … 25 26 #endif 26 27 27 REAL , PARAMETER :: r_v = 461.6 28 REAL , PARAMETER :: r_v = 461.6 ! gas constant for water vapor 28 29 REAL , PARAMETER :: cv = cp-r_d 29 REAL , PARAMETER :: cpv = 4.4*r_v 30 REAL , PARAMETER :: cpv = 4.4*r_v ! cp / R = 4.4 ?? ! rcp = r/cp ~0.256793 30 31 REAL , PARAMETER :: cvv = cpv-r_v 31 32 REAL , PARAMETER :: cvpm = -cv/cp … … 40 41 41 42 REAL , PARAMETER :: p1000mb = 610. 42 REAL , PARAMETER :: t0 = 220. 43 REAL , PARAMETER :: t0 = 220. ! earth : 300 ... remplacer par 220 ? 43 44 REAL , PARAMETER :: p0 = 610. 44 45 REAL , PARAMETER :: cpovcv = cp/(cp-r_d) … … 46 47 REAL , PARAMETER :: rvovrd = r_v/r_d 47 48 48 REAL , PARAMETER :: reradius = 1./3397200. 49 REAL , PARAMETER :: reradius = 1./3397200. ! Rayon de mars (m) ~3397200 m 49 50 50 51 !------------------------------- … … 75 76 !------------------------------- 76 77 77 78 79 78 REAL , PARAMETER :: asselin = .025 80 79 ! REAL , PARAMETER :: asselin = .0 … … 92 91 REAL , PARAMETER :: rhowater = 1000. 93 92 REAL , PARAMETER :: rhosnow = 100. 94 REAL , PARAMETER :: rhoair0 = 0.02 95 93 REAL , PARAMETER :: rhoair0 = 0.02 !Earth Surface density: 1.217 kg/m3 94 !Mars Surface density: 0.020 kg/m3 96 95 REAL , PARAMETER :: DEGRAD = 3.1415926/180. 97 REAL , PARAMETER :: DPD = 360./669. 96 REAL , PARAMETER :: DPD = 360./669. ! longitude solaire d'un jour? 98 97 99 98 REAL , PARAMETER :: SVP1=0.6112 … … 103 102 REAL , PARAMETER :: EP_1=R_v/R_d-1. 104 103 REAL , PARAMETER :: EP_2=R_d/R_v 105 REAL , PARAMETER :: KARMAN=0.4 106 REAL , PARAMETER :: EOMEG=7.0721E-5 107 REAL , PARAMETER :: STBOLT=5.67051E-8 104 REAL , PARAMETER :: KARMAN=0.4 ! von karman constant 105 REAL , PARAMETER :: EOMEG=7.0721E-5 ! angular rotation rate (rad.s-1) 106 REAL , PARAMETER :: STBOLT=5.67051E-8 ! stefan-boltzmann constant 107 108 109 110 108 111 109 112 REAL , PARAMETER :: prandtl = 1./3.0
Note: See TracChangeset
for help on using the changeset viewer.
