Changeset 1585 for trunk/MESOSCALE/LMD_MM_MARS/SRC

Timestamp:

Aug 28, 2016, 2:24:01 PM (8 years ago)

Author:

aslmd

Message:

adapted module_model_constants (and module_diffusion_em) so that it is common between WRFV2 mesoscale and WRFV3 LES versions.

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

Files:

• dyn_em/module_diffusion_em.F (modified) (1 diff)
• share/module_model_constants.F (modified) (6 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_diffusion_em.F

@@ -9 +9 @@
     USE module_model_constants    
     USE module_wrf_error
+
+!!!******MARS MARS
+!!!******MARS MARS
+!!! mixing coefficients are defined here so that
+!!! module_model_constants can be common btw WRFV2&3 
+                                      ! proportionality constants for eddy viscosity coefficient calc
+   REAL    , PARAMETER ::  c_s = .25  ! turbulence parameterization constant, for smagorinsky
+!   REAL    , PARAMETER ::  c_s = .1  ! turbulence parameterization constant, for smagorinsky
+   REAL    , PARAMETER ::  c_k = .15  ! turbulence parameterization constant, for TKE
+!!!******MARS MARS
+!!!******MARS MARS
 
     CONTAINS

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/share/module_model_constants.F

@@ -28 +28 @@
    REAL    , PARAMETER :: r_d          = 191.
    REAL    , PARAMETER :: cp           = 744.5
+
 
    REAL    , PARAMETER :: r_v          = 461.6  ! gas constant for water vapor
@@ -42 +43 @@
    REAL    , PARAMETER :: rovg         = r_d/g
    REAL    , PARAMETER :: c2           = cp * rcv
+   !! added in WRFV3
+   real    , parameter :: mwdry        = 43.49 ! molecular weight of dry air (g/mole)
 
    REAL    , PARAMETER :: p1000mb      = 610.
@@ -82 +85 @@
 !-------------------------------
 
-!   REAL    , PARAMETER :: asselin      = .025
-   REAL    , PARAMETER :: asselin      = .0
+   REAL    , PARAMETER :: asselin      = .025
+!   REAL    , PARAMETER :: asselin      = .0
    REAL    , PARAMETER :: cb           = 25.
 
@@ -112 +115 @@
    REAL    , PARAMETER ::  STBOLT=5.67051E-8 ! stefan-boltzmann constant
 
-                                      ! proportionality constants for eddy viscosity coefficient calc
-   REAL    , PARAMETER ::  c_s = .25  ! turbulence parameterization constant, for smagorinsky
-!   REAL    , PARAMETER ::  c_s = .1  ! turbulence parameterization constant, for smagorinsky
-   REAL    , PARAMETER ::  c_k = .15  ! turbulence parameterization constant, for TKE
+
+
+
+   !! NB: c_s and c_k in WRFV2 moved to module_diffusion_em.F
+   !! NB: no need in WRFV3 because set in namelist
    REAL    , PARAMETER ::  prandtl = 1./3.0
                                          ! constants for w-damping option
@@ -123 +127 @@
    REAL    , PARAMETER ::  w_beta  = 0.8 ! activation cfl number
 
+!! in WRFV3
+!   REAL    , PARAMETER ::  w_alpha = 0.3 ! strength m/s/s
+!   REAL    , PARAMETER ::  w_beta  = 1.0 ! activation cfl number
+!!
        REAL , PARAMETER ::  pq0=379.90516
        REAL , PARAMETER ::  epsq2=0.2
-!try both ?
-!       REAL , PARAMETER ::  epsq2=0.02
+!      REAL , PARAMETER ::  epsq2=0.02
        REAL , PARAMETER ::  a2=17.2693882
        REAL , PARAMETER ::  a3=273.16
@@ -172 +179 @@
 !#endif
 
+!!!! needed in WRFV3
+   !  Earth
+
+   !  The value for P2SI *must* be set to 1.0 for Earth
+   !  Although, now we may not need this declaration here (see above)
+   !REAL    , PARAMETER :: P2SI         = 1.0
+
+   !  Orbital constants:
+
+   INTEGER , PARAMETER :: PLANET_YEAR = 365
+   REAL , PARAMETER :: OBLIQUITY = 23.5
+   REAL , PARAMETER :: ECCENTRICITY = 0.014
+   REAL , PARAMETER :: SEMIMAJORAXIS = 1.0 ! In AU
+   ! Don't know the following values, so we'll fake them for now
+   REAL , PARAMETER :: zero_date = 0.0   ! Time of perihelion passage
+   !  Fraction into the year (from perhelion) of the
+   !  occurrence of the Northern Spring Equinox
+   REAL , PARAMETER :: EQUINOX_FRACTION= 0.0
+!!!! needed in WRFV3
 
  CONTAINS