Changeset 6a678e3 in flexpart.git for src/par_mod.f90

Timestamp:

Jan 25, 2016, 2:17:34 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:

26f6039

Parents:

fddc6ec

Message:

Added option to use double precision for calculating the deposition fields

File:

• src/par_mod.f90 (modified) (2 diffs)

src/par_mod.f90

@@ -41 +41 @@
 
   !****************************************************************
-  ! Parameter defining KIND parameter for double/single precision
-  ! 'cp' sets the precision for some internal calculations (sp or dp),
-  ! sp is default but dp can be used for increased precision.
-  ! WARNING: for now, do not use 'dp'
+  ! Parameters defining KIND parameter for double/single precision
   !****************************************************************
 
   integer,parameter :: dp=selected_real_kind(P=15)
   integer,parameter :: sp=selected_real_kind(6)
-  integer,parameter :: cp=sp ! sp is default
+
+  !****************************************************************
+  ! dep_prec sets the precision for deposition calculations (sp or 
+  ! dp). sp is default, dp can be used for increased precision.
+  !****************************************************************
+
+  integer,parameter :: dep_prec=dp
 
   !***********************************************************
@@ -79 +82 @@
   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]
+  real,parameter :: rho_water=1000. !ZHG 2015 [kg/m3]
   ! karman                  Karman's constant
   ! href [m]                Reference height for dry deposition