Changeset 292 for palm/trunk/SOURCE


Ignore:
Timestamp:
Apr 16, 2009 4:10:51 PM (13 years ago)
Author:
letzel
Message:
  • Typographical errors (netcdf)
  • If the inversion height calculated by the prerun is zero, inflow_damping_height must be explicitly specified (init_3d_model)
  • Small bugfix concerning 3d 64bit netcdf output format (header)
  • Adjustments for lcsgib and lcsgih (subjob)
Location:
palm/trunk/SOURCE
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/CURRENT_MODIFICATIONS

    r291 r292  
    9494on. (time_integration)
    9595
    96 advec_particles, check_parameters, cpu_log, init_particles, modules, time_integration
     96Typographical errors (netcdf)
     97
     98If the inversion height calculated by the prerun is zero, inflow_damping_height
     99must be explicitly specified (init_3d_model)
     100
     101Small bugfix concerning 3d 64bit netcdf output format (header)
     102
     103
     104advec_particles, check_parameters, cpu_log, header, init_3d_model, init_particles, modules, netcdf, time_integration
  • palm/trunk/SOURCE/header.f90

    r291 r292  
    1111! dp_smooth, dpdxy, u_bulk, v_bulk
    1212! topography_grid_convention moved from user_header
     13! small bugfix concerning 3d 64bit netcdf output format
    1314!
    1415! Former revisions:
     
    693694          ENDIF
    694695       ENDIF
    695        WRITE ( io, 345 )  output_format
     696       WRITE ( io, 344 )  output_format
    696697
    697698       DO  i = 1, dopr_n
     
    771772             ENDIF
    772773          ENDIF
    773           WRITE ( io, 345 )  output_format
     774          WRITE ( io, 344 )  output_format
    774775
    775776          IF ( do2d_xy /= ''  .AND.  section(1,1) /= -9999 )  THEN
     
    915916          output_format = ''
    916917          IF ( netcdf_output )  THEN
    917              IF ( netcdf_64bit .AND. netcdf_64bit_3d )  THEN
     918             IF ( netcdf_64bit_3d )  THEN
    918919                output_format = 'netcdf (64 bit offset)'
    919920             ELSE
     
    928929             ENDIF
    929930          ENDIF
    930           WRITE ( io, 345 )  output_format
     931          WRITE ( io, 344 )  output_format
    931932
    932933          IF ( do3d_at_begin )  THEN
     
    10051006          ENDIF
    10061007       ENDIF
    1007        WRITE ( io, 345 )  output_format
     1008       WRITE ( io, 344 )  output_format
    10081009       WRITE ( io, 341 )  dt_dots
    10091010    ENDIF
     
    10611062          ENDIF
    10621063       ENDIF
    1063        WRITE ( io, 345 )  output_format
     1064       WRITE ( io, 344 )  output_format
    10641065       WRITE ( io, 371 )  dt_dosp
    10651066       IF ( skip_time_dosp /= 0.0 )  WRITE ( io, 339 )  skip_time_dosp
     
    13201321             output_format = 'binary'
    13211322          ENDIF
    1322           WRITE ( io, 345 )  output_format
     1323          WRITE ( io, 344 )  output_format
    13231324       ENDIF
    13241325       IF ( dt_dopts /= 9999999.9 )  WRITE ( io, 494 )  dt_dopts
     
    15931594            '       Averaging input every    ',F8.2,' s'/ &
    15941595            '       Upper output limit at    ',F8.2,' m  (GP ',I4,')'/)
    1595 345 FORMAT ('       Output format: ',A/)
     1596344 FORMAT ('       Output format: ',A/)
    15961597#if defined( __dvrp_graphics )
    15971598360 FORMAT ('    Plot-Sequence with dvrp-software:'/ &
  • palm/trunk/SOURCE/init_3d_model.f90

    r274 r292  
    1212! +conserve_volume_flow_mode: 'default', 'initial_profiles', 'inflow_profile'
    1313! and 'bulk_velocity'
     14! If the inversion height calculated by the prerun is zero,
     15! inflow_damping_height must be explicitly specified.
    1416!
    1517! Former revisions:
     
    538540
    539541!
    540 !--       Set initial horizontal velocities at the lowest computational grid levels
    541 !--       to zero in order to avoid too small time steps caused by the diffusion
    542 !--       limit in the initial phase of a run (at k=1, dz/2 occurs in the
    543 !--       limiting formula!). The original values are stored to be later used for
    544 !--       volume flow control.
     542!--       Set initial horizontal velocities at the lowest computational grid
     543!--       levels to zero in order to avoid too small time steps caused by the
     544!--       diffusion limit in the initial phase of a run (at k=1, dz/2 occurs
     545!--       in the limiting formula!). The original values are stored to be later
     546!--       used for volume flow control.
    545547          DO  i = nxl-1, nxr+1
    546548             DO  j = nys-1, nyn+1
     
    968970          IF ( inflow_damping_height == 9999999.9 )  THEN
    969971!
    970 !--          Default: use the inversion height calculated by the prerun
    971              inflow_damping_height = hom_sum(nzb+6,pr_palm,0)
     972!--          Default: use the inversion height calculated by the prerun; if
     973!--          this is zero, inflow_damping_height must be explicitly specified.
     974             IF ( hom_sum(nzb+6,pr_palm,0) /= 0.0 )  THEN
     975                inflow_damping_height = hom_sum(nzb+6,pr_palm,0)
     976             ELSE
     977                WRITE( message_string, * ) 'inflow_damping_height must be ', &
     978                     'explicitly specified because&the inversion height ', &
     979                     'calculated by the prerun is zero.'
     980                CALL message( 'init_3d_model', 'PA0318', 1, 2, 0, 6, 0 )
     981             ENDIF
    972982
    973983          ENDIF
  • palm/trunk/SOURCE/netcdf.f90

    r288 r292  
    1111! Output of NetCDF messages with aid of message handling routine.
    1212! Output of messages replaced by message handling routine.
     13! Typographical errors fixed.
    1314!
    1415!
     
    181182
    182183!
    183 !-- Select the mode to be processed. Possibilities are xy, xz, yz, pr and ts.
     184!-- Select the mode to be processed. Possibilities are 3d, xy, xz, yz,
     185!-- pr and ts.
    184186    SELECT CASE ( mode )
    185187
     
    563565          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    564566             message_string = 'NetCDF file for volume data ' //             & 
    565                               TRIM( var ) // ' from previuos run found,' // &
     567                              TRIM( var ) // ' from previous run found,' // &
    566568                              '&but this file cannot be extended due to' // &
    567569                              ' variable mismatch.' //                      &
     
    588590          IF ( nz_do3d-nzb+1 /= nz_old )  THEN
    589591              message_string = 'NetCDF file for volume data ' //             &
    590                                TRIM( var ) // ' from previuos run found,' // &
     592                               TRIM( var ) // ' from previous run found,' // &
    591593                               '&but this file cannot be extended due to' // &
    592594                               ' mismatch in number of' //                   &
     
    624626          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    625627             message_string = 'NetCDF file for volume data ' // &
    626                               TRIM( var ) // ' from previuos run found,' // &
     628                              TRIM( var ) // ' from previous run found,' // &
    627629                              '&but this file cannot be extended becaus' // &
    628630                              'e the current output time' //                &
     
    648650
    649651!
    650 !--       Change the titel attribute on file
     652!--       Change the title attribute on file
    651653          nc_stat = NF90_PUT_ATT( id_set_3d(av), NF90_GLOBAL, 'title', &
    652654                                  TRIM( run_description_header ) )
     
    11391141          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    11401142             message_string = 'NetCDF file for cross-sections ' //           &
    1141                               TRIM( var ) // ' from previuos run found,' //  &
     1143                              TRIM( var ) // ' from previous run found,' //  &
    11421144                              '& but this file cannot be extended due to' // &
    11431145                              ' variable mismatch.' //                       &
     
    11721174          IF ( ns /= ns_old )  THEN
    11731175             message_string = 'NetCDF file for cross-sections ' //          &
    1174                               TRIM( var ) // ' from previuos run found,' // &
     1176                              TRIM( var ) // ' from previous run found,' // &
    11751177                              '&but this file cannot be extended due to' // &
    11761178                              ' mismatch in number of' //                   &
     
    11931195                IF ( zu(section(i,1)) /= netcdf_data(i) )  THEN
    11941196                   message_string = 'NetCDF file for cross-sections ' //     &
    1195                                TRIM( var ) // ' from previuos run found,' // &
     1197                               TRIM( var ) // ' from previous run found,' // &
    11961198                               '&but this file cannot be extended' //        &
    11971199                               ' due to mismatch in cross' //                &
     
    12061208                IF ( -1.0 /= netcdf_data(i) )  THEN
    12071209                   message_string = 'NetCDF file for cross-sections ' //     &
    1208                                TRIM( var ) // ' from previuos run found,' // &
     1210                               TRIM( var ) // ' from previous run found,' // &
    12091211                               '&but this file cannot be extended' //        &
    12101212                               ' due to mismatch in cross' //                &
     
    12461248          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    12471249             message_string = 'NetCDF file for cross sections ' //          &
    1248                               TRIM( var ) // ' from previuos run found,' // &
     1250                              TRIM( var ) // ' from previous run found,' // &
    12491251                              '&but this file cannot be extended becaus' // &
    12501252                              'e the current output time' //                &
     
    12741276
    12751277!
    1276 !--       Change the titel attribute on file
     1278!--       Change the title attribute on file
    12771279          nc_stat = NF90_PUT_ATT( id_set_xy(av), NF90_GLOBAL, 'title', &
    12781280                                  TRIM( run_description_header ) )
     
    16611663          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    16621664             message_string = 'NetCDF file for cross-sections ' //           &
    1663                               TRIM( var ) // ' from previuos run found,' //  &
     1665                              TRIM( var ) // ' from previous run found,' //  &
    16641666                              '& but this file cannot be extended due to' // &
    16651667                              ' variable mismatch.' //                       &
     
    16941696          IF ( ns /= ns_old )  THEN
    16951697             message_string = 'NetCDF file for cross-sections ' //          &
    1696                               TRIM( var ) // ' from previuos run found,' // &
     1698                              TRIM( var ) // ' from previous run found,' // &
    16971699                              '&but this file cannot be extended due to' // &
    16981700                              ' mismatch in number of' //                   &
     
    17151717                IF ( ( section(i,2) * dy ) /= netcdf_data(i) )  THEN
    17161718                   message_string = 'NetCDF file for cross-sections ' //     &
    1717                                TRIM( var ) // ' from previuos run found,' // &
     1719                               TRIM( var ) // ' from previous run found,' // &
    17181720                               '&but this file cannot be extended' //        &
    17191721                               ' due to mismatch in cross' //                &
     
    17281730                IF ( -1.0 /= netcdf_data(i) )  THEN
    17291731                   message_string = 'NetCDF file for cross-sections ' //     &
    1730                                TRIM( var ) // ' from previuos run found,' // &
     1732                               TRIM( var ) // ' from previous run found,' // &
    17311733                               '&but this file cannot be extended' //        &
    17321734                               ' due to mismatch in cross' //                &
     
    17681770          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    17691771             message_string = 'NetCDF file for cross sections ' //          &
    1770                               TRIM( var ) // ' from previuos run found,' // &
     1772                              TRIM( var ) // ' from previous run found,' // &
    17711773                              '&but this file cannot be extended becaus' // &
    17721774                              'e the current output time' //                &
     
    17961798
    17971799!
    1798 !--       Change the titel attribute on file
     1800!--       Change the title attribute on file
    17991801          nc_stat = NF90_PUT_ATT( id_set_xz(av), NF90_GLOBAL, 'title', &
    18001802                                  TRIM( run_description_header ) )
     
    21812183          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    21822184             message_string = 'NetCDF file for cross-sections ' //           &
    2183                               TRIM( var ) // ' from previuos run found,' //  &
     2185                              TRIM( var ) // ' from previous run found,' //  &
    21842186                              '& but this file cannot be extended due to' // &
    21852187                              ' variable mismatch.' //                       &
     
    22142216          IF ( ns /= ns_old )  THEN
    22152217             message_string = 'NetCDF file for cross-sections ' //          &
    2216                               TRIM( var ) // ' from previuos run found,' // &
     2218                              TRIM( var ) // ' from previous run found,' // &
    22172219                              '&but this file cannot be extended due to' // &
    22182220                              ' mismatch in number of' //                   &
     
    22352237                IF ( ( section(i,3) * dx ) /= netcdf_data(i) )  THEN
    22362238                   message_string = 'NetCDF file for cross-sections ' //    &
    2237                               TRIM( var ) // ' from previuos run found,' // &
     2239                              TRIM( var ) // ' from previous run found,' // &
    22382240                              '&but this file cannot be extended' //        &
    22392241                              ' due to mismatch in cross' //                &
     
    22482250                IF ( -1.0 /= netcdf_data(i) )  THEN
    22492251                   message_string = 'NetCDF file for cross-sections ' //    &
    2250                               TRIM( var ) // ' from previuos run found,' // &
     2252                              TRIM( var ) // ' from previous run found,' // &
    22512253                              '&but this file cannot be extended' //        &
    22522254                              ' due to mismatch in cross' //                &
     
    22882290          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    22892291             message_string = 'NetCDF file for cross sections ' //          &
    2290                               TRIM( var ) // ' from previuos run found,' // &
     2292                              TRIM( var ) // ' from previous run found,' // &
    22912293                              '&but this file cannot be extended becaus' // &
    22922294                              'e the current output time' //                &
     
    23162318
    23172319!
    2318 !--       Change the titel attribute on file
     2320!--       Change the title attribute on file
    23192321          nc_stat = NF90_PUT_ATT( id_set_yz(av), NF90_GLOBAL, 'title', &
    23202322                                  TRIM( run_description_header ) )
     
    25282530          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    25292531             message_string = 'NetCDF file for vertical profiles ' //        &
    2530                               'from previuos run found,' //                  &
     2532                              'from previous run found,' //                  &
    25312533                              '& but this file cannot be extended due to' // &
    25322534                              ' variable mismatch.' //                       &
     
    25622564          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    25632565             message_string = 'NetCDF file for vertical profiles ' //       &
    2564                               'from previuos run found,' //                 &
     2566                              'from previous run found,' //                 &
    25652567                              '&but this file cannot be extended becaus' // &
    25662568                              'e the current output time' //                &
     
    27432745          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    27442746             message_string = 'NetCDF file for time series ' //              &
    2745                               'from previuos run found,' //                  &
     2747                              'from previous run found,' //                  &
    27462748                              '& but this file cannot be extended due to' // &
    27472749                              ' variable mismatch.' //                       &
     
    27772779          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    27782780             message_string = 'NetCDF file for time series ' //             &
    2779                               'from previuos run found,' //                 &
     2781                              'from previous run found,' //                 &
    27802782                              '&but this file cannot be extended becaus' // &
    27812783                              'e the current output time' //                &
     
    30993101          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    31003102             message_string = 'NetCDF file for spectra  ' //                 &
    3101                               'from previuos run found,' //                  &
     3103                              'from previous run found,' //                  &
    31023104                              '& but this file cannot be extended due to' // &
    31033105                              ' variable mismatch.' //                       &
     
    31333135          IF ( ns /= ns_old )  THEN
    31343136             message_string = 'NetCDF file for spectra ' //                 &
    3135                               ' from previuos run found,' //                &
     3137                              ' from previous run found,' //                &
    31363138                              '&but this file cannot be extended due to' // &
    31373139                              ' mismatch in number of' //                   &
     
    31533155             IF ( zu(comp_spectra_level(i)) /= netcdf_data(i) )  THEN
    31543156                message_string = 'NetCDF file for spectra ' //                 &
    3155                                  ' from previuos run found,' //                &
     3157                                 ' from previous run found,' //                &
    31563158                                 '&but this file cannot be extended due to' // &
    31573159                                 ' mismatch in heights of' //                  &
     
    31913193          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    31923194             message_string = 'NetCDF file for spectra ' //                 &
    3193                               'from previuos run found,' //                 &
     3195                              'from previous run found,' //                 &
    31943196                              '&but this file cannot be extended becaus' // &
    31953197                              'e the current output time' //                &
     
    32283230
    32293231!
    3230 !--       Change the titel attribute on file
     3232!--       Change the title attribute on file
    32313233          IF ( averaging_interval_sp /= 0.0 )  THEN
    32323234             WRITE (time_average_text,'('', '',F7.1,'' s average'')') &
     
    33503352          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    33513353             message_string = 'NetCDF file for particles ' //               &
    3352                               'from previuos run found,' //                 &
     3354                              'from previous run found,' //                 &
    33533355                              '&but this file cannot be extended becaus' // &
    33543356                              'e the current output time' //                &
     
    35083510          IF ( TRIM( var_list ) /= TRIM( var_list_old ) )  THEN
    35093511             message_string = 'NetCDF file for particle time series ' //     &
    3510                               'from previuos run found,' //                  &
     3512                              'from previous run found,' //                  &
    35113513                              '& but this file cannot be extended due to' // &
    35123514                              ' variable mismatch.' //                       &
     
    35423544          IF ( last_time_coordinate(1) >= simulated_time )  THEN
    35433545             message_string = 'NetCDF file for particle time series ' //    &
    3544                               'from previuos run found,' //                 &
     3546                              'from previous run found,' //                 &
    35453547                              '&but this file cannot be extended becaus' // &
    35463548                              'e the current output time' //                &
Note: See TracChangeset for help on using the changeset viewer.