Changeset 4352 for LMDZ6/trunk

Timestamp:

Nov 20, 2022, 12:01:18 PM (3 years ago)

Author:

evignon

Message:

option pour conditionner l'activation des params SSO a un critere
sur la "planitude" de la maille afin de ne pas freiner de facon irrealiste
les vents sur les calottes
Travail de Valentin Wiener

Location:

LMDZ6/trunk/libf/phylmd

Files:

• clesphys.h (modified) (2 diffs)
• conf_phys_m.F90 (modified) (4 diffs)
• physiq_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/clesphys.h

@@ -7 +7 @@
 !
 !..include cles_phys.h
-!
+
+       ! threshold on zrel_mount to activate SSO schemes
+       REAL zrel_mount_t
        INTEGER iflag_cycle_diurne
        LOGICAL soil_model,new_oliq,ok_orodr,ok_orolf 
@@ -122 +124 @@
      &     , co2_ppm0                                                   &
 !FC
-     &     , Cd_frein                                                   &
+     &     , Cd_frein, zrel_mount_t                                     &
      &     , ecrit_LES                                                  &
      &     , ecrit_ins, ecrit_hf, ecrit_day                             &

LMDZ6/trunk/libf/phylmd/conf_phys_m.F90

@@ -226 +226 @@
     LOGICAL, SAVE :: ok_lic_cond_omp
     !
+    REAL, SAVE    :: zrel_mount_t_omp
     INTEGER, SAVE :: iflag_cycle_diurne_omp
     LOGICAL, SAVE :: soil_model_omp,new_oliq_omp
@@ -900 +901 @@
     CALL getin('ok_orolf', ok_orolf_omp)
 
+
+    !Config  Key  =  zrel_mount_t
+    !Config  Desc = zrel_mount_t
+    !Config  Def  = 0.
+    !Config  Help = Connais pas !
+    zrel_mount_t_omp = 0.
+    CALL getin('zrel_mount_t', zrel_mount_t_omp)
+
+
     !Config  Key  = ok_limitvrai
     !Config  Desc = Force la lecture de la bonne annee
@@ -2426 +2436 @@
     ok_orodr = ok_orodr_omp
     ok_orolf = ok_orolf_omp
+    zrel_mount_t=zrel_mount_t_omp
     ok_limitvrai = ok_limitvrai_omp
     nbapp_rad = nbapp_rad_omp
@@ -2851 +2862 @@
     WRITE(lunout,*) ' ok_orodr=',ok_orodr
     WRITE(lunout,*) ' ok_orolf=',ok_orolf
+    WRITE(lunout,*) ' zrel_mount_t=',zrel_mount_t
     WRITE(lunout,*) ' ok_limitvrai=',ok_limitvrai
     WRITE(lunout,*) ' nbapp_rad=',nbapp_rad

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -972 +972 @@
 
     REAL picefra(klon,klev)
+    REAL zrel_mount(klon)
     !IM cf. AM 081204 END
     !
@@ -4547 +4548 @@
        DO i=1,klon
           itest(i)=0
-          !        IF ((zstd(i).gt.10.0)) THEN
-          IF (((zpic(i)-zmea(i)).GT.100.).AND.(zstd(i).GT.10.0)) THEN
+          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
+          !    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
              itest(i)=1
              igwd=igwd+1
@@ -4597 +4601 @@
        DO i=1,klon
           itest(i)=0
-          IF ((zpic(i)-zmea(i)).GT.100.) THEN
+          !zrel_mount: relative mountain height wrt relief explained by mean slope
+          ! -> condition on zrel_mount 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
              itest(i)=1
              igwd=igwd+1