Changeset 341f4b7 in flexpart.git for src/readspecies.f90

Timestamp:

Apr 20, 2016, 1:39:51 PM (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:

61df8d9

Parents:

32b49c3

Message:

SPECIES files converted to namelist format. Fixed-format SPCIES files renamed with .oldformat extension. Added 2 wet depo parameters to old-format SPECIES. Renamed internal variables and parameters used for wet deposition. incloud_ratio increased to 2.2

File:

• src/readspecies.f90 (modified) (13 diffs)

src/readspecies.f90

@@ -42 +42 @@
   ! decaytime(maxtable)  half time for radiological decay                      *
   ! specname(maxtable)   names of chemical species, radionuclides              *
-  ! weta, wetb           Parameters for determining below-cloud scavenging     *
-  ! weta_in              Parameter for determining in-cloud scavenging         *
-  ! wetb_in              Parameter for determining in-cloud scavenging         *
+  ! weta_gas, wetb_gas   Parameters for determining below-cloud scavenging     *
+  ! ccn_aero              Parameter for determining in-cloud scavenging         *
+  ! in_aero              Parameter for determining in-cloud scavenging         *
   ! ohcconst             OH reaction rate constant C                           *
   ! ohdconst             OH reaction rate constant D                           *
@@ -66 +66 @@
 
   character(len=16) :: pspecies
-  real :: pdecay, pweta, pwetb, preldiff, phenry, pf0, pdensity, pdquer
+  real :: pdecay, pweta_gas, pwetb_gas, preldiff, phenry, pf0, pdensity, pdquer
   real :: pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst, pspec_ass, pkao
-  real :: pweta_in, pwetb_in, pwetc_in, pwetd_in
+  real :: pcrain_aero, pcsnow_aero, pccn_aero, pin_aero
   integer :: readerror
 
 ! declare namelist
   namelist /species_params/ &
-       pspecies, pdecay, pweta, pwetb, &
-       pweta_in, pwetb_in, pwetc_in, pwetd_in, &
+       pspecies, pdecay, pweta_gas, pwetb_gas, &
+       pcrain_aero, pcsnow_aero, pccn_aero, pin_aero, &
        preldiff, phenry, pf0, pdensity, pdquer, &
        pdsigma, pdryvel, pweightmolar, pohcconst, pohdconst, pohnconst, pspec_ass, pkao
 
-  pspecies=" "
+  pspecies="" ! read failure indicator value
   pdecay=-999.9
-  pweta=-9.9E-09
-  pwetb=0.0
-  pweta_in=-9.9E-09
-  pwetb_in=-9.9E-09
-!  pwetc_in=-9.9E-09
-!  pwetd_in=-9.9E-09
+  pweta_gas=-9.9E-09
+  pwetb_gas=0.0
+  pcrain_aero=-9.9E-09
+  pcsnow_aero=-9.9E-09
+  pccn_aero=-9.9E-09
+  pin_aero=-9.9E-09
   preldiff=-9.9
   phenry=0.0
@@ -98 +98 @@
   pspec_ass=-9
   pkao=-99.99
-  pweightmolar=-789.0 ! read failure indicator value
+  pweightmolar=-999.9
 
 ! Open the SPECIES file and read species names and properties
@@ -113 +113 @@
   close(unitspecies)
 
-  if ((pweightmolar.eq.-789.0).or.(readerror.ne.0)) then ! no namelist found
+!  if ((pweightmolar.eq.-789.0).or.(readerror.ne.0)) then ! no namelist found
+  if ((len(pspecies).eq.0).or.(readerror.ne.0)) then ! no namelist found
 
     readerror=1
@@ -127 +128 @@
     read(unitspecies,'(f18.1)',end=22) decay(pos_spec)
 !  write(*,*) decay(pos_spec)
-    read(unitspecies,'(e18.1)',end=22) weta(pos_spec)
-!  write(*,*) weta(pos_spec)
-    read(unitspecies,'(f18.2)',end=22) wetb(pos_spec)
-!  write(*,*) wetb(pos_spec)
-
+    read(unitspecies,'(e18.1)',end=22) weta_gas(pos_spec)
+!  write(*,*) weta_gas(pos_spec)
+    read(unitspecies,'(f18.2)',end=22) wetb_gas(pos_spec)
+!  write(*,*) wetb_gas(pos_spec)
+    read(unitspecies,'(e18.1)',end=22) crain_aero(pos_spec)
+!  write(*,*) crain_aero(pos_spec)
+    read(unitspecies,'(f18.2)',end=22) csnow_aero(pos_spec)
+!  write(*,*) csnow_aero(pos_spec)
 !*** NIK 31.01.2013: including in-cloud scavening parameters
-    read(unitspecies,'(e18.1)',end=22) weta_in(pos_spec)
-!  write(*,*) weta_in(pos_spec)
-    read(unitspecies,'(f18.2)',end=22) wetb_in(pos_spec)
-!  write(*,*) wetb_in(pos_spec)
-! read(unitspecies,'(f18.2)',end=22) wetc_in(pos_spec)
-!  write(*,*) wetc_in(pos_spec)
-! read(unitspecies,'(f18.2)',end=22) wetd_in(pos_spec)
-!  write(*,*) wetd_in(pos_spec)
-
+    read(unitspecies,'(e18.1)',end=22) ccn_aero(pos_spec)
+!  write(*,*) ccn_aero(pos_spec)
+    read(unitspecies,'(f18.2)',end=22) in_aero(pos_spec)
+!  write(*,*) in_aero(pos_spec)
     read(unitspecies,'(f18.1)',end=22) reldiff(pos_spec)
 !  write(*,*) reldiff(pos_spec)
@@ -171 +170 @@
     pspecies=species(pos_spec)
     pdecay=decay(pos_spec)
-    pweta=weta(pos_spec)
-    pwetb=wetb(pos_spec)
-    pweta_in=weta_in(pos_spec)
-    pwetb_in=wetb_in(pos_spec)
-!    pwetc_in=wetc_in(pos_spec)
-!    pwetd_in=wetd_in(pos_spec)
+    pweta_gas=weta_gas(pos_spec)
+    pwetb_gas=wetb_gas(pos_spec)
+    pcrain_aero=crain_aero(pos_spec)
+    pcsnow_aero=csnow_aero(pos_spec)
+    pccn_aero=ccn_aero(pos_spec)
+    pin_aero=in_aero(pos_spec)
     preldiff=reldiff(pos_spec)
     phenry=henry(pos_spec)
@@ -195 +194 @@
     species(pos_spec)=pspecies
     decay(pos_spec)=pdecay
-    weta(pos_spec)=pweta
-    wetb(pos_spec)=pwetb
-    weta_in(pos_spec)=pweta_in
-    wetb_in(pos_spec)=pwetb_in
-!    wetc_in(pos_spec)=pwetc_in
-!    wetd_in(pos_spec)=pwetd_in
+    weta_gas(pos_spec)=pweta_gas
+    wetb_gas(pos_spec)=pwetb_gas
+    crain_aero=pcrain_aero
+    csnow_aero=pcsnow_aero
+    ccn_aero(pos_spec)=pccn_aero
+    in_aero(pos_spec)=pin_aero
     reldiff(pos_spec)=preldiff
     henry(pos_spec)=phenry
@@ -222 +221 @@
 ! Check scavenging parameters given in SPECIES file
 
-  if (weta(pos_spec).gt.0.0 .or. wetb(pos_spec).gt.0.0 .or. weta_in(pos_spec).gt.0.0 .or. wetb_in(pos_spec).gt.0.0) then
+  if (weta_gas(pos_spec).gt.0.0 .or. wetb_gas(pos_spec).gt.0.0 .or. ccn_aero(pos_spec).gt.0.0 .or. in_aero(pos_spec).gt.0.0) then
 
   if (dquer(pos_spec).gt.0) then !is particle
     if (lroot) then
       write(*,'(a,f5.2)') '  Particle below-cloud scavenging parameter A &
-           &(Rain collection efficiency)  ', weta(pos_spec)
+           &(Rain collection efficiency)  ', crain_aero(pos_spec)
       write(*,'(a,f5.2)') '  Particle below-cloud scavenging parameter B &
-           &(Snow collection efficiency)  ', wetb(pos_spec)
-    end if
-    if (weta(pos_spec).gt.1.0 .or. wetb(pos_spec).gt.1.0) then
+           &(Snow collection efficiency)  ', csnow_aero(pos_spec)
+    end if
+    if (weta_gas(pos_spec).gt.1.0 .or. wetb_gas(pos_spec).gt.1.0) then
       if (lroot) then
         write(*,*) '*******************************************'
@@ -242 +241 @@
     if (lroot) then
       write(*,'(a,f5.2)') '  Particle in-cloud scavenging parameter Ai &
-           &(CCN efficiency)  ', weta_in(pos_spec)
+           &(CCN efficiency)  ', ccn_aero(pos_spec)
       write(*,'(a,f5.2)') '  Particle in-cloud scavenging parameter Bi &
-           &(IN efficiency)  ', wetb_in(pos_spec)
-    end if
-    if (weta_in(pos_spec).gt.1.0 .or. weta_in(pos_spec).lt.0.01) then
+           &(IN efficiency)  ', in_aero(pos_spec)
+    end if
+    if (ccn_aero(pos_spec).gt.1.0 .or. ccn_aero(pos_spec).lt.0.01) then
       if (lroot) then
         write(*,*) '*******************************************'
@@ -254 +253 @@
       end if
     endif
-    if (wetb_in(pos_spec).gt.1.0 .or. wetb_in(pos_spec).lt.0.01) then
+    if (in_aero(pos_spec).gt.1.0 .or. in_aero(pos_spec).lt.0.01) then
       if (lroot) then
         write(*,*) '*******************************************'
@@ -265 +264 @@
   else !is gas
     if (lroot) then
-      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 below-cloud scavenging parameter A  ', weta_gas(pos_spec)
+      write(*,'(a,f5.2)') '  Gas below-cloud scavenging parameter B  ', wetb_gas(pos_spec)
       write(*,*) ' Gas in-cloud scavenging uses default values as in Hertel et al 1995'
     end if
-    if (weta(pos_spec).gt.0.) then !if wet deposition is turned on
-      if (weta(pos_spec).gt.1E-04 .or. weta(pos_spec).lt.1E-09) then
+    if (weta_gas(pos_spec).gt.0.) then !if wet deposition is turned on
+      if (weta_gas(pos_spec).gt.1E-04 .or. weta_gas(pos_spec).lt.1E-09) then
         if (lroot) then
           write(*,*) '*******************************************'
@@ -279 +278 @@
       endif
     end if
-    if (wetb(pos_spec).gt.0.) then !if wet deposition is turned on
-      if (wetb(pos_spec).gt.0.8 .or. wetb(pos_spec).lt.0.6) then
+    if (wetb_gas(pos_spec).gt.0.) then !if wet deposition is turned on
+      if (wetb_gas(pos_spec).gt.0.8 .or. wetb_gas(pos_spec).lt.0.6) then
         if (lroot) then
           write(*,*) '*******************************************'
@@ -292 +291 @@
   endif
 
-  if (((weta(pos_spec).gt.0).or.(wetb(pos_spec).gt.0)).and.(henry(pos_spec).le.0)) then
+  if (((weta_gas(pos_spec).gt.0).or.(wetb_gas(pos_spec).gt.0)).and.(henry(pos_spec).le.0)) then
     if (dquer(pos_spec).le.0) goto 996 ! no particle, no henry set
   endif