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

Timestamp:

Jul 30, 2019 2:47:10 PM (5 years ago)

Author:

suehring

Message:

Merge with branch resler: biomet- output of bio_mrt added; plant_canopy - separate vertical dimension for 3D output (to save disk space); radiation - remove unused plant canopy variables; urban-surface model - do not add anthropogenic heat during wall spin-up

File:

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

palm/trunk/SOURCE/data_output_3d.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Adjustment for top boundary index for plant-canopy model outputs 
+! (merge from branch resler)
+! 
+! 4048 2019-06-21 21:00:21Z knoop
 ! Moved tcm_data_output_3d to module_interface
 ! 
@@ -276 +280 @@
         ONLY:  debug_output_timestep,                                          &
                do3d, do3d_no, do3d_time_count, io_blocks, io_group,            &
-               land_surface, message_string, ntdim_3d, nz_do3d, psolver,       &
-               time_since_reference_point, urban_surface, varnamelength
+               land_surface, message_string, ntdim_3d, nz_do3d, plant_canopy,  &
+               psolver, time_since_reference_point, urban_surface,             &
+               varnamelength
 
     USE cpulog,                                                                &
@@ -313 +318 @@
 
     USE pegrid
+
+    USE plant_canopy_model_mod,                                                &
+        ONLY:  pch_index
 
     USE radiation_model_mod,                                                   &
@@ -410 +418 @@
        found = .FALSE.
        resorted = .FALSE.
+       trimvar = TRIM( do3d(av,ivar) )
+
 !
 !--    Temporary solution to account for data output within the new urban
@@ -415 +425 @@
 !--    Store the array chosen on the temporary array.
        nzb_do   = nzb
-       nzt_do   = nz_do3d
-
-       trimvar = TRIM( do3d(av,ivar) )
+!
+!--    Set top index for 3D output. Note in case of plant-canopy model
+!--    these index is determined by pch_index. 
+       IF ( plant_canopy  .AND.  trimvar(1:4) == 'pcm_' )  THEN
+          nzt_do   = pch_index
+       ELSE
+          nzt_do   = nz_do3d
+       ENDIF
+
 !
 !--    Allocate a temporary array with the desired output dimensions.