Changeset 3632 for palm/trunk/SOURCE/data_output_mask.f90

Timestamp:

Dec 17, 2018 2:52:50 PM (5 years ago)

Author:

kanani

Message:

Fix for masked output, change default chemistry mechanism

File:

• palm/trunk/SOURCE/data_output_mask.f90 (modified) (5 diffs)

palm/trunk/SOURCE/data_output_mask.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Fix output time levels (use time_since_reference_point)
+! 
+! 3589 2018-11-30 15:09:51Z suehring
 ! Move the control parameter "salsa" from salsa_mod.f90 to control_parameters
 ! 
@@ -168 +171 @@
                mask_j, mask_k, mask_size, mask_size_l, mask_start_l,           &
                mask_surface,                                                   &
-               max_masks, message_string, mid, nz_do3d, salsa, simulated_time
+               max_masks, message_string, mid, nz_do3d, salsa,                 &
+               time_since_reference_point
+
     USE cpulog,                                                                &
         ONLY:  cpu_log, log_point
@@ -265 +270 @@
     IF ( myid == 0  .OR.  netcdf_data_format > 4 )  THEN
        nc_stat = NF90_PUT_VAR( id_set_mask(mid,av), id_var_time_mask(mid,av), &
-                               (/ simulated_time /),                          &
+                               (/ time_since_reference_point /),              &
                                start = (/ domask_time_count(mid,av) /),       &
                                count = (/ 1 /) )
@@ -370 +375 @@
           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
@@ -514 +519 @@
           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