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

Timestamp:

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

Author:

jleconte

Message:

Changed nomenclature for planet variables from m_ to p_ in interface V4

Location:

trunk/WRF.COMMON/INTERFACES_V4/dynphy_wrf_generic_lmd

Files:

• callphysiq_mod.F (modified) (1 diff)
• update_inputs_physiq_mod.F (modified) (18 diffs)
• update_outputs_physiq_mod.F (modified) (9 diffs)

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

@@ -40 +40 @@
   ! these two quantities do not even have the same units (kg/kg and kg/m^2)!!!
   ! seems it is the same for some other physics. It is probably a bug.
-  print*,'zdt_split',zdt_split
 
   call allocate_comm_wrf(klon,llm)

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
@@ -93 +92 @@
   endif
   nq=nqtot
-  print*,'in update_inputs_physiq_tracers: noms: ',noms
-  print*,'nq',nq
-  print*,'TRACER_MODE',TRACER_MODE
 
 END SUBROUTINE update_inputs_physiq_tracers
@@ -245 +241 @@
             ips,ipe,jps,jpe,&
             JULYR,TRACER_MODE,&
-            M_ALBEDO,CST_AL,&
-            P_TSURF,M_EMISS,M_CO2ICE,&
-            M_GW,M_Z0,CST_Z0,&
-            M_H2OICE,&
+            P_ALBEDO,CST_AL,&
+            P_TSURF,P_EMISS,P_CO2ICE,&
+            P_GW,P_Z0,CST_Z0,&
+            P_H2OICE,&
             phisfi_val)
 
@@ -262 +258 @@
    REAL, INTENT(IN  ) :: CST_AL, phisfi_val, CST_Z0
    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN)  :: &
-     M_ALBEDO,P_TSURF,M_EMISS,M_CO2ICE,M_H2OICE,M_Z0
-   REAL, DIMENSION( ims:ime, 5, jms:jme ), INTENT(IN   )  :: M_GW  
+     P_ALBEDO,P_TSURF,P_EMISS,P_CO2ICE,P_H2OICE,P_Z0
+   REAL, DIMENSION( ims:ime, 5, jms:jme ), INTENT(IN   )  :: P_GW  
 
    DO j = jps,jpe
@@ -282 +278 @@
      IF (JULYR .le. 8999) THEN 
       IF (CST_AL == 0) THEN
-       albedodat(subs)=M_ALBEDO(i,j)
+       albedodat(subs)=P_ALBEDO(i,j)
       ELSE 
        albedodat(subs)=CST_AL
@@ -311 +307 @@
      !-----------------------------------------------!
      tsurf(subs) = P_TSURF(i,j)
-     emis(subs) = M_EMISS(i,j)     
+     emis(subs) = P_EMISS(i,j)     
      !do i=1,noceanmx
      tslab(subs,:)=tsurf(subs)
@@ -321 +317 @@
      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)=P_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]
@@ -353 +349 @@
             ips,ipe,jps,jpe,&
             JULYR,nsoil,&
-            M_TI,CST_TI,&
-            M_ISOIL,M_DSOIL,&
-            M_TSOIL,P_TSURF)
+            P_TI,CST_TI,&
+            P_ISOIL,P_DSOIL,&
+            P_TSOIL,P_TSURF)
 
    use comsoil_h, only: inertiedat,mlayer,layer,volcapa
@@ -365 +361 @@
    REAL, INTENT(IN  ) :: CST_TI 
    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN)  :: &
-     M_TI, P_TSURF
+     P_TI, P_TSURF
    REAL, DIMENSION( ims:ime, nsoil, jms:jme ), INTENT(IN)  :: &
-     M_TSOIL, M_ISOIL, M_DSOIL
+     P_TSOIL, P_ISOIL, P_DSOIL
    REAL :: inertiedat_val
    REAL :: lay1,alpha
@@ -384 +380 @@
      IF (JULYR .le. 8999) THEN
       IF (CST_TI == 0) THEN
