Changeset 3c1da52 in flexpart.git for src/par_mod.f90

Timestamp:

Jul 15, 2015, 9:45:31 AM (9 years ago)

Author:

Ignacio Pisso <Ignacio.Pisso@…>

Branches:

svn-petra

Children:

c07b09e

Parents:

31674b5

git-author:

Ignacio Pisso <Ignacio.Pisso@…> (07/15/15 09:23:56)

git-committer:

Ignacio Pisso <Ignacio.Pisso@…> (07/15/15 09:45:31)

Message:

from svn branch petra
r37 | pesei | 2015-04-22 17:42:35 +0200 (Wed, 22 Apr 2015) | 1 line
Changed paths:

M /branches/petra/src/FLEXPART.f90
M /branches/petra/src/com_mod.f90
M /branches/petra/src/getfields.f90
M /branches/petra/src/gridcheck.f90
M /branches/petra/src/interpol_rain.f90
M /branches/petra/src/interpol_rain_nests.f90
M /branches/petra/src/par_mod.f90
M /branches/petra/src/readcommand.f90
M /branches/petra/src/readspecies.f90
M /branches/petra/src/readwind.f90
M /branches/petra/src/readwind_nests.f90
M /branches/petra/src/timemanager.f90
A /branches/petra/src/verboutput.f90
M /branches/petra/src/verttransform.f90
M /branches/petra/src/verttransform_nests.f90
M /branches/petra/src/wetdepo.f90
M /branches/petra/src/wetdepokernel.f90
M /branches/petra/src/wetdepokernel_nest.f90

Wet dep quick fix and other small changes. Wet depo quick fix not final yet.

File:

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

src/par_mod.f90

@@ -103 +103 @@
 
   integer,parameter :: idiffnorm=10800, idiffmax=2*idiffnorm, minstep=1
+  integer,parameter :: itagekernmin=10800
 
   ! idiffnorm [s]           normal time interval between two wind fields
   ! idiffmax [s]            maximum time interval between two wind fields
   ! minstep [s]             minimum time step to be used within FLEXPART
-
+  ! itagekernmin [s]        minimum particle age [s] for using output kernel
 
   !*****************************************************************
@@ -265 +266 @@
 
 !******************************************************
-! integer code for missing values, used in wet scavenging (PS, 2012)
+! for use in wet scavenging (PS, 2012, 2015)
 !******************************************************
 
-    !  integer icmv
-    !  parameter(icmv=-9999)
-      integer,parameter ::  icmv=-9999
+  integer,parameter :: icmv=-9999 ! integer code for missing values 
+  ! some cloud heights used when trying to construct reasonable cloud:
+  integer,parameter :: icloudtopconvmin = 6000
+  integer,parameter :: icloudtopmin = 3000
+  integer,parameter :: icloudbot1 = 500, icloudtop1 =  8000
+  integer,parameter :: icloudbot2 = 0,   icloudtop2 = 10000
+  
+  real, parameter :: precmin = 0.002 ! minimum prec in mm/h for cloud diagnostics
+  real, parameter :: rhmininit = 0.90 ! standard condition for presence of clouds
+! PS       note that original by Sabine Eckhart was 80%
+! PS       however, for T<-20 C we consider saturation over ice
+! PS       so I think 90% should be enough          
+  real, parameter :: rhminx = 0.30 ! extreme condition for presence of clouds
 
 ! Parameters for testing