Changeset 4458 for LMDZ6/trunk/libf/phylmdiso

Timestamp:

Mar 8, 2023, 3:54:13 PM (22 months ago)

Author:

evignon

Message:

mise des seuils d'activation des params de SSO sous flag
pour faciliter les tests de sensibilité à venir

File:

• LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -1071 +1071 @@
 
     REAL picefra(klon,klev)
-    REAL zrel_mount(klon)
+    REAL zrel_oro(klon)
     !IM cf. AM 081204 END
     !
@@ -5803 +5803 @@
        DO i=1,klon
           itest(i)=0
-          zrel_mount(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i)))
-          !zrel_mount: relative mountain height wrt relief explained by mean slope
-          ! -> condition on zrel_mount can deactivate the drag on tilted planar terrains
+          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)
-          IF (((zpic(i)-zmea(i)).GT.100.).AND.(zstd(i).GT.10.0).AND.(zrel_mount(i).GE.zrel_mount_t)) THEN
+          ! 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).GE.zrel_oro_t)) THEN
              itest(i)=1
              igwd=igwd+1
@@ -5860 +5862 @@
        DO i=1,klon
           itest(i)=0
-          !zrel_mount: relative mountain height wrt relief explained by mean slope
-          ! -> condition on zrel_mount can deactivate the lifting on tilted planar terrains
+          !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_mount(i)=zstd(i)/(max(zsig(i),1.E-8)*sqrt(cell_area(i)))
-          IF (((zpic(i)-zmea(i)).GT.100.).AND.(zrel_mount(i).GE.zrel_mount_t)) THEN
+          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).GE.zrel_oro_t)) THEN
              itest(i)=1
              igwd=igwd+1
@@ -6080 +6082 @@
 ! 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) THEN
+          IF ((zstd(i).GT.1.0).AND.(zrel_oro(i).GE.zrel_oro_t)) THEN
              itest(i)=1
              igwd=igwd+1
@@ -6092 +6094 @@
        DO i=1,klon
           itest(i)=0
-        IF (((zpic(i)-zmea(i)).GT.100.).AND.(zstd(i).GT.10.0)) THEN
+        IF (((zpic(i)-zmea(i)).GT.zpmm_orodr_t).AND.(zstd(i).GT.zstd_orodr_t).AND.(zrel_oro(i).GE.zrel_oro_t)) THEN
              itest(i)=1
              igwd=igwd+1