Changeset 4246

Timestamp:

Sep 14, 2022, 3:14:16 PM (20 months ago)

Author:

evignon

Message:

controle dans les .def du presnivs limite en dessous duquel on ne guide plus
si guide_BL=false
Etienne, pour valentin

Location:

LMDZ6/trunk/libf

Files:

• dyn3d/guide_mod.F90 (modified) (3 diffs)
• dyn3dmem/guide_loc_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/dyn3d/guide_mod.F90

@@ -39 +39 @@
   REAL, PRIVATE, SAVE     :: lon_min_g,lon_max_g
   REAL, PRIVATE, SAVE     :: tau_lon,tau_lat
+
+  REAL, PRIVATE, SAVE     :: plim_guide_BL
 
   REAL, ALLOCATABLE, DIMENSION(:), PRIVATE, SAVE     :: alpha_u,alpha_v 
@@ -113 +115 @@
     CALL getpar('gamma4',.false.,gamma4,'Zone sans rappel elargie')
     CALL getpar('guide_BL',.true.,guide_BL,'guidage dans C.Lim')
+    CALL getpar('plim_guide_BL',85000.,plim_guide_BL,'BL top presnivs value')
+
 
 ! Sauvegarde du for�age
@@ -399 +403 @@
         else
             do l=1,llm
-                alpha_pcor(l)=(1.+tanh((0.85-presnivs(l)/preff)/0.05))/2.
+               alpha_pcor(l)=(1.+tanh(((plim_guide_BL-presnivs(l))/preff)/0.05))/2.
             enddo
         endif

LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.F90

@@ -40 +40 @@
   REAL, PRIVATE, SAVE     :: lon_min_g,lon_max_g
   REAL, PRIVATE, SAVE     :: tau_lon,tau_lat
+
+  REAL, PRIVATE, SAVE     :: plim_guide_BL
 
   REAL, ALLOCATABLE, DIMENSION(:), PRIVATE, SAVE     :: alpha_u,alpha_v 
@@ -118 +120 @@
     CALL getpar('gamma4',.false.,gamma4,'Zone sans rappel elargie')
     CALL getpar('guide_BL',.true.,guide_BL,'guidage dans C.Lim')
-    
+    CALL getpar('plim_guide_BL',85000.,plim_guide_BL,'BL top presnivs value')
+
 ! Sauvegarde du for�age
     CALL getpar('guide_sav',.false.,guide_sav,'sauvegarde guidage')
@@ -453 +456 @@
         else
             do l=1,llm
-                alpha_pcor(l)=(1.+tanh((0.85-presnivs(l)/preff)/0.05))/2.
+                alpha_pcor(l)=(1.+tanh(((plim_guide_BL-presnivs(l))/preff)/0.05))/2.
             enddo
         endif