Commit 3ea93bb3 authored by Ignacio Pisso's avatar Ignacio Pisso

Merge branch 'release-10' into dev

parents 95a8cb6c 3d7eebf7
This diff is collapsed.
This diff is collapsed.
......@@ -20,12 +20,12 @@
IOUT= 1, ! Output type: [1]mass 2]pptv 3]1&2 4]plume 5]1&4, +8 for NetCDF output
IPOUT= 0, ! Particle position output: 0]no 1]every output 2]only at end 3]time averaged
LSUBGRID= 0, ! Increase of ABL heights due to sub-grid scale orographic variations;[0]off 1]on
LCONVECTION= 1, ! Switch for convection parameterization;0]off [1]on
LAGESPECTRA= 0, ! Switch for calculation of age spectra (needs AGECLASSES);[0]off 1]on
LCONVECTION= 1, ! Switch for convection parameterization;0]off [1]on
LAGESPECTRA= 0, ! Switch for calculation of age spectra (needs AGECLASSES);[0]off 1]on
IPIN= 0, ! Warm start from particle dump (needs previous partposit_end file); [0]no 1]yes
IOUTPUTFOREACHRELEASE= 1, ! Separate output fields for each location in the RELEASE file; [0]no 1]yes
IFLUX= 0, ! Output of mass fluxes through output grid box boundaries
MDOMAINFILL= 0, ! Switch for domain-filling, if limited-area particles generated at boundary
MDOMAINFILL= 0, ! Switch for domain-filling, if limited-area particles generated at boundary
IND_SOURCE= 1, ! Unit to be used at the source ; [1]mass 2]mass mixing ratio
IND_RECEPTOR= 1, ! Unit to be used at the receptor; [1]mass 2]mass mixing ratio 3]wet depo. 4]dry depo.
MQUASILAG= 0, ! Quasi-Lagrangian mode to track individual numbered particles
......
......@@ -18,44 +18,44 @@ WET DEPOSITION
Below-cloud scavenging
- Gases:
A and B are the parameters for below-cloud scavenging defined by Eq. scavenging coefficient = A x precipitation-rate^B
A and B are the parameters for below-cloud scavenging defined by Eq. scavenging coefficient = A x precipitation-rate^B
If you migrate from FLEXPART version 9.2 and below this is equal to A and B as it was.
- Particles:
A "rain collection efficiency"
B "snow collection efficiency" are both defined by Equations 46-47.
Suggested values for particle scavenging are A=1.0,B=1.0.
Using values A=0.2 and B=0.2 will give results close to the old scavenging scheme of FLEXPART version 9.2 and below.
A "rain collection efficiency"
B "snow collection efficiency" are both defined by Equations 46-47.
Suggested values for particle scavenging are A=1.0,B=1.0.
Using values A=0.2 and B=0.2 will give results close to the old scavenging scheme of FLEXPART version 9.2 and below.
In-cloud scavenging
- Gases:
Is switched on if Ai or Bi are given greater than 0
If you migrate from FLEXPART version 9.2 and below this is equal to Ai and Bi should be set to 1
- Particles:
Ai ("Cloud Condensation Nuclei (CCN) efficiency")
Bi ("Ice Nuclei (IN) efficiency") are both from Eq. 51.
Suggested values for particle scavenging are Ai=0.9 and Bi=0.1.
Using values Ai=0.9 and Bi=0.9 should give results cto the old scavenging scheme of FLEXPART version 9.2 and below.
Ai ("Cloud Condensation Nuclei (CCN) efficiency")
Bi ("Ice Nuclei (IN) efficiency") are both from Eq. 51.
Suggested values for particle scavenging are Ai=0.9 and Bi=0.1.
Using values Ai=0.9 and Bi=0.9 should give results cto the old scavenging scheme of FLEXPART version 9.2 and below.
---------------------------------------------------------------------------------------------------------------------------
DRY DEPOSITION
- Gases:
D is the diffusivity of water vapor and
D_i is the diffusivity of the species,
H is the effective Henry’s constant,
f0 varies between 0 and 1 and gives the reactivity of a species relative to that of ozone. For non-reactive species f0 is 0, for slightly reactive it is 0.1 and for highly reactive it is 1.
Dry deposition of gases is switched off by negative D.
D is the diffusivity of water vapor and
D_i is the diffusivity of the species,
H is the effective Henry’s constant,
f0 varies between 0 and 1 and gives the reactivity of a species relative to that of ozone. For non-reactive species f0 is 0, for slightly reactive it is 0.1 and for highly reactive it is 1.
Dry deposition of gases is switched off by negative D.
- Particles:
rho specifies the density of the substance,
dquer its mean diameter
dsig the measure of variation.
If no detailed information for deposition velocity calculation is available, a constant deposition velocity vd [cm s-1] can be used.
Dry deposition of particles is switched off by negative rho.
rho specifies the density of the substance,
dquer its mean diameter
dsig the measure of variation.
If no detailed information for deposition velocity calculation is available, a constant deposition velocity vd [cm s-1] can be used.
Dry deposition of particles is switched off by negative rho.
---------------------------------------------------------------------------------------------------------------------------
OH REACTION:
can be turned on by giving parameter C [cm^3/molecule/sec]], D [K] and N [no unit] positive values, defined by Eq. 39.
OH reaction is switched off by negative C.
can be turned on by giving parameter C [cm^3/molecule/sec]], D [K] and N [no unit] positive values, defined by Eq. 39.
OH reaction is switched off by negative C.
---------------------------------------------------------------------------------------------------------------------------
Optionally emission variation information
can be added at the end of the file, defined as following: Since FLEXPART version 6.0, emission factors can be defined that change the temporal variation of particle releases. This is useful, for instance,
can be added at the end of the file, defined as following: Since FLEXPART version 6.0, emission factors can be defined that change the temporal variation of particle releases. This is useful, for instance,
to simulate the typical daily and weekly cycle of anthropogenic emissions. The emission factors are given in the file of the
corresponding species SPECIES_nnn, where nnn is the species number defined in file RELEASES. If no emission variation
information is given, emission rates for species nnn are taken as constant. Release rates can vary with the hour of the day
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
program flexpart
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
program flexpart
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine advance(itime,nrelpoint,ldt,up,vp,wp, &
usigold,vsigold,wsigold,nstop,xt,yt,zt,prob,icbt)
! i i i/oi/oi/o
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine assignland
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine boundcond_domainfill(itime,loutend)
! i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine boundcond_domainfill(itime,loutend)
! i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine calcfluxes(nage,jpart,xold,yold,zold)
! i i i i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine calcmatrix(lconv,delt,cbmf,metdata_format)
! o i o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine calcpar(n,uuh,vvh,pvh,metdata_format)
! i i i o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine calcpar_nests(n,uuhn,vvhn,pvhn,metdata_format)
! i i i o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine calcpv(n,uuh,vvh,pvh)
! i i i o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine calcpv_nests(l,n,uuhn,vvhn,pvhn)
! i i i i o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine caldate(juldate,yyyymmdd,hhmiss)
! i o o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine cbl(wp,zp,ust,wst,h,rhoa,rhograd,sigmaw,dsigmawdz,tlw,ptot,Q,phi,ath,bth,ol,flagrein)
! i i i i i i i i i i o o o o o i o
......@@ -105,7 +108,7 @@ subroutine cbl(wp,zp,ust,wst,h,rhoa,rhograd,sigmaw,dsigmawdz,tlw,ptot,Q,phi,ath,
radw2=w2**0.5
dradw2=0.5*w2**(-0.5)*dw2
!costluar4=0.66667 ! costante da LHH
fluarw=costluar4*(cuberoot(skew)) !skew**(1./3.)
fluarw=costluar4*(cuberoot(skew)) !skew**(1./3.)
fluarw2=fluarw*fluarw
if (skew.ne.0) then
......@@ -183,7 +186,7 @@ subroutine cbl(wp,zp,ust,wst,h,rhoa,rhograd,sigmaw,dsigmawdz,tlw,ptot,Q,phi,ath,
pa=(usurad2p*(1./sigmawa))*(exp(-(0.5*((deltawa/sigmawa)**2.))))
pb=(usurad2p*(1./sigmawb))*(exp(-(0.5*((deltawb/sigmawb)**2.))))
ptot=dens*aluarw*pa+dens*bluarw*pb
aperfa=deltawa*usurad2/sigmawa
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine centerofmass(xl,yl,n,xcenter,ycenter)
! i i i o o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine clustering(xl,yl,zl,n,xclust,yclust,zclust,fclust,rms, &
rmsclust,zrms)
! i i i i o o o o o
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
! Changes to the routines by A. Stohl
! xi,xi0,eta,eta0 are double precision variables to avoid problems
! at poles
......@@ -228,7 +231,7 @@ real function cgszll (strcmp, xlat,xlong)
ymerc = log((1. + slat) / (1. - slat))/2.
!efact = exp(ymerc)
!cgszll = 2. * strcmp(7) * exp (strcmp(1) * ymerc)
!c / (efact + 1./efact)
!c / (efact + 1./efact)
endif
cgszll = strcmp(7) * cos(radpdg * xlat) * exp(strcmp(1) *ymerc)
return
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine conccalc(itime,weight)
! i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine conccalc(itime,weight)
! i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine concoutput(itime,outnum,gridtotalunc,wetgridtotalunc, &
drygridtotalunc)
! i i o o
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine concoutput_inversion(itime,outnum,gridtotalunc,wetgridtotalunc, &
drygridtotalunc)
! i i o o
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine concoutput_inversion_nest(itime,outnum)
! i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine concoutput(itime,outnum,gridtotalunc,wetgridtotalunc, &
drygridtotalunc)
! i i o o
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine concoutput_nest(itime,outnum)
! i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine concoutput_nest(itime,outnum)
! i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine concoutput_surf(itime,outnum,gridtotalunc,wetgridtotalunc, &
drygridtotalunc)
! i i o o
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine concoutput_surf(itime,outnum,gridtotalunc,wetgridtotalunc, &
drygridtotalunc)
! i i o o
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine concoutput_surf_nest(itime,outnum)
! i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine concoutput_surf_nest(itime,outnum)
! i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
!**************************************************************************
!**** SUBROUTINE CONVECT *****
!**** VERSION 4.3c *****
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine convmix(itime,metdata_format)
! i
!**************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine coordtrafo
!**********************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
integer function detectformat()
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
!-----------------------------------------------------------------------
function distance(rlat1,rlon1,rlat2,rlon2)
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
!-----------------------------------------------------------------------
function distance2(rlat1,rlon1,rlat2,rlon2)
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine drydepokernel(nunc,deposit,x,y,nage,kp)
! i i i i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine drydepokernel_nest(nunc,deposit,x,y,nage,kp)
! i i i i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
! Function calculates dynamic viscosity of air (kg/m/s) as function of
! temperature (K) using Sutherland's formula
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
! To be used, if the non-standard Fortran function erf does not exist on
! your machine
!
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
real function ew(x)
!****************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
module flux_mod
! flux eastward, westward, northward, southward, upward and downward
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine fluxoutput(itime)
! i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine get_settling(itime,xt,yt,zt,nsp,settling)
! i i i i i o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine get_vdep_prob(itime,xt,yt,zt,prob)
! i i i i o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine get_wetscav(itime,ltsample,loutnext,jpart,ks,grfraction,inc_count,blc_count,wetscav)
! i i i i i o o o o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine getfields(itime,nstop,metdata_format)
! i o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine getfields(itime,nstop,memstat,metdata_format)
! i o o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine gethourlyOH(itime)
! i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine getrb(nc,ustar,nyl,diffh2o,reldiff,rb)
! i i i i i o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine getrc(nc,i,j,t,gr,rh,rr,rc)
! i i i i i i i o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine getvdep(n,ix,jy,ust,temp,pa,L,gr,rh,rr,snow,vdepo)
! i i i i i i i i i i i o
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine getvdep_nests(n,ix,jy,ust,temp,pa, &
L,gr,rh,rr,snow,vdepo,lnest)
! i i i i i i i i i i i o i
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine gridcheck_ecmwf
!**********************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine gridcheck_gfs
!**********************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine gridcheck_nests
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine hanna(z)
! i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine hanna1(z)
! i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine hanna_short(z)
! i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine init_domainfill
!
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine init_domainfill
!
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine initial_cond_calc(itime,i)
! i i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine initial_cond_output(itime)
! i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine initial_cond_output_inversion(itime)
! i
!*****************************************************************************
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine initialize(itime,ldt,up,vp,wp, &
usigold,vsigold,wsigold,xt,yt,zt,icbt)
! i i o o o
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt
! SPDX-License-Identifier: GPL-3.0-or-later
subroutine initialize_cbl_vel(idum,zp,ust,wst,h,sigmaw,wp, ol)
! i/o i i i i i o i
......
! SPDX-FileCopyrightText: FLEXPART 1998-2019, see flexpart_license.txt