Changeset 5470 for LMDZ6/trunk/libf/phylmd

Timestamp:

Jan 11, 2025, 9:36:52 PM (6 weeks ago)

Author:

evignon

Message:

raffinement de la condition pour ne pas activer la param d'orographie sous maille
sur les terrains pentus non montagneux tels que les calottes. On raisonne désormais sur un proxy
du nombre de montagnes sous-maille. Travail de Valentin et Etienne

Location:

LMDZ6/trunk/libf/phylmd

Files:

• clesphys_mod_h.f90 (modified) (2 diffs)
• conf_phys_m.f90 (modified) (4 diffs)
• physiq_mod.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/clesphys_mod_h.f90

@@ -22 +22 @@
           , co2_ppm0                                                   &
           , tau_thermals                                               &
-          , Cd_frein, zrel_oro_t, zpmm_orodr_t, zpmm_orolf_t, zstd_orodr_t &
+          , Cd_frein, nm_oro_t, zpmm_orodr_t, zpmm_orolf_t, zstd_orodr_t &
           , ecrit_LES                                                  &
           , ecrit_ins, ecrit_hf, ecrit_day                             &
@@ -55 +55 @@
 
   ! threshold on to activate SSO schemes
-  ! threshold on to activate SSO schemes
-  REAL zrel_oro_t, zpmm_orodr_t, zpmm_orolf_t, zstd_orodr_t
+  REAL nm_oro_t, zpmm_orodr_t, zpmm_orolf_t, zstd_orodr_t
   INTEGER iflag_cycle_diurne
   LOGICAL soil_model, new_oliq, ok_orodr, ok_orolf

LMDZ6/trunk/libf/phylmd/conf_phys_m.f90

@@ -213 +213 @@
     LOGICAL, SAVE :: ok_lic_cond_omp
     !
-    REAL, SAVE    :: zrel_oro_t_omp, zstd_orodr_t_omp
+    REAL, SAVE    :: nm_oro_t_omp, zstd_orodr_t_omp
     REAL, SAVE    :: zpmm_orodr_t_omp, zpmm_orolf_t_omp
     INTEGER, SAVE :: iflag_cycle_diurne_omp
@@ -893 +893 @@
 
 
-    !Config  Key  =  zrel_oro_t
-    !Config  Desc = zrel_oro_t
-    !Config  Def  = 9999.
+    !Config  Key  =  nm_oro_t
+    !Config  Desc = nm_oro_t
+    !Config  Def  = -1
     !Config  Help = Connais pas !
-    zrel_oro_t_omp = 9999.
-    CALL getin('zrel_oro_t', zrel_oro_t_omp)
+    nm_oro_t_omp = -1.
+    CALL getin('nm_oro_t', nm_oro_t_omp)
 
     !Config  Key  =  zstd_orodr_t
@@ -2313 +2313 @@
     ok_orodr = ok_orodr_omp
     ok_orolf = ok_orolf_omp
-    zrel_oro_t=zrel_oro_t_omp
+    nm_oro_t=nm_oro_t_omp
     zstd_orodr_t=zstd_orodr_t_omp
     zpmm_orodr_t=zpmm_orodr_t_omp
@@ -2732 +2732 @@
     WRITE(lunout,*) ' ok_orodr=',ok_orodr
     WRITE(lunout,*) ' ok_orolf=',ok_orolf
-    WRITE(lunout,*) ' zrel_oro_t=',zrel_oro_t
+    WRITE(lunout,*) ' nm_oro_t=',nm_oro_t
     WRITE(lunout,*) ' zstd_orodr_t=',zstd_orodr_t
     WRITE(lunout,*) ' zpmm_orodr_t=',zpmm_orodr_t

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -1019 +1019 @@
 
     REAL picefra(klon,klev)
-    REAL zrel_oro(klon)
+    REAL nm_oro(klon)
     !IM cf. AM 081204 END
     !
@@ -4841 +4841 @@
     ! a l'echelle sous-maille:
     !
+
+    ! calculation of nm_oro
+    DO i=1,klon
+          ! nm_oro is a proxy for the number of subgrid scale mountains
+          ! -> condition on nm_oro can deactivate the lifting on tilted planar terrains
+          !    such as ice sheets (work by V. Wiener)
+          ! in such a case, the SSO scheme should activate only where nm_oro>0 i.e. by setting
+          ! nm_oro_t=0.
+          nm_oro(i)=zsig(i)*sqrt(cell_area(i)*(pctsrf(i,is_ter)+pctsrf(i,is_lic)))/(4.*MAX(zstd(i),1.e-8))-1.
+    ENDDO
+
     IF (prt_level .GE.10) THEN
        print *,' call orography ? ', ok_orodr
@@ -4851 +4862 @@
        DO i=1,klon
           itest(i)=0
-          zrel_oro(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i)))
-          !zrel_oro: relative mountain height wrt relief explained by mean slope
-          ! -> condition on zrel_oro can deactivate the drag on tilted planar terrains
-          !    such as ice sheets (work by V. Wiener)
           ! zpmm_orodr_t and zstd_orodr_t are activation thresholds set by F. Lott to
           ! earn computation time but they are not physical.
-          IF (((zpic(i)-zmea(i)).GT.zpmm_orodr_t).AND.(zstd(i).GT.zstd_orodr_t).AND.(zrel_oro(i).LE.zrel_oro_t)) THEN
+          IF (((zpic(i)-zmea(i)).GT.zpmm_orodr_t).AND.(zstd(i).GT.zstd_orodr_t).AND.(nm_oro(i).GT.nm_oro_t)) THEN
              itest(i)=1
              igwd=igwd+1
@@ -4906 +4913 @@
        DO i=1,klon
           itest(i)=0
-          !zrel_oro: relative mountain height wrt relief explained by mean slope
-          ! -> condition on zrel_oro can deactivate the lifting on tilted planar terrains
-          !    such as ice sheets (work by V. Wiener)
-          zrel_oro(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i)))
-          IF (((zpic(i)-zmea(i)).GT.zpmm_orolf_t).AND.(zrel_oro(i).LE.zrel_oro_t)) THEN
+          IF (((zpic(i)-zmea(i)).GT.zpmm_orolf_t).AND.(nm_oro(i).GT.nm_oro_t)) THEN
              itest(i)=1
              igwd=igwd+1
@@ -5151 +5154 @@
 ! car on peut s'attendre a ce que les petites echelles produisent aussi de la TKE
 ! Mais attention, cela ne va pas dans le sens de la conservation de l'energie!
-          IF ((zstd(i).GT.1.0) .AND.(zrel_oro(i).LE.zrel_oro_t)) THEN
+          IF ((zstd(i).GT.1.0) .AND.(nm_oro(i).GT.nm_oro_t)) THEN
              itest(i)=1
              igwd=igwd+1
@@ -5163 +5166 @@
        DO i=1,klon
           itest(i)=0
-        IF (((zpic(i)-zmea(i)).GT.zpmm_orodr_t).AND.(zstd(i).GT.zstd_orodr_t).AND.(zrel_oro(i).LE.zrel_oro_t)) THEN
+        IF (((zpic(i)-zmea(i)).GT.zpmm_orodr_t).AND.(zstd(i).GT.zstd_orodr_t).AND.(nm_oro(i).GT.nm_oro_t)) THEN
              itest(i)=1
              igwd=igwd+1