Changeset 3661 for trunk/WRF.COMMON

Timestamp:

Feb 27, 2025, 2:53:11 PM (4 months ago)

Author:

emoisan

Message:

Titan CRM:
Add Titan interface in INTERFACES_V4
Adapt module_model_constants.F to Titan
Add new tracer_mode for Titan (CH4 scalar)
Add new communication of variables between LMDZ.TITAN and WRF
Allow microphysics for Mesoscale in physiq_mod.F90
EMo

Location:

trunk/WRF.COMMON/INTERFACES_V4

Files:

• dynphy_wrf_titan_lmd (added)
• dynphy_wrf_titan_lmd/callphysiq_mod.F (added)
• dynphy_wrf_titan_lmd/iniphysiq_mod.F (added)
• dynphy_wrf_titan_lmd/update_inputs_physiq_mod.F (added)
• dynphy_wrf_titan_lmd/update_outputs_physiq_mod.F (added)
• dynphy_wrf_titan_lmd/variables_mod.F (added)
• module_lmd_driver.F (modified) (5 diffs)
• module_model_constants.F (modified) (3 diffs)

trunk/WRF.COMMON/INTERFACES_V4/module_lmd_driver.F

@@ -48 +48 @@
         HFMAX,ZMAX,&
         USTM,HFX,&
-        SLPX,SLPY,RESTART)
+        SLPX,SLPY,RESTART,&
+        DT_COND)
 ! NB: module_lmd_driver_output1.inc : output arguments generated from Registry
 
@@ -114 +115 @@
      RTHPLATEN,RUPLATEN,RVPLATEN, &
      HR_SW,HR_LW,HR_DYN,DT_RAD,&
-     CLOUDFRAC,RH,DQICE,DQVAP,DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF
+     CLOUDFRAC,RH,DQICE,DQVAP,DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF,&
+     DT_COND
 REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(INOUT ) :: &
      P_Q2
@@ -445 +447 @@
     q_prof(:,2) = moist(i,kps:kpe,j,P_QC) / (1.d0 + moist(i,kps:kpe,j,P_QV))
     ! conversion from mass mixing ratio in WRF to specific concentration in Physiq
+ELSE IF (TRACER_MODE == 61) 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_mu_m0as) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4))
+    q_prof(:,2) = SCALAR(i,kps:kpe,j,P_mu_m3as) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4))
+    q_prof(:,3) = SCALAR(i,kps:kpe,j,P_mu_m0af) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4))
+    q_prof(:,4) = SCALAR(i,kps:kpe,j,P_mu_m3af) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4))
+    q_prof(:,5) = SCALAR(i,kps:kpe,j,P_mu_m0n) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4))
+    q_prof(:,6) = SCALAR(i,kps:kpe,j,P_mu_m3n) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4))
+    q_prof(:,7) = SCALAR(i,kps:kpe,j,P_mu_m3CH4) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4))
+    q_prof(:,8) = SCALAR(i,kps:kpe,j,P_CH4) / (1.d0 + SCALAR(i,kps:kpe,j,P_CH4))
 ELSE
     q_prof(:,1:nq) = SCALAR(i,kps:kpe,j,2:nq+1)  !! the names were set above !! one dummy tracer in WRF
@@ -651 +663 @@
             DQICE,DQVAP,REEVAP,SURFRAIN,&
             ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,&
-            FLUXSURF_LW,FLXGRD,DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF,LATENT_HF)
+            FLUXSURF_LW,FLXGRD,DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF,LATENT_HF,&
+            DT_COND)
 !!!
 !print *, '** ',planet_type,'** OUTPUT PHYSICS DONE'
@@ -739 +752 @@
       SCALAR(i,kps:kpe,j,P_MARKER) = SCALAR(i,kps:kpe,j,P_MARKER)*exp(-dt/tau_decay)
       SCALAR(i,1,j,P_MARKER) = 1. !! this tracer is emitted in the surface layer