-       inertiedat_val=M_TI(i,j)
+       inertiedat_val=P_TI(i,j)
       ELSE
        inertiedat_val=CST_TI
@@ -396 +392 @@
      !--pb de dimensions???!!???
      IF (JULYR .le. 8999) THEN
-       inertiedat(subs,:)=M_ISOIL(i,:,j) !! verifier que cest bien hires TI en surface
-       mlayer(0:nsoil-1)=M_DSOIL(i,:,j)
+       inertiedat(subs,:)=P_ISOIL(i,:,j) !! verifier que cest bien hires TI en surface
+       mlayer(0:nsoil-1)=P_DSOIL(i,:,j)
      ELSE
         IF ( nsoil .lt. 18 ) THEN
@@ -425 +421 @@
      ! Deep soil temperatures ! 
      !------------------------!
-     IF (M_TSOIL(i,1,j) .gt. 0. .and. JULYR .le. 8999) THEN
-       tsoil(subs,:)=M_TSOIL(i,:,j)
+     IF (P_TSOIL(i,1,j) .gt. 0. .and. JULYR .le. 8999) THEN
+       tsoil(subs,:)=P_TSOIL(i,:,j)
      ELSE
        IF ( (i == ips) .AND. (j == jps) ) PRINT *,'** Mars ** no tsoil. set it to tsurf.'
@@ -457 +453 @@
             ips,ipe,jps,jpe,&
             RESTART,isles,&
-            M_Q2,M_WSTAR)
+            P_Q2,P_WSTAR)
 
    use turb_mod, only: q2,wstar,turb_resolved
@@ -465 +461 @@
    INTEGER, INTENT(IN) :: ips,ipe,jps,jpe
    INTEGER :: i,j,subs,nlast      
-   REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN)  :: M_WSTAR
-   REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(IN) :: M_Q2
+   REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN)  :: P_WSTAR
+   REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(IN) :: P_Q2
    LOGICAL, INTENT(IN ) :: RESTART,isles
 
@@ -487 +483 @@
       wstar(subs)=0.
      ELSE
-      q2(subs,:)=M_Q2(i,:,j)!
+      q2(subs,:)=P_Q2(i,:,j)!
       !q2(subs,:) = 1.e-3
-      wstar(subs)=M_WSTAR(i,j)
+      wstar(subs)=P_WSTAR(i,j)
      ENDIF
 
@@ -510 +506 @@
             ips,ipe,jps,jpe,&
             RESTART,&
-            M_FLUXRAD)
+            P_FLUXRAD)
 
    !use dimradmars_mod, only: fluxrad
@@ -518 +514 @@
    INTEGER, INTENT(IN) :: ips,ipe,jps,jpe
    INTEGER :: i,j,subs,nlast      
-   REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN)  :: M_FLUXRAD
+   REAL, DIMENSION( ims:ime, jms:jme ), INTENT(IN)  :: P_FLUXRAD
    LOGICAL, INTENT(IN ) :: RESTART
 
@@ -533 +529 @@
       fluxrad(subs)=0.
      ELSE
-      fluxrad(subs)=M_FLUXRAD(i,j)
+      fluxrad(subs)=P_FLUXRAD(i,j)
      ENDIF
      !! et fluxrad_sky ???!???

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

@@ -9 +9 @@
             ips,ipe,jps,jpe,&
             TRACER_MODE,&
-            P_TSURF,M_CO2ICE,&
-            M_H2OICE)
+            P_TSURF,P_CO2ICE,&
+            P_H2OICE)
 
    !use surfdat_h, only: tsurf, co2ice, qsurf
@@ -20 +20 @@
    INTEGER :: i,j,subs
    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)  :: &
-     P_TSURF,M_CO2ICE,M_H2OICE
+     P_TSURF,P_CO2ICE,P_H2OICE
 
    DO j = jps,jpe
@@ -46 +46 @@
             ips,ipe,jps,jpe,&
             nsoil,&
