Changes in src/par_mod.f90 [2eefa58:79e0349] in flexpart.git

File:

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

src/par_mod.f90

@@ -30 +30 @@
 !        Update 15 August 2013 IP                                              *
 !                                                                              *
+!        ESO 2016:                                                             *
+!          GFS specific parameters moved to gfs_mod.f90                        *
+!          ECMWF specific parameters moved to ecmwf_mod.f90                    *
 !                                                                              *
 !*******************************************************************************
@@ -77 +80 @@
   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
-  ! additional constants RLT Aug-2017
-  real,parameter :: rgas=8.31447 
-  real,parameter :: r_water=461.495
 
   ! pi                      number "pi"
@@ -89 +89 @@
   ! kappa                   exponent of formula for potential temperature
   ! vonkarman               von Karman constant
-  ! rgas                    universal gas constant [J/mol/K]
-  ! r_water                 specific gas constant for water vapor [J/kg/K]
 
   real,parameter :: karman=0.40, href=15., convke=2.0
@@ -150 +148 @@
   ! 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=0 ! 1.0 degree 138 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
@@ -156 +153 @@
 
 ! GFS
-   integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138
-   integer :: nxshift=0
+   integer,parameter :: nxmax=361,nymax=181,nuvzmax=138,nwzmax=138,nzmax=138,nxshift=0
 
 
@@ -282 +278 @@
 
   integer,parameter :: unitpath=1, unitcommand=1, unitageclasses=1, unitgrid=1
-  integer,parameter :: unitavailab=1, unitreleases=88, unitpartout=93, unitpartout_average=105
+  integer,parameter :: unitavailab=1, unitreleases=88, unitpartout=93
   integer,parameter :: unitpartin=93, unitflux=98, unitouttraj=96
   integer,parameter :: unitvert=1, unitoro=1, unitpoin=1, unitreceptor=1
@@ -292 +288 @@
   integer,parameter :: unitboundcond=89
   integer,parameter :: unittmp=101
-! RLT
-  integer,parameter :: unitoutfactor=102
 
 !******************************************************