Ignore:
Timestamp:
Dec 18, 2019 11:55:56 AM (23 months ago)
Author:
motisi
Message:

Introduction of wall_flags_total_0, which currently sets bits based on static topography information used in wall_flags_static_0

File:
1 edited

Legend:

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

    r4342 r4346  
    2727! -----------------
    2828! $Id$
     29! Introduction of wall_flags_total_0, which currently sets bits based on static
     30! topography information used in wall_flags_static_0
     31!
     32! 4342 2019-12-16 13:49:14Z Giersch
    2933! Use statements moved to module level, ocean dependency removed, redundant
    3034! variables removed
     
    152156    USE indices,                                                               &
    153157        ONLY:  nbgp, nxl, nxlg, nxlu, nxr, nxrg, nyn, nyng, nys, nysg, nysv,   &
    154                nz, nzb, nzt, topo_top_ind, wall_flags_static_0
     158               nz, nzb, nzt, topo_top_ind, wall_flags_total_0
    155159
    156160    USE kinds
     
    10601064!--                   an informative message is given.
    10611065                      IF ( ANY( lad_s(:,j,i) /= 0.0_wp )  .AND.                &
    1062                            ANY( BTEST( wall_flags_static_0(:,j,i), 6 ) ) .AND.        &
    1063                            ANY( BTEST( wall_flags_static_0(:,j,i), 4 ) ) )  THEN
     1066                           ANY( BTEST( wall_flags_total_0(:,j,i), 6 ) ) .AND.  &
     1067                           ANY( BTEST( wall_flags_total_0(:,j,i), 4 ) ) )  THEN
    10641068                         lad_s(:,j,i) = 0.0_wp
    10651069                         WRITE( message_string, * )                            &
     
    14931497!--                of the potential rooftop grid point, while building_edge_e is
    14941498!--                is set from the perspective of the non-building grid point.
    1495                    building_edge_w = ANY( BTEST( wall_flags_static_0(:,j,i),   6 ) )&
    1496                         .AND.  .NOT. ANY( BTEST( wall_flags_static_0(:,j,i-1), 6 ) )
    1497                    building_edge_e = ANY( BTEST( wall_flags_static_0(:,j,i-1), 6 ) )&
    1498                         .AND.  .NOT. ANY( BTEST( wall_flags_static_0(:,j,i),   6 ) )
     1499                   building_edge_w = ANY( BTEST( wall_flags_total_0(:,j,i),   6 ) )&
     1500                        .AND.  .NOT. ANY( BTEST( wall_flags_total_0(:,j,i-1), 6 ) )
     1501                   building_edge_e = ANY( BTEST( wall_flags_total_0(:,j,i-1), 6 ) )&
     1502                        .AND.  .NOT. ANY( BTEST( wall_flags_total_0(:,j,i),   6 ) )
    14991503!
    15001504!--                Determine topography-top index on u-grid
     
    15731577!--                of the potential rooftop grid point, while building_edge_n is
    15741578!--                is set from the perspective of the non-building grid point.
    1575                    building_edge_s = ANY( BTEST( wall_flags_static_0(:,j,i),   6 ) )&
    1576                         .AND.  .NOT. ANY( BTEST( wall_flags_static_0(:,j-1,i), 6 ) )
    1577                    building_edge_n = ANY( BTEST( wall_flags_static_0(:,j-1,i), 6 ) )&
    1578                         .AND.  .NOT. ANY( BTEST( wall_flags_static_0(:,j,i),   6 ) )
     1579                   building_edge_s = ANY( BTEST( wall_flags_total_0(:,j,i),   6 ) )&
     1580                        .AND.  .NOT. ANY( BTEST( wall_flags_total_0(:,j-1,i), 6 ) )
     1581                   building_edge_n = ANY( BTEST( wall_flags_total_0(:,j-1,i), 6 ) )&
     1582                        .AND.  .NOT. ANY( BTEST( wall_flags_total_0(:,j,i),   6 ) )
    15791583!
    15801584!--                Determine topography-top index on v-grid
     
    18811885!--          of the potential rooftop grid point, while building_edge_e is
    18821886!--          is set from the perspective of the non-building grid point.
    1883              building_edge_w = ANY( BTEST( wall_flags_static_0(:,j,i),   6 ) )  .AND. &
    1884                          .NOT. ANY( BTEST( wall_flags_static_0(:,j,i-1), 6 ) )
    1885              building_edge_e = ANY( BTEST( wall_flags_static_0(:,j,i-1), 6 ) )  .AND. &
    1886                          .NOT. ANY( BTEST( wall_flags_static_0(:,j,i),   6 ) )
     1887             building_edge_w = ANY( BTEST( wall_flags_total_0(:,j,i),   6 ) )  .AND. &
     1888                         .NOT. ANY( BTEST( wall_flags_total_0(:,j,i-1), 6 ) )
     1889             building_edge_e = ANY( BTEST( wall_flags_total_0(:,j,i-1), 6 ) )  .AND. &
     1890                         .NOT. ANY( BTEST( wall_flags_total_0(:,j,i),   6 ) )
    18871891!
    18881892!--          Determine topography-top index on u-grid
     
    19601964!--          of the potential rooftop grid point, while building_edge_n is
    19611965!--          is set from the perspective of the non-building grid point.
    1962              building_edge_s = ANY( BTEST( wall_flags_static_0(:,j,i),   6 ) )  .AND. &
    1963                          .NOT. ANY( BTEST( wall_flags_static_0(:,j-1,i), 6 ) )
    1964              building_edge_n = ANY( BTEST( wall_flags_static_0(:,j-1,i), 6 ) )  .AND. &
    1965                          .NOT. ANY( BTEST( wall_flags_static_0(:,j,i),   6 ) )
     1966             building_edge_s = ANY( BTEST( wall_flags_total_0(:,j,i),   6 ) )  .AND. &
     1967                         .NOT. ANY( BTEST( wall_flags_total_0(:,j-1,i), 6 ) )
     1968             building_edge_n = ANY( BTEST( wall_flags_total_0(:,j-1,i), 6 ) )  .AND. &
     1969                         .NOT. ANY( BTEST( wall_flags_total_0(:,j,i),   6 ) )
    19661970!
    19671971!--          Determine topography-top index on v-grid
Note: See TracChangeset for help on using the changeset viewer.