Home


About us
PALM
PALM-4U
Installation
Release notes
Documentation
Getting Help
Gallery
Developers area

Ignore:
Timestamp:
Jul 20, 2017 5:27:19 PM (7 years ago)
Author:
suehring
Message:

get topograpyh top index via function call

File:
1 edited

Legend:

Unmodified
Added
Removed

  • palm/trunk/SOURCE/microphysics_mod.f90

    r2312 r2317  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Get topography top index via Function call
    2323!
    2424! Former revisions:
    … …  
    445445
    446446       USE indices,                                                            &
    447            ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzb_max, nzt, wall_flags_0
     447           ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzt, wall_flags_0
    448448
    449449       USE kinds
    … …  
    787787
    788788       USE indices,                                                            &
    789            ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzb_max, nzt, wall_flags_0
     789           ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzt, wall_flags_0
    790790
    791791       USE kinds
    … …  
    923923
    924924       USE indices,                                                            &
    925            ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzb_max, nzt, wall_flags_0
     925           ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzt, wall_flags_0
    926926
    927927       USE kinds
     928
     929       USE surface_mod,                                                        &
     930           ONLY:  get_topography_top_index
    928931
    929932
    … …  
    942945!
    943946!--          Determine vertical index of topography top
    944              k_wall = MAXLOC(                                                  &
    945                           MERGE( 1, 0,                                         &
    946                                  BTEST( wall_flags_0(nzb:nzb_max,j,i), 12 )    &
    947                                ), DIM = 1                                      &
    948                             ) - 1
     947             k_wall = get_topography_top_index( j, i, 's' )
    949948             DO  k = nzb+1, nzt
    950949!
    … …  
    994993
    995994       USE indices,                                                            &
    996            ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzb_max, nzt, wall_flags_0
     995           ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzt, wall_flags_0
    997996
    998997       USE kinds
    … …  
    10931092
    10941093       USE indices,                                                            &
    1095            ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzb_max, nzt, wall_flags_0
     1094           ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzt, wall_flags_0
    10961095
    10971096       USE kinds
    … …  
    11701169
    11711170       USE indices,                                                            &
    1172            ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzb_max, nzt, wall_flags_0
     1171           ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzt, wall_flags_0
    11731172
    11741173       USE kinds
    … …  
    13311330
    13321331       USE indices,                                                            &
    1333            ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzb_max, nzt, wall_flags_0
     1332           ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzt, wall_flags_0
    13341333
    13351334       USE kinds
    … …  
    14431442
    14441443       USE indices,                                                            &
    1445            ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzb_max, nzt, wall_flags_0
     1444           ONLY:  nxlg, nxrg, nyng, nysg, nzb, nzt, wall_flags_0
    14461445
    14471446       USE kinds
    … …  
    23292328
    23302329       USE indices,                                                            &
    2331            ONLY:  nzb, nzb_max, nzt, wall_flags_0
     2330           ONLY:  nzb, nzt, wall_flags_0
    23322331
    23332332       USE kinds
     2333
     2334       USE surface_mod,                                                        &
     2335           ONLY:  get_topography_top_index
    23342336
    23352337
    … …  
    23462348!
    23472349!--    Determine vertical index of topography top
    2348        k_wall = MAXLOC(                                                        &
    2349                         MERGE( 1, 0,                                           &
    2350                                BTEST( wall_flags_0(nzb:nzb_max,j,i), 12 )      &
    2351                              ), DIM = 1                                        &
    2352                       ) - 1
     2350       k_wall = get_topography_top_index( j, i, 's' )
    23532351       DO  k = nzb+1, nzt
    23542352!
Note: See TracChangeset for help on using the changeset viewer.
                                                                                                                                                                                                                                                                                                                                                                               
  | Impressum | ©Leibniz Universität Hannover |