Changeset 08a38b5 in flexpart.git for src/par_mod.f90

Timestamp:

Oct 3, 2017, 3:44:21 PM (7 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:

de4c5e9

Parents:

bb579a9

Message:

Added (hardcoded) option to output number of particles per grid cell.

File:

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

src/par_mod.f90

@@ -60 +60 @@
   logical, parameter :: lnokernel=.false.
 
+  !*********************************************************************
+  ! 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
+  ! number of directories/files used for FLEXPART input/output
   !***********************************************************
 
@@ -141 +146 @@
   !*********************************************
   
-!  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 :: nxshift=359  ! for ECMWF
-  integer,parameter :: nxshift=0     ! for GFS
+  ! ECMWF
+! integer,parameter :: nxmax=361,nymax=181,nuvzmax=92,nwzmax=92,nzmax=92,nxshift=359 ! 1.0 degree 92 level
+   integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=359 ! 1.0 degree 138 level
+! integer,parameter :: nxmax=721,nymax=361,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=359  ! 0.5 degree 138 level
+!  integer,parameter :: nxmax=181,nymax=91,nuvzmax=92,nwzmax=92,nzmax=92,nxshift=0  ! CERA 2.0 degree 92 level
+
+! GFS
+!   integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=0
+
 
   !*********************************************
@@ -155 +160 @@
   !*********************************************
 
-  integer,parameter :: maxnests=1,nxmaxn=451,nymaxn=226
+  integer,parameter :: maxnests=0,nxmaxn=451,nymaxn=226
 
   ! nxmax,nymax        maximum dimension of wind fields in x and y
@@ -213 +218 @@
   !**************************************************
 
-  integer,parameter :: maxpart=10000000
+  integer,parameter :: maxpart=3000000
   integer,parameter :: maxspec=4