Changeset 3045


Ignore:
Timestamp:
May 28, 2018 7:55:41 AM (6 years ago)
Author:
Giersch
Message:

Code adjusted according to coding standards, renamed namelists, error messages revised until PA0347, output CASE 108 disabled

Location:
palm/trunk
Files:
42 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/INSTALL/example_cbl_p3d

    r2490 r3045  
    33!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/inipar
    44!-------------------------------------------------------------------------------
    5 &inipar
     5&initialization_parameters
    66!
    77!-- grid parameters
     
    4949!   Documentation: https://palm.muk.uni-hannover.de/trac/wiki/doc/app/d3par
    5050!-------------------------------------------------------------------------------
    51 &d3par
     51&runtime_parameters
    5252!
    5353!-- run steering
  • palm/trunk/INSTALL/test_oceanml_p3d

    r2575 r3045  
    1 &inipar nx = 63, ny = 63, nz=64,
     1&initialization_parameters
     2        nx = 63, ny = 63, nz=64,
    23        dx = 1.25, dy = 1.25, dz = 1.25,
    34
     
    3233        bc_sa_t = 'neumann', /
    3334
    34 &d3par  end_time = 300.0,
     35&runtime_parameters
     36        end_time = 300.0,
    3537        create_disturbances = .T.,
    3638        dt_disturb = 60.0,
  • palm/trunk/SOURCE/check_for_restart.f90

    r2718 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    141144    IF ( terminate_run  .AND.  myid == 0 )  THEN
    142145       WRITE( message_string, * ) 'run will be terminated because it is ',     &
    143                        'running out of job cpu limit & ',                      &
     146                       'running out of job cpu limit, ',                       &
    144147                       'remaining time:         ', remaining_time, ' s',       &
    145148                       'termination time needed:', termination_time_needed, ' s'
     
    184187                                  'forced a job finalization using a flag',    &
    185188                                  'file:',                                     &
    186                                   '&DO_STOP_NOW: ', do_stop_now,               &
    187                                   '&DO_RESTART_NOW: ', do_restart_now
     189                                  ' DO_STOP_NOW: ', do_stop_now,               &
     190                                  ' DO_RESTART_NOW: ', do_restart_now
    188191          CALL message( 'check_for_restart', 'PA0398', 0, 0, 0, 6, 0 )
    189192
     
    247250
    248251          WRITE( message_string, * ) 'run will be terminated due to user ',    &
    249                                   'settings of',                               &
    250                                   '&restart_time / dt_restart',                &
    251                                   '&new restart time is: ', time_restart, ' s'
     252                                  'settings of ',                              &
     253                                  'restart_time / dt_restart, ',               &
     254                                  'new restart time is: ', time_restart, ' s'
    252255          CALL message( 'check_for_restart', 'PA0164', 0, 0, 0, 6, 0 )
    253256 
  • palm/trunk/SOURCE/check_open.f90

    r2964 r3045  
    2525! -----------------
    2626! $Id$
     27! Output case 108 disabled
     28!
     29! 2964 2018-04-12 16:04:03Z Giersch
    2730! Error message moved to radiation_model_mod
    2831!
     
    971974          ENDIF
    972975
    973 
    974        CASE ( 108 )
    975 
    976           IF ( myid_char == '' )  THEN
    977              filename = 'DATA_PRT_NETCDF' // TRIM( coupling_char )
    978           ELSE
    979              filename = 'DATA_PRT_NETCDF' // TRIM( coupling_char ) // '/' //   &
    980                         myid_char
    981           ENDIF
     976!
     977!--     Currently disabled
     978!       CASE ( 108 )
     979
     980!          IF ( myid_char == '' )  THEN
     981!             filename = 'DATA_PRT_NETCDF' // TRIM( coupling_char )
     982!          ELSE
     983!             filename = 'DATA_PRT_NETCDF' // TRIM( coupling_char ) // '/' //   &
     984!                        myid_char
     985!          ENDIF
    982986!
    983987!--       Inquire, if there is a netCDF file from a previuos run. This should
    984988!--       be opened for extension, if its variables match the actual run.
    985           INQUIRE( FILE=filename, EXIST=netcdf_extend )
    986 
    987           IF ( netcdf_extend )  THEN
     989!          INQUIRE( FILE=filename, EXIST=netcdf_extend )
     990
     991!          IF ( netcdf_extend )  THEN
    988992!
    989993!--          Open an existing netCDF file for output
    990              CALL netcdf_open_write_file( filename, id_set_prt, .FALSE., 41 )
     994!             CALL netcdf_open_write_file( filename, id_set_prt, .FALSE., 41 )
    991995!
    992996!--          Read header information and set all ids. If there is a mismatch
    993997!--          between the previuos and the actual run, netcdf_extend is returned
    994998!--          as .FALSE.
    995              CALL netcdf_define_header( 'pt', netcdf_extend, 0 )
     999!             CALL netcdf_define_header( 'pt', netcdf_extend, 0 )
    9961000
    9971001!
    9981002!--          Remove the local file, if it can not be extended
    999              IF ( .NOT. netcdf_extend )  THEN
    1000                 nc_stat = NF90_CLOSE( id_set_prt )
    1001                 CALL netcdf_handle_error( 'check_open', 42 )
    1002                 CALL local_system( 'rm ' // TRIM( filename ) )
    1003              ENDIF
    1004 
    1005           ENDIF         
    1006 
    1007           IF ( .NOT. netcdf_extend )  THEN
     1003!             IF ( .NOT. netcdf_extend )  THEN
     1004!                nc_stat = NF90_CLOSE( id_set_prt )
     1005!                CALL netcdf_handle_error( 'check_open', 42 )
     1006!                CALL local_system( 'rm ' // TRIM( filename ) )
     1007!             ENDIF
     1008
     1009!          ENDIF         
     1010
     1011!          IF ( .NOT. netcdf_extend )  THEN
    10081012
    10091013!
    10101014!--          For runs on multiple processors create the subdirectory
    1011              IF ( myid_char /= '' )  THEN
    1012                 IF ( myid == 0  .AND. .NOT. openfile(file_id)%opened_before )  &
    1013                 THEN    ! needs modification in case of non-extendable sets
    1014                    CALL local_system( 'mkdir  DATA_PRT_NETCDF' //              &
    1015                                        TRIM( coupling_char ) // '/' )
    1016                 ENDIF
     1015!             IF ( myid_char /= '' )  THEN
     1016!                IF ( myid == 0  .AND. .NOT. openfile(file_id)%opened_before )  &
     1017!                THEN    ! needs modification in case of non-extendable sets
     1018!                   CALL local_system( 'mkdir  DATA_PRT_NETCDF' //              &
     1019!                                       TRIM( coupling_char ) // '/' )
     1020!                ENDIF
    10171021#if defined( __parallel )
    10181022!
    10191023!--             Set a barrier in order to allow that all other processors in the
    10201024!--             directory created by PE0 can open their file
    1021                 CALL MPI_BARRIER( comm2d, ierr )
     1025!                CALL MPI_BARRIER( comm2d, ierr )
    10221026#endif
    1023              ENDIF
     1027!             ENDIF
    10241028
    10251029!
    10261030!--          Create a new netCDF output file with requested netCDF format
    1027              CALL netcdf_create_file( filename, id_set_prt, .FALSE., 43 )
     1031!             CALL netcdf_create_file( filename, id_set_prt, .FALSE., 43 )
    10281032
    10291033!
    10301034!--          Define the header
    1031              CALL netcdf_define_header( 'pt', netcdf_extend, 0 )
    1032 
    1033           ENDIF
     1035!             CALL netcdf_define_header( 'pt', netcdf_extend, 0 )
     1036
     1037!          ENDIF
    10341038
    10351039       CASE ( 109 )
  • palm/trunk/SOURCE/check_parameters.f90

    r3035 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 3035 2018-05-24 09:35:20Z schwenkel
    2730! Add option to initialize warm air bubble close to surface
    28 ! 
     31!
    2932! 3034 2018-05-24 08:41:20Z raasch
    3033! bugfix: check that initializing_actions has been set
    31 ! 
     34!
    3235! 2980 2018-04-17 15:19:27Z suehring
    3336! Further improvement for spinup checks.
     
    751754!
    752755!-- Check the coupling mode
    753 !> @todo Check if any queries for other coupling modes (e.g. precursor_ocean) are missing
    754     IF ( coupling_mode /= 'uncoupled'            .AND.  &
    755          coupling_mode /= 'vnested_crse'         .AND.  &
    756          coupling_mode /= 'vnested_fine'         .AND.  &
    757          coupling_mode /= 'atmosphere_to_ocean'  .AND.  &
     756    IF ( coupling_mode /= 'uncoupled'            .AND.                         &
     757         coupling_mode /= 'precursor_atmos'      .AND.                         &
     758         coupling_mode /= 'precursor_ocean'      .AND.                         &
     759         coupling_mode /= 'vnested_crse'         .AND.                         &
     760         coupling_mode /= 'vnested_fine'         .AND.                         &
     761         coupling_mode /= 'atmosphere_to_ocean'  .AND.                         &
    758762         coupling_mode /= 'ocean_to_atmosphere' )  THEN
    759763       message_string = 'illegal coupling mode: ' // TRIM( coupling_mode )
     
    764768!-- Check if humidity is set to TRUE in case of the atmospheric run (for coupled runs)
    765769    IF ( coupling_mode == 'atmosphere_to_ocean' .AND. .NOT. humidity) THEN
    766        message_string = ' Humidity has to be set to .T. in the _p3d file for ' //      &
    767                         'coupled runs between ocean and atmosphere.'
     770       message_string = ' Humidity has to be set to .T. in the _p3d file ' //  &
     771                        'for coupled runs between ocean and atmosphere.'
    768772       CALL message( 'check_parameters', 'PA0476', 1, 2, 0, 6, 0 )
    769773    ENDIF
     
    771775!
    772776!-- Check dt_coupling, restart_time, dt_restart, end_time, dx, dy, nx and ny
    773     IF ( coupling_mode /= 'uncoupled'  .AND.                                   &
    774          coupling_mode(1:8) /= 'vnested_' )  THEN
     777    IF ( coupling_mode /= 'uncoupled'       .AND.                              &
     778         coupling_mode(1:8) /= 'vnested_'   .AND.                              &
     779         coupling_mode /= 'precursor_atmos' .AND.                              &
     780         coupling_mode /= 'precursor_ocean' )  THEN
    775781
    776782       IF ( dt_coupling == 9999999.9_wp )  THEN
     
    914920             WRITE( message_string, * ) 'coupling mode "',                     &
    915921                   TRIM( coupling_mode ),                                      &
    916              '": nx+1 in ocean is not divisible by nx+1 in', &
     922             '": nx+1 in ocean is not divisible by nx+1 in',                   &
    917923             ' atmosphere without remainder'
    918924             CALL message( 'check_parameters', 'PA0339', 1, 2, 0, 6, 0 )
     
    10541060!
    10551061!-- Check approximation
    1056     IF ( TRIM( approximation ) /= 'boussinesq'   .AND.   &
     1062    IF ( TRIM( approximation ) /= 'boussinesq'   .AND.                         &
    10571063         TRIM( approximation ) /= 'anelastic' )  THEN
    1058        message_string = 'unknown approximation: approximation = "' //    &
     1064       message_string = 'unknown approximation: approximation = "' //          &
    10591065                        TRIM( approximation ) // '"'
    10601066       CALL message( 'check_parameters', 'PA0446', 1, 2, 0, 6, 0 )
     
    10631069!
    10641070!-- Check approximation requirements
    1065     IF ( TRIM( approximation ) == 'anelastic'   .AND.   &
     1071    IF ( TRIM( approximation ) == 'anelastic'   .AND.                          &
    10661072         TRIM( momentum_advec ) /= 'ws-scheme' )  THEN
    1067        message_string = 'Anelastic approximation requires: ' //                &
     1073       message_string = 'Anelastic approximation requires ' //                 &
    10681074                        'momentum_advec = "ws-scheme"'
    10691075       CALL message( 'check_parameters', 'PA0447', 1, 2, 0, 6, 0 )
    10701076    ENDIF
    1071     IF ( TRIM( approximation ) == 'anelastic'   .AND.   &
     1077    IF ( TRIM( approximation ) == 'anelastic'   .AND.                          &
    10721078         TRIM( psolver ) == 'multigrid' )  THEN
    1073        message_string = 'Anelastic approximation currently only supports: ' // &
    1074                         'psolver = "poisfft", ' // &
    1075                         'psolver = "sor" and ' // &
     1079       message_string = 'Anelastic approximation currently only supports ' // &
     1080                        'psolver = "poisfft", ' //                             &
     1081                        'psolver = "sor" and ' //                              &
    10761082                        'psolver = "multigrid_noopt"'
    10771083       CALL message( 'check_parameters', 'PA0448', 1, 2, 0, 6, 0 )
    10781084    ENDIF
    1079     IF ( TRIM( approximation ) == 'anelastic'   .AND.   &
     1085    IF ( TRIM( approximation ) == 'anelastic'   .AND.                          &
    10801086         conserve_volume_flow )  THEN
    1081        message_string = 'Anelastic approximation is not allowed with:' //      &
     1087       message_string = 'Anelastic approximation is not allowed with ' //      &
    10821088                        'conserve_volume_flow = .TRUE.'
    10831089       CALL message( 'check_parameters', 'PA0449', 1, 2, 0, 6, 0 )
     
    10861092!
    10871093!-- Check flux input mode
    1088     IF ( TRIM( flux_input_mode ) /= 'dynamic'    .AND.   &
    1089          TRIM( flux_input_mode ) /= 'kinematic'  .AND.   &
     1094    IF ( TRIM( flux_input_mode ) /= 'dynamic'    .AND.                         &
     1095         TRIM( flux_input_mode ) /= 'kinematic'  .AND.                         &
    10901096         TRIM( flux_input_mode ) /= 'approximation-specific' )  THEN
    10911097       message_string = 'unknown flux input mode: flux_input_mode = "' //      &
     
    11041110!
    11051111!-- Check flux output mode
    1106     IF ( TRIM( flux_output_mode ) /= 'dynamic'    .AND.   &
    1107          TRIM( flux_output_mode ) /= 'kinematic'  .AND.   &
     1112    IF ( TRIM( flux_output_mode ) /= 'dynamic'    .AND.                        &
     1113         TRIM( flux_output_mode ) /= 'kinematic'  .AND.                        &
    11081114         TRIM( flux_output_mode ) /= 'approximation-specific' )  THEN
    11091115       message_string = 'unknown flux output mode: flux_output_mode = "' //    &
     
    12481254!
    12491255!-- Advection schemes:
    1250     IF ( momentum_advec /= 'pw-scheme'  .AND.  momentum_advec /= 'ws-scheme' ) &
     1256    IF ( momentum_advec /= 'pw-scheme'  .AND.                                  & 
     1257         momentum_advec /= 'ws-scheme'  .AND.                                  &
     1258         momentum_advec /= 'up-scheme' )                                       &
    12511259    THEN
    12521260       message_string = 'unknown advection scheme: momentum_advec = "' //      &
     
    12591267    THEN
    12601268       message_string = 'momentum_advec or scalar_advec = "'                   &
    1261          // TRIM( momentum_advec ) // '" is not allowed with timestep_scheme = "' // &
    1262          TRIM( timestep_scheme ) // '"'
     1269         // TRIM( momentum_advec ) // '" is not allowed with ' //              &
     1270         'timestep_scheme = "' // TRIM( timestep_scheme ) // '"'
    12631271       CALL message( 'check_parameters', 'PA0023', 1, 2, 0, 6, 0 )
    12641272    ENDIF
    12651273    IF ( scalar_advec /= 'pw-scheme'  .AND.  scalar_advec /= 'ws-scheme' .AND. &
    1266          scalar_advec /= 'bc-scheme' )                                         &
     1274         scalar_advec /= 'bc-scheme' .AND. scalar_advec /= 'up-scheme' )       &
    12671275    THEN
    12681276       message_string = 'unknown advection scheme: scalar_advec = "' //        &
     
    12731281    THEN
    12741282       message_string = 'advection_scheme scalar_advec = "'                    &
    1275          // TRIM( scalar_advec ) // '" not implemented for & loop_optimization = "' // &
    1276          TRIM( loop_optimization ) // '"'
     1283         // TRIM( scalar_advec ) // '" not implemented for ' //                &
     1284         'loop_optimization = "' // TRIM( loop_optimization ) // '"'
    12771285       CALL message( 'check_parameters', 'PA0026', 1, 2, 0, 6, 0 )
    12781286    ENDIF
     
    13251333!-- Check for proper settings for microphysics
    13261334    IF ( cloud_physics  .AND.  cloud_droplets )  THEN
    1327        message_string = 'cloud_physics = .TRUE. is not allowed with ' //  &
     1335       message_string = 'cloud_physics = .TRUE. is not allowed with ' //       &
    13281336                        'cloud_droplets = .TRUE.'
    13291337       CALL message( 'check_parameters', 'PA0442', 1, 2, 0, 6, 0 )
     
    13551363    IF ( TRIM( initializing_actions ) == '' )  THEN
    13561364       message_string = 'no value specified for initializing_actions'
    1357        CALL message( 'check_parameters', 'PA0017', 1, 2, 0, 6, 0 )
     1365       CALL message( 'check_parameters', 'PA0149', 1, 2, 0, 6, 0 )
    13581366    ENDIF
    13591367
     
    13731381
    13741382             CASE DEFAULT
    1375                 message_string = 'initializing_actions = "' //                 &
     1383                message_string = 'initializing_action = "' //                  &
    13761384                                 TRIM( action ) // '" unknown or not allowed'
    13771385                CALL message( 'check_parameters', 'PA0030', 1, 2, 0, 6, 0 )
     
    14721480    IF ( large_scale_forcing  .OR.  nudging )  CALL lsf_nudging_check_parameters
    14731481
    1474 
    1475 
    1476     IF ( .NOT. ( loop_optimization == 'cache'  .OR.                            &
    1477                  loop_optimization == 'vector' )                               &
    1478          .AND.  cloud_physics  .AND.  microphysics_seifert )  THEN
    1479        message_string = 'cloud_scheme = seifert_beheng requires ' //           &
    1480                         'loop_optimization = "cache" or "vector"'
    1481        CALL message( 'check_parameters', 'PA0362', 1, 2, 0, 6, 0 )
    1482     ENDIF
    1483 
    14841482!
    14851483!-- In case of no model continuation run, check initialising parameters and
     
    18251823!-- Check time step and cfl_factor
    18261824    IF ( dt /= -1.0_wp )  THEN
    1827        IF ( dt <= 0.0_wp  .AND.  dt /= -1.0_wp )  THEN
     1825       IF ( dt <= 0.0_wp )  THEN
    18281826          WRITE( message_string, * ) 'dt = ', dt , ' <= 0.0'
    18291827          CALL message( 'check_parameters', 'PA0044', 1, 2, 0, 6, 0 )
     
    18441842       ELSE
    18451843          WRITE( message_string, * ) 'cfl_factor = ', cfl_factor,              &
    1846                  ' out of range & 0.0 < cfl_factor <= 1.0 is required'
     1844                 ' out of range 0.0 < cfl_factor <= 1.0 is required'
    18471845          CALL message( 'check_parameters', 'PA0045', 1, 2, 0, 6, 0 )
    18481846       ENDIF
     
    18871885          CALL message( 'check_parameters', 'PA0047', 1, 2, 0, 6, 0 )
    18881886       ELSE
    1889           message_string = 'variable translation speed used for galilei-' //   &
    1890              'transformation, which may cause & instabilities in stably ' //   &
     1887          message_string = 'variable translation speed used for Galilei-' //   &
     1888             'transformation, which may cause instabilities in stably ' //     &
    18911889             'stratified regions'
    18921890          CALL message( 'check_parameters', 'PA0048', 0, 1, 0, 6, 0 )
     
    20642062         surface_heatflux /= 0.0_wp )  THEN
    20652063       message_string = 'boundary_condition: bc_pt_b = "' // TRIM( bc_pt_b ) //&
    2066                         '& is not allowed with constant_heatflux = .TRUE.'
     2064                        'is not allowed with constant_heatflux = .TRUE.'
    20672065       CALL message( 'check_parameters', 'PA0065', 1, 2, 0, 6, 0 )
    20682066    ENDIF
     
    21542152    IF ( passive_scalar )  THEN
    21552153
    2156        IF ( ANY( wall_scalarflux /= 0.0_wp )  .AND.                        &
     2154       IF ( ANY( wall_scalarflux /= 0.0_wp )  .AND.                            &
    21572155            surface_scalarflux == 9999999.9_wp )  THEN
    2158           message_string = 'wall_scalarflux additionally requires ' //     &
     2156          message_string = 'wall_scalarflux additionally requires ' //         &
    21592157                           'setting of surface_scalarflux'
    21602158          CALL message( 'check_parameters', 'PA0445', 1, 2, 0, 6, 0 )
     
    23262324    IF ( averaging_interval > dt_data_output_av )  THEN
    23272325       WRITE( message_string, * )  'averaging_interval = ',                    &
    2328              averaging_interval, ' must be <= dt_data_output = ', dt_data_output
     2326             averaging_interval, ' must be <= dt_data_output_av = ',           &
     2327             dt_data_output_av
    23292328       CALL message( 'check_parameters', 'PA0085', 1, 2, 0, 6, 0 )
    23302329    ENDIF
     
    25492548          CASE ( 'sa', '#sa' )
    25502549             IF ( .NOT. ocean )  THEN
    2551                 message_string = 'data_output_pr = ' // &
     2550                message_string = 'data_output_pr = ' //                        &
    25522551                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
    25532552                                 'lemented for ocean = .FALSE.'
     
    26222621          CASE ( 'q', '#q' )
    26232622             IF ( .NOT. humidity )  THEN
    2624                 message_string = 'data_output_pr = ' // &
     2623                message_string = 'data_output_pr = ' //                        &
    26252624                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
    26262625                                 'lemented for humidity = .FALSE.'
     
    28072806                message_string = 'data_output_pr = ' //                        &
    28082807                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
    2809                                  'lemented for cloud_physics = .FALSE. an&' // &
     2808                                 'lemented for cloud_physics = .FALSE. an' // &
    28102809                                 'd humidity = .FALSE.'
    28112810                CALL message( 'check_parameters', 'PA0095', 1, 2, 0, 6, 0 )
     
    28252824                message_string = 'data_output_pr = ' //                        &
    28262825                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
    2827                                  'lemented for cloud_physics = .FALSE. an&' // &
     2826                                 'lemented for cloud_physics = .FALSE. an' // &
    28282827                                 'd humidity = .FALSE.'
    28292828                CALL message( 'check_parameters', 'PA0095', 1, 2, 0, 6, 0 )
     
    28422841                message_string = 'data_output_pr = ' //                        &
    28432842                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
    2844                                  'lemented for cloud_physics = .FALSE. an&' // &
     2843                                 'lemented for cloud_physics = .FALSE. an' // &
    28452844                                 'd humidity = .FALSE.'
    28462845                CALL message( 'check_parameters', 'PA0095', 1, 2, 0, 6, 0 )
     
    28512850                message_string = 'data_output_pr = ' //                        &
    28522851                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
    2853                                  'lemented for cloud_physics = .FALSE. or' // &
    2854                                  '&cloud_droplets = .FALSE.'
     2852                                 'lemented for cloud_physics = .FALSE. and' // &
     2853                                 'cloud_droplets = .FALSE.'
    28552854                CALL message( 'check_parameters', 'PA0096', 1, 2, 0, 6, 0 )
    28562855             ELSE
     
    32493248             ELSEIF ( .NOT.  microphysics_morrison )  THEN
    32503249                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    3251                          'res = microphysics morrison '
     3250                         'res = morrison '
    32523251                CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
    32533252             ENDIF
     
    32693268             IF (  .NOT.  particle_advection )  THEN
    32703269                message_string = 'output of "' // TRIM( var ) // '" requir' // &
    3271                    'es a "particle_parameters"-NAMELIST in the parameter file (PARIN)'
     3270                   'es a "particle_parameters"-NAMELIST in the parameter ' //  &
     3271                   'file (PARIN)'
    32723272                CALL message( 'check_parameters', 'PA0104', 1, 2, 0, 6, 0 )
    32733273             ENDIF
     
    32813281                CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
    32823282             ELSEIF ( microphysics_sat_adjust )  THEN
    3283                 message_string = 'output of "' // TRIM( var ) // '" is ' //  &
     3283                message_string = 'output of "' // TRIM( var ) // '" is ' //    &
    32843284                         'not available for cloud_scheme = saturation_adjust'
    32853285                CALL message( 'check_parameters', 'PA0423', 1, 2, 0, 6, 0 )
     
    33783378             IF ( k == 0  .OR.  data_output(i)(ilen-2:ilen) /= '_xy' )  THEN
    33793379                message_string = 'illegal value for data_output: "' //         &
    3380                                  TRIM( var ) // '" & only 2d-horizontal ' //   &
     3380                                 TRIM( var ) // '" only 2d-horizontal ' //     &
    33813381                                 'cross sections are allowed for this value'
    33823382                CALL message( 'check_parameters', 'PA0111', 1, 2, 0, 6, 0 )
     
    34113411             ENDIF
    34123412
    3413           IF ( TRIM( var ) == 'ghf*'  .AND.  .NOT.  land_surface )  THEN
    3414              message_string = 'output of "' // TRIM( var ) // '" requi' //     &
    3415                               'res land_surface = .TRUE.'
    3416              CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
    3417           ENDIF
     3413             IF ( TRIM( var ) == 'ghf*'  .AND.  .NOT.  land_surface )  THEN
     3414                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
     3415                                 'res land_surface = .TRUE.'
     3416                CALL message( 'check_parameters', 'PA0404', 1, 2, 0, 6, 0 )
     3417             ENDIF
    34183418
    34193419             IF ( ( TRIM( var ) == 'r_a*' .OR.  TRIM( var ) == 'ghf*' )        &
     
    35583558       WRITE( message_string, * )  'output of averaged quantity "',            &
    35593559                                   TRIM( doav(1) ), '_av" requires to set a ', &
    3560                                    'non-zero & averaging interval'
     3560                                   'non-zero averaging interval'
    35613561       CALL message( 'check_parameters', 'PA0323', 1, 2, 0, 6, 0 )
    35623562    ENDIF
     
    35793579    section(:,2) = section_xz
    35803580    section(:,3) = section_yz
    3581 
    3582 !
    3583 !-- Upper plot limit for 2D vertical sections
    3584     IF ( z_max_do2d == -1.0_wp )  z_max_do2d = zu(nzt)
    3585     IF ( z_max_do2d < zu(nzb+1)  .OR.  z_max_do2d > zu(nzt) )  THEN
    3586        WRITE( message_string, * )  'z_max_do2d = ', z_max_do2d,                &
    3587                     ' must be >= ', zu(nzb+1), '(zu(nzb+1)) and <= ', zu(nzt), &
    3588                     ' (zu(nzt))'
    3589        CALL message( 'check_parameters', 'PA0116', 1, 2, 0, 6, 0 )
    3590     ENDIF
    35913581
    35923582!
     
    36473637                           '5 (parallel netCDF 4) and ' //                     &
    36483638                           '6 (parallel netCDF 4 Classic model) '//            &
    3649                            '&are currently not supported (not yet tested) ' // &
    3650                            'for masked data.&Using respective non-parallel' // &
     3639                           ' are currently not supported (not yet tested) ' // &
     3640                           'for masked data. Using respective non-parallel' // &
    36513641                           ' output for masked data.'
    36523642          CALL message( 'check_parameters', 'PA0383', 0, 0, 0, 6, 0 )
     
    36633653#else
    36643654       message_string = 'netCDF: netCDF4 format requested but no ' //          &
    3665                         'cpp-directive __netcdf4 given & switch '  //          &
     3655                        'cpp-directive __netcdf4 given  switch '  //           &
    36663656                        'back to 64-bit offset format'
    36673657       CALL message( 'check_parameters', 'PA0171', 0, 1, 0, 6, 0 )
     
    36743664#else
    36753665       message_string = 'netCDF: netCDF4 parallel output requested but no ' // &
    3676                         'cpp-directive __netcdf4_parallel given & switch '  // &
     3666                        'cpp-directive __netcdf4_parallel given, switch '   // &
    36773667                        'back to netCDF4 non-parallel output'
    36783668       CALL message( 'check_parameters', 'PA0099', 0, 1, 0, 6, 0 )
     
    38653855    IF ( disturbance_level_ind_t < disturbance_level_ind_b )  THEN
    38663856       WRITE( message_string, * )  'disturbance_level_ind_t = ',               &
    3867                 disturbance_level_ind_t, ' must be >= disturbance_level_ind_b = ', &
    3868                 disturbance_level_ind_b
     3857                disturbance_level_ind_t, ' must be >= ',                      &
     3858                'disturbance_level_ind_b = ', disturbance_level_ind_b
    38693859       CALL message( 'check_parameters', 'PA0130', 1, 2, 0, 6, 0 )
    38703860    ENDIF
     
    39843974    IF ( turbulent_inflow )  THEN
    39853975       IF ( recycling_width <= dx  .OR.  recycling_width >= nx * dx )  THEN
    3986           WRITE( message_string, * )  'illegal value for recycling_width:', &
    3987                                       ' ', recycling_width
     3976          WRITE( message_string, * )  'illegal value for recycling_width: ',  &
     3977                                      recycling_width
    39883978          CALL message( 'check_parameters', 'PA0134', 1, 2, 0, 6, 0 )
    39893979       ENDIF
     
    40874077       IF ( dp_level_b < zu(nzb)  .OR.  dp_level_b > zu(nzt) )  THEN
    40884078          WRITE( message_string, * )  'dp_level_b = ', dp_level_b, ' is out ', &
    4089                ' of range'
     4079               ' of range [zu(nzb), zu(nzt)]'
    40904080          CALL message( 'check_parameters', 'PA0151', 1, 2, 0, 6, 0 )
    40914081       ENDIF
    40924082       IF ( .NOT. ANY( dpdxy /= 0.0_wp ) )  THEN
    40934083          WRITE( message_string, * )  'dp_external is .TRUE. but dpdxy is ze', &
    4094                'ro, i.e. the external pressure gradient & will not be applied'
     4084               'ro, i.e. the external pressure gradient will not be applied'
    40954085          CALL message( 'check_parameters', 'PA0152', 0, 1, 0, 6, 0 )
    40964086       ENDIF
     
    41474137       IF ( particle_dvrpsize /= 'absw' )  THEN
    41484138          message_string = 'illegal value for parameter particle_dvrpsize:' // &
    4149                            ' ' // TRIM( particle_color)
     4139                           ' ' // TRIM( particle_dvrpsize)
    41504140          CALL message( 'check_parameters', 'PA0314', 1, 2, 0, 6, 0 )
    41514141       ELSE
     
    41834173!-- Check roughness length, which has to be smaller than dz/2
    41844174    IF ( ( constant_flux_layer .OR.  &
    4185            INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) &
     4175           INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 )       &
    41864176         .AND. roughness_length >= 0.5 * dz )  THEN
    41874177       message_string = 'roughness_length must be smaller than dz/2'
     
    42374227          IF ( dt_fixed )  THEN
    42384228             WRITE( message_string, '(A,F9.4,A)' )  'Output at every '  //     &
    4239                     'timestep is wanted (' // dt_do_name // ' = 0.0).&'//      &
    4240                     'Setting the output interval to the fixed timestep '//     &
    4241                     'dt = ', dt, 's.'
     4229                    'timestep is wanted (' // dt_do_name // ' = 0.0). '//      &
     4230                    'The output interval is set to the fixed timestep dt '//   &
     4231                    '= ', dt, 's.'
    42424232             CALL message( 'check_parameters', 'PA0060', 0, 0, 0, 6, 0 )
    42434233             dt_do = dt
  • palm/trunk/SOURCE/chemistry_model_mod.f90

    r3014 r3045  
    2727! -----------------
    2828! $Id$
     29! error messages revised
     30!
     31! 3014 2018-05-09 08:42:38Z maronga
    2932! Bugfix: nzb_do and nzt_do were not used for 3d data output
    3033!
     
    11341137                 TRIM( data_output_pr(var_count) ) // ' is not imp' // &
    11351138                          'lemented for air_chemistry = .FALSE.'
    1136 !                CALL message( 'check_parameters', 'PA0185', 1, 2, 0, 6, 0 )
     1139!                CALL message( 'check_parameters', 'PA0XXX', 1, 2, 0, 6, 0 )
    11371140          ELSE
    11381141             DO lsp = 1, nspec
  • palm/trunk/SOURCE/close_file.f90

    r2718 r3045  
    2525! -----------------
    2626! $Id$
     27! z_max_do2d removed and output case 108 disabled
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    8891
    8992    USE control_parameters,                                                    &
    90         ONLY:  max_masks, mid, nz_do3d, openfile, run_description_header,      &
    91                z_max_do2d
     93        ONLY:  max_masks, mid, nz_do3d, openfile, run_description_header
    9294               
    9395    USE grid_variables,                                                        &
     
    205207                ENDIF
    206208
    207              CASE ( 108 )
    208 
    209                 nc_stat = NF90_CLOSE( id_set_prt )
    210                 CALL netcdf_handle_error( 'close_file', 51 )
     209!
     210!--           Currently disabled
     211!             CASE ( 108 )
     212
     213!                nc_stat = NF90_CLOSE( id_set_prt )
     214!                CALL netcdf_handle_error( 'close_file', 51 )
    211215
    212216             CASE ( 109 )
  • palm/trunk/SOURCE/cpulog_mod.f90

    r3028 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 3028 2018-05-23 10:27:23Z knoop
    2730! renamed d3par to runtime_parameters
    2831!
     
    222225          log_event%place = place
    223226       ELSEIF ( log_event%place /= place )  THEN
    224           WRITE( message_string, * ) 'wrong argument & expected: ',            &
    225                             TRIM(log_event%place), '  given: ',  TRIM( place )
     227          WRITE( message_string, * ) 'wrong argument expected: ',              &
     228                            TRIM(log_event%place), ' given: ',  TRIM( place )
    226229          CALL message( 'cpu_log', 'PA0174', 1, 2, 0, 6, 0 )
    227230       ENDIF
     
    259262          IF ( ( mtime - log_event%mtime ) < 0.0  .AND.  first )  THEN
    260263             WRITE( message_string, * ) 'negative time interval occured',      &
    261                          ' &PE',myid,' L=PAUSE "',TRIM(log_event%place),       &
     264                         ' PE',myid,' L=PAUSE "',TRIM(log_event%place),        &
    262265                         '" new=', mtime,' last=',log_event%mtime
    263266             CALL message( 'cpu_log', 'PA0176', 0, 1, -1, 6, 0 )
     
    270273               first )  THEN
    271274             WRITE( message_string, * ) 'negative time interval occured',      &
    272                        ' &PE',myid,' L=STOP "',TRIM(log_event%place),'" new=', &
     275                       ' PE',myid,' L=STOP "',TRIM(log_event%place),'" new=', &
    273276                       mtime,' last=',log_event%mtime,' isum=',log_event%isum
    274277             CALL message( 'cpu_log', 'PA0177', 0, 1, -1, 6, 0 )
     
    280283          IF ( log_event%sum < 0.0  .AND.  first )  THEN
    281284             WRITE( message_string, * ) 'negative time interval occured',      &
    282                        ' &PE',myid,' L=STOP "',TRIM(log_event%place),'" sum=', &
     285                       ' PE',myid,' L=STOP "',TRIM(log_event%place),'" sum=', &
    283286                       log_event%sum,' mtime=',log_event%mtime
    284287             CALL message( 'cpu_log', 'PA0178', 0, 1, -1, 6, 0 )
  • palm/trunk/SOURCE/data_output_2d.f90

    r3014 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 3014 2018-05-09 08:42:38Z maronga
    2730! Added nzb_do and nzt_do for some modules for 2d output
    2831!
     
    462465            ntdim_2d_xy(av) )  THEN
    463466          WRITE ( message_string, * ) 'Output of xy cross-sections is not ',   &
    464                           'given at t=', simulated_time, '&because the',       &
     467                          'given at t=', simulated_time, ' because the',       &
    465468                          ' maximum number of output time levels is exceeded.'
    466469          CALL message( 'data_output_2d', 'PA0384', 0, 1, 0, 6, 0 )
     
    471474            ntdim_2d_xz(av) )  THEN
    472475          WRITE ( message_string, * ) 'Output of xz cross-sections is not ',   &
    473                           'given at t=', simulated_time, '&because the',       &
     476                          'given at t=', simulated_time, ' because the',       &
    474477                          ' maximum number of output time levels is exceeded.'
    475478          CALL message( 'data_output_2d', 'PA0385', 0, 1, 0, 6, 0 )
     
    480483            ntdim_2d_yz(av) )  THEN
    481484          WRITE ( message_string, * ) 'Output of yz cross-sections is not ',   &
    482                           'given at t=', simulated_time, '&because the',       &
     485                          'given at t=', simulated_time, ' because the',       &
    483486                          ' maximum number of output time levels is exceeded.'
    484487          CALL message( 'data_output_2d', 'PA0386', 0, 1, 0, 6, 0 )
  • palm/trunk/SOURCE/data_output_3d.f90

    r3014 r3045  
    2525! -----------------
    2626! $Id$
     27! Error message revised
     28!
     29! 3014 2018-05-09 08:42:38Z maronga
    2730! Added nzb_do and nzt_do for some modules for 3d output
    2831!
     
    309312       IF ( do3d_time_count(av) + 1 > ntdim_3d(av) )  THEN
    310313          WRITE ( message_string, * ) 'Output of 3d data is not given at t=',  &
    311                                       simulated_time, '&because the maximum ', &
     314                                      simulated_time, ' because the maximum ', &
    312315                                      'number of output time levels is ',      &
    313316                                      'exceeded.'
  • palm/trunk/SOURCE/data_output_dvrp.f90

    r2718 r3045  
    2525! -----------------
    2626! $Id$
     27! Code adjusted according to PALM coding standards
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    255258!--    Select the plot mode (in case of isosurface or slicer only if user has
    256259!--    defined a variable which shall be plotted; otherwise do nothing)
    257        IF ( ( mode_dvrp(m)(1:10) == 'isosurface'  .OR.                     &
     260       IF ( ( mode_dvrp(m)(1:10) == 'isosurface'  .OR.                         &
    258261                  mode_dvrp(m)(1:6)  == 'slicer'           )                   &
    259262                  .AND.  output_variable /= ' ' )  THEN
     
    362365                   ELSE                   
    363366                      message_string = 'if humidity = '    //                  &
    364                             '.FALSE. output of ' // TRIM( output_variable ) //   &
     367                            '.FALSE. output of ' // TRIM( output_variable ) // &
    365368                            'is not provided'
    366369                      CALL message( 'data_output_dvrp', 'PA0183',&
     
    378381                      ENDDO
    379382                   ELSE                     
    380                       message_string = 'if cloud_physics = .FALSE. and '       //    &
     383                      message_string = 'if cloud_physics = .FALSE. and ' //    &
    381384                                  'cloud_droplets = .FALSE. '
    382385                                  'output of ' // TRIM( output_variable) //    &
     
    561564          CALL cpu_log( log_point_s(31), 'dvrp_pathlines', 'start' )
    562565
    563           CALL DVRP_DATA( m-1, local_pfi, 4, nx_dvrp, ny_dvrp, nz_dvrp, &
     566          CALL DVRP_DATA( m-1, local_pfi, 4, nx_dvrp, ny_dvrp, nz_dvrp,        &
    564567                          cyclic_dvrp, cyclic_dvrp, cyclic_dvrp )
    565568          CALL DVRP_VISUALIZE( m-1, 20, dvrp_filecount )
  • palm/trunk/SOURCE/data_output_mask.f90

    r3030 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 3030 2018-05-23 14:37:00Z raasch
    2730! variable if renamed ivar
    2831!
     
    528531
    529532             IF ( .NOT. found )  THEN
    530                 WRITE ( message_string, * ) 'no output available for: ', &
     533                WRITE ( message_string, * ) 'no masked output available for: ',&
    531534                                            TRIM( domask(mid,av,ivar) )
    532535                CALL message( 'data_output_mask', 'PA0327', 0, 0, 0, 6, 0 )
  • palm/trunk/SOURCE/fft_xy_mod.f90

    r2718 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    723726          ENDIF
    724727
    725 #else
    726           message_string = 'no system-specific fft-call available'
    727           CALL message( 'fft_x', 'PA0188', 1, 2, 0, 6, 0 )
    728 #endif
    729 
    730        ELSE
    731 
    732           message_string = 'fft method "' // TRIM( fft_method) // &
    733                            '" not available'
    734           CALL message( 'fft_x', 'PA0189', 1, 2, 0, 6, 0 )
     728#endif
    735729
    736730       ENDIF
     
    945939
    946940          ENDIF
    947 #else
    948           message_string = 'no system-specific fft-call available'
    949           CALL message( 'fft_x_1d', 'PA0188', 1, 2, 0, 6, 0 )
    950 #endif
    951        ELSE
    952           message_string = 'fft method "' // TRIM( fft_method) // &
    953                            '" not available'
    954           CALL message( 'fft_x_1d', 'PA0189', 1, 2, 0, 6, 0 )
     941#endif
    955942
    956943       ENDIF
     
    12881275
    12891276          ENDIF
    1290 #else
    1291           message_string = 'no system-specific fft-call available'
    1292           CALL message( 'fft_y', 'PA0188', 1, 2, 0, 6, 0 )
    1293 #endif
    1294 
    1295        ELSE
    1296 
    1297           message_string = 'fft method "' // TRIM( fft_method) // &
    1298                            '" not available'
    1299           CALL message( 'fft_y', 'PA0189', 1, 2, 0, 6, 0 )
     1277#endif
    13001278
    13011279       ENDIF
     
    15121490
    15131491          ENDIF
    1514 #else
    1515           message_string = 'no system-specific fft-call available'
    1516           CALL message( 'fft_y_1d', 'PA0188', 1, 2, 0, 6, 0 )
    1517 
    1518 #endif
    1519 
    1520        ELSE
    1521 
    1522           message_string = 'fft method "' // TRIM( fft_method) // &
    1523                            '" not available'
    1524           CALL message( 'fft_y_1d', 'PA0189', 1, 2, 0, 6, 0 )
     1492#endif
    15251493
    15261494       ENDIF
     
    16531621
    16541622          DEALLOCATE( work )
    1655 #else
    1656           message_string = 'no system-specific fft-call available'
    1657           CALL message( 'fft_x_m', 'PA0188', 1, 2, 0, 6, 0 )
    1658 #endif
    1659 
    1660        ELSE
    1661 
    1662           message_string = 'fft method "' // TRIM( fft_method) // &
    1663                            '" not available'
    1664           CALL message( 'fft_x_m', 'PA0189', 1, 2, 0, 6, 0 )
     1623#endif
    16651624
    16661625       ENDIF
     
    17951754
    17961755          DEALLOCATE( work )
    1797 #else
    1798           message_string = 'no system-specific fft-call available'
    1799           CALL message( 'fft_y_m', 'PA0188', 1, 2, 0, 6, 0 )
    1800 #endif
    1801 
    1802        ELSE
    1803          
    1804           message_string = 'fft method "' // TRIM( fft_method) // &
    1805                            '" not available'
    1806           CALL message( 'fft_y_m', 'PA0189', 1, 2, 0, 6, 0 )
     1756#endif
    18071757
    18081758       ENDIF
  • palm/trunk/SOURCE/header.f90

    r2967 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 2967 2018-04-13 11:22:08Z raasch
    2730! bugfix: missing parallel cpp-directives added
    2831!
     
    550553    ELSEIF ( INDEX( initializing_actions, 'by_user' ) /=0 )  THEN
    551554       run_classification = 'run initialized by user'
     555    ELSEIF ( INDEX( initializing_actions, 'initialize_vortex' ) /=0 )  THEN
     556       run_classification = 'run additionally initialized by a Rankine-vortex'
     557    ELSEIF ( INDEX( initializing_actions, 'initialize_ptanom' ) /=0 )  THEN
     558       run_classification = 'run additionally initialized by temperature anomaly'
    552559    ELSE
    553560       message_string = ' unknown action(s): ' // TRIM( initializing_actions )
  • palm/trunk/SOURCE/init_3d_model.f90

    r3042 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 3042 2018-05-25 10:44:37Z schwenkel
    2730! Changed the name specific humidity to mixing ratio
    2831!
     
    26112614       WRITE( message_string, * ) 'number of time series quantities exceeds',  &
    26122615                                  ' its maximum of dots_max = ', dots_max,     &
    2613                                   ' &Please increase dots_max in modules.f90.'
     2616                                  ' Please increase dots_max in modules.f90.'
    26142617       CALL message( 'init_3d_model', 'PA0194', 1, 2, 0, 6, 0 )   
    26152618    ENDIF
  • palm/trunk/SOURCE/init_grid.f90

    r2968 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 2968 2018-04-13 11:52:24Z suehring
    2730! Bugfix in initialization in case of elevated model surface. Introduce
    2831! index for minimum topography-top.
     
    651654       IF ( nzb_local_min < 0  .OR.  nzb_local_max  > nz + 1 )  THEN
    652655          WRITE( message_string, * ) 'nzb_local values are outside the',       &
    653                                 'model domain',                                &
    654                                 '&MINVAL( nzb_local ) = ', nzb_local_min,      &
    655                                 '&MAXVAL( nzb_local ) = ', nzb_local_max
     656                                ' model domain',                               &
     657                                ' MINVAL( nzb_local ) = ', nzb_local_min,      &
     658                                ' MAXVAL( nzb_local ) = ', nzb_local_max
    656659          CALL message( 'init_grid', 'PA0210', 1, 2, 0, 6, 0 )
    657660       ENDIF
     
    14971500               ( bys < 1 ) .OR. ( byn > ny-1 ) .OR. ( byn < bys+3 ) )  THEN
    14981501             WRITE( message_string, * ) 'inconsistent building parameters:',   &
    1499                                       '& bxl=', bxl, 'bxr=', bxr, 'bys=', bys, &
     1502                                      ' bxl=', bxl, 'bxr=', bxr, 'bys=', bys, &
    15001503                                      'byn=', byn, 'nx=', nx, 'ny=', ny
    15011504             CALL message( 'init_grid', 'PA0203', 1, 2, 0, 6, 0 )
     
    15891592                  ( ch < 3 ) )  THEN
    15901593                WRITE( message_string, * ) 'inconsistent canyon parameters:',  &
    1591                                            '&cxl=', cxl, 'cxr=', cxr,          &
    1592                                            'cwx=', cwx,                        &
    1593                                            'ch=', ch, 'nx=', nx, 'ny=', ny
     1594                                           ' cxl=', cxl, ' cxr=', cxr,         &
     1595                                           ' cwx=', cwx,                       &
     1596                                           ' ch=', ch, ' nx=', nx, ' ny=', ny
    15941597                CALL message( 'init_grid', 'PA0205', 1, 2, 0, 6, 0 )
    15951598             ENDIF
     
    15981601                  ( ch < 3 ) )  THEN
    15991602                WRITE( message_string, * ) 'inconsistent canyon parameters:',  &
    1600                                            '&cys=', cys, 'cyn=', cyn,          &
    1601                                            'cwy=', cwy,                        &
    1602                                            'ch=', ch, 'nx=', nx, 'ny=', ny
     1603                                           ' cys=', cys, ' cyn=', cyn,         &
     1604                                           ' cwy=', cwy,                       &
     1605                                           ' ch=', ch, ' nx=', nx, ' ny=', ny
    16031606                CALL message( 'init_grid', 'PA0206', 1, 2, 0, 6, 0 )
    16041607             ENDIF
     
    16071610               canyon_width_y /= 9999999.9_wp )  THEN
    16081611             message_string = 'inconsistent canyon parameters:' //             &   
    1609                               '&street canyon can only be oriented' //         &
    1610                               '&either in x- or in y-direction'
     1612                              ' street canyon can only be oriented' //         &
     1613                              ' either in x- or in y-direction'
    16111614             CALL message( 'init_grid', 'PA0207', 1, 2, 0, 6, 0 )
    16121615          ENDIF
     
    18951898               TRIM( topography ) /= 'read_from_file')  THEN
    18961899!--          The default value is not applicable here, because it is only valid
    1897 !--          for the two standard cases 'single_building' and 'read_from_file'
     1900!--          for the four standard cases 'single_building',
     1901!--          'single_street_canyon', 'tunnel' and 'read_from_file'
    18981902!--          defined in init_grid.
    18991903             WRITE( message_string, * )                                        &
    19001904               'The value for "topography_grid_convention" ',                  &
    1901                'is not set. Its default value is & only valid for ',           &
    1902                '"topography" = ''single_building'', ',                         &
    1903                '''single_street_canyon'' & or ''read_from_file''.',            &
    1904                ' & Choose ''cell_edge'' or ''cell_center''.'
     1905               'is not set. Its default value is only valid for ',             &
     1906               '"topography" = ''single_building'', ''tunnel'' ',              &
     1907               '''single_street_canyon'' or ''read_from_file''.',              &
     1908               ' Choose ''cell_edge'' or ''cell_center''.'
    19051909             CALL message( 'init_grid', 'PA0239', 1, 2, 0, 6, 0 )
    19061910          ELSE
     
    19191923          WRITE( message_string, * )                                           &
    19201924            'The value for "topography_grid_convention" is ',                  &
    1921             'not recognized. & Choose ''cell_edge'' or ''cell_center''.'
     1925            'not recognized. Choose ''cell_edge'' or ''cell_center''.'
    19221926          CALL message( 'init_grid', 'PA0240', 1, 2, 0, 6, 0 )
    19231927       ENDIF
  • palm/trunk/SOURCE/init_masks.f90

    r2718 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    192195   
    193196!
    194 !-- Parallel mask output not yet tested
     197!-- Parallel mask output not yet supported. In check_parameters data format
     198!-- is restricted and is switched back to non-parallel output. Therefore the
     199!-- following error can not occur at the moment.
    195200    IF ( netcdf_data_format > 4 )  THEN
    196        message_string = 'netCDF file formats '//                         &
    197                         '3 (netCDF 4) and 4 (netCDF 4 Classic model)'//  &
    198                         '&are currently not supported (not yet tested).'
     201       message_string = 'netCDF file formats '//                               &
     202                        '5 and 6 (with parallel I/O support)'//                &
     203                        ' are currently not supported.'
    199204       CALL message( 'init_masks', 'PA0328', 1, 2, 0, 6, 0 )
    200205    ENDIF
     
    258263          DO  WHILE ( do_mask_user(mid,j) /= ' '  .AND.  j <= 100 )
    259264             IF ( i > 100 )  THEN
    260                 WRITE ( message_string, * ) 'number of output quantitities ', &
    261                      'given by data_output_mask and data_output_mask_user ', &
     265                WRITE ( message_string, * ) 'number of output quantitities ',  &
     266                     'given by data_output_mask and data_output_mask_user ',   &
    262267                     'exceeds the limit of 100'
    263268                CALL message( 'init_masks', 'PA0329', 1, 2, 0, 6, 0 )
     
    290295             CASE ( 'e' )
    291296                IF ( constant_diffusion )  THEN
    292                    WRITE ( message_string, * ) 'output of "', TRIM( var ), &
     297                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
    293298                        '" requires constant_diffusion = .FALSE.'
    294299                   CALL message( 'init_masks', 'PA0103', 1, 2, 0, 6, 0 )
     
    298303             CASE ( 'lpt' )
    299304                IF ( .NOT. cloud_physics )  THEN
    300                    WRITE ( message_string, * ) 'output of "', TRIM( var ), &
     305                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
    301306                        '" requires cloud_physics = .TRUE.'
    302307                   CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
     
    306311             CASE ( 'nc' )
    307312                IF ( .NOT. cloud_physics )  THEN
    308                    WRITE ( message_string, * ) 'output of "', TRIM( var ), &
     313                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
    309314                        '" requires cloud_physics = .TRUE.'
    310315                   CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
    311316                 ELSEIF ( .NOT. microphysics_morrison ) THEN
    312                    message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    313                          'res  = morrison'
     317                   message_string = 'output of "' // TRIM( var ) // '" ' //    &
     318                         'requires  = morrison'
    314319                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
    315320                ENDIF
     
    318323             CASE ( 'nr' )
    319324                IF ( .NOT. cloud_physics )  THEN
    320                    WRITE ( message_string, * ) 'output of "', TRIM( var ), &
     325                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
    321326                        '" requires cloud_physics = .TRUE.'
    322327                   CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
    323328                 ELSEIF ( .NOT. microphysics_seifert ) THEN
    324                    message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    325                          'res cloud_scheme = seifert_beheng'
     329                   message_string = 'output of "' // TRIM( var ) // '"' //     &
     330                         'requires cloud_scheme = seifert_beheng'
    326331                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
    327332                ENDIF
     
    330335             CASE ( 'pc', 'pr' )
    331336                IF ( .NOT. particle_advection )  THEN
    332                    WRITE ( message_string, * ) 'output of "', TRIM( var ), &
    333                         '" requires a "particles_par"-NAMELIST in the ', &
     337                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
     338                        '" requires a "particles_par"-NAMELIST in the ',       &
    334339                        'parameter file (PARIN)'
    335340                   CALL message( 'init_masks', 'PA0104', 1, 2, 0, 6, 0 )
     
    340345             CASE ( 'q', 'vpt' )
    341346                IF ( .NOT. humidity )  THEN
    342                    WRITE ( message_string, * ) 'output of "', TRIM( var ), &
     347                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
    343348                        '" requires humidity = .TRUE.'
    344349                   CALL message( 'init_masks', 'PA0105', 1, 2, 0, 6, 0 )
     
    349354             CASE ( 'qc' )
    350355                IF ( .NOT. cloud_physics )  THEN
    351                    message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    352                             'res cloud_physics = .TRUE.'
     356                   message_string = 'output of "' // TRIM( var ) // '"' //     &
     357                            'requires cloud_physics = .TRUE.'
    353358                   CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
    354359                ENDIF
     
    357362             CASE ( 'ql' )
    358363                IF ( .NOT. ( cloud_physics  .OR.  cloud_droplets ) )  THEN
    359                    WRITE ( message_string, * ) 'output of "', TRIM( var ), &
     364                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
    360365                        '" requires cloud_physics = .TRUE. or cloud_droplets', &
    361366                        ' = .TRUE.'
    362                    CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
     367                   CALL message( 'init_masks', 'PA0106', 1, 2, 0, 6, 0 )
    363368                ENDIF
    364369                unit = 'kg/kg'
     
    366371             CASE ( 'ql_c', 'ql_v', 'ql_vp' )
    367372                IF ( .NOT. cloud_droplets )  THEN
    368                    WRITE ( message_string, * ) 'output of "', TRIM( var ), &
     373                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
    369374                        '" requires cloud_droplets = .TRUE.'
    370375                   CALL message( 'init_masks', 'PA0107', 1, 2, 0, 6, 0 )
     
    376381             CASE ( 'qv' )
    377382                IF ( .NOT. cloud_physics )  THEN
    378                    WRITE ( message_string, * ) 'output of "', TRIM( var ), &
     383                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
    379384                        '" requires cloud_physics = .TRUE.'
    380385                   CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
     
    384389             CASE ( 'qr' )
    385390                IF ( .NOT. cloud_physics )  THEN
    386                    message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    387                             'res cloud_physics = .TRUE.'
     391                   message_string = 'output of "' // TRIM( var ) // '" ' //    &
     392                            'requires cloud_physics = .TRUE.'
    388393                   CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
    389394                ELSEIF ( .NOT. microphysics_seifert ) THEN
    390                    message_string = 'output of "' // TRIM( var ) // '" requi' //  &
    391                             'res cloud_scheme = seifert_beheng'
     395                   message_string = 'output of "' // TRIM( var ) // '" ' //    &
     396                            'requires cloud_scheme = seifert_beheng'
    392397                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
    393398                ENDIF
     
    396401             CASE ( 'rho_ocean' )
    397402                IF ( .NOT. ocean )  THEN
    398                    WRITE ( message_string, * ) 'output of "', TRIM( var ), &
     403                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
    399404                        '" requires ocean = .TRUE.'
    400405                   CALL message( 'init_masks', 'PA0109', 1, 2, 0, 6, 0 )
     
    404409             CASE ( 's' )
    405410                IF ( .NOT. passive_scalar )  THEN
    406                    WRITE ( message_string, * ) 'output of "', TRIM( var ), &
     411                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
    407412                        '" requires passive_scalar = .TRUE.'
    408413                   CALL message( 'init_masks', 'PA0110', 1, 2, 0, 6, 0 )
     
    412417             CASE ( 'sa' )
    413418                IF ( .NOT. ocean )  THEN
    414                    WRITE ( message_string, * ) 'output of "', TRIM( var ), &
     419                   WRITE ( message_string, * ) 'output of "', TRIM( var ),     &
    415420                        '" requires ocean = .TRUE.'
    416421                   CALL message( 'init_masks', 'PA0109', 1, 2, 0, 6, 0 )
     
    419424
    420425             CASE ( 'u*', 't*', 'lwp*', 'pra*', 'prr*', 'z0*', 'z0h*' )
    421                 WRITE ( message_string, * ) 'illegal value for data_',&
    422                      'output: "', TRIM( var ), '" is only allowed', &
     426                WRITE ( message_string, * ) 'illegal value for data_',         &
     427                     'output: "', TRIM( var ), '" is only allowed',            &
    423428                     'for horizontal cross section'
    424429                CALL message( 'init_masks', 'PA0111', 1, 2, 0, 6, 0 )
     
    437442                IF ( unit == 'illegal' )  THEN
    438443                   IF ( do_mask_user(mid,1) /= ' ' )  THEN
    439                       WRITE ( message_string, * ) 'illegal value for data_',&
    440                            'output or data_output_user: "', &
     444                      WRITE ( message_string, * ) 'illegal value for data_',   &
     445                           'output_masks or data_output_masks_user: "',        &
    441446                           TRIM( do_mask(mid,i) ), '"'
    442                       CALL message( 'init_masks', 'PA0114', 1, 2, 0, 6, 0 )
     447                      CALL message( 'init_masks', 'PA0018', 1, 2, 0, 6, 0 )
    443448                   ELSE
    444                       WRITE ( message_string, * ) 'illegal value for data_',&
    445                            'output: "', TRIM( do_mask(mid,i) ), '"'
     449                      WRITE ( message_string, * ) 'illegal value for data_',   &
     450                           ' output_masks : "', TRIM( do_mask(mid,i) ), '"'
    446451                      CALL message( 'init_masks', 'PA0330', 1, 2, 0, 6, 0 )
    447452                   ENDIF
     
    643648             ENDIF
    644649             IF ( m > (nxyz+1) )  THEN
    645                 WRITE ( message_string, '(I3,A,I3,A,I1,3A,I3)' )  &
    646                      m,' in mask ',mid,' along dimension ',dim,  &
     650                WRITE ( message_string, '(I3,A,I3,A,I1,3A,I3)' )               &
     651                     m,' in mask ',mid,' along dimension ',dim,                &
    647652                     ' exceeds (',nxyz_string,'+1) = ',nxyz+1
    648653                CALL message( 'init_masks', 'PA0331', 1, 2, 0, 6, 0 )
     
    681686                  > (nxyz+1) * dxyz / mask_scale(dim) )  THEN
    682687                WRITE ( message_string, '(2(A,I3,A,I1,A,F9.3),5A,I1,A,F9.3)' ) &
    683                      'mask_loop(',mid,',',dim,',1)=',mask_loop(mid,dim,1), &
     688                     'mask_loop(',mid,',',dim,',1)=',mask_loop(mid,dim,1),     &
    684689                     ' and/or mask_loop(',mid,',',dim,',2)=', &
    685                      mask_loop(mid,dim,2),'&exceed (', &
    686                      nxyz_string,'+1)*',dxyz_string,'/mask_scale(',dim,')=', &
     690                     mask_loop(mid,dim,2),' exceed (', &
     691                     nxyz_string,'+1)*',dxyz_string,'/mask_scale(',dim,')=',   &
    687692                     (nxyz+1)*dxyz/mask_scale(dim)
    688693                CALL message( 'init_masks', 'PA0332', 1, 2, 0, 6, 0 )
    689694             ENDIF
    690              loop_begin  = NINT( mask_loop(mid,dim,1) * mask_scale(dim) &
     695             loop_begin  = NINT( mask_loop(mid,dim,1) * mask_scale(dim)        &
    691696                  * ddxyz - 0.5_wp )
    692              loop_end    = NINT( mask_loop(mid,dim,2) * mask_scale(dim) &
     697             loop_end    = NINT( mask_loop(mid,dim,2) * mask_scale(dim)        &
    693698                  * ddxyz - 0.5_wp )
    694              loop_stride = NINT( mask_loop(mid,dim,3) * mask_scale(dim) &
     699             loop_stride = NINT( mask_loop(mid,dim,3) * mask_scale(dim)        &
    695700                  * ddxyz )
    696701             IF ( loop_begin == -1 )  loop_begin = 0  ! avoid negative values
     
    702707             IF ( MAXVAL( mask_loop(mid,dim,1:2) )  &
    703708                  > zu(nz+1) / mask_scale(dim) )  THEN
    704                 WRITE ( message_string, '(2(A,I3,A,I1,A,F9.3),A,I1,A,F9.3)' ) &
    705                      'mask_loop(',mid,',',dim,',1)=',mask_loop(mid,dim,1), &
     709                WRITE ( message_string, '(2(A,I3,A,I1,A,F9.3),A,I1,A,F9.3)' )  &
     710                     'mask_loop(',mid,',',dim,',1)=',mask_loop(mid,dim,1),     &
    706711                     ' and/or mask_loop(',mid,',',dim,',2)=', &
    707                      mask_loop(mid,dim,2),'&exceed zu(nz+1)/mask_scale(',dim, &
     712                     mask_loop(mid,dim,2),' exceed zu(nz+1)/mask_scale(',dim, &
    708713                     ')=',zu(nz+1)/mask_scale(dim)
    709714                CALL message( 'init_masks', 'PA0333', 1, 2, 0, 6, 0 )
     
    724729             IF ( mask_loop(mid,dim,2) * mask_scale(dim) > dz_stretch_level )  &
    725730                  THEN
    726                 WRITE ( message_string, '(A,I3,A,I1,A,F9.3,A,F8.2,3A)' ) &
    727                      'mask_loop(',mid,',',dim,',2)=', mask_loop(mid,dim,2),&
    728                      ' exceeds dz_stretch_level=',dz_stretch_level, &
    729                      '.&Vertical mask locations will not ', &
    730                      'match the desired heights&within the stretching ', &
     731                WRITE ( message_string, '(A,I3,A,I1,A,F9.3,A,F8.2,3A)' )       &
     732                     'mask_loop(',mid,',',dim,',2)=', mask_loop(mid,dim,2),    &
     733                     ' exceeds dz_stretch_level=',dz_stretch_level,            &
     734                     '. Vertical mask locations will not ',                    &
     735                     'match the desired heights within the stretching ',      &
    731736                     'region.'
    732737                CALL message( 'init_masks', 'PA0334', 0, 1, 0, 6, 0 )
  • palm/trunk/SOURCE/init_pegrid.f90

    r2938 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 2938 2018-03-27 15:52:42Z suehring
    2730! - No checks for domain decomposition in case of turbulence generator
    2831!  (is done in stg module)
     
    341344!--    must be equal to the number of PEs available to the job
    342345       IF ( ( npex * npey ) /= numprocs )  THEN
    343           WRITE( message_string, * ) 'number of PEs of the prescribed ',   &
    344               'topology (', npex*npey,') does not match & the number of ', &
     346          WRITE( message_string, * ) 'number of PEs of the prescribed ',       &
     347              'topology (', npex*npey,') does not match the number of ',      &
    345348              'PEs available to the job (', numprocs, ')'
    346349          CALL message( 'init_pegrid', 'PA0221', 1, 2, 0, 6, 0 )
     
    353356!--    If the processor topology is prescribed by the user, the number of
    354357!--    PEs must be given in both directions
    355        message_string = 'if the processor topology is prescribed by th' //  &
    356                 'e user& both values of "npex" and "npey" must be given' // &
     358       message_string = 'if the processor topology is prescribed by th' //     &
     359                'e user both values of "npex" and "npey" must be given' //    &
    357360                ' in the &NAMELIST-parameter file'
    358361       CALL message( 'init_pegrid', 'PA0222', 1, 2, 0, 6, 0 )
     
    450453!
    451454!-- Calculate array bounds along x-direction for every PE.
    452     ALLOCATE( nxlf(0:pdims(1)-1), nxrf(0:pdims(1)-1), nynf(0:pdims(2)-1), &
     455    ALLOCATE( nxlf(0:pdims(1)-1), nxrf(0:pdims(1)-1), nynf(0:pdims(2)-1),      &
    453456              nysf(0:pdims(2)-1) )
    454457
    455458    IF ( MOD( nx+1 , pdims(1) ) /= 0 )  THEN
    456        WRITE( message_string, * ) 'x-direction: gridpoint number (',nx+1,') ',&
    457                                'is not an& integral divisor of the number ',  &
    458                                'processors (', pdims(1),')'
     459       WRITE( message_string, * ) 'x-direction: gridpoint number (',nx+1,') ', &
     460                               'is not an integral divisor of the number ',    &
     461                               'of processors (', pdims(1),')'
    459462       CALL message( 'init_pegrid', 'PA0225', 1, 2, 0, 6, 0 )
    460463    ELSE
    461464       nnx  = ( nx + 1 ) / pdims(1)
    462        IF ( nnx*pdims(1) - ( nx + 1) > nnx )  THEN
    463           WRITE( message_string, * ) 'x-direction: nx does not match the',    &
    464                        'requirements given by the number of PEs &used',       &
    465                        '& please use nx = ', nx - ( pdims(1) - ( nnx*pdims(1) &
    466                                       - ( nx + 1 ) ) ), ' instead of nx =', nx
    467           CALL message( 'init_pegrid', 'PA0226', 1, 2, 0, 6, 0 )
    468        ENDIF
    469465    ENDIF   
    470466
     
    480476    IF ( MOD( ny+1 , pdims(2) ) /= 0 )  THEN
    481477       WRITE( message_string, * ) 'y-direction: gridpoint number (',ny+1,') ', &
    482                            'is not an& integral divisor of the number of',     &
     478                           'is not an integral divisor of the number of',      &
    483479                           'processors (', pdims(2),')'
    484480       CALL message( 'init_pegrid', 'PA0227', 1, 2, 0, 6, 0 )
    485481    ELSE
    486482       nny  = ( ny + 1 ) / pdims(2)
    487        IF ( nny*pdims(2) - ( ny + 1) > nny )  THEN
    488           WRITE( message_string, * ) 'y-direction: ny does not match the',    &
    489                        'requirements given by the number of PEs &used ',      &
    490                        '& please use ny = ', ny - ( pdims(2) - ( nnx*pdims(2) &
    491                                      - ( ny + 1 ) ) ), ' instead of ny =', ny
    492           CALL message( 'init_pegrid', 'PA0228', 1, 2, 0, 6, 0 )
    493        ENDIF
    494483    ENDIF   
    495484
     
    533522          IF ( MOD( nz , pdims(1) ) /= 0 )  THEN
    534523             WRITE( message_string, * ) 'transposition z --> x:',              &
    535                        '&nz=',nz,' is not an integral divisior of pdims(1)=',  &
     524                       ' nz=',nz,' is not an integral divisior of pdims(1)=',  &
    536525                                                                   pdims(1)
    537526             CALL message( 'init_pegrid', 'PA0230', 1, 2, 0, 6, 0 )
     
    555544       IF ( MOD( nx+1 , pdims(2) ) /= 0 )  THEN
    556545          WRITE( message_string, * ) 'transposition x --> y:',                 &
    557                             '&nx+1=',nx+1,' is not an integral divisor of ',   &
     546                            ' nx+1=',nx+1,' is not an integral divisor of ',   &
    558547                            'pdims(2)=',pdims(2)
    559548          CALL message( 'init_pegrid', 'PA0231', 1, 2, 0, 6, 0 )
     
    584573          IF ( MOD( ny+1 , pdims(1) ) /= 0 )  THEN
    585574             WRITE( message_string, * ) 'transposition y --> z:',              &
    586                                '& ny+1=',ny+1,' is not an integral divisor of',&
     575                               ' ny+1=',ny+1,' is not an integral divisor of', &
    587576                               ' pdims(1)=',pdims(1)
    588577             CALL message( 'init_pegrid', 'PA0232', 1, 2, 0, 6, 0 )
     
    595584          IF ( MOD( ny+1 , pdims(1) ) /= 0 )  THEN
    596585             WRITE( message_string, * ) 'transposition x --> y:',              &
    597                                '& ny+1=',ny+1,' is not an integral divisor of',&
     586                               ' ny+1=',ny+1,' is not an integral divisor of', &
    598587                               ' pdims(1)=',pdims(1)
    599588             CALL message( 'init_pegrid', 'PA0233', 1, 2, 0, 6, 0 )
     
    609598       IF ( MOD( nz, pdims(2) ) /= 0 )  THEN
    610599          WRITE( message_string, * ) 'direct transposition z --> y (needed ',  &
    611                     'for spectra):& nz=',nz,' is not an integral divisor of ', &
     600                    'for spectra): nz=',nz,' is not an integral divisor of ', &
    612601                    'pdims(2)=',pdims(2)
    613602          CALL message( 'init_pegrid', 'PA0234', 1, 2, 0, 6, 0 )
     
    974963!
    975964!--          Check pre-defined value and reset to default, if neccessary
    976              IF ( mg_switch_to_pe0_level < mg_switch_to_pe0_level_l  .OR.  &
     965             IF ( mg_switch_to_pe0_level < mg_switch_to_pe0_level_l  .OR.      &
    977966                  mg_switch_to_pe0_level >= maximum_grid_level_l )  THEN
    978                 message_string = 'mg_switch_to_pe0_level ' // &
     967                message_string = 'mg_switch_to_pe0_level ' //                  &
    979968                                 'out of range and reset to 0'
    980969                CALL message( 'init_pegrid', 'PA0235', 0, 1, 0, 6, 0 )
     
    10381027             subdomain_size = ( nxr - nxl + 2 * nbgp + 1 ) * &
    10391028                              ( nyn - nys + 2 * nbgp + 1 ) * ( nzt - nzb + 2 )
    1040              gathered_size  = ( nxr_l - nxl_l + 3 ) * ( nyn_l - nys_l + 3 ) * &
     1029             gathered_size  = ( nxr_l - nxl_l + 3 ) * ( nyn_l - nys_l + 3 ) *  &
    10411030                              ( nzt_l - nzb + 2 )
    10421031
    10431032#else
    1044              message_string = 'multigrid gather/scatter impossible ' // &
     1033             message_string = 'multigrid gather/scatter impossible ' //        &
    10451034                          'in non parallel mode'
    10461035             CALL message( 'init_pegrid', 'PA0237', 1, 2, 0, 6, 0 )
     
    10631052
    10641053!
    1065 !--    Temporary problem: Currently calculation of maxerror iin routine poismg crashes
     1054!--    Temporary problem: Currently calculation of maxerror in routine poismg crashes
    10661055!--    if grid data are collected on PE0 already on the finest grid level.
    10671056!--    To be solved later.
  • palm/trunk/SOURCE/land_surface_model_mod.f90

    r3026 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised and added
     28!
     29! 3026 2018-05-22 10:30:53Z schwenkel
    2730! Changed the name specific humidity to mixing ratio, since we are computing
    2831! mixing ratios.
     
    431434    USE radiation_model_mod,                                                   &
    432435        ONLY:  albedo, albedo_type, emissivity, force_radiation_call,          &
    433                radiation_scheme, unscheduled_radiation_calls
     436               radiation, radiation_scheme, unscheduled_radiation_calls
    434437       
    435438    USE statistics,                                                            &
     
    12741277         TRIM( surface_type ) /= 'water'       .AND.                           &
    12751278         TRIM( surface_type ) /= 'netcdf' )  THEN 
    1276        message_string = 'unknown surface type surface_type = "' //             &
     1279       message_string = 'unknown surface type: surface_type = "' //            &
    12771280                        TRIM( surface_type ) // '"'
    12781281       CALL message( 'lsm_check_parameters', 'PA0019', 1, 2, 0, 6, 0 )
     
    12951298       CALL message( 'lsm_check_parameters', 'PA0400', 1, 2, 0, 6, 0 )
    12961299    ENDIF
     1300   
     1301    IF (  .NOT.  radiation )  THEN
     1302       message_string = 'lsm requires '//                                      &
     1303                        'the radiation model to be switched on'
     1304       CALL message( 'lsm_check_parameters', 'PA0400', 1, 2, 0, 6, 0 )
     1305    ENDIF
    12971306
    12981307    IF ( TRIM( surface_type ) == 'vegetation' )  THEN
     
    13801389                                     ' is not allowed in combination with ',   &
    13811390                                     'most_method = ', most_method
    1382           CALL message( 'lsm_check_parameters', 'PA0417', 1, 2, 0, 6, 0 )
     1391          CALL message( 'lsm_check_parameters', 'PA0414', 1, 2, 0, 6, 0 )
    13831392       ENDIF
    13841393
     
    14261435             CALL message( 'lsm_check_parameters', 'PA0352', 1, 2, 0, 6, 0 )
    14271436          ENDIF
    1428 
     1437         
    14291438          IF ( z0h_pavement == 9999999.9_wp )  THEN
    14301439             message_string = 'pavement_type = 0 (user_defined)'//             &
     
    16301639    DO  k = nzb_soil, nzt_soil-1
    16311640       dz_soil_center(k) = zs(k+1) - zs(k)
    1632        IF ( dz_soil_center(k) == 0.0_wp )  THEN
     1641       IF ( dz_soil_center(k) <= 0.0_wp )  THEN
    16331642          message_string = 'invalid soil layer configuration found ' //        &
    1634                            '(dz_soil_center(k) = 0.0)'
     1643                           '(dz_soil_center(k) <= 0.0)'
    16351644          CALL message( 'lsm_rrd_local', 'PA0140', 1, 2, 0, 6, 0 )
    16361645       ENDIF 
  • palm/trunk/SOURCE/large_scale_forcing_nudging_mod.f90

    r3026 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 3026 2018-05-22 10:30:53Z schwenkel
    2730! Changed the name specific humidity to mixing ratio, since we are computing
    2831! mixing ratios.
     
    732735!--    Check nudging and large scale forcing from external file
    733736       IF ( nudging  .AND.  (  .NOT.  large_scale_forcing ) )  THEN
    734           message_string = 'Nudging requires large_scale_forcing = .T.. &'//   &
    735                         'Surface fluxes and geostrophic wind should be &'//    &
     737          message_string = 'Nudging requires large_scale_forcing = .T.. '//    &
     738                        'Surface fluxes and geostrophic wind should be '//     &
    736739                        'prescribed in file LSF_DATA'
    737740          CALL message( 'check_parameters', 'PA0374', 1, 2, 0, 6, 0 )
     
    740743       IF ( large_scale_forcing  .AND.  ( bc_lr /= 'cyclic'  .OR.              &
    741744                                          bc_ns /= 'cyclic' ) )  THEN
    742           message_string = 'Non-cyclic lateral boundaries do not allow for &'//&
     745          message_string = 'Non-cyclic lateral boundaries do not allow for  '//&
    743746                        'the usage of large scale forcing from external file.'
    744747          CALL message( 'check_parameters', 'PA0375', 1, 2, 0, 6, 0 )
     
    746749
    747750       IF ( large_scale_forcing  .AND.  (  .NOT.  humidity ) )  THEN
    748           message_string = 'The usage of large scale forcing from external &'//&
     751          message_string = 'The usage of large scale forcing from external  '//&
    749752                        'file LSF_DATA requires humidity = .T..'
    750753          CALL message( 'check_parameters', 'PA0376', 1, 2, 0, 6, 0 )
     
    752755
    753756       IF ( large_scale_forcing  .AND.  passive_scalar )  THEN
    754           message_string = 'The usage of large scale forcing from external &'//&
     757          message_string = 'The usage of large scale forcing from external '// &
    755758                        'file LSF_DATA is not implemented for passive scalars'
    756759          CALL message( 'check_parameters', 'PA0440', 1, 2, 0, 6, 0 )
     
    759762       IF ( large_scale_forcing  .AND.  topography /= 'flat'                   &
    760763                              .AND.  .NOT.  lsf_exception )  THEN
    761           message_string = 'The usage of large scale forcing from external &'//&
     764          message_string = 'The usage of large scale forcing from external  '//&
    762765                        'file LSF_DATA is not implemented for non-flat topography'
    763766          CALL message( 'check_parameters', 'PA0377', 1, 2, 0, 6, 0 )
     
    765768
    766769       IF ( large_scale_forcing  .AND.  ocean )  THEN
    767           message_string = 'The usage of large scale forcing from external &'//&
     770          message_string = 'The usage of large scale forcing from external '// &
    768771                        'file LSF_DATA is not implemented for ocean runs'
    769772          CALL message( 'check_parameters', 'PA0378', 1, 2, 0, 6, 0 )
     
    858861                                 TRIM( data_output_pr(var_count) ) //          &
    859862                                 ' is not implemented for ' //                 &
    860                                  'large_scale_forcing = .FALSE.'
    861                 CALL message( 'lsf_nudging_check_data_output_pr', 'PA0393',    &
     863                                 'nudging = .FALSE.'
     864                CALL message( 'lsf_nudging_check_data_output_pr', 'PA0394',    &
    862865                               1, 2, 0, 6, 0 )
    863866             ELSE
     
    873876                                 TRIM( data_output_pr(var_count) ) //          &
    874877                                 ' is not implemented for ' //                 &
    875                                  'large_scale_forcing = .FALSE.'
    876                 CALL message( 'lsf_nudging_check_data_output_pr', 'PA0393',    &
     878                                 'nudging = .FALSE.'
     879                CALL message( 'lsf_nudging_check_data_output_pr', 'PA0394',    &
    877880                               1, 2, 0, 6, 0 )
    878881             ELSE
     
    888891                                 TRIM( data_output_pr(var_count) ) //          &
    889892                                 ' is not implemented for ' //                 &
    890                                  'large_scale_forcing = .FALSE.'
    891                 CALL message( 'lsf_nudging_check_data_output_pr', 'PA0393',    &
     893                                 'nudging = .FALSE.'
     894                CALL message( 'lsf_nudging_check_data_output_pr', 'PA0394',    &
    892895                               1, 2, 0, 6, 0 )
    893896             ELSE
     
    903906                                 TRIM( data_output_pr(var_count) ) //          &
    904907                                 ' is not implemented for ' //                 &
    905                                  'large_scale_forcing = .FALSE.'
    906                 CALL message( 'lsf_nudging_check_data_output_pr', 'PA0393',    &
     908                                 'nudging = .FALSE.'
     909                CALL message( 'lsf_nudging_check_data_output_pr', 'PA0394',    &
    907910                               1, 2, 0, 6, 0 )
    908911             ELSE
     
    12801283
    12811284             IF ( ierrn /= 0 )  THEN
    1282                WRITE ( message_string, * ) 'No time dependent surface variables ' //&
    1283                                  'in&LSF_DATA for end of run found'
     1285               WRITE ( message_string, * ) 'No time dependent surface ' //     &
     1286                                 'variables in LSF_DATA for end of run found'
    12841287
    12851288                CALL message( 'ls_forcing', 'PA0363', 1, 2, 0, 6, 0 )
     
    12891292          IF ( time_surf(1) > end_time )  THEN
    12901293             WRITE ( message_string, * ) 'Time dependent surface variables in ' // &
    1291                                          '&LSF_DATA set in after end of ' ,        &
     1294                                         'LSF_DATA set in after end of ' ,         &
    12921295                                         'simulation - lsf_surf is set to FALSE'
    12931296             CALL message( 'ls_forcing', 'PA0371', 0, 0, 0, 6, 0 )
     
    13161319                IF ( ierrn < 0 )  THEN
    13171320                   WRITE( message_string, * ) 'No time dependent vertical profiles',&
    1318                                     ' in&LSF_DATA for end of run found'
     1321                                    ' in LSF_DATA for end of run found'
    13191322                   CALL message( 'ls_forcing', 'PA0372', 1, 2, 0, 6, 0 )
    13201323                ENDIF
     
    13621365
    13631366                   IF ( ierrn /= 0 )  THEN
    1364                       WRITE( message_string, * ) 'zu(',k,') = ', zu(k), 'm ', &
    1365                            'is higher than the maximum height in LSF_DATA which ',&
    1366                            'is ', lowheight, 'm. Interpolation on PALM ',         &
     1367                      WRITE( message_string, * ) 'zu(',k,') = ', zu(k), 'm ',  &
     1368                           'is higher than the maximum height in LSF_DATA ',   &
     1369                           'which is ', lowheight, 'm. Interpolation on PALM ',&
    13671370                           'grid is not possible.'
    13681371                      CALL message( 'ls_forcing', 'PA0395', 1, 2, 0, 6, 0 )
     
    18321835 
    18331836                ierrn = 0
    1834                 READ ( finput, *, IOSTAT=ierrn )  highheight , hightnudge , &
    1835                                                   highunudge , highvnudge , &
    1836                                                   highwnudge , highptnudge, &
     1837                READ ( finput, *, IOSTAT=ierrn )  highheight , hightnudge ,    &
     1838                                                  highunudge , highvnudge ,    &
     1839                                                  highwnudge , highptnudge,    &
    18371840                                                  highqnudge
    18381841                IF (ierrn /= 0 )  THEN
    1839                    WRITE( message_string, * ) 'zu(',k,') = ', zu(k), 'm is ',&
    1840                         'higher than the maximum height in NUDING_DATA which ',  &
    1841                         'is ', lowheight, 'm. Interpolation on PALM ',           &
     1842                   WRITE( message_string, * ) 'zu(',k,') = ', zu(k), 'm is ',  &
     1843                        'higher than the maximum height in NUDING_DATA which ',&
     1844                        'is ', lowheight, 'm. Interpolation on PALM ',         &
    18421845                        'grid is not possible.'
    18431846                   CALL message( 'nudging', 'PA0364', 1, 2, 0, 6, 0 )
     
    18721875          ENDIF
    18731876
    1874           WRITE( message_string, * ) 'Initial profiles of u, v and ',          &
    1875                                      'scalars from NUDGING_DATA are used.'
     1877          WRITE( message_string, * ) 'Initial profiles of u, v, pt and q ',    &
     1878                                     'from NUDGING_DATA are used.'
    18761879          CALL message( 'large_scale_forcing_nudging', 'PA0370', 0, 0, 0, 6, 0 )
    18771880       ENDIF
  • palm/trunk/SOURCE/lpm_droplet_condensation.f90

    r3039 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 3039 2018-05-24 13:13:11Z schwenkel
    2730! bugfix for lcm with grid stretching
    2831!
     
    380383!--    the model timestep might be too long.
    381384       IF ( ql_c(kp,jp,ip) > 100.0_wp )  THEN
    382           WRITE( message_string, * ) 'k=',kp,' j=',jp,' i=',ip,      &
    383                        ' ql_c=',ql_c(kp,jp,ip), ' &part(',n,')%wf=', &
     385          WRITE( message_string, * ) 'k=',kp,' j=',jp,' i=',ip,                &
     386                       ' ql_c=',ql_c(kp,jp,ip), ' part(',n,')%wf=',            &
    384387                       particles(n)%weight_factor,' delta_r=',delta_r
    385388          CALL message( 'lpm_droplet_condensation', 'PA0143', 2, 2, -1, 6, 1 )
     
    390393       delta_r = new_r(n) - particles(n)%radius
    391394       IF ( delta_r < 0.0_wp  .AND. new_r(n) < 0.0_wp )  THEN
    392           WRITE( message_string, * ) '#1 k=',kp,' j=',jp,' i=',ip,    &
    393                        ' e_s=',e_s, ' e_a=',e_a,' t_int=',t_int,      &
    394                        ' &delta_r=',delta_r,                          &
     395          WRITE( message_string, * ) '#1 k=',kp,' j=',jp,' i=',ip,             &
     396                       ' e_s=',e_s, ' e_a=',e_a,' t_int=',t_int,               &
     397                       ' delta_r=',delta_r,                                    &
    395398                       ' particle_radius=',particles(n)%radius
    396399          CALL message( 'lpm_droplet_condensation', 'PA0144', 2, 2, -1, 6, 1 )
  • palm/trunk/SOURCE/lpm_init.f90

    r3039 r3045  
    2525! -----------------
    2626! $Id$
     27! Error message revised
     28!
     29! 3039 2018-05-24 13:13:11Z schwenkel
    2730! bugfix for lcm with grid stretching
    2831!
     
    333336!-- Check the number of particle groups.
    334337    IF ( number_of_particle_groups > max_number_of_particle_groups )  THEN
    335        WRITE( message_string, * ) 'max_number_of_particle_groups =',      &
    336                                   max_number_of_particle_groups ,         &
    337                                   '&number_of_particle_groups reset to ', &
     338       WRITE( message_string, * ) 'max_number_of_particle_groups =',           &
     339                                  max_number_of_particle_groups ,              &
     340                                  ' number_of_particle_groups reset to ',      &
    338341                                  max_number_of_particle_groups
    339342       CALL message( 'lpm_init', 'PA0213', 0, 1, 0, 6, 0 )
     
    345348!-- propably (not realized so far).
    346349    IF ( surf_def_h(1)%ns >= 1 )  THEN
    347        WRITE( message_string, * ) 'Overhanging topography do not work '// &
     350       WRITE( message_string, * ) 'Overhanging topography do not work '//      &
    348351                                  'with particles'
    349352       CALL message( 'lpm_init', 'PA0212', 0, 1, 0, 6, 0 )
     
    466469
    467470       CASE DEFAULT
    468           WRITE( message_string, * )  'unknown boundary condition ',   &
     471          WRITE( message_string, * )  'unknown boundary condition ',           &
    469472                                       'bc_par_b = "', TRIM( bc_par_b ), '"'
    470473          CALL message( 'lpm_init', 'PA0217', 1, 2, 0, 6, 0 )
     
    483486
    484487       CASE DEFAULT
    485           WRITE( message_string, * ) 'unknown boundary condition ',   &
     488          WRITE( message_string, * ) 'unknown boundary condition ',            &
    486489                                     'bc_par_t = "', TRIM( bc_par_t ), '"'
    487490          CALL message( 'lpm_init', 'PA0218', 1, 2, 0, 6, 0 )
     
    540543
    541544       CASE DEFAULT
    542           WRITE( message_string, * )  'unknown splitting condition ',   &
    543                                        'splitting_mode = "', TRIM( splitting_mode ), '"'
     545          WRITE( message_string, * )  'unknown splitting_mode = "',            &
     546                                      TRIM( splitting_mode ), '"'
    544547          CALL message( 'lpm_init', 'PA0146', 1, 2, 0, 6, 0 )
    545548
     
    557560
    558561       CASE DEFAULT
    559           WRITE( message_string, * )  'unknown splitting function ',   &
    560                                        'splitting_function = "', TRIM( splitting_function ), '"'
     562          WRITE( message_string, * )  'unknown splitting function = "',        &
     563                                       TRIM( splitting_function ), '"'
    561564          CALL message( 'lpm_init', 'PA0147', 1, 2, 0, 6, 0 )
    562565
     
    619622       DO  i = 1, number_of_particle_groups
    620623          IF ( density_ratio(i) /= 0.0_wp  .AND.  radius(i) == 0 )  THEN
    621              WRITE( message_string, * ) 'particle group #', i, 'has a', &
     624             WRITE( message_string, * ) 'particle group #', i, ' has a',      &
    622625                                        'density ratio /= 0 but radius = 0'
    623626             CALL message( 'lpm_init', 'PA0215', 1, 2, 0, 6, 0 )
  • palm/trunk/SOURCE/lpm_read_restart_file.f90

    r2718 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised, code adjusted according to PALMs coding standard
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    114117!-- First open the input unit.
    115118    IF ( myid_char == '' )  THEN
    116        OPEN ( 90, FILE='PARTICLE_RESTART_DATA_IN'//myid_char, &
     119       OPEN ( 90, FILE='PARTICLE_RESTART_DATA_IN'//myid_char,                  &
    117120                  FORM='UNFORMATTED' )
    118121    ELSE
    119        OPEN ( 90, FILE='PARTICLE_RESTART_DATA_IN/'//myid_char, &
     122       OPEN ( 90, FILE='PARTICLE_RESTART_DATA_IN/'//myid_char,                 &
    120123                  FORM='UNFORMATTED' )
    121124    ENDIF
     
    126129    particle_binary_version = '4.0'
    127130    IF ( TRIM( version_on_file ) /= TRIM( particle_binary_version ) )  THEN
    128        message_string = 'version mismatch concerning data from prior ' // &
    129                         'run &version on file    = "' //                  &
    130                                       TRIM( version_on_file ) //          &
    131                         '&version in program = "' //                      &
     131       message_string = 'version mismatch concerning data from prior ' //      &
     132                        'run version on file = "' //                           &
     133                                      TRIM( version_on_file ) //               &
     134                        ' version in program = "' //                           &
    132135                                      TRIM( particle_binary_version ) // '"'
    133136       CALL message( 'lpm_read_restart_file', 'PA0214', 1, 2, 0, 6, 0 )
     
    138141!-- min_nr_particle, the remainder is initialized by zero_particle to avoid
    139142!-- errors.
    140     zero_particle = particle_type( 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp,  &
    141                                    0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp,  &
    142                                    0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp,  &
    143                                    0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp,  &
     143    zero_particle = particle_type( 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp,     &
     144                                   0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp,     &
     145                                   0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp,     &
     146                                   0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp, 0.0_wp,     &
    144147                                   0, 0, 0_idp, .FALSE., -1 )
    145148!
    146149!-- Read some particle parameters and the size of the particle arrays,
    147150!-- allocate them and read their contents.
    148     READ ( 90 )  bc_par_b, bc_par_lr, bc_par_ns, bc_par_t,                  &
    149                  number_of_particle_groups, particle_groups, time_prel,     &
     151    READ ( 90 )  bc_par_b, bc_par_lr, bc_par_ns, bc_par_t,                     &
     152                 number_of_particle_groups, particle_groups, time_prel,        &
    150153                 time_write_particle_data
    151154
    152     ALLOCATE( prt_count(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                     &
     155    ALLOCATE( prt_count(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                        &
    153156              grid_particles(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    154157
  • palm/trunk/SOURCE/model_1d_mod.f90

    r2965 r3045  
    2525! -----------------
    2626! $Id$
     27! Error message revised
     28!
     29! 2965 2018-04-13 07:37:25Z scharf
    2730! adjusted format string for 1D run control output
    2831!
     
    10321035       stop_dt_1d = .TRUE.
    10331036
    1034        WRITE( message_string, * ) 'timestep has exceeded the lower limit &', &
     1037       WRITE( message_string, * ) 'timestep has exceeded the lower limit ',    &
    10351038                                  'dt_1d = ',dt_1d,' s   simulation stopped!'
    10361039       CALL message( 'timestep_1d', 'PA0192', 1, 2, 0, 6, 0 )
  • palm/trunk/SOURCE/modules.f90

    r3026 r3045  
    2525! -----------------
    2626! $Id$
     27! z_max_do2d removed
     28!
     29! 3026 2018-05-22 10:30:53Z schwenkel
    2730! Changed the name specific humidity to mixing ratio, since we are computing
    2831! mixing ratios.
     
    15191522    REAL(wp) ::  zeta_max = 20.0_wp                            !< namelist parameter
    15201523    REAL(wp) ::  zeta_min = -20.0_wp                           !< namelist parameter
    1521     REAL(wp) ::  z_max_do2d = -1.0_wp                          !< namelist parameter
    15221524    REAL(wp) ::  z0h_factor = 1.0_wp                           !< namelist parameter
    15231525
  • palm/trunk/SOURCE/netcdf_data_input_mod.f90

    r3041 r3045  
    2020! Current revisions:
    2121! -----------------
    22 !
    23 !
     22! 
     23! 
    2424! Former revisions:
    2525! -----------------
    2626! $Id$
     27! Error message revised
     28!
     29! 3041 2018-05-25 10:39:54Z gronemeier
    2730! Add data type for global file attributes
    2831! Add read of global attributes of static driver
     
    20032006                GOTO 12
    20042007
    2005  10             message_string = 'file TOPOGRAPHY'//TRIM( coupling_char )//    &
    2006                                  ' does not exist'
     2008 10             message_string = 'file TOPOGRAPHY_DATA'//                      &
     2009                                 TRIM( coupling_char )// ' does not exist'
    20072010                CALL message( 'netcdf_data_input_mod', 'PA0208', 1, 2, 0, 6, 0 )
    20082011
  • palm/trunk/SOURCE/netcdf_interface_mod.f90

    r3004 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised, code adjusted to PALMs coding standards, CASE pt_ext
     28! pt_new disabled, comment revised
     29!
     30! 3004 2018-04-27 12:33:25Z Giersch
    2731! .NOT. found in if-query added to account for variables found in tcm
    2832!
     
    708712       IF ( netcdf_deflate < 0  .OR.  netcdf_deflate > 9 )  THEN
    709713          WRITE ( message_string, '(A,I3,A)' ) 'netcdf_deflate out of ' //     &
    710                                       'range & given value: ', netcdf_deflate, &
     714                                      'range given value: ', netcdf_deflate,  &
    711715                                      ', allowed range: 0-9'
    712716          CALL message( 'netcdf_define_header', 'PA0355', 2, 2, 0, 6, 0 )
     
    734738
    735739!
    736 !-- Select the mode to be processed. Possibilities are 3d, mask, xy, xz, yz,
    737 !-- pr and ts.
     740!-- Select the mode to be processed. Possibilities are 3d, ma (mask), xy, xz,
     741!-- yz, pr (profiles), ps (particle timeseries), fl (flight data), ts
     742!-- (timeseries) or sp (spectra)
    738743    SELECT CASE ( mode )
    739744
     
    11851190
    11861191          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    1187              WRITE ( message_string, * ) 'netCDF file for ', TRIM( var ), &
    1188                   ' data for mask', mid, ' from previous run found,', &
    1189                   '&but this file cannot be extended due to variable ', &
    1190                   'mismatch.&New file is created instead.'
     1192             WRITE ( message_string, * ) 'netCDF file for ', TRIM( var ),       &
     1193                  ' data for mask', mid, ' from previous run found,',           &
     1194                  ' but this file cannot be extended due to variable ',        &
     1195                  'mismatch. New file is created instead.'
    11911196             CALL message( 'define_netcdf_header', 'PA0335', 0, 1, 0, 6, 0 )
    11921197             extend = .FALSE.
     
    12061211          id_dim_zu_mask(mid,av) = id_dim_zu_mask_old(1)
    12071212
    1208           nc_stat = NF90_INQUIRE_DIMENSION( id_set_mask(mid,av),    &
    1209                                             id_dim_zu_mask(mid,av), &
     1213          nc_stat = NF90_INQUIRE_DIMENSION( id_set_mask(mid,av),               &
     1214                                            id_dim_zu_mask(mid,av),            &
    12101215                                            len = nz_old )
    12111216          CALL netcdf_handle_error( 'netcdf_define_header', 510 )
    12121217
    12131218          IF ( mask_size(mid,3) /= nz_old )  THEN
    1214              WRITE ( message_string, * ) 'netCDF file for ', TRIM( var ), &
    1215                   ' data for mask', mid, ' from previous run found,', &
    1216                   '&but this file cannot be extended due to mismatch in ', &
    1217                   ' number of&vertical grid points.', &
    1218                   '&New file is created instead.'
     1219             WRITE ( message_string, * ) 'netCDF file for ', TRIM( var ),      &
     1220                  ' data for mask', mid, ' from previous run found,',          &
     1221                  ' but this file cannot be extended due to mismatch in ',    &
     1222                  ' number of vertical grid points.',                          &
     1223                  ' New file is created instead.'
    12191224             CALL message( 'define_netcdf_header', 'PA0336', 0, 1, 0, 6, 0 )
    12201225             extend = .FALSE.
     
    12271232!--       The current time must be larger than the last output time
    12281233!--       on the file.
    1229           nc_stat = NF90_INQ_VARID( id_set_mask(mid,av), 'time', &
     1234          nc_stat = NF90_INQ_VARID( id_set_mask(mid,av), 'time',               &
    12301235                                    id_var_time_mask(mid,av) )
    12311236          CALL netcdf_handle_error( 'netcdf_define_header', 511 )
    12321237
    1233           nc_stat = NF90_INQUIRE_VARIABLE( id_set_mask(mid,av), &
    1234                                            id_var_time_mask(mid,av), &
     1238          nc_stat = NF90_INQUIRE_VARIABLE( id_set_mask(mid,av),                &
     1239                                           id_var_time_mask(mid,av),           &
    12351240                                           dimids = id_dim_time_old )
    12361241          CALL netcdf_handle_error( 'netcdf_define_header', 512 )
    12371242          id_dim_time_mask(mid,av) = id_dim_time_old(1)
    12381243
    1239           nc_stat = NF90_INQUIRE_DIMENSION( id_set_mask(mid,av), &
    1240                                             id_dim_time_mask(mid,av), &
     1244          nc_stat = NF90_INQUIRE_DIMENSION( id_set_mask(mid,av),               &
     1245                                            id_dim_time_mask(mid,av),          &
    12411246                                            len = domask_time_count(mid,av) )
    12421247          CALL netcdf_handle_error( 'netcdf_define_header', 513 )
    12431248
    1244           nc_stat = NF90_GET_VAR( id_set_mask(mid,av), &
    1245                                   id_var_time_mask(mid,av), &
    1246                                   last_time_coordinate,              &
    1247                                   start = (/ domask_time_count(mid,av) /), &
     1249          nc_stat = NF90_GET_VAR( id_set_mask(mid,av),                         &
     1250                                  id_var_time_mask(mid,av),                    &
     1251                                  last_time_coordinate,                        &
     1252                                  start = (/ domask_time_count(mid,av) /),     &
    12481253                                  count = (/ 1 /) )
    12491254          CALL netcdf_handle_error( 'netcdf_define_header', 514 )
    12501255
    12511256          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    1252              WRITE ( message_string, * ) 'netCDF file for ', TRIM( var ), &
    1253                   ' data for mask', mid, ' from previous run found,', &
    1254                   '&but this file cannot be extended because the current ', &
    1255                   'output time&is less or equal than the last output time ', &
    1256                   'on this file.&New file is created instead.'
     1257             WRITE ( message_string, * ) 'netCDF file for ', TRIM( var ),      &
     1258                  ' data for mask', mid, ' from previous run found,',          &
     1259                  ' but this file cannot be extended because the current ',    &
     1260                  'output time is less or equal than the last output time ',  &
     1261                  'on this file. New file is created instead.'
    12571262             CALL message( 'define_netcdf_header', 'PA0337', 0, 1, 0, 6, 0 )
    12581263             domask_time_count(mid,av) = 0
     
    12831288             time_average_text = ' '
    12841289          ELSE
    1285              WRITE (time_average_text, '('', '',F7.1,'' s average'')') &
     1290             WRITE (time_average_text, '('', '',F7.1,'' s average'')')         &
    12861291                                                            averaging_interval
    12871292          ENDIF
    12881293          nc_stat = NF90_REDEF( id_set_mask(mid,av) )
    12891294          CALL netcdf_handle_error( 'netcdf_define_header', 516 )
    1290           nc_stat = NF90_PUT_ATT( id_set_mask(mid,av), NF90_GLOBAL, 'title', &
    1291                                   TRIM( run_description_header ) //    &
     1295          nc_stat = NF90_PUT_ATT( id_set_mask(mid,av), NF90_GLOBAL, 'title',   &
     1296                                  TRIM( run_description_header ) //            &
    12921297                                  TRIM( time_average_text ) )
    12931298          CALL netcdf_handle_error( 'netcdf_define_header', 517 )
    12941299          nc_stat = NF90_ENDDEF( id_set_mask(mid,av) )
    12951300          CALL netcdf_handle_error( 'netcdf_define_header', 518 )
    1296           WRITE ( message_string, * ) 'netCDF file for ', TRIM( var ), &
    1297                ' data for mask', mid, ' from previous run found.', &
    1298                '&This file will be extended.'
     1301          WRITE ( message_string, * ) 'netCDF file for ', TRIM( var ),         &
     1302               ' data for mask', mid, ' from previous run found.',             &
     1303               ' This file will be extended.'
    12991304          CALL message( 'define_netcdf_header', 'PA0338', 0, 0, 0, 6, 0 )
    13001305!
     
    13071312!
    13081313!--       Define some global attributes of the dataset
    1309           nc_stat = NF90_PUT_ATT( id_set_3d(av), NF90_GLOBAL, 'Conventions', &
     1314          nc_stat = NF90_PUT_ATT( id_set_3d(av), NF90_GLOBAL, 'Conventions',   &
    13101315                                  'COARDS' )
    13111316          CALL netcdf_handle_error( 'netcdf_define_header', 62 )
     
    17561761             message_string = 'netCDF file for volume data ' //             &
    17571762                              TRIM( var ) // ' from previous run found,' // &
    1758                               '&but this file cannot be extended due to' // &
     1763                              ' but this file cannot be extended due to' // &
    17591764                              ' variable mismatch.' //                      &
    1760                               '&New file is created instead.'
     1765                              ' New file is created instead.'
    17611766             CALL message( 'define_netcdf_header', 'PA0245', 0, 1, 0, 6, 0 )
    17621767             extend = .FALSE.
     
    17811786              message_string = 'netCDF file for volume data ' //             &
    17821787                               TRIM( var ) // ' from previous run found,' // &
    1783                                '&but this file cannot be extended due to' // &
     1788                               ' but this file cannot be extended due to' // &
    17841789                               ' mismatch in number of' //                   &
    1785                                '&vertical grid points (nz_do3d).' //         &
    1786                                '&New file is created instead.'
     1790                               ' vertical grid points (nz_do3d).' //         &
     1791                               ' New file is created instead.'
    17871792             CALL message( 'define_netcdf_header', 'PA0246', 0, 1, 0, 6, 0 )
    17881793             extend = .FALSE.
     
    18221827
    18231828          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    1824              message_string = 'netCDF file for volume data ' // &
     1829             message_string = 'netCDF file for volume data ' //             &
    18251830                              TRIM( var ) // ' from previous run found,' // &
    1826                               '&but this file cannot be extended becaus' // &
     1831                              ' but this file cannot be extended becaus' // &
    18271832                              'e the current output time' //                &
    1828                               '&is less or equal than the last output t' // &
     1833                              ' is less or equal than the last output t' // &
    18291834                              'ime on this file.' //                        &
    1830                               '&New file is created instead.'
     1835                              ' New file is created instead.'
    18311836             CALL message( 'define_netcdf_header', 'PA0247', 0, 1, 0, 6, 0 )
    18321837             do3d_time_count(av) = 0
     
    18421847                message_string = 'netCDF file for volume data ' // &
    18431848                                 TRIM( var ) // ' from previous run found,' // &
    1844                                  '&but this file cannot be extended becaus' // &
    1845                                  'e the number of output time levels&has b' // &
     1849                                 ' but this file cannot be extended becaus' // &
     1850                                 'e the number of output time levels has b' // &
    18461851                                 'een increased compared to the previous s' // &
    18471852                                 'imulation.' //                               &
    1848                                  '&New file is created instead.'
     1853                                 ' New file is created instead.'
    18491854                CALL message( 'define_netcdf_header', 'PA0388', 0, 1, 0, 6, 0 )
    18501855                do3d_time_count(av) = 0
     
    19121917          message_string = 'netCDF file for volume data ' //             &
    19131918                           TRIM( var ) // ' from previous run found.' // &
    1914                            '&This file will be extended.'
     1919                           ' This file will be extended.'
    19151920          CALL message( 'define_netcdf_header', 'PA0248', 0, 0, 0, 6, 0 )
    19161921
     
    24702475             message_string = 'netCDF file for cross-sections ' //           &
    24712476                              TRIM( var ) // ' from previous run found,' //  &
    2472                               '& but this file cannot be extended due to' // &
     2477                              ' but this file cannot be extended due to' // &
    24732478                              ' variable mismatch.' //                       &
    2474                               '&New file is created instead.'
     2479                              ' New file is created instead.'
    24752480             CALL message( 'define_netcdf_header', 'PA0249', 0, 1, 0, 6, 0 )
    24762481             extend = .FALSE.
     
    25032508             message_string = 'netCDF file for cross-sections ' //          &
    25042509                              TRIM( var ) // ' from previous run found,' // &
    2505                               '&but this file cannot be extended due to' // &
     2510                              ' but this file cannot be extended due to' // &
    25062511                              ' mismatch in number of' //                   &
    2507                               '&cross sections.' //                         &
    2508                               '&New file is created instead.'
     2512                              ' cross sections.' //                         &
     2513                              ' New file is created instead.'
    25092514             CALL message( 'define_netcdf_header', 'PA0250', 0, 1, 0, 6, 0 )
    25102515             extend = .FALSE.
     
    25222527             IF ( section(i,1) /= -1 )  THEN
    25232528                IF ( zu(section(i,1)) /= netcdf_data(i) )  THEN
    2524                    message_string = 'netCDF file for cross-sections ' //     &
    2525                                TRIM( var ) // ' from previous run found,' // &
    2526                                '&but this file cannot be extended' //        &
    2527                                ' due to mismatch in cross' //                &
    2528                                '&section levels.' //                         &
    2529                                '&New file is created instead.'
    2530                    CALL message( 'define_netcdf_header', 'PA0251', &
     2529                   message_string = 'netCDF file for cross-sections ' //       &
     2530                               TRIM( var ) // ' from previous run found,' //   &
     2531                               ' but this file cannot be extended' //          &
     2532                               ' due to mismatch in cross' //                  &
     2533                               ' section levels.' //                           &
     2534                               ' New file is created instead.'
     2535                   CALL message( 'define_netcdf_header', 'PA0251',             &
    25312536                                                                 0, 1, 0, 6, 0 )
    25322537                   extend = .FALSE.
     
    25352540             ELSE
    25362541                IF ( -1.0_wp /= netcdf_data(i) )  THEN
    2537                    message_string = 'netCDF file for cross-sections ' //     &
    2538                                TRIM( var ) // ' from previous run found,' // &
    2539                                '&but this file cannot be extended' //        &
    2540                                ' due to mismatch in cross' //                &
    2541                                '&section levels.' //                         &
    2542                                '&New file is created instead.'
    2543                    CALL message( 'define_netcdf_header', 'PA0251',&
     2542                   message_string = 'netCDF file for cross-sections ' //       &
     2543                               TRIM( var ) // ' from previous run found,' //   &
     2544                               ' but this file cannot be extended' //          &
     2545                               ' due to mismatch in cross' //                  &
     2546                               ' section levels.' //                           &
     2547                               ' New file is created instead.'
     2548                   CALL message( 'define_netcdf_header', 'PA0251',             &
    25442549                                                                 0, 1, 0, 6, 0 )
    25452550                   extend = .FALSE.
     
    25752580          IF ( netcdf_data_format < 5 ) do2d_xy_time_count(av) = ntime_count
    25762581
    2577           nc_stat = NF90_GET_VAR( id_set_xy(av), id_var_time_xy(av),    &
    2578                                   last_time_coordinate,                 &
    2579                                   start = (/ do2d_xy_time_count(av) /), &
     2582          nc_stat = NF90_GET_VAR( id_set_xy(av), id_var_time_xy(av),           &
     2583                                  last_time_coordinate,                        &
     2584                                  start = (/ do2d_xy_time_count(av) /),        &
    25802585                                  count = (/ 1 /) )
    25812586          CALL netcdf_handle_error( 'netcdf_define_header', 137 )
    25822587
    25832588          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    2584              message_string = 'netCDF file for cross sections ' //          &
    2585                               TRIM( var ) // ' from previous run found,' // &
    2586                               '&but this file cannot be extended becaus' // &
    2587                               'e the current output time' //                &
    2588                               '&is less or equal than the last output t' // &
    2589                               'ime on this file.' //                        &
    2590                               '&New file is created instead.'
     2589             message_string = 'netCDF file for cross sections ' //             &
     2590                              TRIM( var ) // ' from previous run found,' //    &
     2591                              ' but this file cannot be extended becaus' //    &
     2592                              'e the current output time' //                   &
     2593                              ' is less or equal than the last output t' //    &
     2594                              'ime on this file.' //                           &
     2595                              ' New file is created instead.'
    25912596             CALL message( 'define_netcdf_header', 'PA0252', 0, 1, 0, 6, 0 )
    25922597             do2d_xy_time_count(av) = 0
     
    26022607                message_string = 'netCDF file for cross sections ' //          &
    26032608                                 TRIM( var ) // ' from previous run found,' // &
    2604                                  '&but this file cannot be extended becaus' // &
    2605                                  'e the number of output time levels&has b' // &
     2609                                 ' but this file cannot be extended becaus' // &
     2610                                 'e the number of output time levels has b' // &
    26062611                                 'een increased compared to the previous s' // &
    26072612                                 'imulation.' //                               &
    2608                                  '&New file is created instead.'
     2613                                 ' New file is created instead.'
    26092614                CALL message( 'define_netcdf_header', 'PA0389', 0, 1, 0, 6, 0 )
    26102615                do2d_xy_time_count(av) = 0
     
    26662671          nc_stat = NF90_REDEF( id_set_xy(av) )
    26672672          CALL netcdf_handle_error( 'netcdf_define_header', 431 )
    2668           nc_stat = NF90_PUT_ATT( id_set_xy(av), NF90_GLOBAL, 'title', &
    2669                                   TRIM( run_description_header ) //    &
     2673          nc_stat = NF90_PUT_ATT( id_set_xy(av), NF90_GLOBAL, 'title',         &
     2674                                  TRIM( run_description_header ) //            &
    26702675                                  TRIM( time_average_text ) )
    26712676          CALL netcdf_handle_error( 'netcdf_define_header', 139 )
    26722677          nc_stat = NF90_ENDDEF( id_set_xy(av) )
    26732678          CALL netcdf_handle_error( 'netcdf_define_header', 432 )
    2674           message_string = 'netCDF file for cross-sections ' //           &
    2675                             TRIM( var ) // ' from previous run found.' // &
    2676                            '&This file will be extended.'
     2679          message_string = 'netCDF file for cross-sections ' //                &
     2680                            TRIM( var ) // ' from previous run found.' //      &
     2681                           ' This file will be extended.'
    26772682          CALL message( 'define_netcdf_header', 'PA0253', 0, 0, 0, 6, 0 )
    26782683         
     
    26822687!
    26832688!--       Define some global attributes of the dataset
    2684           nc_stat = NF90_PUT_ATT( id_set_xz(av), NF90_GLOBAL, 'Conventions', &
     2689          nc_stat = NF90_PUT_ATT( id_set_xz(av), NF90_GLOBAL, 'Conventions',   &
    26852690                                  'COARDS' )
    26862691          CALL netcdf_handle_error( 'netcdf_define_header', 140 )
     
    26892694             time_average_text = ' '
    26902695          ELSE
    2691              WRITE (time_average_text, '('', '',F7.1,'' s average'')') &
     2696             WRITE (time_average_text, '('', '',F7.1,'' s average'')')         &
    26922697                                                            averaging_interval
    26932698          ENDIF
    2694           nc_stat = NF90_PUT_ATT( id_set_xz(av), NF90_GLOBAL, 'title', &
    2695                                   TRIM( run_description_header )  //   &
     2699          nc_stat = NF90_PUT_ATT( id_set_xz(av), NF90_GLOBAL, 'title',         &
     2700                                  TRIM( run_description_header )  //           &
    26962701                                  TRIM( time_average_text ) )
    26972702          CALL netcdf_handle_error( 'netcdf_define_header', 141 )
     
    28792884
    28802885                      IF ( .NOT. found )  THEN
    2881                          WRITE ( message_string, * ) 'no grid defined for', &
     2886                         WRITE ( message_string, * ) 'no grid defined for',    &
    28822887                                                ' variable ', TRIM( do2d(av,i) )
    2883                          CALL message( 'define_netcdf_header', 'PA0244', &
     2888                         CALL message( 'define_netcdf_header', 'PA0244',       &
    28842889                                       0, 1, 0, 6, 0 )
    28852890                      ENDIF
     
    31133118             message_string = 'netCDF file for cross-sections ' //           &
    31143119                              TRIM( var ) // ' from previous run found,' //  &
    3115                               '& but this file cannot be extended due to' // &
     3120                              ' but this file cannot be extended due to' // &
    31163121                              ' variable mismatch.' //                       &
    3117                               '&New file is created instead.'
     3122                              ' New file is created instead.'
    31183123             CALL message( 'define_netcdf_header', 'PA0249', 0, 1, 0, 6, 0 )
    31193124             extend = .FALSE.
     
    31463151             message_string = 'netCDF file for cross-sections ' //          &
    31473152                              TRIM( var ) // ' from previous run found,' // &
    3148                               '&but this file cannot be extended due to' // &
     3153                              ' but this file cannot be extended due to' // &
    31493154                              ' mismatch in number of' //                   &
    3150                               '&cross sections.' //                         &
    3151                               '&New file is created instead.'
     3155                              ' cross sections.' //                         &
     3156                              ' New file is created instead.'
    31523157             CALL message( 'define_netcdf_header', 'PA0250', 0, 1, 0, 6, 0 )
    31533158             extend = .FALSE.
     
    31653170             IF ( section(i,2) /= -1 )  THEN
    31663171                IF ( ( ( section(i,2) + 0.5 ) * dy ) /= netcdf_data(i) )  THEN
    3167                    message_string = 'netCDF file for cross-sections ' //     &
    3168                                TRIM( var ) // ' from previous run found,' // &
    3169                                '&but this file cannot be extended' //        &
    3170                                ' due to mismatch in cross' //                &
    3171                                '&section levels.' //                         &
    3172                                 '&New file is created instead.'
    3173                    CALL message( 'define_netcdf_header', 'PA0251', &
     3172                   message_string = 'netCDF file for cross-sections ' //       &
     3173                               TRIM( var ) // ' from previous run found,' //   &
     3174                               ' but this file cannot be extended' //          &
     3175                               ' due to mismatch in cross' //                  &
     3176                               ' section levels.' //                           &
     3177                               ' New file is created instead.'
     3178                   CALL message( 'define_netcdf_header', 'PA0251',             &
    31743179                                                                 0, 1, 0, 6, 0 )
    31753180                   extend = .FALSE.
     
    31783183             ELSE
    31793184                IF ( -1.0_wp /= netcdf_data(i) )  THEN
    3180                    message_string = 'netCDF file for cross-sections ' //     &
    3181                                TRIM( var ) // ' from previous run found,' // &
    3182                                '&but this file cannot be extended' //        &
    3183                                ' due to mismatch in cross' //                &
    3184                                '&section levels.' //                         &
    3185                                '&New file is created instead.'
    3186                    CALL message( 'define_netcdf_header', 'PA0251', &
     3185                   message_string = 'netCDF file for cross-sections ' //       &
     3186                               TRIM( var ) // ' from previous run found,' //   &
     3187                               ' but this file cannot be extended' //          &
     3188                               ' due to mismatch in cross' //                  &
     3189                               ' section levels.' //                           &
     3190                               ' New file is created instead.'
     3191                   CALL message( 'define_netcdf_header', 'PA0251',             &
    31873192                                                                 0, 1, 0, 6, 0 )
    31883193                   extend = .FALSE.
     
    32183223          IF ( netcdf_data_format < 5 ) do2d_xz_time_count(av) = ntime_count
    32193224
    3220           nc_stat = NF90_GET_VAR( id_set_xz(av), id_var_time_xz(av),    &
    3221                                   last_time_coordinate,                 &
    3222                                   start = (/ do2d_xz_time_count(av) /), &
     3225          nc_stat = NF90_GET_VAR( id_set_xz(av), id_var_time_xz(av),           &
     3226                                  last_time_coordinate,                        &
     3227                                  start = (/ do2d_xz_time_count(av) /),        &
    32233228                                  count = (/ 1 /) )
    32243229          CALL netcdf_handle_error( 'netcdf_define_header', 176 )
    32253230
    32263231          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    3227              message_string = 'netCDF file for cross sections ' //          &
    3228                               TRIM( var ) // ' from previous run found,' // &
    3229                               '&but this file cannot be extended becaus' // &
    3230                               'e the current output time' //                &
    3231                               '&is less or equal than the last output t' // &
    3232                               'ime on this file.' //                        &
    3233                               '&New file is created instead.'
     3232             message_string = 'netCDF file for cross sections ' //             &
     3233                              TRIM( var ) // ' from previous run found,' //    &
     3234                              ' but this file cannot be extended becaus' //    &
     3235                              'e the current output time' //                   &
     3236                              ' is less or equal than the last output t' //    &
     3237                              'ime on this file.' //                           &
     3238                              ' New file is created instead.'
    32343239             CALL message( 'define_netcdf_header', 'PA0252', 0, 1, 0, 6, 0 )
    32353240             do2d_xz_time_count(av) = 0
     
    32453250                message_string = 'netCDF file for cross sections ' // &
    32463251                                 TRIM( var ) // ' from previous run found,' // &
    3247                                  '&but this file cannot be extended becaus' // &
    3248                                  'e the number of output time levels&has b' // &
     3252                                 ' but this file cannot be extended becaus' // &
     3253                                 'e the number of output time levels has b' // &
    32493254                                 'een increased compared to the previous s' // &
    32503255                                 'imulation.' //                               &
    3251                                  '&New file is created instead.'
     3256                                 ' New file is created instead.'
    32523257                CALL message( 'define_netcdf_header', 'PA0390', 0, 1, 0, 6, 0 )
    32533258                do2d_xz_time_count(av) = 0
     
    33253330          nc_stat = NF90_REDEF( id_set_xz(av) )
    33263331          CALL netcdf_handle_error( 'netcdf_define_header', 433 )
    3327           nc_stat = NF90_PUT_ATT( id_set_xz(av), NF90_GLOBAL, 'title', &
    3328                                   TRIM( run_description_header ) //    &
     3332          nc_stat = NF90_PUT_ATT( id_set_xz(av), NF90_GLOBAL, 'title',         &
     3333                                  TRIM( run_description_header ) //            &
    33293334                                  TRIM( time_average_text ) )
    33303335          CALL netcdf_handle_error( 'netcdf_define_header', 178 )
    33313336          nc_stat = NF90_ENDDEF( id_set_xz(av) )
    33323337          CALL netcdf_handle_error( 'netcdf_define_header', 434 )
    3333           message_string = 'netCDF file for cross-sections ' //           &
    3334                             TRIM( var ) // ' from previous run found.' // &
    3335                            '&This file will be extended.'
     3338          message_string = 'netCDF file for cross-sections ' //                &
     3339                            TRIM( var ) // ' from previous run found.' //      &
     3340                           ' This file will be extended.'
    33363341          CALL message( 'define_netcdf_header', 'PA0253', 0, 0, 0, 6, 0 )
    33373342
     
    33413346!
    33423347!--       Define some global attributes of the dataset
    3343           nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'Conventions', &
     3348          nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'Conventions',   &
    33443349                                  'COARDS' )
    33453350          CALL netcdf_handle_error( 'netcdf_define_header', 179 )
     
    33483353             time_average_text = ' '
    33493354          ELSE
    3350              WRITE (time_average_text, '('', '',F7.1,'' s average'')') &
     3355             WRITE (time_average_text, '('', '',F7.1,'' s average'')')         &
    33513356                                                            averaging_interval
    33523357          ENDIF
    3353           nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'title', &
    3354                                   TRIM( run_description_header ) //    &
     3358          nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'title',         &
     3359                                  TRIM( run_description_header ) //            &
    33553360                                  TRIM( time_average_text ) )
    33563361          CALL netcdf_handle_error( 'netcdf_define_header', 180 )
    33573362          IF ( av == 1 )  THEN
    33583363             WRITE ( time_average_text,'(F7.1,'' s avg'')' )  averaging_interval
    3359              nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'time_avg', &
     3364             nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'time_avg',   &
    33603365                                     TRIM( time_average_text ) )
    33613366             CALL netcdf_handle_error( 'netcdf_define_header', 180 )
     
    37593764             message_string = 'netCDF file for cross-sections ' //           &
    37603765                              TRIM( var ) // ' from previous run found,' //  &
    3761                               '& but this file cannot be extended due to' // &
     3766                              ' but this file cannot be extended due to' // &
    37623767                              ' variable mismatch.' //                       &
    3763                               '&New file is created instead.'
     3768                              ' New file is created instead.'
    37643769             CALL message( 'define_netcdf_header', 'PA0249', 0, 1, 0, 6, 0 )
    37653770             extend = .FALSE.
     
    37923797             message_string = 'netCDF file for cross-sections ' //          &
    37933798                              TRIM( var ) // ' from previous run found,' // &
    3794                               '&but this file cannot be extended due to' // &
     3799                              ' but this file cannot be extended due to' // &
    37953800                              ' mismatch in number of' //                   &
    3796                               '&cross sections.' //                         &
    3797                               '&New file is created instead.'
     3801                              ' cross sections.' //                         &
     3802                              ' New file is created instead.'
    37983803             CALL message( 'define_netcdf_header', 'PA0250', 0, 1, 0, 6, 0 )
    37993804             extend = .FALSE.
     
    38113816             IF ( section(i,3) /= -1 )  THEN
    38123817                IF ( ( ( section(i,3) + 0.5 ) * dx ) /= netcdf_data(i) )  THEN
    3813                    message_string = 'netCDF file for cross-sections ' //    &
    3814                               TRIM( var ) // ' from previous run found,' // &
    3815                               '&but this file cannot be extended' //        &
    3816                               ' due to mismatch in cross' //                &
    3817                               '&section levels.' //                         &
    3818                              '&New file is created instead.'
    3819                    CALL message( 'define_netcdf_header', 'PA0251', &
     3818                   message_string = 'netCDF file for cross-sections ' //       &
     3819                              TRIM( var ) // ' from previous run found,' //    &
     3820                              ' but this file cannot be extended' //           &
     3821                              ' due to mismatch in cross' //                   &
     3822                              ' section levels.' //                            &
     3823                              ' New file is created instead.'
     3824                   CALL message( 'define_netcdf_header', 'PA0251',             &
    38203825                                                                 0, 1, 0, 6, 0 )
    38213826                   extend = .FALSE.
     
    38243829             ELSE
    38253830                IF ( -1.0_wp /= netcdf_data(i) )  THEN
    3826                    message_string = 'netCDF file for cross-sections ' //    &
    3827                               TRIM( var ) // ' from previous run found,' // &
    3828                               '&but this file cannot be extended' //        &
    3829                               ' due to mismatch in cross' //                &
    3830                               '&section levels.' //                         &
    3831                               '&New file is created instead.'
    3832                    CALL message( 'define_netcdf_header', 'PA0251', &
     3831                   message_string = 'netCDF file for cross-sections ' //       &
     3832                              TRIM( var ) // ' from previous run found,' //    &
     3833                              ' but this file cannot be extended' //           &
     3834                              ' due to mismatch in cross' //                   &
     3835                              ' section levels.' //                            &
     3836                              ' New file is created instead.'
     3837                   CALL message( 'define_netcdf_header', 'PA0251',             &
    38333838                                                                 0, 1, 0, 6, 0 )
    38343839                   extend = .FALSE.
     
    38643869          IF ( netcdf_data_format < 5 ) do2d_yz_time_count(av) = ntime_count
    38653870
    3866           nc_stat = NF90_GET_VAR( id_set_yz(av), id_var_time_yz(av),    &
    3867                                   last_time_coordinate,                 &
    3868                                   start = (/ do2d_yz_time_count(av) /), &
     3871          nc_stat = NF90_GET_VAR( id_set_yz(av), id_var_time_yz(av),           &
     3872                                  last_time_coordinate,                        &
     3873                                  start = (/ do2d_yz_time_count(av) /),        &
    38693874                                  count = (/ 1 /) )
    38703875          CALL netcdf_handle_error( 'netcdf_define_header', 215 )
    38713876
    38723877          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    3873              message_string = 'netCDF file for cross sections ' //          &
    3874                               TRIM( var ) // ' from previous run found,' // &
    3875                               '&but this file cannot be extended becaus' // &
    3876                               'e the current output time' //                &
    3877                               '&is less or equal than the last output t' // &
    3878                               'ime on this file.' //                        &
    3879                               '&New file is created instead.'
     3878             message_string = 'netCDF file for cross sections ' //             &
     3879                              TRIM( var ) // ' from previous run found,' //    &
     3880                              ' but this file cannot be extended becaus' //    &
     3881                              'e the current output time' //                   &
     3882                              ' is less or equal than the last output t' //    &
     3883                              'ime on this file.' //                           &
     3884                              ' New file is created instead.'
    38803885             CALL message( 'define_netcdf_header', 'PA0252', 0, 1, 0, 6, 0 )
    38813886             do2d_yz_time_count(av) = 0
     
    38913896                message_string = 'netCDF file for cross sections ' //          &
    38923897                                 TRIM( var ) // ' from previous run found,' // &
    3893                                  '&but this file cannot be extended becaus' // &
    3894                                  'e the number of output time levels&has b' // &
     3898                                 ' but this file cannot be extended becaus' // &
     3899                                 'e the number of output time levels has b' // &
    38953900                                 'een increased compared to the previous s' // &
    38963901                                 'imulation.' //                               &
    3897                                  '&New file is created instead.'
     3902                                 ' New file is created instead.'
    38983903                CALL message( 'define_netcdf_header', 'PA0391', 0, 1, 0, 6, 0 )
    38993904                do2d_yz_time_count(av) = 0
     
    39713976          nc_stat = NF90_REDEF( id_set_yz(av) )
    39723977          CALL netcdf_handle_error( 'netcdf_define_header', 435 )
    3973           nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'title', &
    3974                                   TRIM( run_description_header ) //    &
     3978          nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'title',         &
     3979                                  TRIM( run_description_header ) //            &
    39753980                                  TRIM( time_average_text ) )
    39763981          CALL netcdf_handle_error( 'netcdf_define_header', 217 )
    39773982          nc_stat = NF90_ENDDEF( id_set_yz(av) )
    39783983          CALL netcdf_handle_error( 'netcdf_define_header', 436 )
    3979           message_string = 'netCDF file for cross-sections ' //           &
    3980                             TRIM( var ) // ' from previous run found.' // &
    3981                            '&This file will be extended.'
     3984          message_string = 'netCDF file for cross-sections ' //                &
     3985                            TRIM( var ) // ' from previous run found.' //      &
     3986                           ' This file will be extended.'
    39823987          CALL message( 'define_netcdf_header', 'PA0253', 0, 0, 0, 6, 0 )
    39833988
     
    39883993!--       Define some global attributes of the dataset
    39893994          IF ( averaging_interval_pr /= 0.0_wp )  THEN
    3990              WRITE (time_average_text,'('', '',F7.1,'' s average'')') &
     3995             WRITE (time_average_text,'('', '',F7.1,'' s average'')')          &
    39913996                                                            averaging_interval_pr
    3992              nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'title',   &
    3993                                      TRIM( run_description_header ) //  &
     3997             nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'title',          &
     3998                                     TRIM( run_description_header ) //         &
    39943999                                     TRIM( time_average_text ) )
    39954000             CALL netcdf_handle_error( 'netcdf_define_header', 218 )
    39964001
    39974002             WRITE ( time_average_text,'(F7.1,'' s avg'')' ) averaging_interval_pr
    3998              nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'time_avg', &
     4003             nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'time_avg',       &
    39994004                                     TRIM( time_average_text ) )
    40004005          ELSE
    4001              nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'title', &
     4006             nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'title',          &
    40024007                                     TRIM( run_description_header ) )
    40034008          ENDIF
     
    40924097             DO  j = 1, cross_profiles_count
    40934098
    4094                 IF ( TRIM(cross_profiles_char(j)) == TRIM(data_output_pr(i))) &
     4099                IF ( TRIM(cross_profiles_char(j)) == TRIM(data_output_pr(i)))  &
    40954100                THEN
    40964101                   EXIT
     
    41014106                   cross_profiles_count = cross_profiles_count + 1
    41024107                   cross_profiles_maxi  = cross_profiles_maxi  + 1
    4103                    cross_profiles_char(MIN( crmax, cross_profiles_count )) =  &
     4108                   cross_profiles_char(MIN( crmax, cross_profiles_count )) =   &
    41044109                                                      TRIM( data_output_pr(i) )
    4105                    cross_profiles_numb(MIN( crmax, cross_profiles_count )) =  &
     4110                   cross_profiles_numb(MIN( crmax, cross_profiles_count )) =   &
    41064111                                                      cross_profiles_maxi
    41074112                ENDIF
     
    41114116
    41124117          IF ( cross_profiles_count >= crmax )  THEN
    4113              message_string = 'It is not allowed to arrange more than '     &
    4114                               // '100 profiles with & cross_profiles. Apart'&
    4115                               // ' from that, all profiles are saved & to ' &
     4118             message_string = 'It is not allowed to arrange more than '        &
     4119                              // '100 profiles with cross_profiles. Apart'     &
     4120                              // ' from that, all profiles are saved to '      &
    41164121                              // 'the netCDF file.'
    41174122             CALL message( 'define_netcdf_header', 'PA0354', 0, 0, 0, 6, 0 )
     
    42854290
    42864291          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    4287              message_string = 'netCDF file for vertical profiles ' //        &
    4288                               'from previous run found,' //                  &
    4289                               '& but this file cannot be extended due to' // &
    4290                               ' variable mismatch.' //                       &
    4291                               '&New file is created instead.'
     4292             message_string = 'netCDF file for vertical profiles ' //          &
     4293                              'from previous run found,' //                    &
     4294                              ' but this file cannot be extended due to' //    &
     4295                              ' variable mismatch.' //                         &
     4296                              ' New file is created instead.'
    42924297             CALL message( 'define_netcdf_header', 'PA0254', 0, 1, 0, 6, 0 )
    42934298             extend = .FALSE.
     
    43194324
    43204325          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    4321              message_string = 'netCDF file for vertical profiles ' //       &
    4322                               'from previous run found,' //                 &
    4323                               '&but this file cannot be extended becaus' // &
    4324                               'e the current output time' //                &
    4325                               '&is less or equal than the last output t' // &
    4326                               'ime on this file.' //                        &
    4327                               '&New file is created instead.'
     4326             message_string = 'netCDF file for vertical profiles ' //          &
     4327                              'from previous run found,' //                    &
     4328                              ' but this file cannot be extended becaus' //    &
     4329                              'e the current output time' //                   &
     4330                              ' is less or equal than the last output t' //    &
     4331                              'ime on this file.' //                           &
     4332                              ' New file is created instead.'
    43284333             CALL message( 'define_netcdf_header', 'PA0255', 0, 1, 0, 6, 0 )
    43294334             dopr_time_count = 0
     
    43394344 
    43404345             IF ( statistic_regions == 0 )  THEN
    4341                 nc_stat = NF90_INQ_VARID( id_set_pr, data_output_pr(i), &
     4346                nc_stat = NF90_INQ_VARID( id_set_pr, data_output_pr(i),        &
    43424347                                          id_var_dopr(i,0) )
    43434348                CALL netcdf_handle_error( 'netcdf_define_header', 245 )
     
    43464351                   WRITE ( suffix, '(''_'',I2.2)' )  j
    43474352                   netcdf_var_name = TRIM( data_output_pr(i) ) // suffix
    4348                    nc_stat = NF90_INQ_VARID( id_set_pr, netcdf_var_name, &
     4353                   nc_stat = NF90_INQ_VARID( id_set_pr, netcdf_var_name,       &
    43494354                                             id_var_dopr(i,j) )
    43504355                   CALL netcdf_handle_error( 'netcdf_define_header', 246 )
     
    43784383          nc_stat = NF90_REDEF( id_set_pr )
    43794384          CALL netcdf_handle_error( 'netcdf_define_header', 437 )
    4380           nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'title', &
    4381                                   TRIM( run_description_header ) //    &
     4385          nc_stat = NF90_PUT_ATT( id_set_pr, NF90_GLOBAL, 'title',             &
     4386                                  TRIM( run_description_header ) //            &
    43824387                                  TRIM( time_average_text ) )
    43834388          CALL netcdf_handle_error( 'netcdf_define_header', 248 )
     
    43854390          nc_stat = NF90_ENDDEF( id_set_pr )
    43864391          CALL netcdf_handle_error( 'netcdf_define_header', 438 )
    4387           message_string = 'netCDF file for vertical profiles ' // &
    4388                            'from previous run found.' //           &
    4389                            '&This file will be extended.'
     4392          message_string = 'netCDF file for vertical profiles ' //             &
     4393                           'from previous run found.' //                       &
     4394                           ' This file will be extended.'
    43904395          CALL message( 'define_netcdf_header', 'PA0256', 0, 0, 0, 6, 0 )
    43914396
     
    43954400!
    43964401!--       Define some global attributes of the dataset
    4397           nc_stat = NF90_PUT_ATT( id_set_ts, NF90_GLOBAL, 'title', &
     4402          nc_stat = NF90_PUT_ATT( id_set_ts, NF90_GLOBAL, 'title',             &
    43984403                                  TRIM( run_description_header ) )
    43994404          CALL netcdf_handle_error( 'netcdf_define_header', 249 )
     
    44764481                DO  j = 0, statistic_regions
    44774482                   WRITE ( suffix, '(''_'',I2.2)' )  j
    4478                    var_list = TRIM( var_list ) // TRIM( dots_label(i) ) // &
     4483                   var_list = TRIM( var_list ) // TRIM( dots_label(i) ) //     &
    44794484                              suffix // ';'
    44804485                ENDDO
     
    44844489
    44854490          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    4486              message_string = 'netCDF file for time series ' //              &
    4487                               'from previous run found,' //                  &
    4488                               '& but this file cannot be extended due to' // &
    4489                               ' variable mismatch.' //                       &
    4490                               '&New file is created instead.'
     4491             message_string = 'netCDF file for time series ' //                &
     4492                              'from previous run found,' //                    &
     4493                              ' but this file cannot be extended due to' //    &
     4494                              ' variable mismatch.' //                         &
     4495                              ' New file is created instead.'
    44914496             CALL message( 'define_netcdf_header', 'PA0257', 0, 1, 0, 6, 0 )
    44924497             extend = .FALSE.
     
    45024507          CALL netcdf_handle_error( 'netcdf_define_header', 261 )
    45034508
    4504           nc_stat = NF90_INQUIRE_VARIABLE( id_set_ts, id_var_time_ts, &
     4509          nc_stat = NF90_INQUIRE_VARIABLE( id_set_ts, id_var_time_ts,          &
    45054510                                           dimids = id_dim_time_old )
    45064511          CALL netcdf_handle_error( 'netcdf_define_header', 262 )
    45074512          id_dim_time_ts = id_dim_time_old(1)
    45084513
    4509           nc_stat = NF90_INQUIRE_DIMENSION( id_set_ts, id_dim_time_ts, &
     4514          nc_stat = NF90_INQUIRE_DIMENSION( id_set_ts, id_dim_time_ts,         &
    45104515                                            len = dots_time_count )
    45114516          CALL netcdf_handle_error( 'netcdf_define_header', 263 )
    45124517
    4513           nc_stat = NF90_GET_VAR( id_set_ts, id_var_time_ts,        &
    4514                                   last_time_coordinate,             &
    4515                                   start = (/ dots_time_count /), &
     4518          nc_stat = NF90_GET_VAR( id_set_ts, id_var_time_ts,                   &
     4519                                  last_time_coordinate,                        &
     4520                                  start = (/ dots_time_count /),               &
    45164521                                  count = (/ 1 /) )
    45174522          CALL netcdf_handle_error( 'netcdf_define_header', 264 )
    45184523
    45194524          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    4520              message_string = 'netCDF file for time series ' //             &
    4521                               'from previous run found,' //                 &
    4522                               '&but this file cannot be extended becaus' // &
    4523                               'e the current output time' //                &
    4524                               '&is less or equal than the last output t' // &
    4525                               'ime on this file.' //                        &
    4526                               '&New file is created instead.'
     4525             message_string = 'netCDF file for time series ' //                &
     4526                              'from previous run found,' //                    &
     4527                              ' but this file cannot be extended becaus' //    &
     4528                              'e the current output time' //                   &
     4529                              ' is less or equal than the last output t' //    &
     4530                              'ime on this file.' //                           &
     4531                              ' New file is created instead.'
    45274532             CALL message( 'define_netcdf_header', 'PA0258', 0, 1, 0, 6, 0 )
    45284533             dots_time_count = 0
     
    45624567          nc_stat = NF90_REDEF( id_set_ts )
    45634568          CALL netcdf_handle_error( 'netcdf_define_header', 439 )
    4564           nc_stat = NF90_PUT_ATT( id_set_ts, NF90_GLOBAL, 'title', &
     4569          nc_stat = NF90_PUT_ATT( id_set_ts, NF90_GLOBAL, 'title',             &
    45654570                                  TRIM( run_description_header ) )
    45664571          CALL netcdf_handle_error( 'netcdf_define_header', 267 )
    45674572          nc_stat = NF90_ENDDEF( id_set_ts )
    45684573          CALL netcdf_handle_error( 'netcdf_define_header', 440 )
    4569           message_string = 'netCDF file for time series ' // &
    4570                            'from previous run found.' //     &
    4571                            '&This file will be extended.'
     4574          message_string = 'netCDF file for time series ' //                   &
     4575                           'from previous run found.' //                       &
     4576                           ' This file will be extended.'
    45724577          CALL message( 'define_netcdf_header', 'PA0259', 0, 0, 0, 6, 0 )
    45734578
     
    45784583!--       Define some global attributes of the dataset
    45794584          IF ( averaging_interval_sp /= 0.0_wp )  THEN
    4580              WRITE (time_average_text,'('', '',F7.1,'' s average'')') &
     4585             WRITE (time_average_text,'('', '',F7.1,'' s average'')')          &
    45814586                                                            averaging_interval_sp
    4582              nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'title',  &
    4583                                      TRIM( run_description_header ) // &
     4587             nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'title',          &
     4588                                     TRIM( run_description_header ) //         &
    45844589                                     TRIM( time_average_text ) )
    45854590             CALL netcdf_handle_error( 'netcdf_define_header', 268 )
     
    48014806
    48024807          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    4803              message_string = 'netCDF file for spectra  ' //                 &
    4804                               'from previous run found,' //                  &
    4805                               '& but this file cannot be extended due to' // &
    4806                               ' variable mismatch.' //                       &
    4807                               '&New file is created instead.'
     4808             message_string = 'netCDF file for spectra  ' //                   &
     4809                              'from previous run found,' //                    &
     4810                              ' but this file cannot be extended due to' //    &
     4811                              ' variable mismatch.' //                         &
     4812                              ' New file is created instead.'
    48084813             CALL message( 'define_netcdf_header', 'PA0260', 0, 1, 0, 6, 0 )
    48094814             extend = .FALSE.
     
    48354840
    48364841          IF ( ns /= ns_old )  THEN
    4837              message_string = 'netCDF file for spectra ' //                 &
    4838                               ' from previous run found,' //                &
    4839                               '&but this file cannot be extended due to' // &
    4840                               ' mismatch in number of' //                   &
    4841                               '&vertical levels.' //                        &
    4842                               '&New file is created instead.'
     4842             message_string = 'netCDF file for spectra ' //                    &
     4843                              ' from previous run found,' //                   &
     4844                              ' but this file cannot be extended due to' //    &
     4845                              ' mismatch in number of' //                      &
     4846                              ' vertical levels.' //                           &
     4847                              ' New file is created instead.'
    48434848             CALL message( 'define_netcdf_header', 'PA0261', 0, 1, 0, 6, 0 )
    48444849             extend = .FALSE.
     
    48574862                message_string = 'netCDF file for spectra ' //                 &
    48584863                                 ' from previous run found,' //                &
    4859                                  '&but this file cannot be extended due to' // &
     4864                                 ' but this file cannot be extended due to' // &
    48604865                                 ' mismatch in heights of' //                  &
    4861                                  '&vertical levels.' //                        &
    4862                                  '&New file is created instead.'
     4866                                 ' vertical levels.' //                        &
     4867                                 ' New file is created instead.'
    48634868                CALL message( 'define_netcdf_header', 'PA0262', 0, 1, 0, 6, 0 )
    48644869                extend = .FALSE.
     
    48934898
    48944899          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    4895              message_string = 'netCDF file for spectra ' //                 &
    4896                               'from previous run found,' //                 &
    4897                               '&but this file cannot be extended becaus' // &
    4898                               'e the current output time' //                &
    4899                               '&is less or equal than the last output t' // &
    4900                               'ime on this file.' //                        &
    4901                               '&New file is created instead.'
     4900             message_string = 'netCDF file for spectra ' //                    &
     4901                              'from previous run found,' //                    &
     4902                              ' but this file cannot be extended becaus' //    &
     4903                              'e the current output time' //                   &
     4904                              ' is less or equal than the last output t' //    &
     4905                              'ime on this file.' //                           &
     4906                              ' New file is created instead.'
    49024907             CALL message( 'define_netcdf_header', 'PA0263', 0, 1, 0, 6, 0 )
    49034908             dosp_time_count = 0
     
    49484953
    49494954             WRITE ( time_average_text,'(F7.1,'' s avg'')' )  averaging_interval_sp
    4950              nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'time_avg', &
     4955             nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'time_avg',       &
    49514956                                     TRIM( time_average_text ) )
    49524957          ELSE
    4953              nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'title', &
     4958             nc_stat = NF90_PUT_ATT( id_set_sp, NF90_GLOBAL, 'title',          &
    49544959                                     TRIM( run_description_header ) )
    49554960          ENDIF
     
    49574962          nc_stat = NF90_ENDDEF( id_set_sp )
    49584963          CALL netcdf_handle_error( 'netcdf_define_header', 442 )
    4959           message_string = 'netCDF file for spectra ' //     &
    4960                            'from previous run found.' //     &
    4961                            '&This file will be extended.'
     4964          message_string = 'netCDF file for spectra ' //                       &
     4965                           'from previous run found.' //                       &
     4966                           ' This file will be extended.'
    49624967          CALL message( 'define_netcdf_header', 'PA0264', 0, 0, 0, 6, 0 )
    49634968
    4964 
    4965        CASE ( 'pt_new' )
     4969!
     4970!--     Currently disabled (DATA_PRT_NETCDF)
     4971!       CASE ( 'pt_new' )
    49664972
    49674973!
    49684974!--       Define some global attributes of the dataset
    4969           nc_stat = NF90_PUT_ATT( id_set_prt, NF90_GLOBAL, 'title', &
    4970                                   TRIM( run_description_header ) )
    4971           CALL netcdf_handle_error( 'netcdf_define_header', 310 )
     4975!          nc_stat = NF90_PUT_ATT( id_set_prt, NF90_GLOBAL, 'title',            &
     4976!                                  TRIM( run_description_header ) )
     4977!          CALL netcdf_handle_error( 'netcdf_define_header', 310 )
    49724978
    49734979!
    49744980!--       Define time coordinate for particles (unlimited dimension)
    4975           CALL netcdf_create_dim( id_set_prt, 'time', NF90_UNLIMITED,          &
    4976                                   id_dim_time_prt, 311 )
    4977           CALL netcdf_create_var( id_set_prt, (/ id_dim_time_prt /), 'time',   &
    4978                                   NF90_DOUBLE, id_var_time_prt, 'seconds', '', &
    4979                                   312, 313, 000 )
     4981!          CALL netcdf_create_dim( id_set_prt, 'time', NF90_UNLIMITED,          &
     4982!                                  id_dim_time_prt, 311 )
     4983!          CALL netcdf_create_var( id_set_prt, (/ id_dim_time_prt /), 'time',   &
     4984!                                  NF90_DOUBLE, id_var_time_prt, 'seconds', '', &
     4985!                                  312, 313, 000 )
    49804986!
    49814987!--       netCDF4 allows more than one unlimited dimension
    4982           CALL netcdf_create_dim( id_set_prt, 'particle_number',            &
    4983                                   NF90_UNLIMITED, id_dim_prtnum, 314 )
    4984 
    4985           CALL netcdf_create_var( id_set_prt, (/ id_dim_prtnum /),             &
    4986                                   'particle_number', NF90_DOUBLE,              &
    4987                                   id_var_prtnum, 'particle number', '', 315,   &
    4988                                   316, 000 )
     4988!          CALL netcdf_create_dim( id_set_prt, 'particle_number',            &
     4989!                                  NF90_UNLIMITED, id_dim_prtnum, 314 )
     4990
     4991!          CALL netcdf_create_var( id_set_prt, (/ id_dim_prtnum /),             &
     4992!                                  'particle_number', NF90_DOUBLE,              &
     4993!                                  id_var_prtnum, 'particle number', '', 315,   &
     4994!                                  316, 000 )
    49894995!
    49904996!--       Define variable which contains the real number of particles in use
    4991           CALL netcdf_create_var( id_set_prt, (/ id_dim_time_prt /),           &
    4992                                   'real_num_of_prt', NF90_DOUBLE,              &
    4993                                   id_var_rnop_prt, 'particle number', '', 317, &
    4994                                   318, 000 )
     4997!          CALL netcdf_create_var( id_set_prt, (/ id_dim_time_prt /),           &
     4998!                                  'real_num_of_prt', NF90_DOUBLE,              &
     4999!                                  id_var_rnop_prt, 'particle number', '', 317, &
     5000!                                  318, 000 )
    49955001!
    49965002!--       Define the variables
    4997           DO  i = 1, 17
    4998 
    4999              CALL netcdf_create_var( id_set_prt, (/ id_dim_prtnum,             &
    5000                                      id_dim_time_prt /), prt_var_names(i),     &
    5001                                      nc_precision(8), id_var_prt(i),           &
    5002                                      TRIM( prt_var_units(i) ),                 &
    5003                                      TRIM( prt_var_names(i) ), 319, 320, 321 )
    5004 
    5005           ENDDO
     5003!          DO  i = 1, 17
     5004
     5005!             CALL netcdf_create_var( id_set_prt, (/ id_dim_prtnum,             &
     5006!                                     id_dim_time_prt /), prt_var_names(i),     &
     5007!                                     nc_precision(8), id_var_prt(i),           &
     5008!                                     TRIM( prt_var_units(i) ),                 &
     5009!                                     TRIM( prt_var_names(i) ), 319, 320, 321 )
     5010
     5011!          ENDDO
    50065012
    50075013!
    50085014!--       Leave netCDF define mode
    5009           nc_stat = NF90_ENDDEF( id_set_prt )
    5010           CALL netcdf_handle_error( 'netcdf_define_header', 322 )
    5011 
    5012 
    5013        CASE ( 'pt_ext' )
     5015!          nc_stat = NF90_ENDDEF( id_set_prt )
     5016!          CALL netcdf_handle_error( 'netcdf_define_header', 322 )
     5017
     5018!
     5019!--     Currently disabled (DATA_PRT_NETCDF)
     5020!       CASE ( 'pt_ext' )
    50145021
    50155022!
     
    50185025!--       The current time must be larger than the last output time
    50195026!--       on the file.
    5020           nc_stat = NF90_INQ_VARID( id_set_prt, 'time', id_var_time_prt )
    5021           CALL netcdf_handle_error( 'netcdf_define_header', 323 )
    5022 
    5023           nc_stat = NF90_INQUIRE_VARIABLE( id_set_prt, id_var_time_prt, &
    5024                                            dimids = id_dim_time_old )
    5025           CALL netcdf_handle_error( 'netcdf_define_header', 324 )
    5026           id_dim_time_prt = id_dim_time_old(1)
    5027 
    5028           nc_stat = NF90_INQUIRE_DIMENSION( id_set_prt, id_dim_time_prt, &
    5029                                             len = prt_time_count )
    5030           CALL netcdf_handle_error( 'netcdf_define_header', 325 )
    5031 
    5032           nc_stat = NF90_GET_VAR( id_set_prt, id_var_time_prt,  &
    5033                                   last_time_coordinate,         &
    5034                                   start = (/ prt_time_count /), &
    5035                                   count = (/ 1 /) )
    5036           CALL netcdf_handle_error( 'netcdf_define_header', 326 )
    5037 
    5038           IF ( last_time_coordinate(1) >= simulated_time )  THEN
    5039              message_string = 'netCDF file for particles ' //               &
    5040                               'from previous run found,' //                 &
    5041                               '&but this file cannot be extended becaus' // &
    5042                               'e the current output time' //                &
    5043                               '&is less or equal than the last output t' // &
    5044                               'ime on this file.' //                        &
    5045                               '&New file is created instead.'
    5046              CALL message( 'define_netcdf_header', 'PA0265', 0, 1, 0, 6, 0 )
    5047              prt_time_count = 0
    5048              extend = .FALSE.
    5049              RETURN
    5050           ENDIF
     5027!          nc_stat = NF90_INQ_VARID( id_set_prt, 'time', id_var_time_prt )
     5028!          CALL netcdf_handle_error( 'netcdf_define_header', 323 )
     5029
     5030!          nc_stat = NF90_INQUIRE_VARIABLE( id_set_prt, id_var_time_prt, &
     5031!                                           dimids = id_dim_time_old )
     5032!          CALL netcdf_handle_error( 'netcdf_define_header', 324 )
     5033!          id_dim_time_prt = id_dim_time_old(1)
     5034
     5035!          nc_stat = NF90_INQUIRE_DIMENSION( id_set_prt, id_dim_time_prt, &
     5036!                                            len = prt_time_count )
     5037!          CALL netcdf_handle_error( 'netcdf_define_header', 325 )
     5038
     5039!          nc_stat = NF90_GET_VAR( id_set_prt, id_var_time_prt,  &
     5040!                                  last_time_coordinate,         &
     5041!                                  start = (/ prt_time_count /), &
     5042!                                  count = (/ 1 /) )
     5043!          CALL netcdf_handle_error( 'netcdf_define_header', 326 )
     5044
     5045!          IF ( last_time_coordinate(1) >= simulated_time )  THEN
     5046!             message_string = 'netCDF file for particles ' //                  &
     5047!                              'from previous run found,' //                    &
     5048!                              ' but this file cannot be extended becaus' //    &
     5049!                              'e the current output time' //                   &
     5050!                              ' is less or equal than the last output t' //    &
     5051!                              'ime on this file.' //                           &
     5052!                              ' New file is created instead.'
     5053!             CALL message( 'define_netcdf_header', 'PA0265', 0, 1, 0, 6, 0 )
     5054!             prt_time_count = 0
     5055!             extend = .FALSE.
     5056!             RETURN
     5057!          ENDIF
    50515058
    50525059!
    50535060!--       Dataset seems to be extendable.
    50545061!--       Now get the variable ids.
    5055           nc_stat = NF90_INQ_VARID( id_set_prt, 'real_num_of_prt', &
    5056                                     id_var_rnop_prt )
    5057           CALL netcdf_handle_error( 'netcdf_define_header', 327 )
    5058 
    5059           DO  i = 1, 17
    5060 
    5061              nc_stat = NF90_INQ_VARID( id_set_prt, prt_var_names(i), &
    5062                                        id_var_prt(i) )
    5063              CALL netcdf_handle_error( 'netcdf_define_header', 328 )
    5064 
    5065           ENDDO
    5066 
    5067           message_string = 'netCDF file for particles ' // &
    5068                            'from previous run found.' //   &
    5069                            '&This file will be extended.'
    5070           CALL message( 'define_netcdf_header', 'PA0266', 0, 0, 0, 6, 0 )
     5062!         nc_stat = NF90_INQ_VARID( id_set_prt, 'real_num_of_prt',            &
     5063!                                   id_var_rnop_prt )
     5064!         CALL netcdf_handle_error( 'netcdf_define_header', 327 )
     5065
     5066!          DO  i = 1, 17
     5067
     5068!             nc_stat = NF90_INQ_VARID( id_set_prt, prt_var_names(i),          &
     5069!                                       id_var_prt(i) )
     5070!             CALL netcdf_handle_error( 'netcdf_define_header', 328 )
     5071
     5072!          ENDDO
     5073
     5074!          message_string = 'netCDF file for particles ' //                    &
     5075!                           'from previous run found.' //                       &
     5076!                           ' This file will be extended.'
     5077!          CALL message( 'define_netcdf_header', 'PA0266', 0, 0, 0, 6, 0 )
    50715078       
    50725079
     
    50765083!
    50775084!--       Define some global attributes of the dataset
    5078           nc_stat = NF90_PUT_ATT( id_set_pts, NF90_GLOBAL, 'title', &
     5085          nc_stat = NF90_PUT_ATT( id_set_pts, NF90_GLOBAL, 'title',            &
    50795086                                  TRIM( run_description_header ) )
    50805087          CALL netcdf_handle_error( 'netcdf_define_header', 396 )
     
    51735180
    51745181          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    5175              message_string = 'netCDF file for particle time series ' //     &
    5176                               'from previous run found,' //                  &
    5177                               '& but this file cannot be extended due to' // &
    5178                               ' variable mismatch.' //                       &
    5179                               '&New file is created instead.'
     5182             message_string = 'netCDF file for particle time series ' //       &
     5183                              'from previous run found,' //                    &
     5184                              ' but this file cannot be extended due to' //    &
     5185                              ' variable mismatch.' //                         &
     5186                              ' New file is created instead.'
    51805187             CALL message( 'define_netcdf_header', 'PA0267', 0, 1, 0, 6, 0 )
    51815188             extend = .FALSE.
     
    52075214
    52085215          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    5209              message_string = 'netCDF file for particle time series ' //    &
    5210                               'from previous run found,' //                 &
    5211                               '&but this file cannot be extended becaus' // &
    5212                               'e the current output time' //                &
    5213                               '&is less or equal than the last output t' // &
    5214                               'ime on this file.' //                        &
    5215                               '&New file is created instead.'
     5216             message_string = 'netCDF file for particle time series ' //       &
     5217                              'from previous run found,' //                    &
     5218                              ' but this file cannot be extended becaus' //    &
     5219                              'e the current output time' //                   &
     5220                              ' is less or equal than the last output t' //    &
     5221                              'ime on this file.' //                           &
     5222                              ' New file is created instead.'
    52165223             CALL message( 'define_netcdf_header', 'PA0268', 0, 1, 0, 6, 0 )
    52175224             dopts_time_count = 0
     
    52555262          nc_stat = NF90_REDEF( id_set_pts )
    52565263          CALL netcdf_handle_error( 'netcdf_define_header', 443 )
    5257           nc_stat = NF90_PUT_ATT( id_set_pts, NF90_GLOBAL, 'title', &
     5264          nc_stat = NF90_PUT_ATT( id_set_pts, NF90_GLOBAL, 'title',            &
    52585265                                  TRIM( run_description_header ) )
    52595266          CALL netcdf_handle_error( 'netcdf_define_header', 411 )
    52605267          nc_stat = NF90_ENDDEF( id_set_pts )
    52615268          CALL netcdf_handle_error( 'netcdf_define_header', 444 )
    5262           message_string = 'netCDF file for particle time series ' // &
    5263                            'from previous run found.' //              &
    5264                            '&This file will be extended.'
     5269          message_string = 'netCDF file for particle time series ' //          &
     5270                           'from previous run found.' //                       &
     5271                           ' This file will be extended.'
    52655272          CALL message( 'netcdf_define_header', 'PA0269', 0, 0, 0, 6, 0 )
    52665273
     
    53555362             message_string = 'netCDF file for flight time series ' //         &
    53565363                              'from previous run found,' //                    &
    5357                               '& but this file cannot be extended due to' //   &
     5364                              ' but this file cannot be extended due to' //    &
    53585365                              ' variable mismatch.' //                         &
    5359                               '&New file is created instead.'
     5366                              ' New file is created instead.'
    53605367             CALL message( 'define_netcdf_header', 'PA0257', 0, 1, 0, 6, 0 )
    53615368             extend = .FALSE.
     
    53715378          CALL netcdf_handle_error( 'netcdf_define_header', 261 )
    53725379
    5373           nc_stat = NF90_INQUIRE_VARIABLE( id_set_fl, id_var_time_fl, &
     5380          nc_stat = NF90_INQUIRE_VARIABLE( id_set_fl, id_var_time_fl,          &
    53745381                                           dimids = id_dim_time_old )
    53755382          CALL netcdf_handle_error( 'netcdf_define_header', 262 )
    53765383          id_dim_time_fl = id_dim_time_old(1)
    53775384
    5378           nc_stat = NF90_INQUIRE_DIMENSION( id_set_fl, id_dim_time_fl, &
     5385          nc_stat = NF90_INQUIRE_DIMENSION( id_set_fl, id_dim_time_fl,         &
    53795386                                            len = dofl_time_count )
    53805387          CALL netcdf_handle_error( 'netcdf_define_header', 263 )
     
    53875394
    53885395          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    5389              message_string = 'netCDF file for flight-time series ' //      &
    5390                               'from previous run found,' //                 &
    5391                               '&but this file cannot be extended becaus' // &
    5392                               'e the current output time' //                &
    5393                               '&is less or equal than the last output t' // &
    5394                               'ime on this file.' //                        &
    5395                               '&New file is created instead.'
     5396             message_string = 'netCDF file for flight-time series ' //         &
     5397                              'from previous run found,' //                    &
     5398                              ' but this file cannot be extended becaus' //    &
     5399                              'e the current output time' //                   &
     5400                              ' is less or equal than the last output t' //    &
     5401                              'ime on this file.' //                           &
     5402                              ' New file is created instead.'
    53965403             CALL message( 'define_netcdf_header', 'PA0258', 0, 1, 0, 6, 0 )
    53975404             dofl_time_count = 0
     
    54415448          message_string = 'netCDF file for flight-time series ' //            &
    54425449                           'from previous run found.' //                       &
    5443                            '&This file will be extended.'
     5450                           ' This file will be extended.'
    54445451          CALL message( 'define_netcdf_header', 'PA0259', 0, 0, 0, 6, 0 )
    54455452
  • palm/trunk/SOURCE/parin.f90

    r2995 r3045  
    2525! -----------------
    2626! $Id$
     27! z_max_do2d removed, error messages revised
     28!
     29! 2995 2018-04-19 12:13:16Z Giersch
    2730! time_since_reference_point must be calculated/initialized before the first 
    2831! call of functions related to the radiation model which occur in
     
    639642             skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz,          &
    640643             skip_time_do3d, skip_time_domask, synchronous_exchange,           &
    641              termination_time_needed, vnest_start_time, z_max_do2d
     644             termination_time_needed, vnest_start_time
    642645
    643646    NAMELIST /runtime_parameters/  averaging_interval, averaging_interval_pr,  &
     
    660663             skip_time_do2d_xy, skip_time_do2d_xz, skip_time_do2d_yz,          &
    661664             skip_time_do3d, skip_time_domask, synchronous_exchange,           &
    662              termination_time_needed, vnest_start_time, z_max_do2d
     665             termination_time_needed, vnest_start_time
    663666
    664667    NAMELIST /envpar/  batch_job, host, local_dvrserver_running,               &
     
    669672!
    670673!-- First read values of environment variables (this NAMELIST file is
    671 !-- generated by mrun)
     674!-- generated by palmrun)
    672675    CALL location_message( 'reading environment parameters from ENVPAR', .FALSE. )
    673676
     
    676679    IF ( ioerr /= 0 )  THEN
    677680       message_string = 'local file ENVPAR not found' //                       &
    678                         '&some variables for steering may not be properly set'
     681                        ' some variables for steering may not be properly set'
    679682       CALL message( 'parin', 'PA0276', 0, 1, 0, 6, 0 )
    680683    ELSE
     
    682685       IF ( ioerr < 0 )  THEN
    683686          message_string = 'no envpar-NAMELIST found in local file '  //       &
    684                            'ENVPAR&some variables for steering may '  //       &
     687                           'ENVPAR or some variables for steering may '  //    &
    685688                           'not be properly set'
    686689          CALL message( 'parin', 'PA0278', 0, 1, 0, 6, 0 )
    687690       ELSEIF ( ioerr > 0 )  THEN
    688691          message_string = 'errors in local file ENVPAR' //                    &
    689                            '&some variables for steering may not be properly set'
     692                           ' some variables for steering may not be properly set'
    690693          CALL message( 'parin', 'PA0277', 0, 1, 0, 6, 0 )
    691694       ENDIF
     
    742745          GOTO 12
    743746
    744  10       message_string = 'errors in initialization_parameters &or no ' //    &
    745                            'initialization_parameters-namelist ' //            &
     747 10       message_string = 'errors in &initialization_parameters or no ' //    &
     748                           '&initialization_parameters-namelist ' //           &
    746749                           'found (CRAY-machines only)'
    747750          CALL message( 'parin', 'PA0271', 1, 2, 0, 6, 0 )
     
    750753          READ ( 11, inipar, ERR=13, END=14 )
    751754 
    752           message_string = 'namelist inipar is deprecated and will be ' //     &
    753                           'removed in near future. Please &use namelist ' //   &
    754                           'initialization_parameters instead'
    755           CALL message( 'parin', 'PA0272', 0, 1, 0, 6, 0 )
     755          message_string = 'namelist &inipar is deprecated and will be ' //    &
     756                          'removed in near future. Please use namelist ' //    &
     757                          '&initialization_parameters instead'
     758          CALL message( 'parin', 'PA0017', 0, 1, 0, 6, 0 )
    756759 
    757760          GOTO 12
    758761 
    759  13       message_string = 'errors in inipar &or no inipar-namelist ' //       &
     762 13       message_string = 'errors in inipar  or no inipar-namelist ' //       &
    760763                           'found (CRAY-machines only)'
    761764          CALL message( 'parin', 'PA0271', 1, 2, 0, 6, 0 )
    762765         
    763  14       message_string = 'no initialization_parameters-namelist found'
     766 14       message_string = 'no &initialization_parameters-namelist found'
    764767          CALL message( 'parin', 'PA0272', 1, 2, 0, 6, 0 )
    765768
  • palm/trunk/SOURCE/plant_canopy_model_mod.f90

    r3022 r3045  
    2525! -----------------
    2626! $Id$
     27! Error message revised
     28!
     29! 3022 2018-05-18 11:12:35Z suehring
    2730! Bugfix in allocation of transpiration rate
    2831!
     
    359362       IF ( canopy_drag_coeff == 0.0_wp )  THEN
    360363          message_string = 'plant_canopy = .TRUE. requires a non-zero drag '// &
    361                            'coefficient & given value is canopy_drag_coeff = 0.0'
     364                           'coefficient, given value is canopy_drag_coeff = 0.0'
    362365          CALL message( 'pcm_check_parameters', 'PA0041', 1, 2, 0, 6, 0 )
    363366       ENDIF
    364367   
    365        IF ( ( alpha_lad /= 9999999.9_wp  .AND.  beta_lad == 9999999.9_wp )  .OR.&
     368       IF ( ( alpha_lad /= 9999999.9_wp  .AND.  beta_lad == 9999999.9_wp ) .OR.&
    366369              beta_lad /= 9999999.9_wp   .AND.  alpha_lad == 9999999.9_wp )  THEN
    367370          message_string = 'using the beta function for the construction ' //  &
  • palm/trunk/SOURCE/pmc_interface_mod.f90

    r3022 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 3022 2018-05-18 11:12:35Z suehring
    2730! Minor fix - working precision added to real number
    2831!
     
    41074110       IF ( end_time /= end_time_root )  THEN
    41084111          WRITE( message_string, * )  'mismatch between root model and ',      &
    4109                'child settings &   end_time(root) = ', end_time_root,          &
    4110                ' &   end_time(child) = ', end_time, ' & child value is set',   &
     4112               'child settings: end_time(root) = ', end_time_root,             &
     4113               ' end_time(child) = ', end_time, ' child value is set',         &
    41114114               ' to root value'
    41124115          CALL message( 'pmci_check_setting_mismatches', 'PA0419', 0, 1, 0, 6, &
     
    41234126       IF ( restart_time /= restart_time_root )  THEN
    41244127          WRITE( message_string, * )  'mismatch between root model and ',      &
    4125                'child settings &   restart_time(root) = ', restart_time_root,  &
    4126                ' &   restart_time(child) = ', restart_time, ' & child ',       &
     4128               'child settings: restart_time(root) = ', restart_time_root,     &
     4129               ' restart_time(child) = ', restart_time, ' child ',             &
    41274130               'value is set to root value'
    41284131          CALL message( 'pmci_check_setting_mismatches', 'PA0419', 0, 1, 0, 6, &
     
    41394142       IF ( dt_restart /= dt_restart_root )  THEN
    41404143          WRITE( message_string, * )  'mismatch between root model and ',      &
    4141                'child settings &   dt_restart(root) = ', dt_restart_root,      &
    4142                ' &   dt_restart(child) = ', dt_restart, ' & child ',           &
     4144               'child settings: dt_restart(root) = ', dt_restart_root,          &
     4145               ' dt_restart(child) = ', dt_restart, ' child ',                 &
    41434146               'value is set to root value'
    41444147          CALL message( 'pmci_check_setting_mismatches', 'PA0419', 0, 1, 0, 6, &
     
    41554158       IF ( time_restart /= time_restart_root )  THEN
    41564159          WRITE( message_string, * )  'mismatch between root model and ',      &
    4157                'child settings &  time_restart(root) = ', time_restart_root,  &
    4158                ' &   time_restart(child) = ', time_restart, ' & child ',       &
     4160               'child settings: time_restart(root) = ', time_restart_root,  &
     4161               ' time_restart(child) = ', time_restart, ' child ',             &
    41594162               'value is set to root value'
    41604163          CALL message( 'pmci_check_setting_mismatches', 'PA0419', 0, 1, 0, 6, &
  • palm/trunk/SOURCE/radiation_model_mod.f90

    r3026 r3045  
    2828! -----------------
    2929! $Id$
     30! Error message revised
     31!
     32! 3026 2018-05-22 10:30:53Z schwenkel
    3033! Changed the name specific humidity to mixing ratio, since we are computing
    3134! mixing ratios.
     
    13321335          IF ( albedo_type == 0  .AND.  albedo == 9999999.9_wp  .AND.          &
    13331336               radiation_scheme == 'clear-sky')  THEN
    1334              message_string = 'radiation_scheme = "clear-sky" in combination' //&
    1335                               'with albedo_type = 0 requires setting of albedo'//&
    1336                               ' /= 9999999.9'
     1337             message_string = 'radiation_scheme = "clear-sky" in combination'//&
     1338                              'with albedo_type = 0 requires setting of'//     &
     1339                              'albedo /= 9999999.9'
    13371340             CALL message( 'check_parameters', 'PA0410', 1, 2, 0, 6, 0 )
    13381341          ENDIF
  • palm/trunk/SOURCE/read_restart_data_mod.f90

    r3004 r3045  
    2525! -----------------
    2626! $Id$
     27! Error messages revised
     28!
     29! 3004 2018-04-27 12:33:25Z Giersch
    2730! precipitation_rate_av removed
    2831!
     
    164167       binary_version_global = '4.7'
    165168       IF ( TRIM( version_on_file ) /= TRIM( binary_version_global ) )  THEN
    166           WRITE( message_string, * ) 'version mismatch concerning control ',   &
    167                                      'variables',                              &
    168                                      '&version on file    = "',                &
     169          WRITE( message_string, * ) 'version mismatch concerning ',           &
     170                                     'binary_version_global:',                 &
     171                                     ' version on file    = "',                &
    169172                                     TRIM( version_on_file ), '"',             &
    170                                      '&version on program = "',                &
     173                                     ' version on program = "',                &
    171174                                     TRIM( binary_version_global ), '"'
    172175          CALL message( 'rrd_global', 'PA0296', 1, 2, 0, 6, 0 )
     
    887890       IF ( nz_on_file /= nz )  THEN
    888891          WRITE( message_string, * ) 'mismatch concerning number of ',         &
    889                                      'gridpoints along z',                     &
    890                                      '&nz on file    = "', nz_on_file, '"',    &
    891                                      '&nz from run   = "', nz, '"'
     892                                     'gridpoints along z:',                    &
     893                                     ' nz on file    = "', nz_on_file, '"',    &
     894                                     ' nz from run   = "', nz, '"'
    892895          CALL message( 'rrd_read_parts_of_global', 'PA0304', 1, 2, 0, 6, 0 )
    893896       ENDIF
     
    904907          WRITE( message_string, * ) 'number of user profiles on res',         &
    905908                                     'tart data file differs from the ',       &
    906                                      'current run&max_pr_user on file    = "', &
     909                                     'current run: max_pr_user on file    = "',&
    907910                                     max_pr_user_on_file, '"',                 &
    908                                      '&max_pr_user from run   = "',            &
     911                                     ' max_pr_user from run   = "',            &
    909912                                     max_pr_user, '"'
    910913          CALL message( 'rrd_read_parts_of_global', 'PA0306', 0, 0, 0, 6, 0 )
     
    924927       IF ( statistic_regions_on_file /= statistic_regions )  THEN
    925928          WRITE( message_string, * ) 'statistic regions on restart data file ',&
    926                                      'differ from the current run',            &
    927                                      '&statistic regions on file    = "',      &
     929                                     'differ from the current run:',           &
     930                                     ' statistic regions on file    = "',      &
    928931                                     statistic_regions_on_file, '"',           &
    929                                      '&statistic regions from run   = "',      &
     932                                     ' statistic regions from run   = "',      &
    930933                                      statistic_regions, '"',                  &
    931                                      '&statistic data may be lost!'
     934                                     ' statistic data may be lost!'
    932935          CALL message( 'rrd_read_parts_of_global', 'PA0308', 0, 1, 0, 6, 0 )
    933936          tmp_sr = MIN( statistic_regions_on_file, statistic_regions )
     
    949952                IF ( average_count_pr /= 0 )  THEN
    950953                   WRITE( message_string, * ) 'inflow profiles not ',          &
    951                                   'temporally averaged. &Averaging will be ',  &
     954                                  'temporally averaged.  Averaging will be ',  &
    952955                                  'done now using', average_count_pr,          &
    953956                                  ' samples.'
     
    12421245    THEN
    12431246       WRITE( message_string, * ) 'number of PEs or virtual PE-grid changed ', &
    1244                         'in restart run&  PE', myid, ' will read from files ', &
     1247                        'in restart run  PE', myid, ' will read from files ', &
    12451248                         file_list(1:files_to_be_opened)
    12461249       CALL message( 'rrd_local', 'PA0285', 0, 0, 0, 6, 0 )
     
    12711274       binary_version_local = '4.7'
    12721275       IF ( TRIM( version_on_file ) /= TRIM( binary_version_local ) )  THEN
    1273           WRITE( message_string, * ) 'version mismatch concerning data ',      &
    1274                       'from prior run',                                        &
    1275                       '&version on file    = "', TRIM( version_on_file ), '"', &
    1276                       '&version in program = "', TRIM( binary_version_local ), '"'
     1276          WRITE( message_string, * ) 'version mismatch concerning ',           &
     1277                      'binary_version_local:',                                 &
     1278                      ' version on file    = "', TRIM( version_on_file ), '"', &
     1279                      ' version in program = "', TRIM( binary_version_local ), '"'
    12771280          CALL message( 'rrd_local', 'PA0286', 1, 2, 0, 6, 0 )
    12781281       ENDIF
     
    12871290          WRITE( message_string, * ) 'problem with index bound nxl on ',       &
    12881291                            'restart file "', myid_char, '"',                  &
    1289                             '&nxl = ', nxl_on_file, ' but it should be',       &
    1290                             '&= ', hor_index_bounds_previous_run(1,j),         &
    1291                             '&from the index bound information array'
     1292                            ' nxl = ', nxl_on_file, ' but it should be',       &
     1293                            ' = ', hor_index_bounds_previous_run(1,j),         &
     1294                            ' from the index bound information array'
    12921295          CALL message( 'rrd_local', 'PA0287', 2, 2, -1, 6, 1 )
    12931296       ENDIF
     
    12961299           WRITE( message_string, * ) 'problem with index bound nxr on ',      &
    12971300                               'restart file "', myid_char, '"'  ,             &
    1298                                '&nxr = ', nxr_on_file, ' but it should be',    &
    1299                                '&= ', hor_index_bounds_previous_run(2,j),      &
    1300                                '&from the index bound information array'
     1301                               ' nxr = ', nxr_on_file, ' but it should be',    &
     1302                               ' = ', hor_index_bounds_previous_run(2,j),      &
     1303                               ' from the index bound information array'
    13011304          CALL message( 'rrd_local', 'PA0288', 2, 2, -1, 6, 1 )
    13021305
     
    13061309          WRITE( message_string, * ) 'problem with index bound nys on ',       &
    13071310                                 'restart file "', myid_char, '"',             &
    1308                                  '&nys = ', nys_on_file, ' but it should be',  &
    1309                                  '&= ', hor_index_bounds_previous_run(3,j),    &
    1310                                      '&from the index bound information array'
     1311                                 ' nys = ', nys_on_file, ' but it should be',  &
     1312                                 ' = ', hor_index_bounds_previous_run(3,j),    &
     1313                                 ' from the index bound information array'
    13111314          CALL message( 'rrd_local', 'PA0289', 2, 2, -1, 6, 1 )
    13121315       ENDIF
     
    13151318          WRITE( message_string, * ) 'problem with index bound nyn on ',       &
    13161319                               'restart file "', myid_char, '"',               &
    1317                                '&nyn = ', nyn_on_file, ' but it should be',    &
    1318                                '&= ', hor_index_bounds_previous_run(4,j),      &
    1319                                '&from the index bound information array'
     1320                               ' nyn = ', nyn_on_file, ' but it should be',    &
     1321                               ' = ', hor_index_bounds_previous_run(4,j),      &
     1322                               ' from the index bound information array'
    13201323          CALL message( 'rrd_local', 'PA0290', 2, 2, -1, 6, 1 )
    13211324       ENDIF
     
    13231326       IF ( nzb_on_file /= nzb )  THEN
    13241327          WRITE( message_string, * ) 'mismatch between actual data and data ', &
    1325                                      '&from prior run on PE ', myid,           &
    1326                                      '&nzb on file = ', nzb_on_file,           &
    1327                                      '&nzb         = ', nzb
     1328                                     'from prior run on PE ', myid,            &
     1329                                     ' nzb on file = ', nzb_on_file,           &
     1330                                     ' nzb         = ', nzb
    13281331          CALL message( 'rrd_local', 'PA0291', 1, 2, 0, 6, 0 ) 
    13291332       ENDIF
     
    13311334       IF ( nzt_on_file /= nzt )  THEN
    13321335          WRITE( message_string, * ) 'mismatch between actual data and data ', &
    1333                                      '&from prior run on PE ', myid,           &
    1334                                      '&nzt on file = ', nzt_on_file,           &
    1335                                      '&nzt         = ', nzt
     1336                                     'from prior run on PE ', myid,            &
     1337                                     ' nzt on file = ', nzt_on_file,           &
     1338                                     ' nzt         = ', nzt
    13361339          CALL message( 'rrd_local', 'PA0292', 1, 2, 0, 6, 0 ) 
    13371340       ENDIF
  • palm/trunk/SOURCE/spectra_mod.f90

    r2956 r3045  
    2525! -----------------
    2626! $Id$
     27! Error message revised
     28!
     29! 2956 2018-04-10 11:01:03Z Giersch
    2730! spectrum_x and spectrum_y will only be allocated if they are not allocated
    2831! before (e.g. in case of restart runs)
     
    472475
    473476                message_string = 'non-cyclic lateral boundaries along x do'//  &
    474                                  ' not & allow calculation of spectra along x'
     477                                 ' not allow calculation of spectra along x'
    475478                CALL message( 'calc_spectra', 'PA0160', 1, 2, 0, 6, 0 )
    476479             ENDIF
     
    488491#else
    489492             message_string = 'sorry, calculation of spectra in non paral' //  &
    490                               'lel mode& is still not realized'
     493                              'lel mode is still not realized'
    491494             CALL message( 'calc_spectra', 'PA0161', 1, 2, 0, 6, 0 )
    492495#endif
     
    504507                IF ( myid == 0 )  THEN
    505508                   message_string = 'non-cyclic lateral boundaries along y' // &
    506                                     ' do not & allow calculation of spectr' // &
    507                                     'a along y'
     509                                    ' do not allow calculation of spectra' // &
     510                                    ' along y'
    508511                   CALL message( 'calc_spectra', 'PA0162', 1, 2, 0, 6, 0 )
    509512                ENDIF
     
    518521#else
    519522             message_string = 'sorry, calculation of spectra in non paral' //  &
    520                               'lel mode& is still not realized'
     523                              'lel mode is still not realized'
    521524             CALL message( 'calc_spectra', 'PA0161', 1, 2, 0, 6, 0 )
    522525#endif
  • palm/trunk/SOURCE/subsidence_mod.f90

    r2718 r3045  
    2525! -----------------
    2626! $Id$
     27! Error message revised
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    145148
    146149       IF ( ocean )  THEN
    147           message_string = 'Applying large scale vertical motion is not ' // &
     150          message_string = 'Applying large scale vertical motion is not ' //   &
    148151                           'allowed for ocean runs'
    149152          CALL message( 'init_w_subsidence', 'PA0324', 2, 2, 0, 6, 0 )
  • palm/trunk/SOURCE/surface_coupler.f90

    r2718 r3045  
    2525! -----------------
    2626! $Id$
     27! Error message revised
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    167170
    168171    IF ( terminate_coupled_remote > 0 )  THEN
    169        WRITE( message_string, * ) 'remote model "',                         &
    170                                   TRIM( coupling_mode_remote ),             &
    171                                   '" terminated',                           &
    172                                   '&with terminate_coupled_remote = ',      &
    173                                   terminate_coupled_remote,                 &
    174                                   '&local model  "', TRIM( coupling_mode ), &
    175                                   '" has',                                  &
    176                                   '&terminate_coupled = ',                  &
     172       WRITE( message_string, * ) 'remote model "',                            &
     173                                  TRIM( coupling_mode_remote ),                &
     174                                  '" terminated',                              &
     175                                  ' with terminate_coupled_remote = ',         &
     176                                  terminate_coupled_remote,                    &
     177                                  ' local model  "', TRIM( coupling_mode ),    &
     178                                  '" has',                                     &
     179                                  ' terminate_coupled = ',                     &
    177180                                   terminate_coupled
    178181       CALL message( 'surface_coupler', 'PA0310', 1, 2, 0, 6, 0 )
  • palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

    r2766 r3045  
    2626! -----------------
    2727! $Id$
     28! Error message revised
     29!
     30! 2766 2018-01-22 17:17:47Z kanani
    2831! Removed preprocessor directive __chem
    2932!
     
    833836                              'combination with a prescribed roughness '  //   &
    834837                              'heterogeneity'
    835              CALL message( 'surface_layer_fluxes', 'PA0417', 1, 2, 0, 6, 0 )
     838             CALL message( 'surface_layer_fluxes', 'PA0116', 1, 2, 0, 6, 0 )
    836839          ENDIF
    837840
  • palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90

    r3044 r3045  
    2525! -----------------
    2626! $Id$
     27! Error message revised
     28!
     29! 3044 2018-05-25 10:59:41Z gronemeier
    2730! Add missing variable descriptions
    2831!
     
    343346          IF ( bc_lr /= 'dirichlet/radiation' )  THEN
    344347             message_string = 'Using synthetic turbulence generator ' //       &
    345                               'requires &bc_lr = "dirichlet/radiation"'
     348                              'requires bc_lr = "dirichlet/radiation"'
    346349             CALL message( 'stg_check_parameters', 'PA0035', 1, 2, 0, 6, 0 )
    347350          ENDIF
    348351          IF ( bc_ns /= 'cyclic' )  THEN
    349352             message_string = 'Using synthetic turbulence generator ' //       &
    350                               'requires &bc_ns = "cyclic"'
     353                              'requires bc_ns = "cyclic"'
    351354             CALL message( 'stg_check_parameters', 'PA0037', 1, 2, 0, 6, 0 )
    352355          ENDIF
  • palm/trunk/SOURCE/temperton_fft_mod.f90

    r2300 r3045  
    99! -----------------
    1010! $Id$
     11! Error message revised
     12!
     13! 2300 2017-06-29 13:31:14Z raasch
    1114! NEC related CPP directives removed
    1215!
     
    21612164
    21622165!    WRITE (nout,'(A,I4,A)') ' n =',n,' - Contains illegal factors'
    2163     message_string = 'number of gridpoints along x or/and y ' // &
    2164                      'contain illegal  factors' //               &
    2165                      '&only factors 8,6,5,4,3,2 are allowed'
     2166    message_string = 'number of gridpoints along x or/and y ' //               &
     2167                     'contain illegal  factors' //                             &
     2168                     ' only factors 2, 3, 5 are allowed'
    21662169    CALL message( 'temperton_fft', 'PA0311', 1, 2, 0, 6, 0 )
    21672170
  • palm/trunk/SOURCE/timestep.f90

    r2718 r3045  
    2525! -----------------
    2626! $Id$
     27! Error message revised
     28!
     29! 2718 2018-01-02 08:49:38Z maronga
    2730! Corrected "Former revisions" section
    2831!
     
    330333
    331334          WRITE( message_string, * ) 'Time step has reached minimum limit.',   &
    332                '&dt              = ', dt_3d, ' s  Simulation is terminated.',  &
    333                '&old_dt          = ', old_dt, ' s',                            &
    334                '&dt_u            = ', dt_u, ' s',                              &
    335                '&dt_v            = ', dt_v, ' s',                              &
    336                '&dt_w            = ', dt_w, ' s',                              &
    337                '&dt_diff         = ', dt_diff, ' s',                           &
    338                '&u_max           = ', u_max, ' m/s   k=', u_max_ijk(1),        &
     335               ' dt              = ', dt_3d, ' s  Simulation is terminated.',  &
     336               ' old_dt          = ', old_dt, ' s',                            &
     337               ' dt_u            = ', dt_u, ' s',                              &
     338               ' dt_v            = ', dt_v, ' s',                              &
     339               ' dt_w            = ', dt_w, ' s',                              &
     340               ' dt_diff         = ', dt_diff, ' s',                           &
     341               ' u_max           = ', u_max, ' m/s   k=', u_max_ijk(1),        &
    339342               '  j=', u_max_ijk(2), '  i=', u_max_ijk(3),                     &
    340                '&v_max           = ', v_max, ' m/s   k=', v_max_ijk(1),        &
     343               ' v_max           = ', v_max, ' m/s   k=', v_max_ijk(1),        &
    341344               '  j=', v_max_ijk(2), '  i=', v_max_ijk(3),                     &
    342                '&w_max           = ', w_max, ' m/s   k=', w_max_ijk(1),        &
     345               ' w_max           = ', w_max, ' m/s   k=', w_max_ijk(1),        &
    343346               '  j=', w_max_ijk(2), '  i=', w_max_ijk(3)
    344347          CALL message( 'timestep', 'PA0312', 0, 1, 0, 6, 0 )
  • palm/trunk/SOURCE/turbulence_closure_mod.f90

    r3014 r3045  
    2525! -----------------
    2626! $Id$
     27! Error message revised
     28!
     29! 3014 2018-05-09 08:42:38Z maronga
    2730! Bugfix: nzb_do and nzt_do were not used for 3d data output
    2831!
     
    12121215
    12131216       DO  k = 1, nzt
    1214           IF ( l_grid(k) > 1.5_wp * dx * wall_adjustment_factor .OR.  &
     1217          IF ( l_grid(k) > 1.5_wp * dx * wall_adjustment_factor .OR.           &
    12151218               l_grid(k) > 1.5_wp * dy * wall_adjustment_factor )  THEN
    1216              WRITE( message_string, * ) 'grid anisotropy exceeds ', &
    1217                                         'threshold given by only local', &
    1218                                         ' &horizontal reduction of near_wall ', &
    1219                                         'mixing length l_wall', &
    1220                                         ' &starting from height level k = ', k, '.'
     1219             WRITE( message_string, * ) 'grid anisotropy exceeds ',            &
     1220                                        'threshold given by only local',       &
     1221                                        ' horizontal reduction of near_wall ', &
     1222                                        'mixing length l_wall',                &
     1223                                        ' starting from height level k = ', k, &
     1224                                        '.'
    12211225             CALL message( 'init_grid', 'PA0202', 0, 1, 0, 6, 0 )
    12221226             EXIT
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r3029 r3045  
    2828! -----------------
    2929! $Id$
     30! Error message added
     31!
     32! 3029 2018-05-23 12:19:17Z raasch
    3033! bugfix: close unit 151 instead of 90
    3134!
     
    320323    USE radiation_model_mod,                                                   &
    321324        ONLY:  albedo_type, radiation_interaction, calc_zenith, zenith,        &
    322                rad_sw_in, rad_lw_in, rad_sw_out, rad_lw_out,                   &
     325               radiation, rad_sw_in, rad_lw_in, rad_sw_out, rad_lw_out,        &
    323326               sigma_sb, solar_constant, sun_direction, sun_dir_lat,           &
    324327               sun_dir_lon,                                                    &
     
    23032306                           'bc_pt_b = "dirichlet" and '//                      &
    23042307                           'bc_q_b  = "dirichlet"'
    2305           CALL message( 'check_parameters', 'PA0590', 1, 2, 0, 6, 0 )
     2308          CALL message( 'usm_check_parameters', 'PA0590', 1, 2, 0, 6, 0 )
    23062309       ENDIF
    23072310
     
    23092312          message_string = 'urban surface model requires '//                   &
    23102313                           'constant_flux_layer = .T.'
    2311           CALL message( 'check_parameters', 'PA0591', 1, 2, 0, 6, 0 )
     2314          CALL message( 'usm_check_parameters', 'PA0084', 1, 2, 0, 6, 0 )
     2315       ENDIF
     2316
     2317       IF (  .NOT.  radiation )  THEN
     2318          message_string = 'urban surface model requires '//                   &
     2319                           'the radiation model to be switched on'
     2320          CALL message( 'usm_check_parameters', 'PA0084', 1, 2, 0, 6, 0 )
    23122321       ENDIF
    23132322!       
  • palm/trunk/SOURCE/vertical_nesting_mod.f90

    r2718 r3045  
    2626! -----------------
    2727! $Id$
     28! Error message revised
     29!
     30! 2718 2018-01-02 08:49:38Z maronga
    2831! Corrected "Former revisions" section
    2932!
     
    267270   
    268271        IF ( myid == 0) THEN
    269             CALL MPI_SENDRECV( terminate_coupled,        1, MPI_INTEGER, &
    270                 target_id, 0,                             &
    271                 terminate_coupled_remote, 1, MPI_INTEGER, &
    272                 target_id, 0,                             &
     272            CALL MPI_SENDRECV( terminate_coupled,        1, MPI_INTEGER,       &
     273                target_id, 0,                                                  &
     274                terminate_coupled_remote, 1, MPI_INTEGER,                      &
     275                target_id, 0,                                                  &
    273276                comm_inter, status, ierr )
    274277        ENDIF
     
    277280   
    278281        IF ( terminate_coupled_remote > 0 )  THEN
    279             WRITE( message_string, * ) 'remote model "',                         &
    280                 TRIM( coupling_mode_remote ),             &
    281                 '" terminated',                           &
    282                 '&with terminate_coupled_remote = ',      &
    283                 terminate_coupled_remote,                 &
    284                 '&local model  "', TRIM( coupling_mode ), &
    285                 '" has',                                  &
    286                 '&terminate_coupled = ',                  &
     282            WRITE( message_string, * ) 'remote model "',                       &
     283                TRIM( coupling_mode_remote ),                                  &
     284                '" terminated',                                                &
     285                ' with terminate_coupled_remote = ',                           &
     286                terminate_coupled_remote,                                      &
     287                ' local model  "', TRIM( coupling_mode ),                      &
     288                '" has',                                                       &
     289                ' terminate_coupled = ',                                       &
    287290                terminate_coupled
    288291            CALL message( 'vnest_init_fine', 'PA0310', 1, 2, 0, 6, 0 )
     
    27772780   
    27782781        IF ( terminate_coupled_remote > 0 )  THEN
    2779             WRITE( message_string, * ) 'remote model "',                         &
    2780                 TRIM( coupling_mode_remote ),             &
    2781                 '" terminated',                           &
    2782                 '&with terminate_coupled_remote = ',      &
    2783                 terminate_coupled_remote,                 &
    2784                 '&local model  "', TRIM( coupling_mode ), &
    2785                 '" has',                                  &
    2786                 '&terminate_coupled = ',                  &
     2782            WRITE( message_string, * ) 'remote model "',                       &
     2783                TRIM( coupling_mode_remote ),                                  &
     2784                '" terminated',                                                &
     2785                ' with terminate_coupled_remote = ',                           &
     2786                terminate_coupled_remote,                                      &
     2787                ' local model  "', TRIM( coupling_mode ),                      &
     2788                '" has',                                                       &
     2789                ' terminate_coupled = ',                                       &
    27872790                terminate_coupled
    27882791            CALL message( 'vnest_anterpolate', 'PA0310', 1, 2, 0, 6, 0 )
     
    33153318   
    33163319        IF ( myid == 0) THEN
    3317             CALL MPI_SENDRECV( terminate_coupled,        1, MPI_INTEGER, &
    3318                 target_id, 0,                             &
    3319                 terminate_coupled_remote, 1, MPI_INTEGER, &
    3320                 target_id, 0,                             &
     3320            CALL MPI_SENDRECV( terminate_coupled,        1, MPI_INTEGER,       &
     3321                target_id, 0,                                                  &
     3322                terminate_coupled_remote, 1, MPI_INTEGER,                      &
     3323                target_id, 0,                                                  &
    33213324                comm_inter, status, ierr )
    33223325        ENDIF
    3323         CALL MPI_BCAST( terminate_coupled_remote, 1, MPI_INTEGER, 0, comm2d, &
     3326        CALL MPI_BCAST( terminate_coupled_remote, 1, MPI_INTEGER, 0, comm2d,   &
    33243327            ierr )
    33253328   
    33263329        IF ( terminate_coupled_remote > 0 )  THEN
    3327             WRITE( message_string, * ) 'remote model "',                         &
    3328                 TRIM( coupling_mode_remote ),             &
    3329                 '" terminated',                           &
    3330                 '&with terminate_coupled_remote = ',      &
    3331                 terminate_coupled_remote,                 &
    3332                 '&local model  "', TRIM( coupling_mode ), &
    3333                 '" has',                                  &
    3334                 '&terminate_coupled = ',                  &
     3330            WRITE( message_string, * ) 'remote model "',                       &
     3331                TRIM( coupling_mode_remote ),                                  &
     3332                '" terminated',                                                &
     3333                ' with terminate_coupled_remote = ',                           &
     3334                terminate_coupled_remote,                                      &
     3335                ' local model  "', TRIM( coupling_mode ),                      &
     3336                '" has',                                                       &
     3337                ' terminate_coupled = ',                                       &
    33353338                terminate_coupled
    33363339            CALL message( 'vnest_anterpolate_e', 'PA0310', 1, 2, 0, 6, 0 )
Note: See TracChangeset for help on using the changeset viewer.