Changeset 2869

Timestamp:

Jan 16, 2023, 4:47:09 PM (2 years ago)

Author:

jleconte

Message:

Use new tendencies with symbol PLA for planet

File:

• trunk/WRF.COMMON/INTERFACES_V4/module_lmd_driver.F (modified) (8 diffs)

trunk/WRF.COMMON/INTERFACES_V4/module_lmd_driver.F

@@ -30 +30 @@
         RADT, &
         TSK,PSFC, &
-        RTHBLTEN,RUBLTEN,RVBLTEN, &
+        RTHPLATEN,RUPLATEN,RVPLATEN, &
         num_3d_s,SCALAR, &
         num_3d_m,moist, &
@@ -117 +117 @@
      dz8w,p8w,p,exner,t,t8w,rho,u,v,z,th,p_hyd,p_hyd_w
 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(OUT ) :: &
-     RTHBLTEN,RUBLTEN,RVBLTEN, &
+     RTHPLATEN,RUPLATEN,RVPLATEN, &
      HR_SW,HR_LW,HR_DYN,DDT,DT_RAD,DT_VDF,RDUST,VMR_ICE,RICE,&
      CLOUDFRAC,RH,DQICE,DQVAP,DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF
@@ -442 +442 @@
 IF (TRACER_MODE .EQ. 0) THEN
     q_prof(:,1)=0.95
+ELSE IF (TRACER_MODE .GE. 1) THEN
+    ! to be clean we should have an automatized process that makes sure that moist is sent to igcm_h2o_vap and etc.
+    q_prof(:,1) = SCALAR(i,kps:kpe,j,P_QH2O) / (1.d0 + SCALAR(i,kps:kpe,j,P_QH2O)) !! P_xxx is the index for variable xxx.
+    q_prof(:,2) = SCALAR(i,kps:kpe,j,P_QH2O_ICE) / (1.d0 + SCALAR(i,kps:kpe,j,P_QH2O))
+    ! conversion from mass mixing ratio in WRF to specific concentration in Physiq
 ELSE IF (TRACER_MODE .GE. 42) THEN
     ! to be clean we should have an automatized process that makes sure that moist is sent to igcm_h2o_vap and etc.
@@ -684 +689 @@
 !! DEDUCE TENDENCIES FOR WRF !!
 !!***************************!!
-RTHBLTEN(ims:ime,kms:kme,jms:jme)=0.
-RUBLTEN(ims:ime,kms:kme,jms:jme)=0.
-RVBLTEN(ims:ime,kms:kme,jms:jme)=0.
+RTHPLATEN(ims:ime,kms:kme,jms:jme)=0.
+RUPLATEN(ims:ime,kms:kme,jms:jme)=0.
+RVPLATEN(ims:ime,kms:kme,jms:jme)=0.
 PSFC(ims:ime,jms:jme)=p8w(ims:ime,kms,jms:jme) ! was done in surface driver in regular WRF 
 !------------------------------------------------------------------!
@@ -702 +707 @@
 
     ! zonal wind
-  RUBLTEN(i,kps:kpe,j) = zdufi_omp(subs,kps:kpe)
+  RUPLATEN(i,kps:kpe,j) = zdufi_omp(subs,kps:kpe)
     ! meridional wind
-  RVBLTEN(i,kps:kpe,j) = zdvfi_omp(subs,kps:kpe) 
+  RVPLATEN(i,kps:kpe,j) = zdvfi_omp(subs,kps:kpe) 
     ! potential temperature
     ! (dT = dtheta * exner for isobaric coordinates or if pressure variations are negligible)
-  RTHBLTEN(i,kps:kpe,j) = zdtfi_omp(subs,kps:kpe) / exner(i,kps:kpe,j)
+  RTHPLATEN(i,kps:kpe,j) = zdtfi_omp(subs,kps:kpe) / exner(i,kps:kpe,j)
     ! update surface pressure (cf CO2 cycle in physics)
     ! here dt is needed
@@ -722 +727 @@
       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(1)
+      scalar(i,kps:kpe,j,P_QH2O)=scalar(i,kps:kpe,j,P_QH2O) &
+          +zdqfi_omp(subs,kps:kpe,1)*dt * (1.d0+scalar(i,kps:kpe,j,P_QH2O))
+      scalar(i,kps:kpe,j,P_QH2O_ICE)=scalar(i,kps:kpe,j,P_QH2O_ICE) &
+           +zdqfi_omp(subs,kps:kpe,2)*dt * (1.d0+scalar(i,kps:kpe,j,P_QH2O))
+       ! if you want to use this mode, RTHPLATEN should be corrected as below. 
+       ! we keep it like that for the moment for testing.
     CASE(42)
       moist(i,kps:kpe,j,P_QV)=moist(i,kps:kpe,j,P_QV) &
@@ -727 +739 @@
       scalar(i,kps:kpe,j,P_QH2O_ICE)=scalar(i,kps:kpe,j,P_QH2O_ICE) &
            +zdqfi_omp(subs,kps:kpe,2)*dt * (1.d0+moist(i,kps:kpe,j,P_QV))
-       ! if you want to use this mode, RTHBLTEN should be corrected as below. 
+       ! if you want to use this mode, RTHPLATEN should be corrected as below. 
        ! we keep it like that for the moment for testing.
     CASE(43)
@@ -734 +746 @@
       scalar(i,kps:kpe,j,P_QH2O_ICE)=scalar(i,kps:kpe,j,P_QH2O_ICE) &
            +zdqfi_omp(subs,kps:kpe,2)*dt * (1.d0+moist(i,kps:kpe,j,P_QV))
-      RTHBLTEN(i,kps:kpe,j) = RTHBLTEN(i,kps:kpe,j) &
+      RTHPLATEN(i,kps:kpe,j) = RTHPLATEN(i,kps:kpe,j) &
         * (1.d0+moist(i,kps:kpe,j,P_QV))/(1.d0+rvovrd*moist(i,kps:kpe,j,P_QV))
         ! correct dT/dt assuming a constant molar heat capacity.