Changeset 4745 for LMDZ6/trunk/libf/phylmdiso

Timestamp:

Nov 7, 2023, 8:27:07 AM (21 months ago)

Author:

evignon

Message:

nettoyage et corrections dans les routines atke pour utilisation en 3D (terre + mars)

Location:

LMDZ6/trunk/libf/phylmdiso

Files:

• pbl_surface_mod.F90 (modified) (3 diffs)
• physiq_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmdiso/pbl_surface_mod.F90

@@ -1977 +1977 @@
 
         IF (iflag_pbl>=50) THEN
-
-        CALL call_atke(dtime,knon,klev,ycdragm, ycdragh,yus0,yvs0,yts,yu,yv,yt,yq, &
-             ypplay,ypaprs,ytke,ycoefm, ycoefh)
+        CALL call_atke(dtime,knon,klev,ycdragm(1:knon), ycdragh(1:knon),yus0(1:knon),yvs0(1:knon),yts(1:knon), &
+                  yu(1:knon,:),yv(1:knon,:),yt(1:knon,:),yq(1:knon,:),ypplay(1:knon,:),ypaprs(1:knon,:),       &
+                  ytke(1:knon,:),ycoefm(1:knon,:), ycoefh(1:knon,:))
 
         ELSE
@@ -2020 +2020 @@
 
         IF (iflag_pbl>=50) THEN
-
-        CALL call_atke(dtime,knon,klev,ycdragm_x,ycdragh_x,yus0,yvs0,yts_x,yu_x,yv_x,yt_x,yq_x, &
-             ypplay,ypaprs,ytke_x,ycoefm_x, ycoefh_x)
+     
+        CALL call_atke(dtime,knon,klev,ycdragm_x(1:knon),ycdragh_x(1:knon),yus0(1:knon),yvs0(1:knon),yts_x(1:knon),    &
+                       yu_x(1:knon,:),yv_x(1:knon,:),yt_x(1:knon,:),yq_x(1:knon,:),ypplay(1:knon,:),ypaprs(1:knon,:),  &
+                       ytke_x(1:knon,:),ycoefm_x(1:knon,:), ycoefh_x(1:knon,:))
 
         ELSE
-
 
         CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
@@ -2060 +2060 @@
 
         IF (iflag_pbl>=50) THEN
-
-        CALL call_atke(dtime,knon,klev,ycdragm_w,ycdragh_w,yus0,yvs0,yts_w,yu_w,yv_w,yt_w,yq_w, &
-             ypplay,ypaprs,ytke_w,ycoefm_w, ycoefh_w)
+        
+        CALL call_atke(dtime,knon,klev,ycdragm_w(1:knon),ycdragh_w(1:knon),yus0(1:knon),yvs0(1:knon),yts_w(1:knon), &
+                yu_w(1:knon,:),yv_w(1:knon,:),yt_w(1:knon,:),yq_w(1:knon,:),ypplay(1:knon,:),ypaprs(1:knon,:),      &
+                ytke_w(1:knon,:),ycoefm_w(1:knon,:), ycoefh_w(1:knon,:))
 
         ELSE

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -1342 +1342 @@
 
     REAL pi
+    REAL viscom, viscoh
     INTEGER ieru
 
@@ -1932 +1933 @@
        CALL wake_ini(rg,rd,rv,prt_level)
        CALL yamada_ini(klon,lunout,prt_level)
-       CALL atke_ini(prt_level, lunout, RG, RD, RPI, RCPD, RV)
+       viscom=1.46E-5
+       viscoh=2.06E-5
+       CALL atke_ini(RG, RD, RPI, RCPD, RV, viscom, viscoh)
        CALL thermcell_ini(iflag_thermals,prt_level,tau_thermals,lunout, &
    &    RG,RD,RCPD,RKAPPA,RLVTT,RETV)