Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/lsc_scav_spl.F90

Timestamp:

Jul 22, 2024, 9:29:09 PM (2 months ago)

Author:

abarral

Message:

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/lsc_scav_spl.F90 (modified) (9 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/lsc_scav_spl.F90

@@ -12 +12 @@
   USE traclmdz_mod
   USE infotrac,ONLY : nbtr
-!  USE comgeomphy
   USE iophy
+  USE lmdz_yomcst
   IMPLICIT NONE 
 !=====================================================================
@@ -25 +25 @@
   include "dimensions.h"
   include "chem.h"
-  include "YOMCST.h"
-  include "YOECUMF.h" 
+  include "YOECUMF.h"
 
   REAL,INTENT(IN)                        :: pdtime ! time step (s)
@@ -56 +55 @@
  LOGICAL,SAVE :: debut=.true.
 !$OMP THREADPRIVATE(debut)
-!
+
 !JE  REAL,PARAMETER :: henry=1.4  ! constante de Henry en mol/l/atm ~1.4 for gases
   REAL,DIMENSION(nbtr) :: henry  ! constante de Henry en mol/l/atm ~1.4 for gases
@@ -99 +98 @@
   REAL           :: pr, ps, ice, water
   real :: conserv
-!
+
 !!!!!!!!!!!!!!!!!!!! choix lessivage !!!!!!!!!!!!!!!!!!!!!!!!
 !!  logical,save  :: inscav_fisrt
 !!! $OMP THREADPRIVATE(inscav_first)
-!
+
 !!!!!!!!!!!!!!!!!!!!!!!!!!!
   IF (debut) THEN
-!
+
 !  inscav_fisrt=.true.
 !  call getin('inscav_fisrt',inscav_fisrt)
@@ -114 +113 @@
 !   print*,'beta from Reddy and Bocuher 2004 (original version), inscav_fisrt=',inscav_fisrt
 !  endif
-!
+
 !JE      alpha_r=0.001        !  coefficient d'impaction pour la pluie
 !JE      alpha_s=0.01         !  coefficient d'impaction pour la neige  
@@ -123 +122 @@
 !     frac_aer=0.5 ~ droplet size shrinks by evap
       frac_aer=0.5
-!
 
 !JE to speed up, commented 20140219
-!
+
 !      OPEN(99,file='lsc_scav_param.data',status='old', &
 !                form='formatted',err=9999)
@@ -139 +137 @@
 !      CLOSE(99)
 !9999  Continue
-!
+
 !   print*,'JE alpha_r',alpha_r
 !   print*,'JE alpha_s',alpha_s
@@ -147 +145 @@
 !   print*,'frac_coar_scav',frac_coar_scav
 !   print*,'frac_aer ev',frac_aer
-!
+
 ! JE endcomment
-!
+
   ENDIF !(debut)
 !!!!!!!!!!!!!!!!!!!!!!!!!!!
-!
+
 ! initialization
   dxin=0.
@@ -223 +221 @@
    endif ! (iflag_lscav .eq. 4)
    beta_v1(i,k)=beta    !! for output
-!
+
       dxin=tr_seri(i,k,it)*(exp(-scav(i,k)*beta*pdtime)-1.)
 !      his_dh(i)=his_dh(i)-dxin*zrho(i,k)*zdz(i,k)/pdtime !  kg/m2/s