Changeset 2259 for trunk/LMDZ.VENUS/libf

Timestamp:

Mar 13, 2020, 3:23:08 PM (5 years ago)

Author:

emillour

Message:

Venus GCM:
Some cleanup to be able to compile with modern gfortran compiler (5.5+) which
complained about the need for an explicit interface in sed_and_prod_mad.F90
for get_weff().
Turned it into a module in the process.
EM

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• cloudvenus/sed_and_prod_mad.F90 (modified) (14 diffs)
• physiq_mod.F (modified) (1 diff)

trunk/LMDZ.VENUS/libf/phyvenus/cloudvenus/sed_and_prod_mad.F90

@@ +1 @@
+MODULE SED_AND_PROD_MAD
+
+IMPLICIT NONE
+
 !
 ! All parameter variable should be modified regarding the sudied system.
@@ -30 +34 @@
 ! ft = Theoretical Settling velocity at each vertical levels with Fiadero correction
 
+CONTAINS
 
 !============================================================================
@@ -36 +41 @@
 FUNCTION get_r0_1(M0,M3) RESULT(ret)
 
-  use free_param
-  use donnees
+  use free_param, only: sigma1
 
   IMPLICIT NONE
@@ -57 +61 @@
 FUNCTION get_r0_2(M0,M3) RESULT(ret)
 
-  use free_param
-  use donnees
+  use free_param, only: sigma2
 
   IMPLICIT NONE
@@ -81 +84 @@
   !! Compute effective gravitational acceleration.
 
-  use free_param
-  use donnees
+  use donnees, only: g0, rpla
 
   IMPLICIT NONE
@@ -98 +100 @@
 SUBROUTINE aer_production(dt,nlay,plev,zlev,zlay,dm0_1,dm3_1)
 
-  use free_param
-  use donnees
+  use free_param, only: sigz, p_aer, r_aer, rho_aer, sig_aer, tx_prod
+  use donnees, only: pi
 
   IMPLICIT NONE
@@ -166 +168 @@
   !! from ancient times. It is based on \cite{toon1988b,fiadeiro1977,turco1979} and is an update of the algorithm
   !! originally implemented in the LMDZ-Titan 2D GCM.
-
-  use free_param
-  use donnees
 
   IMPLICIT NONE
@@ -240 +239 @@
   !! be computed.
 
-  use free_param
-  use donnees
+  use free_param, only: rho_aer
+  use donnees, only: akn, df, fiadero_max, fiadero_min, no_fiadero_w, rmono
   
   IMPLICIT NONE
@@ -267 +266 @@
   REAL                       :: af1,af2,ar1,ar2
   REAL                       :: csto,cslf,ratio,wdt,dzb
-  REAL                       :: rb2ra, VISAIR, FPLAIR, alpha_k, effg
+  REAL                       :: rb2ra, VISAIR, FPLAIR, alpha_k
   REAL, DIMENSION(nlay+1) :: zcorf
   ! ------------------
@@ -327 +326 @@
   !! It is assumed that the aerosols size-distribution is the same as the mode 1 of cloud drops.
   
-  use free_param
-  use donnees
+  use free_param, only: rho_aer, sig_aer
+  use donnees, only: pi, wsed_m0, wsed_m3
   
   IMPLICIT NONE
@@ -354 +353 @@
   REAL, DIMENSION(nlay+1):: ft,fdcor,wth
   REAL, DIMENSION(nlay)  :: m0vsed,m3vsed,r0
-  REAL                   :: m,n,p,get_r0_1
+  REAL                   :: m,n,p
   REAL, PARAMETER        :: fac = 4.D0/3.D0 * pi
   INTEGER                :: nla,nle,i
@@ -436 +435 @@
   !!     cld_flux = fac *ft * (m3ccn * rho_aer + SUM(rho_liq(:) + m3liq(:)))
   
-  use free_param
-  use donnees
+  use free_param, only: sigma1, sigma2
+  use donnees, only: pi, wsed_m0, wsed_m3
   
   IMPLICIT NONE
@@ -467 +466 @@
   REAL, DIMENSION(nlay+1):: ft, fdcor, wth
   REAL, DIMENSION(nlay)  :: r0, m3sum
-  REAL                   :: m, n, p, sig, get_r0_1, get_r0_2
+  REAL                   :: m, n, p, sig
   REAL, PARAMETER        :: fac = 4.D0/3.D0 * pi
   INTEGER                :: i, nla, nle, nc
@@ -540 +539 @@
 END SUBROUTINE drop_sedimentation
 
-!END MODULE SED_AND_PROD
+END MODULE SED_AND_PROD_MAD

trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F

@@ -76 +76 @@
       use vertical_layers_mod, only: pseudoalt
       use turb_mod, only : sens, turb_resolved
+      use sed_and_prod_mad, only: aer_sedimentation, drop_sedimentation
 #ifdef CPP_XIOS      
       use xios_output_mod, only: initialize_xios_output,