Changeset 3605 for LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/micphy_tstep.F90

Timestamp:

Nov 21, 2019, 4:43:45 PM (5 years ago)

Author:

lguez

Message:

Merge revisions 3427:3600 of trunk into branch Ocean_skin

Location:

LMDZ6/branches/Ocean_skin

Files:

• . (modified) (1 prop)
• libf/phylmd/StratAer/micphy_tstep.F90 (modified) (4 diffs, 1 prop)

LMDZ6/branches/Ocean_skin/libf/phylmd/StratAer/micphy_tstep.F90

@@ +1 @@
+!
+! $Id$
+!
 SUBROUTINE micphy_tstep(pdtphys,tr_seri,t_seri,pplay,paprs,rh,is_strato)
 
+  USE geometry_mod, ONLY : latitude_deg !NL- latitude corr. to local domain
   USE dimphy, ONLY : klon,klev
   USE aerophys
@@ -9 +13 @@
   USE sulfate_aer_mod, ONLY : STRAACT
   USE YOMCST, ONLY : RPI, RD, RG
-
+  USE print_control_mod, ONLY: lunout
+  USE strataer_mod
+  
   IMPLICIT NONE
 
@@ -89 +95 @@
       ! compute nucleation rate in kg(H2SO4)/kgA/s
       CALL nucleation_rate(rhoa,t_seri(ilon,ilev),pplay(ilon,ilev),rh(ilon,ilev), &
-             & a_xm,b_xm,c_xm,nucl_rate,ntot,x)
+           & a_xm,b_xm,c_xm,nucl_rate,ntot,x)
+      !NL - add nucleation box (if flag on)
+      IF (flag_nuc_rate_box) THEN
+         IF (latitude_deg(ilon).LE.nuclat_min .OR. latitude_deg(ilon).GE.nuclat_max &
+              .OR. pplay(ilon,ilev).GE.nucpres_max .AND. pplay(ilon,ilev).LE.nucpres_min) THEN
+            nucl_rate=0.0
+         ENDIF
+      ENDIF
       ! compute cond/evap rate in kg(H2SO4)/kgA/s
       CALL condens_evapor_rate(rhoa,t_seri(ilon,ilev),pplay(ilon,ilev), &
@@ -160 +173 @@
     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
+         WRITE(lunout,*) 'micphy_tstep: negative concentration', tr_seri(ilon,ilev,it), ilon, ilev, it
       ENDIF
     ENDDO