Changeset 3098 for LMDZ6/trunk/libf

Timestamp:

Dec 3, 2017, 4:08:33 PM (7 years ago)

Author:

oboucher

Message:

Cosmetic changes to the routine

File:

• LMDZ6/trunk/libf/phylmd/StratAer/micphy_tstep.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/StratAer/micphy_tstep.F90

@@ -25 +25 @@
   ! local variables in coagulation routine
   INTEGER, PARAMETER        :: nbtstep=4  ! Max number of time steps in microphysics per time step in physics
-  INTEGER                   :: it,ilon,ilev,IK,count_tstep
+  INTEGER                   :: it,ilon,ilev,count_tstep
   REAL                      :: rhoa !H2SO4 number density [molecules/cm3]
   REAL                      :: ntot !total number of molecules in the critical cluster (ntot>4)
@@ -43 +43 @@
   REAL H2SO4_sat(nbtr_bin)
 
-  DO IK=1,nbtr_bin
-    Vbin(IK)=4.0*RPI*((mdw(IK)/2.)**3)/3.0
+  DO it=1,nbtr_bin
+    Vbin(it)=4.0*RPI*((mdw(it)/2.)**3)/3.0
   ENDDO
 
@@ -95 +95 @@
              & RRSI,Vbin,FL,ASO4,DNDR)
       ! consider only condensation (positive FL)
-      DO IK=1,nbtr_bin
-        FL(IK)=MAX(FL(IK),0.)
+      DO it=1,nbtr_bin
+        FL(it)=MAX(FL(it),0.)
       ENDDO
       ! compute total H2SO4 cond flux for all particles
       cond_evap_rate=0.0
-      DO IK=1, nbtr_bin
-        cond_evap_rate=cond_evap_rate+tr_seri(ilon,ilev,IK+nbtr_sulgas)*FL(IK)*mH2SO4mol
+      DO it=1, nbtr_bin
+        cond_evap_rate=cond_evap_rate+tr_seri(ilon,ilev,it+nbtr_sulgas)*FL(it)*mH2SO4mol
       ENDDO
       ! determine appropriate time step
@@ -131 +131 @@
            & RRSI,Vbin,FL,ASO4,DNDR)
     ! limit evaporation (negative FL) over one physics time step to H2SO4 content of the droplet
-    DO IK=1,nbtr_bin
-      FL(IK)=MAX(FL(IK)*pdtphys,0.-ASO4(IK))/pdtphys
+    DO it=1,nbtr_bin
+      FL(it)=MAX(FL(it)*pdtphys,0.-ASO4(it))/pdtphys
       ! consider only evap (negative FL)
-      FL(IK)=MIN(FL(IK),0.)
+      FL(it)=MIN(FL(it),0.)
     ENDDO
     ! compute total H2SO4 evap flux for all particles
     evap_rate=0.0
-    DO IK=1, nbtr_bin
-      evap_rate=evap_rate+tr_seri(ilon,ilev,IK+nbtr_sulgas)*FL(IK)*mH2SO4mol
+    DO it=1, nbtr_bin
+      evap_rate=evap_rate+tr_seri(ilon,ilev,it+nbtr_sulgas)*FL(it)*mH2SO4mol
     ENDDO
     ! update H2SO4 concentration after evap
@@ -158 +158 @@
     DO ilon=1, klon
     DO ilev=1, klev    
-    DO IK=1, nbtr
-      IF (tr_seri(ilon,ilev,IK).LT.0.0) THEN
-        PRINT *, 'micphy_tstep: negative concentration', tr_seri(ilon,ilev,IK), ilon, ilev, IK
+    DO it=1, nbtr
+      IF (tr_seri(ilon,ilev,it).LT.0.0) THEN
+        PRINT *, 'micphy_tstep: negative concentration', tr_seri(ilon,ilev,it), ilon, ilev, it
       ENDIF
     ENDDO