Changeset 4426 for palm/trunk


Ignore:
Timestamp:
Feb 27, 2020 10:02:19 AM (20 months ago)
Author:
oliver.maas
Message:

define time as unlimited dimension so that no maximum number of time steps has to be given for wtm_data_output

File:
1 edited

Legend:

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

    r4423 r4426  
    2626! -----------------
    2727! $Id$
     28! define time as unlimited dimension so that no maximum number
     29! of time steps has to be given for wtm_data_output
     30!
     31! 4423 2020-02-25 07:17:11Z maronga
    2832! Switched to serial output as data is aggerated before anyway.
    2933!
     
    173177    LOGICAL ::  initial_write_coordinates = .FALSE.
    174178   
    175     REAL(wp) ::  dt_data_output_wtm = 1.0_wp       !< data output interval
     179    REAL(wp) ::  dt_data_output_wtm = 0.0_wp       !< data output interval
    176180    REAL(wp) ::  time_wtm           = 0.0_wp       !< time since last data output
    177181   
     
    654658       ENDIF
    655659
    656        IF ( dt_data_output_wtm <= 0.0 )  THEN
    657           message_string = 'dt_data_output_wtm must be > 0.0'
    658        CALL message( 'check_parameters', 'PA0709', 1, 2, 0, 6, 0 )
    659     ENDIF
    660        
    661  
    662660    END SUBROUTINE wtm_check_parameters
    663661!     
     
    12881286   
    12891287   
    1290     INTEGER(iwp) ::  ntimesteps               !< number of timesteps defined in NetCDF output file
    1291     INTEGER(iwp) ::  ntimesteps_max = 80000   !< number of maximum timesteps defined in NetCDF output file
     1288!    INTEGER(iwp) ::  ntimesteps               !< number of timesteps defined in NetCDF output file
     1289!    INTEGER(iwp) ::  ntimesteps_max = 80000   !< number of maximum timesteps defined in NetCDF output file
    12921290    INTEGER(iwp) ::  return_value             !< returned status value of called function
    12931291   
     
    13011299    nc_filename = 'DATA_1D_TS_WTM_NETCDF' // TRIM( coupling_char )     
    13021300    return_value = dom_def_file( nc_filename, 'netcdf4-serial' )
    1303                                                                          
    1304     ntimesteps = MIN( CEILING(                                                      &
    1305                     ( end_time - MAX( time_turbine_on, time_since_reference_point ) &
    1306                     ) / MAX( 0.1_wp, dt_data_output_wtm ) ), ntimesteps_max )
    1307 
    1308      
     1301
    13091302    IF ( myid == 0 )  THEN
    13101303!
     
    13211314       DEALLOCATE( ndim )
    13221315
    1323        
    13241316!
    13251317!--    time
    13261318       return_value = dom_def_dim( nc_filename,                                &
    1327                                    dimension_name = 'time',                   &
    1328                                    output_type = 'real32',                      &
    1329                                    bounds = (/1_iwp, ntimesteps/),             &
    1330                                    values_realwp = (/0.0_wp/) )   
    1331        
     1319                                   dimension_name = 'time',                    &
     1320                                   output_type = 'real32',                     &
     1321                                   bounds = (/1_iwp/),                         &
     1322                                   values_realwp = (/0.0_wp/) )
    13321323 
    13331324!
Note: See TracChangeset for help on using the changeset viewer.