Commit 44746a80 authored by Espen Sollum's avatar Espen Sollum
Browse files

CTBTO version (2nd NILU delivery)

parent 876598a7
# ibdate,ibtime, iedate, ietime, flexversion
20130625 0 20130625 20000 Version 10.2beta (2017-08-01)
20130625 0 20130625 20000 Version 10.1beta (2016-11-02)
# interval, averaging time, sampling time
-1800 -1800 -900
# information on grid setup
......
# ibdate,ibtime, iedate, ietime, flexversion
20130625 0 20130625 20000 Ver. 10.2beta MPI (2017-08-01)
20130625 0 20130625 20000 Ver. 10.1beta MPI (2016-11-02)
# interval, averaging time, sampling time
-1800 -1800 -900
# information on grid setup
......
......@@ -39,6 +39,12 @@ module par_mod
implicit none
! Set to T to disable use of kernel for concentrations/deposition
!****************************************************************
logical, parameter :: lusekerneloutput=.false.
!****************************************************************
! Parameters defining KIND parameter for double/single precision
!****************************************************************
......@@ -53,17 +59,6 @@ module par_mod
integer,parameter :: dep_prec=sp
!****************************************************************
! Set to T to disable use of kernel for concentrations/deposition
!****************************************************************
logical, parameter :: lusekerneloutput=.true.
!*********************************************************************
! Set to T to change output units to number of particles per grid cell
!*********************************************************************
LOGICAL, parameter :: lparticlecountoutput=.false.
!***********************************************************
! Number of directories/files used for FLEXPART input/output
!***********************************************************
......@@ -73,12 +68,20 @@ module par_mod
! numpath Number of different pathnames for input/output files
!*********************************************************************
! Set to T to change output units to number of particles per grid cell
!*********************************************************************
logical, parameter :: lparticlecountoutput=.false.
!*****************************
! Physical and other constants
!*****************************
real,parameter :: pi=3.14159265, r_earth=6.371e6, r_air=287.05, ga=9.81
real,parameter :: cpa=1004.6, kappa=0.286, pi180=pi/180., vonkarman=0.4
REAL,PARAMETER :: cpa=1004.6, kappa=0.286, pi180=pi/180., vonkarman=0.4
! additional constants RLT Aug-2017
REAL,PARAMETER :: rgas=8.31447
REAL,PARAMETER :: r_water=461.495
! pi number "pi"
! pi180 pi/180.
......@@ -141,25 +144,29 @@ module par_mod
! switchsouth use polar stereographic grid south of switchsouth
!*********************************************
!*********************************************
! Maximum dimensions of the input mother grids
!*********************************************
! integer,parameter :: nxmax=361,nymax=181,nuvzmax=92,nwzmax=92,nzmax=92 !ECMWF new
integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138 !ECMWF new
! integer,parameter :: nxmax=181,nymax=91,nuvzmax=138,nwzmax=138,nzmax=138 !ECMWF new
! INTEGER,PARAMETER :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138 !NCEP data
!integer,parameter :: nxmax=361,nymax=181,nuvzmax=92,nwzmax=92,nzmax=92 !ECMWF new
!integer,parameter :: nxmax=15,nymax=15,nuvzmax=140,nwzmax=140,nzmax=140 !ECMWF new
! integer,parameter :: nxmax=361,nymax=181,nuvzmax=140,nwzmax=140,nzmax=140 !NCEP data
! integer,parameter :: nxshift=359 ! for ECMWF
integer,parameter :: nxshift=0 ! for GFS
!*********************************************
! integer,parameter :: nxshift=0 ! for GFS
! integer,parameter :: nxmax=15,nymax=15,nuvzmax=140,nwzmax=140,nzmax=140
! integer,parameter :: nxmax=361,nymax=181,nuvzmax=140,nwzmax=140,nzmax=140
! integer,parameter :: nxmax=721,nymax=361,nuvzmax=140,nwzmax=140,nzmax=140
!integer,parameter :: nxshift=359 ! for ECMWF
!integer,parameter :: nxshift=0 ! for GFS
integer,parameter :: nxmax=361,nymax=181,nuvzmax=140,nwzmax=140,nzmax=140
integer,parameter :: nxshift=0
!*********************************************
! Maximum dimensions of the nested input grids
!*********************************************
integer,parameter :: maxnests=1,nxmaxn=451,nymaxn=226
integer,parameter :: maxnests=0,nxmaxn=451,nymaxn=226
! integer,parameter :: maxnests=0
! nxmax,nymax maximum dimension of wind fields in x and y
! direction, respectively
......@@ -171,10 +178,9 @@ module par_mod
! nxshift grid points, in order to accomodate nested
! grids, and output grids overlapping the domain "boundary"
! nxshift must not be negative; "normal" setting would be 0
integer,parameter :: nconvlevmax = nuvzmax-1
integer,parameter :: na = nconvlevmax+1
integer,parameter :: nconvlevmax = nuvzmax-1
integer,parameter :: na = nconvlevmax+1
! ntracermax maximum number of tracer species in convection
! nconvlevmax maximum number of levels for convection
......@@ -185,7 +191,7 @@ module par_mod
! Parmaters for GRIB file decoding
!*********************************
integer,parameter :: jpack=4*nxmax*nymax, jpunp=4*jpack
integer,parameter :: jpack=4*nxmax*nymax, jpunp=4*jpack
! jpack,jpunp maximum dimensions needed for GRIB file decoding
......@@ -217,9 +223,8 @@ module par_mod
! Maximum number of particles, species, and similar
!**************************************************
integer,parameter :: maxpart=10000000
integer,parameter :: maxspec=4
integer,parameter :: maxpart=2000000
integer,parameter :: maxspec=2
real,parameter :: minmass=0.0001
! maxpart Maximum number of particles
......@@ -261,7 +266,7 @@ module par_mod
! Dimension of random number field
!*********************************
integer,parameter :: maxrand=200000000
integer,parameter :: maxrand=2000000
! maxrand number of random numbers used
......@@ -286,7 +291,9 @@ module par_mod
integer,parameter :: unitOH=1
integer,parameter :: unitdates=94, unitheader=90,unitheader_txt=100, unitshortpart=95, unitprecip=101
integer,parameter :: unitboundcond=89
integer,parameter :: unittmp=101
INTEGER,PARAMETER :: unittmp=101
! RLT
INTEGER,PARAMETER :: unitoutfactor=102
!******************************************************
! integer code for missing values, used in wet scavenging (PS, 2012)
......@@ -294,5 +301,8 @@ module par_mod
integer,parameter :: icmv=-9999
! Parameters for testing
!*******************************************
! integer :: verbosity=0
end module par_mod
+++++++++++++ HEADER +++++++++++++++++
+++++++++++++ HEADER +++++++++++++++++
+++++++++++++ HEADER +++++++++++++++++
+++++++++++++ HEADER +++++++++++++++++
+++++++++++++ HEADER +++++++++++++++++
+++++++++++++ HEADER +++++++++++++++++
+++++++++++++ HEADER +++++++++++++++++
-1
20130625 000000
20130625 020000
1800
1800
900
9999999
900 SYNC
-5.0 CTL
4 IFINE
1 IOUT
0 IPOUT
1 LSUBGRID
0 LCONVECTION
0 LAGESPECTRA
0 IPIN
1 IOFR
0 IFLUX
0 MDOMAINFILL
1 IND_SOURCE
2 IND_RECEPTOR
0 MQUASILAG
0 NESTED_OUTPUT
2 LINIT_COND INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT
0 SURF_ONLY
&COMMAND
LDIRECT= -1,
IBDATE= 20130625,
IBTIME= 000000,
IEDATE= 20130625,
IETIME= 020000,
LOUTSTEP= 1800,
LOUTAVER= 1800,
LOUTSAMPLE= 900,
ITSPLIT= 99999999,
LSYNCTIME= 900,
CTL= -5.0000000 ,
IFINE= 4,
IOUT= 1,
IPOUT= 0,
LSUBGRID= 1,
LCONVECTION= 0,
LAGESPECTRA= 0,
IPIN= 0,
IOUTPUTFOREACHRELEASE= 1,
IFLUX= 0,
MDOMAINFILL= 0,
IND_SOURCE= 1,
IND_RECEPTOR= 2,
MQUASILAG= 0,
NESTED_OUTPUT= 0,
LINIT_COND= 2,
SURF_ONLY= 0,
CBLFLAG= 0,
OHFIELDS_PATH= "../../flexin/",
/
......@@ -27,11 +27,7 @@
! *
! 1997 *
! *
! Update 15 August 2013 IP *
! *
! ESO 2016: *
! GFS specific parameters moved to gfs_mod.f90 *
! ECMWF specific parameters moved to ecmwf_mod.f90 *
! Last update 10 August 2000 *
! *
!*******************************************************************************
......@@ -39,25 +35,18 @@ module par_mod
implicit none
! Set to T to disable use of kernel for concentrations/deposition
!****************************************************************
! Parameters defining KIND parameter for double/single precision
!****************************************************************
integer,parameter :: dp=selected_real_kind(P=15)
integer,parameter :: sp=selected_real_kind(6)
!****************************************************************
! dep_prec sets the precision for deposition calculations (sp or
! dp). sp is default, dp can be used for increased precision.
!****************************************************************
logical, parameter :: lnokernel=.false.
integer,parameter :: dep_prec=sp
!****************************************************************
! Set to T to disable use of kernel for concentrations/deposition
! Parameter defining KIND parameter for "double precision"
!****************************************************************
logical, parameter :: lnokernel=.false.
integer,parameter :: dp=selected_real_kind(P=15)
!***********************************************************
! Number of directories/files used for FLEXPART input/output
......@@ -87,9 +76,6 @@ module par_mod
real,parameter :: karman=0.40, href=15., convke=2.0
real,parameter :: hmixmin=100., hmixmax=4500., turbmesoscale=0.16
real,parameter :: d_trop=50., d_strat=0.1
real,parameter :: rho_water=1000. !ZHG 2015 [kg/m3]
!ZHG MAR2016
real,parameter :: incloud_ratio=6.2
! karman Karman's constant
! href [m] Reference height for dry deposition
......@@ -139,25 +125,18 @@ module par_mod
!*********************************************
! Maximum dimensions of the input mother grids
!*********************************************
! integer,parameter :: nxmax=361,nymax=181,nuvzmax=92,nwzmax=92,nzmax=92 !ECMWF new
! integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138 !ECMWF new
integer,parameter :: nxmax=181,nymax=91,nuvzmax=138,nwzmax=138,nzmax=138 !ECMWF new
! INTEGER,PARAMETER :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138 !NCEP data
!integer,parameter :: nxmax=361,nymax=181,nuvzmax=92,nwzmax=92,nzmax=92
!integer,parameter :: nxmax=361,nymax=181,nuvzmax=61,nwzmax=61,nzmax=61
!integer,parameter :: nxmax=721,nymax=361,nuvzmax=64,nwzmax=64,nzmax=64
!integer,parameter :: nxshift=359 ! for ECMWF
!integer,parameter :: nxshift=0 ! for GFS
integer,parameter :: nxmax=721,nymax=361,nuvzmax=138,nwzmax=138,nzmax=138
integer,parameter :: nxshift=0
! !integer,parameter :: nxshift=359 ! for ECMWF
! integer,parameter :: nxshift=0 ! for GFS
! integer,parameter :: nxmax=15,nymax=15,nuvzmax=140,nwzmax=140,nzmax=140
! integer,parameter :: nxshift=359 ! for ECMWF
integer,parameter :: nxshift=0 ! for GFS
!*********************************************
! Maximum dimensions of the nested input grids
!*********************************************
integer,parameter :: nconvlevmax = nuvzmax-1
integer,parameter :: na = nconvlevmax+1
integer,parameter :: maxnests=0,nxmaxn=451,nymaxn=226
! nxmax,nymax maximum dimension of wind fields in x and y
! direction, respectively
......@@ -169,16 +148,23 @@ module par_mod
! nxshift grid points, in order to accomodate nested
! grids, and output grids overlapping the domain "boundary"
! nxshift must not be negative; "normal" setting would be 0
integer,parameter :: nconvlevmax = nuvzmax-1
integer,parameter :: na = nconvlevmax+1
! ntracermax maximum number of tracer species in convection
! nconvlevmax maximum number of levels for convection
! na parameter used in Emanuel's convect subroutine
!*********************************************
! Maximum dimensions of the nested input grids
!*********************************************
integer,parameter :: maxnests=0, nxmaxn=0, nymaxn=0
!integer,parameter :: maxnests=1,nxmaxn=251,nymaxn=151
! maxnests maximum number of nested grids
! nxmaxn,nymaxn maximum dimension of nested wind fields in
! x and y direction, respectively
!*********************************
! Parmaters for GRIB file decoding
!*********************************
......@@ -215,14 +201,12 @@ module par_mod
! Maximum number of particles, species, and similar
!**************************************************
integer,parameter :: maxpart=10000000
integer,parameter :: maxspec=4
integer,parameter :: maxpart=30000000
integer,parameter :: maxspec=1
real,parameter :: minmass=0.0001
! maxpart Maximum number of particles
! maxspec Maximum number of chemical species per release
! minmass Terminate particles carrying less mass
! maxpoint is also set dynamically during runtime
! maxpoint Maximum number of release locations
......@@ -231,16 +215,12 @@ module par_mod
! Sabine Eckhardt: change of landuse inventary numclass=13
! ---------
integer,parameter :: maxwf=50000, maxtable=1000, numclass=13, ni=11
integer,parameter :: numwfmem=2 ! Serial version/MPI with 2 fields
!integer,parameter :: numwfmem=3 ! MPI with 3 fields
! maxwf maximum number of wind fields to be used for simulation
! maxtable Maximum number of chemical species that can be
! tabulated for FLEXPART
! numclass Number of landuse classes available to FLEXPART
! ni Number of diameter classes of particles
! numwfmem Number of windfields kept in memory. 2 for serial
! version, 2 or 3 for MPI version
!**************************************************************************
! dimension of the OH field
......@@ -259,10 +239,10 @@ module par_mod
! Dimension of random number field
!*********************************
integer,parameter :: maxrand=200000000
integer,parameter :: maxrand=2000000
! maxrand number of random numbers used
!*****************************************************
! Number of clusters to be used for plume trajectories
......@@ -282,18 +262,7 @@ module par_mod
integer,parameter :: unitspecies=1, unitoutrecept=91, unitoutreceptppt=92
integer,parameter :: unitlsm=1, unitsurfdata=1, unitland=1, unitwesely=1
integer,parameter :: unitOH=1
integer,parameter :: unitdates=94, unitheader=90,unitheader_txt=100, unitshortpart=95, unitprecip=101
integer,parameter :: unitdates=94, unitheader=90, unitshortpart=95
integer,parameter :: unitboundcond=89
integer,parameter :: unittmp=101
!******************************************************
! integer code for missing values, used in wet scavenging (PS, 2012)
!******************************************************
integer,parameter :: icmv=-9999
! Parameters for testing
!*******************************************
! integer :: verbosity=0
end module par_mod
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment