Changeset 2872 for trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd

Timestamp:

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

Author:

jleconte

Message:

Cleaning up in interface v4

Location:

trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd

Files:

• update_inputs_physiq_mod.F (modified) (5 diffs)
• update_outputs_physiq_mod.F (modified) (2 diffs)

trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd/update_inputs_physiq_mod.F

@@ -26 +26 @@
   REAL,INTENT(OUT) :: MY
   REAL :: sec,nsec
-  
+  print*,'GMT, elaps/3600=',GMT,elaps/3600.
   !IF (JULYR .le. 8999) THEN
-  !  if (tlocked .eqv. .false.) THEN
+    if (tlocked .eqv. .false.) THEN
       JH_cur_split = (GMT + elaps/3600.) !! universal time (0<JH_cur_split<1): JH_cur_split=0.5 at 12:00 UT 
       JH_cur_split = MODULO(JH_cur_split,24.)   !! the two arguments of MODULO must be of the same type 
@@ -51 +51 @@
   !    !ptime = lct_input - lon_input / 15. + elaps/3600.
   !    !pday =  INT((3600*(lct_input - lon_input / 15.) + elaps)/86400)
-  !  ELSE
-  !    JH_cur_split = lct_input - lon_input / 15. !+ elaps/1500.0
+  ELSE
+       JH_cur_split = 0.
+       JD_cur = 0.
   !    !pday =  INT((sec*(lct_input - lon_input / 15.)+ elaps)/36000)
   !    JD_cur =  INT((sec*(lct_input - lon_input / 15.))/3600)
@@ -63 +64 @@
   !    MY = 2024
   !  ENDIF
-  !ENDIF
+  ENDIF
 
 
@@ -86 +87 @@
     noms(1)="h2o_vap"
     noms(2)="h2o_ice"
-  else if (TRACER_MODE .eq. 22)  then 
-    nqtot=2
-    IF (.not.ALLOCATED(noms)) ALLOCATE(noms(nqtot)) !! est fait dans initracer normalement
-    noms(1)="co2_vap"
-    noms(2)="co2_ice"
   else
     nq=1
@@ -325 +321 @@
      SELECT CASE (TRACER_MODE)
        CASE(1)
-       qsurf(subs,2)=M_H2OICE(i,j)    !! logique avec noms(2) = 'h2o_ice' defini ci-dessus
+         qsurf(subs,2)=M_H2OICE(i,j)  !! logique avec noms(2) = 'h2o_ice' defini ci-dessus
                                       !! ----- retrocompatible ancienne physique
                                       !! ----- [H2O ice is last tracer in qsurf in LMD physics]
-       CASE(2)   
-       qsurf(subs,1)=0.                !! not coupled with lifting for the moment [non remobilise]
-       !CASE(3)
-       !qsurf(subs,1)=q_prof(1,1)                !!! temporaire, a definir       
-       !qsurf(subs,2)=q_prof(1,2)     
-       CASE(11)
-       qsurf(subs,2)=M_H2OICE(i,j)    !! logique avec noms(2) = 'h2o_ice' defini ci-dessus
-       qsurf(subs,3)=0.                !! not coupled with lifting for the moment [non remobilise]
-       CASE(12)
-       qsurf(subs,2)=M_H2OICE(i,j)    !! logique avec noms(2) = 'h2o_ice' defini ci-dessus
-       qsurf(subs,3)=0.                !! not coupled with lifting for the moment [non remobilise]
      END SELECT
 

trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd/update_outputs_physiq_mod.F

@@ -156 +156 @@
             ims,ime,jms,jme,kms,kme,&
             ips,ipe,jps,jpe,kps,kpe,&
-            SWDOWNZ,TAU_DUST,QSURFDUST,&
-            MTOT,ICETOT,TAU_ICE,&
-            HR_SW,HR_LW,HR_DYN,DT,DT_RAD,&
-            RDUST,VMR_ICE,RICE,&
+            HR_SW,HR_LW,HR_DYN,DT_RAD,&
             CLOUDFRAC,TOTCLOUDFRAC,RH,&
             DQICE,DQVAP,REEVAP,SURFRAIN,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,&
@@ -169 +166 @@
    INTEGER, INTENT(IN) :: ips,ipe,jps,jpe,kps,kpe
    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT) :: &
-     SWDOWNZ,TAU_DUST,QSURFDUST,&
-     MTOT,ICETOT,TAU_ICE,&
      TOTCLOUDFRAC,ALBEQ,FLUXTOP_DN,FLUXABS_SW,FLUXTOP_LW,FLUXSURF_SW,&
      FLUXSURF_LW,FLXGRD,LATENT_HF,REEVAP,SURFRAIN
    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT ) :: &
-     HR_SW,HR_LW,RDUST,VMR_ICE,RICE,CLOUDFRAC,HR_DYN,DT,DT_RAD,RH,DQICE,DQVAP,&
+     HR_SW,HR_LW,CLOUDFRAC,HR_DYN,DT_RAD,RH,DQICE,DQVAP,&
      DTLSC,DTRAIN,DT_MOIST,H2OICE_REFF
    INTEGER :: i,j,subs