Changeset 3182 for palm/trunk/SOURCE/check_parameters.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/check_parameters.f90 (modified) (12 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Rename boundary conditions in offline nesting
 ! 
 ! Former revisions:
@@ -1069 +1069 @@
 !-- Check turbulence closure setup
     CALL tcm_check_parameters
+
 !
 !-- Check approximation
@@ -1437 +1438 @@
        CALL MPI_ALLREDUCE( dt_spinup,   dt_spinup_max,   1, MPI_REAL,          &
                            MPI_MAX, MPI_COMM_WORLD, ierr )
+
        IF ( spinup_time /= spinup_time_max  .OR.  dt_spinup /= dt_spinup_max ) &
        THEN
@@ -1781 +1783 @@
 !
 !-- Overwrite latitude if necessary and compute Coriolis parameter. 
-!-- To do - move initialization of f and fs to coriolis_mod.
+!-- @todo - move initialization of f and fs to coriolis_mod.
     IF ( input_pids_static )  THEN
        latitude  = init_model%latitude
@@ -1976 +1978 @@
 !-- TO_DO: later set bc_p_t to neumann before, in case of nested domain
     ELSEIF ( bc_p_t == 'neumann' .OR. bc_p_t == 'nested'  .OR.                 &
-             bc_p_t == 'forcing'  )  THEN
+             bc_p_t == 'nesting_offline'  )  THEN
        ibc_p_t = 1
     ELSE
@@ -2006 +2008 @@
     ELSEIF ( bc_pt_t == 'initial_gradient' )  THEN
        ibc_pt_t = 2
-    ELSEIF ( bc_pt_t == 'nested'  .OR.  bc_pt_t == 'forcing' )  THEN
+    ELSEIF ( bc_pt_t == 'nested'  .OR.  bc_pt_t == 'nesting_offline' )  THEN
        ibc_pt_t = 3
     ELSE
@@ -2234 +2236 @@
        ELSEIF ( bc_uv_t == 'neumann' )  THEN
           ibc_uv_t = 1
-       ELSEIF ( bc_uv_t == 'nested'  .OR.  bc_uv_t == 'forcing' )  THEN
+       ELSEIF ( bc_uv_t == 'nested'  .OR.  bc_uv_t == 'nesting_offline' )  THEN
           ibc_uv_t = 3
        ELSE
@@ -3926 +3928 @@
           dist_nxl    = MAX( inflow_disturbance_begin, nxl )
           dist_nxr(1) = MIN( inflow_disturbance_end, nxr )
-       ELSEIF ( bc_lr == 'nested'  .OR.  bc_lr == 'forcing' )  THEN
+       ELSEIF ( bc_lr == 'nested'  .OR.  bc_lr == 'nesting_offline' )  THEN
           dist_nxl    = MAX( inflow_disturbance_begin, nxl )
           dist_nxr    = MIN( nx - inflow_disturbance_begin, nxr )
@@ -3936 +3938 @@
           dist_nys    = MAX( inflow_disturbance_begin, nys )
           dist_nyn(1) = MIN( inflow_disturbance_end, nyn )
-       ELSEIF ( bc_ns == 'nested'  .OR.  bc_ns == 'forcing' )  THEN
+       ELSEIF ( bc_ns == 'nested'  .OR.  bc_ns == 'nesting_offline' )  THEN
           dist_nys    = MAX( inflow_disturbance_begin, nys )
           dist_nyn    = MIN( ny - inflow_disturbance_begin, nyn )
@@ -3949 +3951 @@
           dist_nxl    = inflow_disturbance_begin
           dist_nxr(1) = inflow_disturbance_end
-       ELSEIF ( bc_lr == 'nested'  .OR.  bc_lr == 'forcing' )  THEN
+       ELSEIF ( bc_lr == 'nested'  .OR.  bc_lr == 'nesting_offline' )  THEN
           dist_nxr    = nx - inflow_disturbance_begin
           dist_nxl    = inflow_disturbance_begin
@@ -3959 +3961 @@
           dist_nys    = inflow_disturbance_begin
           dist_nyn(1) = inflow_disturbance_end
-       ELSEIF ( bc_ns == 'nested'  .OR.  bc_ns == 'forcing' )  THEN
+       ELSEIF ( bc_ns == 'nested'  .OR.  bc_ns == 'nesting_offline' )  THEN
           dist_nyn    = ny - inflow_disturbance_begin
           dist_nys    = inflow_disturbance_begin
@@ -4403 +4405 @@
        ELSEIF ( bc_t == 'initial_gradient' )  THEN
           ibc_t = 2
-       ELSEIF ( bc_t == 'nested'  .OR.  bc_t == 'forcing' )  THEN
+       ELSEIF ( bc_t == 'nested'  .OR.  bc_t == 'nesting_offline' )  THEN
           ibc_t = 3
        ELSE