Ignore:
Timestamp:
Aug 16, 2019 1:50:17 PM (5 years ago)
Author:
suehring
Message:

Replace get_topography_top_index functions by pre-calculated arrays in order to save computational resources

File:
1 edited

Legend:

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

    r4167 r4168  
    2525! -----------------
    2626! $Id$
     27! Remove variable grid
     28!
     29! 4167 2019-08-16 11:01:48Z suehring
    2730! Changed behaviour of masked output over surface to follow terrain and ignore
    2831! buildings (J.Resler, T.Gronemeier)
     
    225228
    226229    IMPLICIT NONE
    227 
    228     CHARACTER(LEN=5) ::  grid !< flag to distinquish between staggered grids
    229230
    230231    INTEGER(iwp) ::  av                      !< flag for (non-)average output
     
    310311       ENDIF
    311312!
    312 !--    Set default grid for terrain-following output
    313        grid = 's'
    314 !
    315313!--    Store the variable chosen.
    316314       resorted = .FALSE.
     
    373371                         im = mask_i(mid,i)
    374372                         jm = mask_j(mid,j)
    375                          ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_0(:,jm,im), 5 )), DIM = 1 ) - 1
     373                         ktt = MINLOC( MERGE( 1, 0, BTEST( wall_flags_0(:,jm,im), 5 )),&
     374                                       DIM = 1 ) - 1
    376375                         DO  k = 1, mask_size_l(mid,3)
    377376                            kk = MIN( ktt+mask_k(mid,k), nzt+1 )
     
    685684
    686685          CASE ( 'w' )
    687              grid = 'w'
    688686             IF ( av == 0 )  THEN
    689687                to_be_resorted => w
Note: See TracChangeset for help on using the changeset viewer.