Ignore:
Timestamp:
Dec 14, 2017 6:46:24 PM (4 years ago)
Author:
suehring
Message:

Particle reflections at downward-facing walls; revision of particle speed interpolations at walls; bugfixes in get_topography_index and in date constants

File:
1 edited

Legend:

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

    r2696 r2698  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Bugfix in get_topography_top_index
    2323!
    2424! Former revisions:
     
    151151
    152152    USE surface_mod,                                                           &
    153         ONLY:  get_topography_top_index
     153        ONLY:  get_topography_top_index_ji
    154154
    155155
     
    400400               DO  j = nys, nyn
    401401                  IF ( pch_index_ji(j,i) /= 0 )  THEN
    402                      k_topo = get_topography_top_index( j, i, 's' )
     402                     k_topo = get_topography_top_index_ji( j, i, 's' )
    403403                     DO  k = k_topo, k_topo + pch_index_ji(j,i)
    404404                        local_pf(i,j,k) = pc_heating_rate(k-k_topo,j,i)
     
    415415               DO  j = nys, nyn
    416416                  IF ( pch_index_ji(j,i) /= 0 )  THEN
    417                      k_topo = get_topography_top_index( j, i, 's' )
     417                     k_topo = get_topography_top_index_ji( j, i, 's' )
    418418                     DO  k = k_topo, k_topo + pch_index_ji(j,i)
    419419                        local_pf(i,j,k) = lad_s(k-k_topo,j,i)
     
    803803!--          Check whether topography and local vegetation on top exceed
    804804!--          height of the model domain.
    805              k = get_topography_top_index( j, i, 's' )
     805             k = get_topography_top_index_ji( j, i, 's' )
    806806             IF ( k + pch_index_ji(j,i) >= nzt + 1 )  THEN
    807807                message_string =  'Local vegetation height on top of ' //      &
     
    11101110!
    11111111!--                Determine topography-top index on u-grid
    1112                    k_wall = get_topography_top_index( j, i, 'u' )
     1112                   k_wall = get_topography_top_index_ji( j, i, 'u' )
    11131113                   DO  k = k_wall+1, k_wall + pch_index_ji(j,i)
    11141114
     
    11741174!
    11751175!--                Determine topography-top index on v-grid
    1176                    k_wall = get_topography_top_index( j, i, 'v' )
     1176                   k_wall = get_topography_top_index_ji( j, i, 'v' )
    11771177
    11781178                   DO  k = k_wall+1, k_wall + pch_index_ji(j,i)
     
    12391239!
    12401240!--                Determine topography-top index on w-grid
    1241                    k_wall = get_topography_top_index( j, i, 'w' )
     1241                   k_wall = get_topography_top_index_ji( j, i, 'w' )
    12421242
    12431243                   DO  k = k_wall+1, k_wall + pch_index_ji(j,i) - 1
     
    12911291!
    12921292!--                Determine topography-top index on scalar-grid
    1293                    k_wall = get_topography_top_index( j, i, 's' )
     1293                   k_wall = get_topography_top_index_ji( j, i, 's' )
    12941294
    12951295                   DO  k = k_wall+1, k_wall + pch_index_ji(j,i)
     
    13081308!
    13091309!--                Determine topography-top index on scalar-grid
    1310                    k_wall = get_topography_top_index( j, i, 's' )
     1310                   k_wall = get_topography_top_index_ji( j, i, 's' )
    13111311
    13121312                   DO  k = k_wall+1, k_wall + pch_index_ji(j,i)
     
    13381338!
    13391339!--                Determine topography-top index on scalar-grid
    1340                    k_wall = get_topography_top_index( j, i, 's' )
     1340                   k_wall = get_topography_top_index_ji( j, i, 's' )
    13411341
    13421342                   DO  k = k_wall+1, k_wall + pch_index_ji(j,i)
     
    13671367!
    13681368!--                Determine topography-top index on scalar-grid
    1369                    k_wall = get_topography_top_index( j, i, 's' )
     1369                   k_wall = get_topography_top_index_ji( j, i, 's' )
    13701370
    13711371                   DO  k = k_wall+1, k_wall + pch_index_ji(j,i)
     
    14611461!
    14621462!--          Determine topography-top index on u-grid
    1463              k_wall = get_topography_top_index( j, i, 'u' )
     1463             k_wall = get_topography_top_index_ji( j, i, 'u' )
    14641464             DO  k = k_wall + 1, k_wall + pch_index_ji(j,i)
    14651465
     
    15211521!
    15221522!--          Determine topography-top index on v-grid
    1523              k_wall = get_topography_top_index( j, i, 'v' )
     1523             k_wall = get_topography_top_index_ji( j, i, 'v' )
    15241524
    15251525             DO  k = k_wall + 1, k_wall + pch_index_ji(j,i)
     
    15811581!
    15821582!--          Determine topography-top index on w-grid
    1583              k_wall = get_topography_top_index( j, i, 'w' )
     1583             k_wall = get_topography_top_index_ji( j, i, 'w' )
    15841584
    15851585             DO  k = k_wall + 1, k_wall + pch_index_ji(j,i) - 1
     
    16281628!
    16291629!--          Determine topography-top index on scalar grid
    1630              k_wall = get_topography_top_index( j, i, 's' )
     1630             k_wall = get_topography_top_index_ji( j, i, 's' )
    16311631
    16321632             DO  k = k_wall + 1, k_wall + pch_index_ji(j,i)
     
    16411641!
    16421642!--          Determine topography-top index on scalar grid
    1643              k_wall = get_topography_top_index( j, i, 's' )
     1643             k_wall = get_topography_top_index_ji( j, i, 's' )
    16441644
    16451645             DO  k = k_wall + 1, k_wall + pch_index_ji(j,i)
     
    16671667!
    16681668!--          Determine topography-top index on scalar grid
    1669              k_wall = get_topography_top_index( j, i, 's' )
     1669             k_wall = get_topography_top_index_ji( j, i, 's' )
    16701670
    16711671             DO  k = k_wall + 1, k_wall + pch_index_ji(j,i)
     
    16931693!
    16941694!--          Determine topography-top index on scalar grid
    1695              k_wall = get_topography_top_index( j, i, 's' )
     1695             k_wall = get_topography_top_index_ji( j, i, 's' )
    16961696
    16971697             DO  k = k_wall + 1, k_wall + pch_index_ji(j,i)
Note: See TracChangeset for help on using the changeset viewer.