Changeset 4519 for LMDZ6/trunk/libf/phylmd

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/phylmd

Files:

• calcratqs.F90 (modified) (4 diffs)
• physiq_mod.F90 (modified) (1 diff)

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)