+    CASE(61) !emoisan tobechecked
+      scalar(i,kps:kpe,j,P_mu_m0as)=scalar(i,kps:kpe,j,P_mu_m0as) &
+           +zdqfi_omp(subs,kps:kpe,1)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4))
+      scalar(i,kps:kpe,j,P_mu_m3as)=scalar(i,kps:kpe,j,P_mu_m3as) &
+           +zdqfi_omp(subs,kps:kpe,2)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4))
+      scalar(i,kps:kpe,j,P_mu_m0af)=scalar(i,kps:kpe,j,P_mu_m0af) &
+           +zdqfi_omp(subs,kps:kpe,3)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4))
+      scalar(i,kps:kpe,j,P_mu_m3af)=scalar(i,kps:kpe,j,P_mu_m3af) &
+           +zdqfi_omp(subs,kps:kpe,4)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4))
+      scalar(i,kps:kpe,j,P_mu_m0n)=scalar(i,kps:kpe,j,P_mu_m0n) &
+           +zdqfi_omp(subs,kps:kpe,5)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4))
+      scalar(i,kps:kpe,j,P_mu_m3n)=scalar(i,kps:kpe,j,P_mu_m3n) &
+           +zdqfi_omp(subs,kps:kpe,6)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4))
+      scalar(i,kps:kpe,j,P_mu_m3CH4)=scalar(i,kps:kpe,j,P_mu_m3CH4) &
+           +zdqfi_omp(subs,kps:kpe,7)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4))
+      scalar(i,kps:kpe,j,P_CH4)=scalar(i,kps:kpe,j,P_CH4) &
+           +zdqfi_omp(subs,kps:kpe,8)*dt * (1.d0+scalar(i,kps:kpe,j,P_CH4))
     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

trunk/WRF.COMMON/INTERFACES_V4/module_model_constants.F

@@ -26 +26 @@
    REAL :: cp !          = 7.*r_d/2.  ! 
 
-   REAL    , PARAMETER :: r_v          = 461.6      ! gas constant for water vapor (J deg^-1 kg^-1)
+   REAL :: r_v !         = 461.6      ! gas constant for water vapor (J deg^-1 kg^-1)
    REAL :: cv !          = cp-r_d     ! Specific heat of air at contant volume (J deg^-1 kg^-1)
    REAL :: cpv !         = 4.*r_v
    REAL :: cvv !         = cpv-r_v    ! 
    REAL :: cvpm !        = -cv/cp
-   REAL    , PARAMETER :: cliq         = 4190.      ! specific heat of liquid water at 0^oC
-   REAL    , PARAMETER :: cice         = 2106.      ! specific heat of ice at 0^oC
-   REAL    , PARAMETER :: psat         = 610.78
+   REAL :: cliq !        = 4190.      ! specific heat of liquid water at 0^oC
+   REAL :: cice !        = 2106.      ! specific heat of ice at 0^oC
+   REAL :: psat !        = 610.78
    REAL :: rcv !         = r_d/cv     ! 
    REAL :: rcp !         = r_d/cp
@@ -49 +49 @@
    REAL :: reradius !     = 1./6370.0e03  ! reciprocal of earth radius (m^-1) 
 
-   REAL    , PARAMETER :: asselin      = .025
+   REAL    , PARAMETER :: asselin      = .025 ! for asselin filter?
 !   REAL    , PARAMETER :: asselin      = .0
-   REAL    , PARAMETER :: cb           = 25.
+   REAL :: cb !          = 25.
 
    REAL    , PARAMETER :: XLV0         = 3.15E6       !  constant defined for calculation of latent heating 
@@ -266 +266 @@
      g            = 1.35
      r_d          = 298.734319568
+     r_v          = 518.          ! gas constant for methane vapor at Titan (J.K-1.kg-1)
      cp           = 1038.72627727
      t0           = 94. ! earth : 300