Changeset 598e9da in flexpart.git for src

Timestamp:

May 13, 2015, 2:54:59 PM (9 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:

b724ea3

Parents:

5f42c27

Message:

Updated SPECIES files

Location:

src

Files:

• readreleases.f90 (modified) (1 diff)
• readspecies.f90 (modified) (1 diff)
• verttransform.f90 (modified) (1 diff)
• wetdepo.f90 (modified) (5 diffs)

src/readreleases.f90

@@ -376 +376 @@
       if (lroot) then
         write (*,*) 'In-cloud scavenging: ON'
-        write (*,*) 'In-cloud scavenging coefficients: ',&
-           &weta_in(i),wetb_in(i),i !,wetc_in(i), wetd_in(i),i
+!        write (*,*) 'In-cloud scavenging coefficients: ',&
+!           &weta_in(i),wetb_in(i),i !,wetc_in(i), wetd_in(i),i
       end if
     else

src/readspecies.f90

@@ -262 +262 @@
   else !is gas
     if (lroot) then
-      write(*,'(a,f5.2)') '  Gas below-cloud scavenging parameter A  ', weta(pos_spec)
+      write(*,*) '  Gas below-cloud scavenging parameter A  ', weta(pos_spec)
       write(*,'(a,f5.2)') '  Gas below-cloud scavenging parameter B  ', wetb(pos_spec)
       write(*,*) ' Gas in-cloud scavenging uses default values as in Hertel et al 1995'

src/verttransform.f90

@@ -510 +510 @@
 
   if (readclouds) write(*,*) 'using cloud water from ECMWF'
-  if (.not.readclouds) write(*,*) 'using cloud water from parameterization'
+!  if (.not.readclouds) write(*,*) 'using cloud water from parameterization'
 
   rcw(:,:)=0

src/wetdepo.f90

@@ -264 +264 @@
       if (clouds_v.ge.4) then !below cloud
 
-        if (weta(ks).gt.0. .or. wetb(ks).gt.0) then !if positive below-cloud parameters given in SPECIES file (either A or B)
+        if (weta(ks).gt.0. .or. wetb(ks).gt.0.) then !if positive below-cloud parameters given in SPECIES file (either A or B)
           blc_count=blc_count+1
+
 
 !GAS
@@ -279 +280 @@
            
   !ZHG snow or rain removal is applied based on the temperature.
-            if (act_temp .ge. 273)  then !Rain 
+            if (act_temp .ge. 273 .and. weta(ks).gt.0.)  then !Rain 
 
   !Particle RAIN scavenging coefficient based on Laakso et al 2003, the below-cloud scavenging (rain efficienty) parameter A (=weta) from SPECIES file
@@ -285 +286 @@
                    (log10(dquer_m))**(-2))+ (bclr(5)*(log10(dquer_m))**(-1))+ bclr(6)* (prec(1))**(0.5))
 
-            elseif (act_temp .lt. 273)  then !snow
+            elseif (act_temp .lt. 273 .and. wetb(ks).gt.0.)  then !snow
 
   !Particle SNOW scavenging coefficient based on Kyro et al 2009, the below-cloud scavenging (Snow efficiency) parameter B (=wetb) from SPECIES file
@@ -305 +306 @@
       if (clouds_v.lt.4) then !in-cloud
 
-        inc_count=inc_count+1
+! NIK 13 may 2015: only do incloud if positive in-cloud scavenging parameters are given in species file
+        if (weta_in(ks).gt.0. .or. wetb_in(ks).gt.0.) then !if positive in-cloud parameters given in SPECIES file (either Ai or Bi)
+
+! if negative coefficients (turned off) set to zero for use in equation
+          if (weta_in(ks).lt.0.) weta_in(ks)=0.
+          if (wetb_in(ks).lt.0.) wetb_in(ks)=0.
+
+          inc_count=inc_count+1
 
   !ZHG liquid water parameterization (CLWC+CIWC) 
-        if (readclouds) then !get cloud water clwc & ciwc units Kg/Kg
-          cl=clwc(ix,jy,hz,n)+ciwc(ix,jy,hz,n)
-        else !parameterize cloudwater
-          cl=2E-7*prec(1)**0.36
-        endif
-
-  ! AEROSOL
-        if (dquer(ks).gt. 0.) then ! is particle
-          if (act_temp .le. 253) then
-            liq_frac=0
-          else if (act_temp .ge. 273) then
-            liq_frac=1
-          else
-            liq_frac =((act_temp-273)/(273-253))**2
+          if (readclouds) then !get cloud water clwc & ciwc units Kg/Kg
+            cl=clwc(ix,jy,hz,n)+ciwc(ix,jy,hz,n)
+          else !parameterize cloudwater
+            cl=2E-7*prec(1)**0.36
           endif
+
+            if (act_temp .le. 253) then
+              liq_frac=0
+            else if (act_temp .ge. 273) then
+              liq_frac=1
+            else
+              liq_frac =((act_temp-273)/(273-253))**2
+            endif
 
 ! ZHG  calculate the activated fraction based on the In-cloud scavenging parameters Ai (=weta_in) and Bi (=wetb_in) from SPECIES file
@@ -328 +334 @@
 ! The default values are 0.9 for CCN and 0.1 IN
 ! This parameterization is based on Verheggen et al. (2007) & Cozich et al. (2006)
-          frac_act = liq_frac*weta_in(ks) +(1-liq_frac)*wetb_in(ks)
+            frac_act = liq_frac*weta_in(ks) +(1-liq_frac)*wetb_in(ks)
  
   !ZHG Use the activated fraction and the liqid water to calculate the washout
-          S_i= frac_act/cl
+
+  ! AEROSOL
+          if (dquer(ks).gt. 0.) then ! is particle
+
+            S_i= frac_act/cl
 
   ! GAS
-        else ! is gas
-          cle=(1-cl)/(henry(ks)*(r_air/3500.)*act_temp)+cl
-          S_i=1/cle
-        endif
+          else ! is gas
+               
+            cle=(1-cl)/(henry(ks)*(r_air/3500.)*act_temp)+cl
+            S_i=frac_act/cle
+
+          endif ! gas or particle
 
   ! scavenging coefficient based on Hertel et al 1995 - using the S_i for either gas or aerosol
-        wetscav=S_i*(prec(1)/3.6E6)/clouds_h
+          wetscav=S_i*(prec(1)/3.6E6)/clouds_h
 
 !          write(*,*) 'in-cloud, act_temp=',act_temp,',prec=',prec(1),',wetscav=',wetscav,',jpart=',jpart,',clouds_h=,', &
 !          clouds_h,',cl=',cl, 'diff to old scheme=', cl-2E-7*prec(1)**0.36
-
+        endif ! positive in-cloud scavenging parameters given in Species file
       endif !incloud