Changeset 2123

Timestamp:

Apr 15, 2019, 4:11:11 PM (6 years ago)

Author:

jleconte

Message:

Changed the routing of the tracers so that the mean molecular weight of water vapor be taken into account in the dynamics.

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC

Files:

• LES/modif_mars/Registry.EM.newphys (modified) (4 diffs)
• LES/modif_mars/module_first_rk_step_part1.F (modified) (1 diff)
• WRFV2/phys/module_lmd_driver.F.new (modified) (4 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys

@@ -113 +113 @@
 state  real  FLUXSURF_SW ij misc  1  -  - "FLUXSURF_SW" "SW Surface flux" "W.m-2"
 state  real  FLUXSURF_LW ij misc  1  -  -  "FLUXSURF_LW" "LW Surface flux" "W.m-2"
-state  real  FLXGRD ij misc  1  -  - "FLXGRD" "heat flux from ground"  "W.m-2"
+state  real  FLXGRD ij misc  1  -  h "FLXGRD" "heat flux from ground"  "W.m-2"
 state  real  LSCEZ ikj misc  1  -  - "LSCEZ" "heat from largescale" "W.m-2"
 state  real  H2OICE_REFF ikj misc  1  -  -  "H2OICE_REFF" "Radii of H2O ice" "m"
@@ -133 +133 @@
 ####
 state  real  -         ikjftb  scalar  1  -  -                         -
-state  real  QH2O      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt) "QH2O"       "Water vapor mixing ratio"   "kg kg-1"
-state  real  QH2O_ICE  ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt) "QH2O_ICE"   "Water ice mixing ratio"     "kg kg-1"
+state  real  QH2O      ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt) "QH2O"       "Water vapor mixing ratio"   "kg kg-1"
+state  real  QH2O_ICE  ikjftb  scalar  1  -  i01rusdf=(bdy_interp:dt) "QH2O_ICE"   "Water ice mixing ratio"     "kg kg-1"
 state  real  qdust     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt) "QDUST"      "Dust mixing ratio"          "kg kg-1"
 state  real  qdustn    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt) "QDUSTN"     "Dust number density"        " m-3 "
@@ -1047 +1047 @@
 state    real  CDA              ij      misc        1         -      r        "CDA"                   "DRAG COEFF AT LOWEST MODEL LVL"                      ""
 state    real  UST              ij      misc        1         -      r        "UST"                   "U* IN SIMILARITY THEORY"                      "m s-1"    
-state    real  USTM             ij      misc        1         -      rh       "USTM"                  "U* IN SIMILARITY THEORY WITHOUT VCONV"        "m s-1"    #SAVEMARS2 ustar
+state    real  USTM             ij      misc        1         -      r       "USTM"                  "U* IN SIMILARITY THEORY WITHOUT VCONV"        "m s-1"    #SAVEMARS2 ustar
 i1       real  HOL              ij      misc        1         -      -        "HOL"                   "PBL HEIGHT OVER MONIN-OBUKHOV LENGTH"         ""
 state    real  RMOL             ij      misc        1         -      r        "RMOL"                  "1./Monin Ob. Length"                      ""
@@ -1544 +1544 @@
 ##
 package   default      mars==0                      -              -
-package   water        mars==1                      -              scalar:qh2o,qh2o_ice
+package   water        mars==1                      -              scalar:qh2o,qh2o_ice;moist:qv,qc
 package   dust1        mars==2                      -              scalar:qdust
 package   dust2eq      mars==3                      -              scalar:qdust,qdustn

trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_first_rk_step_part1.F

@@ -613 +613 @@
            ! scalars
 !     &        ,NUM_3D_S=num_3d_s,SCALAR=grid%scalar                       &
-     &        ,NUM_3D_S=num_scalar,SCALAR=grid%scalar                      &
+     &        ,NUM_3D_S=num_scalar,SCALAR=grid%scalar                       &
+     &        ,NUM_3D_M=num_moist,MOIST=grid%moist                      &
            ! mode
      &        ,MARS_MODE=config_flags%mars   &

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new

@@ -30 +30 @@
         RTHBLTEN,RUBLTEN,RVBLTEN, &
         num_3d_s,SCALAR, &
+        num_3d_m,moist, &
         MARS_MODE, &
         planet_type, &
@@ -124 +125 @@
      M_GW
 ! 4D arrays
-INTEGER, INTENT(IN ) :: num_3d_s
-REAL, DIMENSION( ims:ime, kms:kme, jms:jme, 1:num_3d_s ), INTENT(INOUT ) :: &
-     scalar
+INTEGER, INTENT(IN ) :: num_3d_s,num_3d_m
+REAL, DIMENSION( ims:ime, kms:kme, jms:jme, 1:num_3d_s ), INTENT(INOUT ) :: scalar
+REAL, DIMENSION( ims:ime, kms:kme, jms:jme, 1:num_3d_m ), INTENT(INOUT ) ::  moist
 ! Logical
 LOGICAL, INTENT(IN ) :: restart
@@ -495 +496 @@
     q_prof(:,1)=0.95
 ELSE
-    q_prof(:,1:nq) = SCALAR(i,kps:kpe,j,2:nq+1)  !! the names were set above !! one dummy tracer in WRF
+    q_prof(:,1) = moist(i,kps:kpe,j,P_QV) !! the names were set above !! one dummy tracer in WRF
+    q_prof(:,2) = SCALAR(i,kps:kpe,j,3) !! the names were set above !! one dummy tracer in WRF  
 ENDIF
 
@@ -804 +806 @@
       SCALAR(i,1,j,2) = SCALAR(i,1,j,2) + 1. !! this tracer is emitted in the surface layer
     CASE DEFAULT
-      SCALAR(i,kps:kpe,j,2:nq+1)=SCALAR(i,kps:kpe,j,2:nq+1)+zdqfi_omp(subs,kps:kpe,1:nq)*dt  !!! here dt is needed
+      moist(i,kps:kpe,j,P_QV)=moist(i,kps:kpe,j,P_QV)+zdqfi_omp(subs,kps:kpe,1)*dt  !!! here dt is needed
+      scalar(i,kps:kpe,j,3)=scalar(i,kps:kpe,j,3)+zdqfi_omp(subs,kps:kpe,2)*dt  !!! here dt is needed
   END SELECT