Ignore:
Timestamp:
Jul 30, 2019 2:47:10 PM (2 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:
1 edited

Legend:

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

    r4048 r4127  
    2525! -----------------
    2626! $Id$
     27! Adjustment for top boundary index for plant-canopy model outputs
     28! (merge from branch resler)
     29!
     30! 4048 2019-06-21 21:00:21Z knoop
    2731! Moved tcm_data_output_3d to module_interface
    2832!
     
    276280        ONLY:  debug_output_timestep,                                          &
    277281               do3d, do3d_no, do3d_time_count, io_blocks, io_group,            &
    278                land_surface, message_string, ntdim_3d, nz_do3d, psolver,       &
    279                time_since_reference_point, urban_surface, varnamelength
     282               land_surface, message_string, ntdim_3d, nz_do3d, plant_canopy,  &
     283               psolver, time_since_reference_point, urban_surface,             &
     284               varnamelength
    280285
    281286    USE cpulog,                                                                &
     
    313318
    314319    USE pegrid
     320
     321    USE plant_canopy_model_mod,                                                &
     322        ONLY:  pch_index
    315323
    316324    USE radiation_model_mod,                                                   &
     
    410418       found = .FALSE.
    411419       resorted = .FALSE.
     420       trimvar = TRIM( do3d(av,ivar) )
     421
    412422!
    413423!--    Temporary solution to account for data output within the new urban
     
    415425!--    Store the array chosen on the temporary array.
    416426       nzb_do   = nzb
    417        nzt_do   = nz_do3d
    418 
    419        trimvar = TRIM( do3d(av,ivar) )
     427!
     428!--    Set top index for 3D output. Note in case of plant-canopy model
     429!--    these index is determined by pch_index.
     430       IF ( plant_canopy  .AND.  trimvar(1:4) == 'pcm_' )  THEN
     431          nzt_do   = pch_index
     432       ELSE
     433          nzt_do   = nz_do3d
     434       ENDIF
     435
    420436!
    421437!--    Allocate a temporary array with the desired output dimensions.
Note: See TracChangeset for help on using the changeset viewer.