Changeset 3646 for palm/trunk/SOURCE/data_output_3d.f90

Timestamp:

Dec 28, 2018 5:58:49 PM (5 years ago)

Author:

kanani

Message:

Bugfix: replace simulated_time by time_since_reference_point where required

File:

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

palm/trunk/SOURCE/data_output_3d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix: use time_since_reference_point instead of simulated_time (relevant
+! when using wall/soil spinup)
+! 
+! 3637 2018-12-20 01:51:36Z knoop
 ! Implementation of the PALM module interface
 ! 
@@ -255 +259 @@
                io_blocks, io_group, land_surface, message_string,              &
                ntdim_3d, nz_do3d,  ocean_mode, plant_canopy,                   &
-               psolver, salsa, simulated_time, time_since_reference_point,     &
+               psolver, salsa, time_since_reference_point,                     &
                urban_surface, varnamelength
 
@@ -355 +359 @@
     IF ( netcdf_data_format > 4 )  THEN
        IF ( do3d_time_count(av) + 1 > ntdim_3d(av) )  THEN
-          WRITE ( message_string, * ) 'Output of 3d data is not given at t=',  &
-                                      simulated_time, '&because the maximum ', & 
-                                      'number of output time levels is ',      &
+          WRITE ( message_string, * ) 'Output of 3d data is not given at t=',               &
+                                      time_since_reference_point, 's because the maximum ', & 
+                                      'number of output time levels is ',                   &
                                       'exceeded.'
           CALL message( 'data_output_3d', 'PA0387', 0, 1, 0, 6, 0 )
@@ -447 +451 @@
           CASE ( 'pc' )  ! particle concentration (requires ghostpoint exchange)
              IF ( av == 0 )  THEN
-                IF ( simulated_time >= particle_advection_start )  THEN
+                IF ( time_since_reference_point >= particle_advection_start )  THEN
                    tend = prt_count
                 ELSE
@@ -470 +474 @@
           CASE ( 'pr' )  ! mean particle radius (effective radius)
              IF ( av == 0 )  THEN
-                IF ( simulated_time >= particle_advection_start )  THEN
+                IF ( time_since_reference_point >= particle_advection_start )  THEN
                    DO  i = nxl, nxr
                       DO  j = nys, nyn
@@ -585 +589 @@
           CASE ( 'ql_vp' )
              IF ( av == 0 )  THEN
-                IF ( simulated_time >= particle_advection_start )  THEN
+                IF ( time_since_reference_point >= particle_advection_start )  THEN
                    DO  i = nxl, nxr
                       DO  j = nys, nyn