Ignore:
Timestamp:
Apr 10, 2019 9:31:22 AM (2 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/surface_data_output_mod.f90

    r3817 r3881  
    2525! -----------------
    2626! $Id$
     27! Check for zero output timestep (not allowed in parallel NetCDF output mode)
     28!
     29! 3817 2019-03-26 13:53:57Z suehring
    2730! Correct output coordinates of vertical surface elements
    2831!
     
    125128   USE netcdf_interface,                                                       &
    126129       ONLY:  netcdf_create_att, netcdf_create_dim, netcdf_create_file,        &
    127               netcdf_create_global_atts, netcdf_create_var, netcdf_handle_error
     130              netcdf_create_global_atts, netcdf_create_var, netcdf_data_format,&
     131              netcdf_handle_error
    128132
    129133   USE pegrid
     
    43164320#endif
    43174321!
     4322!--    In case of parallel NetCDF output the output timestep must not be zero.
     4323!--    This is because the number of requiered output timesteps is
     4324!--    pre-calculated, which is not possible with zero output timestep.
     4325       IF ( netcdf_data_format > 4 )  THEN
     4326          IF ( dt_dosurf == 0.0_wp )  THEN
     4327             message_string = 'dt_dosurf = 0.0 while using a ' //             &
     4328                              'variable timestep and parallel netCDF4 ' //    &
     4329                              'is not allowed.'
     4330             CALL message( 'surface_data_output_check_parameters', 'PA0081',  &
     4331                           1, 2, 0, 6, 0 )
     4332          ENDIF
     4333         
     4334          IF ( dt_dosurf_av == 0.0_wp )  THEN
     4335             message_string = 'dt_dosurf_av = 0.0 while using a ' //          &
     4336                              'variable timestep and parallel netCDF4 ' //    &
     4337                              'is not allowed.'
     4338             CALL message( 'surface_data_output_check_parameters', 'PA0081',  &
     4339                           1, 2, 0, 6, 0 )
     4340          ENDIF
     4341       ENDIF
     4342
     4343!
    43184344!--   In case of restart runs, check it the number of cores has been changed.
    43194345!--   With surface output this is not allowed.
Note: See TracChangeset for help on using the changeset viewer.