-            M_TSOIL)
+            P_TSOIL)
 
    !use comsoil_h, only: tsoil
@@ -55 +55 @@
    INTEGER :: i,j,subs
    REAL, DIMENSION( ims:ime, nsoil, jms:jme ), INTENT(INOUT)  :: &
-     M_TSOIL
-
-   DO j = jps,jpe
-   DO i = ips,ipe
-
-     !-----------------------------------!
-     ! 1D subscript for physics "cursor" !
-     !-----------------------------------!
-     subs = (j-jps)*(ipe-ips+1)+(i-ips+1)
-
-     !-------------------------------------------------------!
-     ! Save key variables for restart and output and nesting !  
-     !-------------------------------------------------------!
-     M_TSOIL(i,:,j) = tsoil(subs,:)
+     P_TSOIL
+
+   DO j = jps,jpe
+   DO i = ips,ipe
+
+     !-----------------------------------!
+     ! 1D subscript for physics "cursor" !
+     !-----------------------------------!
+     subs = (j-jps)*(ipe-ips+1)+(i-ips+1)
+
+     !-------------------------------------------------------!
+     ! Save key variables for restart and output and nesting !  
+     !-------------------------------------------------------!
+     P_TSOIL(i,:,j) = tsoil(subs,:)
 
    ENDDO
@@ -80 +80 @@
             ims,ime,jms,jme,&
             ips,ipe,jps,jpe,&
-            M_FLUXRAD)
+            P_FLUXRAD)
 
    !use dimradmars_mod, only: fluxrad
@@ -88 +88 @@
    INTEGER, INTENT(IN) :: ips,ipe,jps,jpe
    INTEGER :: i,j,subs
-   REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)  :: M_FLUXRAD
-
-   DO j = jps,jpe
-   DO i = ips,ipe
-
-     !-----------------------------------!
-     ! 1D subscript for physics "cursor" !
-     !-----------------------------------!
-     subs = (j-jps)*(ipe-ips+1)+(i-ips+1)
-
-     !-------------------------------------------------------!
-     ! Save key variables for restart and output and nesting !  
-     !-------------------------------------------------------!
-     M_FLUXRAD(i,j) = fluxrad(subs)
+   REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)  :: P_FLUXRAD
+
+   DO j = jps,jpe
+   DO i = ips,ipe
+
+     !-----------------------------------!
+     ! 1D subscript for physics "cursor" !
+     !-----------------------------------!
+     subs = (j-jps)*(ipe-ips+1)+(i-ips+1)
+
+     !-------------------------------------------------------!
+     ! Save key variables for restart and output and nesting !  
+     !-------------------------------------------------------!
+     P_FLUXRAD(i,j) = fluxrad(subs)
 
    ENDDO
@@ -113 +113 @@
             ims,ime,jms,jme,kms,kme,&
             ips,ipe,jps,jpe,kps,kpe,&
-            M_Q2,M_WSTAR,&
+            P_Q2,P_WSTAR,&
             HFMAX,ZMAX,USTM,HFX)
 
@@ -124 +124 @@
    INTEGER :: i,j,subs    
    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)  :: &
-     M_WSTAR,HFMAX,ZMAX,USTM,HFX
-   REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(INOUT) :: M_Q2
+     P_WSTAR,HFMAX,ZMAX,USTM,HFX
+   REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(INOUT) :: P_Q2
   
    DO j = jps,jpe
@@ -138 +138 @@
      ! Save key variables for restart and output and nesting !  
      !-------------------------------------------------------!
-     M_Q2(i,kps:kpe+1,j) = q2(subs,:)
-     M_WSTAR(i,j) = wstar(subs)
+     P_Q2(i,kps:kpe+1,j) = q2(subs,:)
+     P_WSTAR(i,j) = wstar(subs)
      !! output only (arrays already in phys modules)
      !HFMAX(i,j) = HFMAX_TH(subs)