Changeset 3182 for palm/trunk/SOURCE/init_grid.f90

Timestamp:

Jul 27, 2018 1:36:03 PM (6 years ago)

Author:

suehring

Message:

New Inifor features: grid stretching, improved command-interface, support start dates in different formats in both YYYYMMDD and YYYYMMDDHH, Ability to manually control input file prefixes (--radiation-prefix, --soil-preifx, --flow-prefix, --soilmoisture-prefix) for compatiblity with DWD forcast naming scheme; GNU-style short and long option; Prepared output of large-scale forcing profiles (no computation yet); Added preprocessor flag netcdf4 to switch output format between netCDF 3 and 4; Updated netCDF variable names and attributes to comply with PIDS v1.9; Inifor bugfixes: Improved compatibility with older Intel Intel compilers by avoiding implicit array allocation; Added origin_lon/_lat values and correct reference time in dynamic driver global attributes; corresponding PALM changes: adjustments to revised Inifor; variables names in dynamic driver adjusted; enable geostrophic forcing also in offline nested mode; variable names in LES-LES and COSMO offline nesting changed; lateral boundary flags for nesting, in- and outflow conditions renamed

File:

• palm/trunk/SOURCE/init_grid.f90 (modified) (7 diffs)

palm/trunk/SOURCE/init_grid.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename variables in mesoscale-offline nesting mode
 ! 
 ! Former revisions:
@@ -342 +342 @@
         
     USE control_parameters,                                                    &
-        ONLY:  bc_lr_cyc, bc_ns_cyc, building_height, building_length_x,       &
+        ONLY:  bc_dirichlet_l, bc_dirichlet_n, bc_dirichlet_r,                 &
+               bc_dirichlet_s, bc_lr_cyc, bc_ns_cyc, bc_radiation_l,           &
+               bc_radiation_n, bc_radiation_r, bc_radiation_s,                 &
+               building_height, building_length_x,                             &
                building_length_y, building_wall_left, building_wall_south,     &
                canyon_height, canyon_wall_left, canyon_wall_south,             &
@@ -350 +353 @@
                dz_stretch_level_end_index, dz_stretch_level_start_index,       &
                dz_stretch_level_start, grid_level,                             &
-               force_bound_l, force_bound_r, force_bound_n, force_bound_s,     &
-               ibc_uv_b, inflow_l, inflow_n, inflow_r, inflow_s,               &
-               masking_method, maximum_grid_level, message_string,             &
-               momentum_advec, nest_domain, nest_bound_l,                      &
-               nest_bound_n, nest_bound_r, nest_bound_s,                       &
-               number_stretch_level_end, number_stretch_level_start, ocean,    &
-               outflow_l, outflow_n, outflow_r, outflow_s, psolver,            & 
-               scalar_advec, topography, topography_grid_convention,           &
-               tunnel_height, tunnel_length, tunnel_width_x, tunnel_width_y,   &
+               ibc_uv_b, masking_method, maximum_grid_level, message_string,   &
+               momentum_advec, number_stretch_level_end,                       &
+               number_stretch_level_start,ocean, psolver, scalar_advec,        &
+               topography, topography_grid_convention, tunnel_height,          &
+               tunnel_length, tunnel_width_x, tunnel_width_y,                  &
                tunnel_wall_depth, use_surface_fluxes, use_top_fluxes,          &
                wall_adjustment_factor
@@ -386 +385 @@
     IMPLICIT NONE
 
-    INTEGER(iwp) ::  i                           !< index variable along x 
-    INTEGER(iwp) ::  j                           !< index variable along y
-    INTEGER(iwp) ::  k                           !< index variable along z
-    INTEGER(iwp) ::  k_top                       !< topography top index on local PE
-    INTEGER(iwp) ::  n                           !< loop variable for stretching
-    INTEGER(iwp) ::  number_dz                   !< number of user-specified dz values       
-    INTEGER(iwp) ::  nzb_local_max               !< vertical grid index of maximum topography height
-    INTEGER(iwp) ::  nzb_local_min               !< vertical grid index of minimum topography height
+    INTEGER(iwp) ::  i             !< index variable along x 
+    INTEGER(iwp) ::  j             !< index variable along y
+    INTEGER(iwp) ::  k             !< index variable along z
+    INTEGER(iwp) ::  k_top         !< topography top index on local PE
+    INTEGER(iwp) ::  l             !< loop variable 
+    INTEGER(iwp) ::  n             !< loop variable for stretching
+    INTEGER(iwp) ::  number_dz     !< number of user-specified dz values       
+    INTEGER(iwp) ::  nzb_local_max !< vertical grid index of maximum topography height
+    INTEGER(iwp) ::  nzb_local_min !< vertical grid index of minimum topography height
                                       
     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  nzb_local  !< index for topography top at cell-center
@@ -834 +834 @@
     nzb_max = k_top + 1
 #endif
-    IF ( inflow_l  .OR.  outflow_l  .OR.  force_bound_l  .OR.  nest_bound_l  .OR.&
-         inflow_r  .OR.  outflow_r  .OR.  force_bound_r  .OR.  nest_bound_r  .OR.&
-         inflow_n  .OR.  outflow_n  .OR.  force_bound_n  .OR.  nest_bound_n  .OR.&
-         inflow_s  .OR.  outflow_s  .OR.  force_bound_s  .OR.  nest_bound_s )    &
+    IF ( bc_dirichlet_l  .OR.  bc_radiation_l     .OR.                         &
+         bc_dirichlet_r  .OR.  bc_radiation_r     .OR.                         &
+         bc_dirichlet_n  .OR.  bc_radiation_n     .OR.                         &
+         bc_dirichlet_s  .OR.  bc_radiation_s )                                &
          nzb_max = nzt
 !   
 !-- Finally, if topography extents up to the model top, limit nzb_max to nzt.
-    nzb_max = MIN( nzb_max, nzt )
+    nzb_max = MIN( nzb_max, nzt ) 
 !
 !-- Determine minimum index of topography. Usually, this will be nzb. In case
@@ -857 +857 @@
 !-- Initialize boundary conditions via surface type
     CALL init_bc
+
 !
 !-- Allocate and set topography height arrays required for data output
@@ -926 +927 @@
 #if defined( __parallel )
        CALL MPI_ALLREDUCE( MAXVAL( get_topography_top_index( 's' ) ),          &
-                           nzb_local_max, 1, MPI_INTEGER, MPI_MAX, comm2d, ierr )             
+                           nzb_local_max, 1, MPI_INTEGER, MPI_MAX, comm2d, ierr )               
 #else
        nzb_local_max = MAXVAL( get_topography_top_index( 's' ) )