Changeset 2968 for palm/trunk/SOURCE/modules.f90

Timestamp:

Apr 13, 2018 11:52:24 AM (6 years ago)

Author:

suehring

Message:

Bugfixes in initalization of land-surface model as well as timeseries data output in case of elevated model surfaces.

File:

• palm/trunk/SOURCE/modules.f90 (modified) (2 diffs)

palm/trunk/SOURCE/modules.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! +topo_min_level
+! 
+! 2964 2018-04-12 16:04:03Z raasch
 ! *_time_count variables are all initialized with zero now
 ! 
@@ -1719 +1722 @@
     USE kinds
 
-    INTEGER(iwp) ::  nbgp = 3     !< number of boundary ghost points
-    INTEGER(iwp) ::  ngp_sums     !< number of vertical profile grid points time number of output profiles - used for allreduce statements in MPI calls
-    INTEGER(iwp) ::  ngp_sums_ls  !< number of vertical profile grid points time number of large-scale forcing profiles - used for allreduce statements in MPI calls
-    INTEGER(iwp) ::  nnx          !< number of subdomain grid points in x-direction
-    INTEGER(iwp) ::  nx = 0       !< nx+1 = total number of grid points in x-direction
-    INTEGER(iwp) ::  nx_a         !< in coupled atmosphere-ocean runs: total number of grid points along x (atmosphere)
-    INTEGER(iwp) ::  nx_o         !< in coupled atmosphere-ocean runs: total number of grid points along x (ocean)
-    INTEGER(iwp) ::  nxl          !< left-most grid index of subdomain (excluding ghost points)
-    INTEGER(iwp) ::  nxlg         !< left-most grid index of subdomain (including ghost points)
-    INTEGER(iwp) ::  nxlu         !< =nxl+1 (at left domain boundary with inflow from left), else =nxl (used for u-velocity component)
-    INTEGER(iwp) ::  nxr          !< right-most grid index of subdomain (excluding ghost points)
-    INTEGER(iwp) ::  nxrg         !< right-most grid index of subdomain (including ghost points)
-    INTEGER(iwp) ::  nx_on_file   !< nx of previous run in job chain
-    INTEGER(iwp) ::  nny          !< number of subdomain grid points in y-direction
-    INTEGER(iwp) ::  ny = 0       !< ny+1 = total number of grid points in y-direction
-    INTEGER(iwp) ::  ny_a         !< in coupled atmosphere-ocean runs: total number of grid points along y (atmosphere)
-    INTEGER(iwp) ::  ny_o         !< in coupled atmosphere-ocean runs: total number of grid points along y (ocean)
-    INTEGER(iwp) ::  nyn          !< north-most grid index of subdomain (excluding ghost points)
-    INTEGER(iwp) ::  nyng         !< north-most grid index of subdomain (including ghost points)
-    INTEGER(iwp) ::  nys          !< south-most grid index of subdomain (excluding ghost points)
-    INTEGER(iwp) ::  nysg         !< south-most grid index of subdomain (including ghost points)
-    INTEGER(iwp) ::  nysv         !< =nys+1 (at south domain boundary with inflow from south), else =nys (used for v-velocity component)
-    INTEGER(iwp) ::  ny_on_file   !< ny of previous run in job chain
-    INTEGER(iwp) ::  nnz          !< number of subdomain grid points in z-direction
-    INTEGER(iwp) ::  nz = 0       !< total number of grid points in z-direction
-    INTEGER(iwp) ::  nzb          !< bottom grid index of computational domain
-    INTEGER(iwp) ::  nzb_diff     !< will be removed
-    INTEGER(iwp) ::  nzb_max      !< vertical index of topography top
-    INTEGER(iwp) ::  nzt          !< nzt+1 = top grid index of computational domain
+    INTEGER(iwp) ::  nbgp = 3       !< number of boundary ghost points
+    INTEGER(iwp) ::  ngp_sums       !< number of vertical profile grid points time number of output profiles - used for allreduce statements in MPI calls
+    INTEGER(iwp) ::  ngp_sums_ls    !< number of vertical profile grid points time number of large-scale forcing profiles - used for allreduce statements in MPI calls
+    INTEGER(iwp) ::  nnx            !< number of subdomain grid points in x-direction
+    INTEGER(iwp) ::  nx = 0         !< nx+1 = total number of grid points in x-direction
+    INTEGER(iwp) ::  nx_a           !< in coupled atmosphere-ocean runs: total number of grid points along x (atmosphere)
+    INTEGER(iwp) ::  nx_o           !< in coupled atmosphere-ocean runs: total number of grid points along x (ocean)
+    INTEGER(iwp) ::  nxl            !< left-most grid index of subdomain (excluding ghost points)
+    INTEGER(iwp) ::  nxlg           !< left-most grid index of subdomain (including ghost points)
+    INTEGER(iwp) ::  nxlu           !< =nxl+1 (at left domain boundary with inflow from left), else =nxl (used for u-velocity component)
+    INTEGER(iwp) ::  nxr            !< right-most grid index of subdomain (excluding ghost points)
+    INTEGER(iwp) ::  nxrg           !< right-most grid index of subdomain (including ghost points)
+    INTEGER(iwp) ::  nx_on_file     !< nx of previous run in job chain
+    INTEGER(iwp) ::  nny            !< number of subdomain grid points in y-direction
+    INTEGER(iwp) ::  ny = 0         !< ny+1 = total number of grid points in y-direction
+    INTEGER(iwp) ::  ny_a           !< in coupled atmosphere-ocean runs: total number of grid points along y (atmosphere)
+    INTEGER(iwp) ::  ny_o           !< in coupled atmosphere-ocean runs: total number of grid points along y (ocean)
+    INTEGER(iwp) ::  nyn            !< north-most grid index of subdomain (excluding ghost points)
+    INTEGER(iwp) ::  nyng           !< north-most grid index of subdomain (including ghost points)
+    INTEGER(iwp) ::  nys            !< south-most grid index of subdomain (excluding ghost points)
+    INTEGER(iwp) ::  nysg           !< south-most grid index of subdomain (including ghost points)
+    INTEGER(iwp) ::  nysv           !< =nys+1 (at south domain boundary with inflow from south), else =nys (used for v-velocity component)
+    INTEGER(iwp) ::  ny_on_file     !< ny of previous run in job chain
+    INTEGER(iwp) ::  nnz            !< number of subdomain grid points in z-direction
+    INTEGER(iwp) ::  nz = 0         !< total number of grid points in z-direction
+    INTEGER(iwp) ::  nzb            !< bottom grid index of computational domain
+    INTEGER(iwp) ::  nzb_diff       !< will be removed
+    INTEGER(iwp) ::  nzb_max        !< vertical index of topography top
+    INTEGER(iwp) ::  nzt            !< nzt+1 = top grid index of computational domain
+    INTEGER(iwp) ::  topo_min_level !< minimum topography-top index (usually equal to nzb)
 
     INTEGER(idp), DIMENSION(:), ALLOCATABLE ::  ngp_3d        !< number of grid points of the total domain