Changeset 4426

Timestamp:

Feb 27, 2020 10:02:19 AM (5 years 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:

• palm/trunk/SOURCE/wind_turbine_model_mod.f90 (modified) (6 diffs)

palm/trunk/SOURCE/wind_turbine_model_mod.f90

@@ -26 +26 @@
 ! -----------------
 ! $Id$
+! define time as unlimited dimension so that no maximum number
+! of time steps has to be given for wtm_data_output
+! 
+! 4423 2020-02-25 07:17:11Z maronga
 ! Switched to serial output as data is aggerated before anyway.
 ! 
@@ -173 +177 @@
     LOGICAL ::  initial_write_coordinates = .FALSE. 
     
-    REAL(wp) ::  dt_data_output_wtm = 1.0_wp       !< data output interval
+    REAL(wp) ::  dt_data_output_wtm = 0.0_wp       !< data output interval
     REAL(wp) ::  time_wtm           = 0.0_wp       !< time since last data output
     
@@ -654 +658 @@
        ENDIF
 
-       IF ( dt_data_output_wtm <= 0.0 )  THEN
-          message_string = 'dt_data_output_wtm must be > 0.0'
-       CALL message( 'check_parameters', 'PA0709', 1, 2, 0, 6, 0 )
-    ENDIF
-       
- 
     END SUBROUTINE wtm_check_parameters 
 !     
@@ -1288 +1286 @@
     
     
-    INTEGER(iwp) ::  ntimesteps               !< number of timesteps defined in NetCDF output file
-    INTEGER(iwp) ::  ntimesteps_max = 80000   !< number of maximum timesteps defined in NetCDF output file
+!    INTEGER(iwp) ::  ntimesteps               !< number of timesteps defined in NetCDF output file
+!    INTEGER(iwp) ::  ntimesteps_max = 80000   !< number of maximum timesteps defined in NetCDF output file
     INTEGER(iwp) ::  return_value             !< returned status value of called function
     
@@ -1301 +1299 @@
     nc_filename = 'DATA_1D_TS_WTM_NETCDF' // TRIM( coupling_char )      
     return_value = dom_def_file( nc_filename, 'netcdf4-serial' )
-                                                                          
-    ntimesteps = MIN( CEILING(                                                      &
-                    ( end_time - MAX( time_turbine_on, time_since_reference_point ) &
-                    ) / MAX( 0.1_wp, dt_data_output_wtm ) ), ntimesteps_max )
-
-      
+
     IF ( myid == 0 )  THEN
 !
@@ -1321 +1314 @@
        DEALLOCATE( ndim )
 
-       
 !
 !--    time
        return_value = dom_def_dim( nc_filename,                                &
-                                   dimension_name = 'time',                   &
-                                   output_type = 'real32',                      &
-                                   bounds = (/1_iwp, ntimesteps/),             &
-                                   values_realwp = (/0.0_wp/) )   
-       
+                                   dimension_name = 'time',                    &
+                                   output_type = 'real32',                     &
+                                   bounds = (/1_iwp/),                         &
+                                   values_realwp = (/0.0_wp/) )
  
 !