Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/coagulate.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/StratAer/coagulate.F90 (modified) (6 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/coagulate.F90

@@ -1 @@
-!
+
 ! $Id$
-!
+
 SUBROUTINE COAGULATE(pdtcoag,mdw,tr_seri,t_seri,pplay,dens_aer,is_strato)
   !     -----------------------------------------------------------------------
-  !
+
   !     Author : Christoph Kleinschmitt (with Olivier Boucher)
   !     ------
-  !
+
   !     purpose
   !     -------
-  !
+
   !     interface
   !     ---------
@@ -17 +17 @@
   !       tr_seri        tracer mixing ratios               [kg/kg]
   !       mdw             # or mass median diameter          [m]
-  !
+
   !     method
   !     ------
-  !
+
   !     -----------------------------------------------------------------------
 
@@ -204 +204 @@
   DO j=1, nbtr_bin
   DO i=1, nbtr_bin
-!
+
      num=4.*RPI*(radiuswet(i)+radiuswet(j))*(Di(i)+Di(j))
      denom=(radiuswet(i)+radiuswet(j))/(radiuswet(i)+radiuswet(j)+sqrt(delta(i)**2.+delta(j)**2.))+ &
    4.*(Di(i)+Di(j))/(sqrt(thvelpar(i)**2.+thvelpar(j)**2.)*(radiuswet(i)+radiuswet(j)))
      beta(i,j)=num/denom
-!
+
 !--compute enhancement factor due to van der Waals forces
    IF (ok_vdw == 0) THEN      !--no enhancement factor
@@ -222 +222 @@
       EvdW = 1. + SQRT(AvdWi/3.)/(1.+bvdW0*SQRT(AvdWi)) + bvdW1*xvdW + bvdW3*xvdW**3.
    ENDIF
-!
+
    beta(i,j)=beta(i,j)*EvdW
 
@@ -247 +247 @@
        numi=0.0
        DO i=1, k-1
-!           
+
 !           see Jacobson: " In order to conserve volume and volume concentration (which
 !           coagulation physically does) while giving up some accuracy in number concentration" 
-!
+
 !           Coagulation of i and j giving k
 !           with V(i) and then V(j) because it considers i,j and j,i with the double loop
-!
+
 !           BUT WHY WET VOLUME V(i) in old STRATAER? tracers are already dry aerosols and coagulation 
 !           kernel beta(i,j) accounts for wet aerosols -> reply below
-!
+
 !             numi=numi+ff(i,j,k)*beta(i,j)*V(i)*tr_tp1(ilon,ilev,i)*tr_t(ilon,ilev,j)
             numi=numi+ff(i,j,k)*beta(i,j)*Vdry(i)*tr_tp1(ilon,ilev,i)*tr_t(ilon,ilev,j)
@@ -266 +266 @@
 !      tr_tp1(ilon,ilev,k)=(V(k)*tr_t(ilon,ilev,k)+pdtcoag*num)/( (1.+pdtcoag*denom)*V(k) )
     tr_tp1(ilon,ilev,k)=(Vdry(k)*tr_t(ilon,ilev,k)+pdtcoag*num)/( (1.+pdtcoag*denom)*Vdry(k) )
-!
+
 !       In constant composition (no dependency on aerosol size because no kelvin effect)
 !       V(l)= (f_r_wet(ilon,ilev)**3)*((mdw(l)/2.)**3) = (f_r_wet(ilon,ilev)**3)*Vdry(i)