Changeset 2123 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

Timestamp:

Apr 15, 2019, 4:11:11 PM (7 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.

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new (modified) (4 diffs)

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