Changeset 4161

Timestamp:

May 20, 2022, 2:54:07 PM (2 years ago)

Author:

fhourdin

Message:

For replay_param.sh

Location:

LMDZ6/trunk/libf/phylmd

Files:

• thermcell_plume.F90 (modified) (3 diffs)
• thermcell_plume_6A.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/thermcell_plume.F90

@@ -63 +63 @@
 
 
-      REAL wa_moy(ngrid,nlay+1)
-      REAL entr(ngrid,nlay),detr(ngrid,nlay)
-      REAL ztv_est(ngrid,nlay)
-      REAL zqla_est(ngrid,nlay)
-      REAL zta_est(ngrid,nlay)
-      REAL ztemp(ngrid),zqsat(ngrid)
+      REAL,dimension(ngrid,nlay+1) :: wa_moy
+      REAL,dimension(ngrid,nlay) :: entr,detr
+      REAL,dimension(ngrid,nlay) :: ztv_est
+      REAL,dimension(ngrid,nlay) :: zqla_est
+      REAL,dimension(ngrid,nlay) :: zta_est
+      REAL,dimension(ngrid) :: ztemp,zqsat
       REAL zdw2,zdw2bis
       REAL zw2modif
       REAL zw2fact,zw2factbis
-      REAL zeps(ngrid,nlay)
-
-      REAL    wmaxa(ngrid)
+      REAL,dimension(ngrid,nlay) :: zeps
+
+      REAL,dimension(ngrid) :: wmaxa
 
       INTEGER ig,l,k,lt,it,lm,nbpb
 
-      real zdz,zbuoy(ngrid,nlay),zalpha,gamma(ngrid,nlay),zdqt(ngrid,nlay),zw2m
-      real zbuoyjam(ngrid,nlay),zdqtjam(ngrid,nlay)
+      real,dimension(ngrid,nlay) :: zbuoy,gamma,zdqt
+      real zdz,zalpha,zw2m
+      real,dimension(ngrid,nlay) :: zbuoyjam,zdqtjam
       real zdz2,zdz3,lmel,entrbis,zdzbis
-      real d_temp(ngrid)
+      real,dimension(ngrid) :: d_temp
       real ztv1,ztv2,factinv,zinv,zlmel
       real zlmelup,zlmeldwn,zlt,zltdwn,zltup
@@ -90 +91 @@
       real eps
       logical Zsat
-      LOGICAL active(ngrid),activetmp(ngrid)
+      LOGICAL,dimension(ngrid) :: active,activetmp
       REAL fact_gamma,fact_gamma2,fact_epsilon2
 
 
-      REAL c2(ngrid,nlay)
+      REAL,dimension(ngrid,nlay) :: c2
 
       if (ngrid==1) print*,'THERMCELL PLUME MODIFIE 2014/07/11'
@@ -443 +444 @@
 
 
-     return 
+   include "dump_param.h" ! replay automatic include 
+ RETURN
      end

LMDZ6/trunk/libf/phylmd/thermcell_plume_6A.F90

@@ -55 +55 @@
       REAL zw2modif
       REAL zw2fact,zw2factbis
-      REAL zeps(ngrid,nlay)
-
-      REAL    wmaxa(ngrid)
+      REAL,dimension(ngrid,nlay) :: zeps
+
+      REAL, dimension(ngrid) ::    wmaxa(ngrid)
 
       INTEGER ig,l,k,lt,it,lm
@@ -70 +70 @@
       real,dimension(ngrid,nlay) :: zta_est
 
-      real zdz,zbuoy(ngrid,nlay),zalpha,gamma(ngrid,nlay),zdqt(ngrid,nlay),zw2m
-      real zbuoyjam(ngrid,nlay),zdqtjam(ngrid,nlay)
+      real,dimension(ngrid,nlay) :: zbuoy,gamma,zdqt
+      real zdz,zalpha,zw2m
+      real,dimension(ngrid,nlay) :: zbuoyjam,zdqtjam
       real zbuoybis,zdz2,zdz3,lmel,entrbis,zdzbis
-      real d_temp(ngrid)
+      real, dimension(ngrid) :: d_temp
       real ztv1,ztv2,factinv,zinv,zlmel
       real zlmelup,zlmeldwn,zlt,zltdwn,zltup
@@ -82 +83 @@
       real eps
       logical Zsat
-      LOGICAL active(ngrid),activetmp(ngrid)
+      LOGICAL,dimension(ngrid) :: active,activetmp
       REAL fact_gamma,fact_gamma2,fact_epsilon2
       REAL coefc
-      REAL c2(ngrid,nlay)
+      REAL,dimension(ngrid,nlay) :: c2
 
       if (ngrid==1) print*,'THERMCELL PLUME MODIFIE 2014/07/11'