Changeset 3856

Timestamp:

Feb 25, 2021, 7:02:16 PM (4 years ago)

Author:

fhourdin

Message:

Introduction d'un ratqs (largeur relative de la distribution sous maille de l'eau)
interactif (travail de Louis D'Alençon, Frédéric Hourdin et Catherine Rio).
Contrôlé par le flag iflag_ratqs=10 (précédemment 4).
N'affecte pas les calculs si iflag_ratqs=4.
Fredho

Location:

LMDZ6/trunk/libf/phylmd

Files:

• calcratqs.F90 (modified) (4 diffs)
• phyetat0.F90 (modified) (2 diffs)
• phyredem.F90 (modified) (2 diffs)
• phys_state_var_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/calcratqs.F90

@@ -2 +2 @@
            iflag_ratqs,iflag_con,iflag_cld_th,pdtphys, &
            ratqsbas,ratqshaut,ratqsp0,ratqsdp, &
-           tau_ratqs,fact_cldcon,   &
+           tau_ratqs,fact_cldcon,wake_s, wake_deltaq,   &
            ptconv,ptconvth,clwcon0th, rnebcon0th,      &
            paprs,pplay,q_seri,zqsat,fm_therm,          &
-           ratqs,ratqsc)
+           ratqs,ratqsc,ratqs_inter)
 
 implicit none
@@ -27 +27 @@
 logical, dimension(klon,klev),intent(in) :: ptconv
 real, dimension(klon,klev),intent(in) :: rnebcon0th,clwcon0th
-
+real, dimension(klon,klev),intent(in) :: wake_deltaq,wake_s
 ! Output
-real, dimension(klon,klev),intent(inout) :: ratqs,ratqsc
+real, dimension(klon,klev),intent(inout) :: ratqs,ratqsc,ratqs_inter
 logical, dimension(klon,klev),intent(inout) :: ptconvth
 
@@ -124 +124 @@
          enddo
 
-      else if (iflag_ratqs==4) then
+      else if (iflag_ratqs==4) then 
          do k=1,klev
            ratqss(:,k)=ratqsbas+0.5*(ratqshaut-ratqsbas) &
@@ -131 +131 @@
          enddo
 
+      else if (iflag_ratqs==10) then ! ratqs interactif dépendant de la présence de poches froides
+         call calcratqs_inter(klon,klev,pdtphys,ratqsbas,wake_deltaq,wake_s,q_seri,ratqs_inter)
+         do k=1,klev
+            do i=1, klon
+                ratqss(i,k)=ratqs_inter(i,k)+0.5*(ratqshaut-ratqs_inter(i,k)) &
+                *( tanh( (ratqsp0-pplay(i,k))/ratqsdp) + 1.)
+            enddo
+         enddo
+
       endif
-
-
 
 

LMDZ6/trunk/libf/phylmd/phyetat0.F90

@@ -19 +19 @@
        wake_s, wake_dens, zgam, zmax0, zmea, zpic, zsig, &
        zstd, zthe, zval, ale_bl, ale_bl_trig, alp_bl, u10m, v10m, treedrg, &
-       ale_wake, ale_bl_stat, ds_ns, dt_ns, delta_sst, delta_sal
+       ale_wake, ale_bl_stat, ds_ns, dt_ns, delta_sst, delta_sal, ratqs_inter
 !FC
   USE geometry_mod, ONLY : longitude_deg, latitude_deg
@@ -435 +435 @@
   found=phyetat0_get(1,ale_bl_stat,"ALE_BL_STAT","ALE_BL_STAT",0.)
 
+! fisrtilp/Clouds
+  found=phyetat0_get(klev,ratqs_inter,"RATQS_INTER","Relative width of the lsc sugrid scale water",ratqs_inter)
+
 !===========================================
   ! Read and send field trs to traclmdz

LMDZ6/trunk/libf/phylmd/phyredem.F90

@@ -28 +28 @@
                                 du_gwd_rando, du_gwd_front, u10m, v10m, &
                                 treedrg, solswfdiff, delta_sal, ds_ns, dt_ns, &
-                                delta_sst
+                                delta_sst, ratqs_inter
 
   USE geometry_mod, ONLY : longitude_deg, latitude_deg
@@ -303 +303 @@
 
     CALL put_field(pass,"ALE_BL_STAT", "ALE_BL_STAT", ale_bl_stat)
+
+
+    ! fisrtilp/clouds
+    CALL put_field(pass,"RATQS_INTER","Relative width of the lsc sugrid scale water",ratqs_inter)
 
 

LMDZ6/trunk/libf/phylmd/phys_state_var_mod.F90

@@ -418 +418 @@
 !$OMP THREADPRIVATE(ccm)
 
-!!! nrlmd le 10/04/2012
       REAL,SAVE,ALLOCATABLE :: ale_bl_trig(:)
 !$OMP THREADPRIVATE(ale_bl_trig)
-!!! fin nrlmd le 10/04/2012
+
+      REAL,SAVE,ALLOCATABLE :: ratqs_inter(:,:)
+!$OMP THREADPRIVATE(ratqs_inter)
 
       REAL, ALLOCATABLE, SAVE:: du_gwd_rando(:, :), du_gwd_front(:, :)
@@ -648 +649 @@
       ALLOCATE(cg_aero_lw_rrtm(klon,klev,2,nbands_lw_rrtm))
       ALLOCATE(ccm(klon,klev,nbands))
-
-!!! nrlmd le 10/04/2012
       ALLOCATE(ale_bl_trig(klon))
-!!! fin nrlmd le 10/04/2012
+      ALLOCATE(ratqs_inter(klon,klev))
       IF (ok_gwd_rando) THEN
         ALLOCATE(du_gwd_rando(klon, klev))
@@ -794 +793 @@
       if (ok_gwd_rando) DEALLOCATE(du_gwd_rando)
       if (.not. ok_hines .and. ok_gwd_rando) DEALLOCATE(du_gwd_front)
-       
-!!! nrlmd le 10/04/2012
       DEALLOCATE(ale_bl_trig)
-!!! fin nrlmd le 10/04/2012
+      DEALLOCATE(ratqs_inter)
 
       if (activate_ocean_skin >= 1) deALLOCATE(delta_sal, ds_ns, dt_ns, &

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -3447 +3447 @@
          iflag_ratqs,iflag_con,iflag_cld_th,pdtphys,  &
          ratqsbas,ratqshaut,ratqsp0, ratqsdp, &
-         tau_ratqs,fact_cldcon,   &
+         tau_ratqs,fact_cldcon,wake_s, wake_deltaq,   &
          ptconv,ptconvth,clwcon0th, rnebcon0th,     &
          paprs,pplay,q_seri,zqsat,fm_therm, &
-         ratqs,ratqsc)
-
+         ratqs,ratqsc,ratqs_inter)
 
     !