Changeset 05cf28d in flexpart.git for src/readreleases.f90

Timestamp:

Apr 22, 2016, 9:35:28 AM (8 years ago)

Author:

Espen Sollum ATMOS <eso@…>

Branches:

master, 10.4.1_pesei, GFS_025, bugfixes+enhancements, dev, release-10, release-10.4.1, scaling-bug, univie

Children:

356dcb9

Parents:

61df8d9

Message:

Updated checks and warning messages for wet deposition parameters given in SPECIES files

File:

• src/readreleases.f90 (modified) (1 diff)

src/readreleases.f90

@@ -362 +362 @@
   ! Check if wet deposition or OH reaction shall be calculated
   !***********************************************************
-! NIK 15.02.2015, with new wet dep scheme either A or B parameters can be assigned a positive value which switches on wet dep
-    if (weta_gas(i).gt.0. .or. wetb_gas(i).gt.0.)  then ! :TODO: or particle
+
+  ! ESO 04.2016 check for below-cloud scavenging (gas or aerosol)
+    if ((dquer(i).le.0..and.(weta_gas(i).gt.0. .or. wetb_gas(i).gt.0.)) .or. &
+         &(dquer(i).gt.0. .and. (crain_aero(i) .gt. 0. .or. csnow_aero(i).gt.0.)))  then
       WETDEP=.true.
       if (lroot) then
-        write (*,*) 'Below-cloud scavenging: ON'
-      !  write (*,*) 'Below-cloud scavenging coefficients: ',weta(i),i
+        write (*,*) '  Below-cloud scavenging: ON'
+!  write (*,*) 'Below-cloud scavenging coefficients: ',weta(i),i
       end if
     else
-      if (lroot) write (*,*) 'Below-cloud scavenging: OFF'
+      if (lroot) write (*,*) '  Below-cloud scavenging: OFF'
     endif
     
 ! NIK 31.01.2013 + 10.12.2013 + 15.02.2015
-    if (ccn_aero(i).gt.0. .or. in_aero(i).gt.0.)  then
+    if (dquer(i).gt.0..and.(ccn_aero(i).gt.0. .or. in_aero(i).gt.0.))  then
       WETDEP=.true.
       if (lroot) then
-        write (*,*) 'In-cloud scavenging: ON'
+        write (*,*) '  In-cloud scavenging: ON'
 !        write (*,*) 'In-cloud scavenging coefficients: ',&
 !           &ccn_aero(i),in_aero(i),i !,wetc_in(i), wetd_in(i),i
       end if
     else
-      if (lroot) write (*,*) 'In-cloud scavenging: OFF' 
+      if (lroot) write (*,*) '  In-cloud scavenging: OFF' 
     endif
 
     if (ohcconst(i).gt.0.) then
       OHREA=.true.
-      if (lroot) write (*,*) 'OHreaction switched on: ',ohcconst(i),i
+      if (lroot) write (*,*) '  OHreaction switched on: ',ohcconst(i),i
     endif