Changeset 4519 for LMDZ6/trunk

Timestamp:

Apr 24, 2023, 6:11:11 PM (19 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)

Location:

LMDZ6/trunk/libf

Files:

• dyn3d/conf_gcm.F90 (modified) (2 diffs)
• dyn3d_common/comconst_mod.F90 (modified) (1 diff)
• dyn3dmem/conf_gcm.F90 (modified) (2 diffs)
• filtrez/filtreg_mod.F90 (modified) (3 diffs)
• phylmd/calcratqs.F90 (modified) (4 diffs)
• phylmd/physiq_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/dyn3d/conf_gcm.F90

@@ -15 +15 @@
   USE comconst_mod, ONLY: dissip_deltaz, dissip_factz, dissip_zref, &
                           iflag_top_bound, mode_top_bound, tau_top_bound, &
-                          ngroup
+                          ngroup, maxlatfilter
   USE logic_mod, ONLY: fxyhypb, iflag_phys, ok_etat0, ok_gradsfile, &
                        ok_guide, ok_limit, ok_strato, purmats, read_start, &
@@ -310 +310 @@
   CALL getin('dissip_deltaz',dissip_deltaz )
   CALL getin('dissip_zref',dissip_zref )
+
+  ! maxlatfilter
+  maxlatfilter=-1.0
+  CALL getin('maxlatfilter',maxlatfilter)
+  if (maxlatfilter > 90) &
+       call abort_gcm("conf_gcm", 'maxlatfilter should be <=90', 1)
+
 
   ! ngroup

LMDZ6/trunk/libf/dyn3d_common/comconst_mod.F90

@@ -26 +26 @@
       INTEGER iflag_top_bound ! sponge type
       INTEGER ngroup ! parameter to group points (along longitude) near poles
+      REAL maxlatfilter ! maximum latitude (in degrees) above which filter is active
       INTEGER mode_top_bound  ! sponge mode
       REAL tau_top_bound ! inverse of sponge characteristic time scale (Hz)

LMDZ6/trunk/libf/dyn3dmem/conf_gcm.F90

@@ -19 +19 @@
   USE comconst_mod, ONLY: dissip_deltaz, dissip_factz, dissip_zref, &
                           iflag_top_bound, mode_top_bound, tau_top_bound, &
-                          ngroup
+                          ngroup, maxlatfilter
   USE logic_mod, ONLY: fxyhypb, iflag_phys, ok_etat0, ok_gradsfile, &
                        ok_guide, ok_limit, ok_strato, purmats, read_start, &
@@ -340 +340 @@
   CALL getin('dissip_deltaz',dissip_deltaz )
   CALL getin('dissip_zref',dissip_zref )
+
+
+  !maxlatfilter
+  maxlatfilter = -1.0
+  CALL getin('maxlatfilter',maxlatfilter)
+  if (maxlatfilter > 90) &
+       call abort_gcm("conf_gcm", 'maxlatfilter should be <=90', 1)
+  
 
   ! ngroup

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

LMDZ6/trunk/libf/phylmd/calcratqs.F90

@@ -7 +7 @@
            qtc_cv, sigt_cv, zqsat,             &
            tke,tke_dissip,lmix,wprime, &
-           t2m,q2m,fm_therm, &
+           t2m,q2m,fm_therm,cell_area,&
            ratqs,ratqsc,ratqs_inter)
 
@@ -37 +37 @@
 real, dimension(klon,klev),intent(in) :: wake_deltaq,wake_s
 real, dimension(klon,nbsrf),intent(in) :: t2m,q2m
+real, dimension(klon), intent(in) :: cell_area
 ! Output
 real, dimension(klon,klev),intent(inout) :: ratqs,ratqsc,ratqs_inter
@@ -47 +48 @@
 real facteur,zfratqs1,zfratqs2
 real, dimension(klon,klev) :: ratqs_hetero,ratqs_oro,ratqs_tke
-
+real resol,resolmax,fact
 
 !-------------------------------------------------------------------------
@@ -142 +143 @@
            *( tanh( (ratqsp0-pplay(:,k))/ratqsdp) + 1.)
          enddo
+
+
+      else if (iflag_ratqs==5) then
+! Dependency of ratqs on model resolution
+! Audran, Meryl, Lea, Gwendal and Etienne
+! April 2023 
+        resolmax=sqrt(maxval(cell_area))
+         do k=1,klev
+            do i=1,klon
+              resol=sqrt(cell_area(i))
+              fact=sqrt(resol/resolmax)
+              ratqss(i,k)=ratqsbas*fact+0.5*(ratqshaut-ratqsbas)*fact &
+              *( tanh( (ratqsp0-pplay(i,k))/ratqsdp) + 1.)
+           enddo
+         enddo
+
 
        else if (iflag_ratqs .GT. 9) then

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -3624 +3624 @@
          ptconv,ptconvth,clwcon0th, rnebcon0th,     &
          paprs,pplay,t_seri,q_seri, qtc_cv, sigt_cv, zqsat, &
-         pbl_tke(:,:,is_ave),tke_dissip_ave,l_mix_ave,wprime_ave,t2m,q2m,fm_therm, &
+         pbl_tke(:,:,is_ave),tke_dissip_ave,l_mix_ave,wprime_ave,t2m,q2m,fm_therm,cell_area, &
          ratqs,ratqsc,ratqs_inter)