Changeset 3881 for palm/trunk/SOURCE/surface_data_output_mod.f90

Timestamp:

Apr 10, 2019 9:31:22 AM (5 years ago)

Author:

suehring

Message:

Make level 3 initialization of urban-surfaces consistent to input data standard; revise flagging of ground-floor level surfaces at buidlings; bugfixes in level 3 initialization of albedo; bugfix in OpenMP directive; check for zero output timestep in surface output; assure that Obukhov lenght does not become zero

File:

• palm/trunk/SOURCE/surface_data_output_mod.f90 (modified) (3 diffs)

palm/trunk/SOURCE/surface_data_output_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Check for zero output timestep (not allowed in parallel NetCDF output mode) 
+! 
+! 3817 2019-03-26 13:53:57Z suehring
 ! Correct output coordinates of vertical surface elements
 ! 
@@ -125 +128 @@
    USE netcdf_interface,                                                       &
        ONLY:  netcdf_create_att, netcdf_create_dim, netcdf_create_file,        &
-              netcdf_create_global_atts, netcdf_create_var, netcdf_handle_error
+              netcdf_create_global_atts, netcdf_create_var, netcdf_data_format,&
+              netcdf_handle_error
 
    USE pegrid
@@ -4316 +4320 @@
 #endif
 !
+!--    In case of parallel NetCDF output the output timestep must not be zero.
+!--    This is because the number of requiered output timesteps is 
+!--    pre-calculated, which is not possible with zero output timestep. 
+       IF ( netcdf_data_format > 4 )  THEN
+          IF ( dt_dosurf == 0.0_wp )  THEN
+             message_string = 'dt_dosurf = 0.0 while using a ' //             &
+                              'variable timestep and parallel netCDF4 ' //    &
+                              'is not allowed.'
+             CALL message( 'surface_data_output_check_parameters', 'PA0081',  &
+                           1, 2, 0, 6, 0 )
+          ENDIF
+          
+          IF ( dt_dosurf_av == 0.0_wp )  THEN
+             message_string = 'dt_dosurf_av = 0.0 while using a ' //          &
+                              'variable timestep and parallel netCDF4 ' //    &
+                              'is not allowed.'
+             CALL message( 'surface_data_output_check_parameters', 'PA0081',  &
+                           1, 2, 0, 6, 0 )
+          ENDIF
+       ENDIF
+
+!
 !--   In case of restart runs, check it the number of cores has been changed.
 !--   With surface output this is not allowed.