Changeset 4829

Timestamp:

Feb 21, 2024, 4:54:27 PM (3 months ago)

Author:

lebasn

Message:

StratAer?: remove unused code containing an uninitialized variable.

File:

• LMDZ6/trunk/libf/phylmd/StratAer/nucleation_tstep_mod.F90 (modified) (7 diffs)

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

@@ -32 +32 @@
   LOGICAL kinetic_n ! true if kinetic neutral nucleation (NOT IN USE)
   LOGICAL kinetic_i ! true if kinetic ion-induced nucleation (NOT IN USE)
-  REAL rhoatres ! threshold concentration of H2SO4 (1/cm^3) for neutral kinetic nucleation (NOT IN USE)
   REAL VH2SO4mol
   REAL ntot_i, x_i, jnuc_i, rc_i, na_i, n_i ! quantities for charged nucleation (NOT IN USE)
@@ -61 +60 @@
 !   airn=pplay/t_seri/RD/1.E3*RNAVO/RMD ! molec cm-3 (for future use, to be confirmed)
     CALL newbinapara(t_seri,rh,rhoa,csi,airn,ipr,jnuc_n,ntot_n,jnuc_i,ntot_i, &
-                   & x_n,x_i,na_n,na_i,rc_n,rc_i,n_i,kinetic_n,kinetic_i,rhoatres)
+                   & x_n,x_i,na_n,na_i,rc_n,rc_i,n_i,kinetic_n,kinetic_i)
   ENDIF
 
@@ -296 +295 @@
 SUBROUTINE newbinapara(t,satrat,rhoa,csi,airn,ipr,jnuc_n_real,ntot_n_real,jnuc_i_real,ntot_i_real,        &
                    &   x_n_real,x_i_real,na_n_real,na_i_real,rc_n_real,rc_i_real,n_i_real,                & 
-                   &   kinetic_n,kinetic_i,rhoatres_real)
+                   &   kinetic_n,kinetic_i)
 
   !    Fortran 90 subroutine newbinapara
@@ -355 +354 @@
   REAL,INTENT(OUT) :: rc_i_real     ! radius of the charged critical cluster in nm 
   REAL,INTENT(OUT) :: n_i_real      ! number of ion pairs in air (cm-3) 
-  REAL,INTENT(OUT) :: rhoatres_real ! threshold concentration of H2SO4 (1/cm^3) for neutral kinetic nucleation
   LOGICAL,INTENT(OUT)  :: kinetic_n        ! true if kinetic neutral nucleation
   LOGICAL,INTENT(OUT)  :: kinetic_i        ! true if kinetic ion-induced nucleation
@@ -372 +370 @@
   DOUBLE PRECISION :: rc_i        ! radius of the charged critical cluster in nm 
   DOUBLE PRECISION :: n_i         ! number of ion pairs in air (cm-3) 
-  DOUBLE PRECISION :: rhoatres    ! threshold concentration of H2SO4 (1/cm^3) for neutral kinetic nucleation
   DOUBLE PRECISION :: x           ! mole fraction of H2SO4 in the critical cluster 
   DOUBLE PRECISION :: satratln    ! bounded water saturation ratio for neutral case (between 5.E-6 - 1.0)
@@ -592 +589 @@
   IF (jnuc_n.LT.1.E-7) THEN
      jnuc_n=0.0
-  ENDIF
-  
-  !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-  
-  ! Threshold neutral nucleation rate (j > 1/cm3s) parameterization (can be commented out if not needed)
-  IF (tln .GE. 310.) THEN
-     rhoatres=EXP(-2.8220714121794250 + 1.1492362322651116E+1*satratln -3.3034839106184218E+3/tln &
-          & + (-7.1828571490168133E+2*satratln)/tln + 1.4649510835204091E-1*tln                   &
-          & -3.0442736551916524E-2*satratln*tln -9.3258567137451497E-5*tln**2                     &
-          & -1.1583992506895649E+1*LOG(satratln) + (1.5184848765906165E+3*LOG(satratln))/tln      &
-          & + 1.8144983916747057E-2*tln*LOG(satratln)) !1/cm3
-  ENDIF
-
-  IF (tln .GT. 190. .AND. tln .LT. 310.) THEN
-     rhoatres=EXP(-3.1820396091231999E+2 + 7.2451289153199676*satratln + 2.6729355170089486E+4/tln &
-          & + (-7.1492506076423069E+2*satratln)/tln + 1.2617291148391978*tln                       &
-          & - 1.6438112080468487E-2*satratln*tln -1.4185518234553220E-3*tln**2                     &
-          & -9.2864597847386694*LOG(satratln) + (1.2607421852455602E+3*LOG(satratln))/tln          &
-          & + 1.3324434472218746E-2*tln*LOG(satratln)) !1/cm3
-  ENDIF
-
-  IF (tln .LT. 185. .AND. tln .GT. 155.) THEN
-     rhoatres=1.1788859232398459E+5 - 1.0244255702550814E+4*satratln +  &
-          & 4.6815029684321962E+3*satratln**2 -1.6755952338499657E+2*tln
   ENDIF
   
@@ -832 +805 @@
   rc_i_real     = rc_i
   n_i_real      = n_i
-  rhoatres_real = rhoatres
   
 END SUBROUTINE newbinapara