Changeset 2208 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F.new

Timestamp:

Jan 7, 2020, 10:32:52 PM (5 years ago)

Author:

aslmd

Message:

MESOSCALE: adapted what was proposed in r2123. created a mode mars=42 for water vapour used in dynamics (what was proposed in r2123). moved back the default behaviour for Mars that was removed in r2123.

File:

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

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

@@ -495 +495 @@
 IF (MARS_MODE .EQ. 0) THEN
     q_prof(:,1)=0.95
-ELSE
+ELSE IF (MARS_MODE .EQ. 42) THEN
     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  
+ELSE
+    q_prof(:,1:nq) = SCALAR(i,kps:kpe,j,2:nq+1)  !! the names were set above !! one dummy tracer in WRF
 ENDIF
 
@@ -512 +514 @@
   ENDIF
 ENDIF
-
 
 IF (firstcall .EQV. .true.) THEN
@@ -805 +806 @@
       SCALAR(i,kps:kpe,j,2) = SCALAR(i,kps:kpe,j,2)*exp(-dt/tau_decay)
       SCALAR(i,1,j,2) = SCALAR(i,1,j,2) + 1. !! this tracer is emitted in the surface layer
-    CASE DEFAULT
+    CASE(42)
       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
+    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
   END SELECT
-    
+
 ENDDO
 ENDDO