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

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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
Note: See TracChangeset for help on using the changeset viewer.