Changeset 4519 for LMDZ6/trunk/libf/filtrez/filtreg_mod.F90

Timestamp:

Apr 24, 2023, 6:11:11 PM (17 months ago)

Author:

evignon

Message:

modif for polar runs: control of the maximum latitude at which the polar filter is active (new key in .def files)
+ dependency of ratqs upon resolution (cell area)

File:

• LMDZ6/trunk/libf/filtrez/filtreg_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/filtrez/filtreg_mod.F90

@@ -16 +16 @@
   USE serre_mod, ONLY: alphax
   USE logic_mod, ONLY: fxyhypb, ysinus
+  USE comconst_mod, ONLY: maxlatfilter
+
     !    ... H. Upadhyaya, O.Sharma   ...
     !
@@ -125 +127 @@
     ! as length dx becomes of the same size as dy 
     !
+    ! if maxlatfilter >0, prescribe the colat0 value from the .def files
+   
+    IF (maxlatfilter .LT. 0.) THEN
+
     colat0  =  MIN( 0.5, dymin/dxmin )
     ! colat0  =  1.
@@ -133 +139 @@
        alphax = 0.
     ENDIF
+
+    ELSE
+
+    colat0=(90.0-maxlatfilter)/180.0*pi        
+
+    ENDIF
+    
+    
+    
     !
     PRINT 50, colat0,alphax