Changeset 3547 for palm/trunk/SOURCE


Ignore:
Timestamp:
Nov 21, 2018 1:21:24 PM (6 years ago)
Author:
suehring
Message:

variable description added some routines

Location:
palm/trunk/SOURCE
Files:
18 edited

Legend:

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

    r3538 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3538 2018-11-20 10:55:41Z suehring
    2730! Remove unnecessary double-masking of topography
    2831!
     
    126129       IMPLICIT NONE
    127130
    128        INTEGER(iwp) ::  i !<
    129        INTEGER(iwp) ::  j !<
    130        INTEGER(iwp) ::  k !<
     131       INTEGER(iwp) ::  i !< grid index along x-direction
     132       INTEGER(iwp) ::  j !< grid index along y-direction
     133       INTEGER(iwp) ::  k !< grid index along z-direction
    131134
    132135       REAL(wp)     ::  gu  !< local additional advective velocity
     
    197200       IMPLICIT NONE
    198201
    199        INTEGER(iwp) ::  i !<
    200        INTEGER(iwp) ::  j !<
    201        INTEGER(iwp) ::  k !<
     202       INTEGER(iwp) ::  i !< grid index along x-direction
     203       INTEGER(iwp) ::  j !< grid index along y-direction
     204       INTEGER(iwp) ::  k !< grid index along z-direction
    202205
    203206       REAL(wp)     ::  gu  !< local additional advective velocity
  • palm/trunk/SOURCE/advec_s_up.f90

    r3538 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3538 2018-11-20 10:55:41Z suehring
    2730! Remove unnecessary double-masking of topography
    2831!
     
    123126       IMPLICIT NONE
    124127
    125        INTEGER(iwp) ::  i !<
    126        INTEGER(iwp) ::  j !<
    127        INTEGER(iwp) ::  k !<
    128 
    129        REAL(wp) ::  ukomp !<
    130        REAL(wp) ::  vkomp !<
    131        REAL(wp) ::  wkomp !<
     128       INTEGER(iwp) ::  i !< grid index along x-direction
     129       INTEGER(iwp) ::  j !< grid index along y-direction
     130       INTEGER(iwp) ::  k !< grid index along z-direction
     131
     132       REAL(wp) ::  ukomp !< advection velocity along x-direction
     133       REAL(wp) ::  vkomp !< advection velocity along y-direction
     134       REAL(wp) ::  wkomp !< advection velocity along z-direction
    132135#if defined( __nopointer )
    133        REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  sk !<
     136       REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  sk !< treated scalar
    134137#else
    135        REAL(wp), DIMENSION(:,:,:), POINTER ::  sk
     138       REAL(wp), DIMENSION(:,:,:), POINTER ::  sk !< treated scalar
    136139#endif
    137140
     
    202205       IMPLICIT NONE
    203206
    204        INTEGER(iwp) ::  i !<
    205        INTEGER(iwp) ::  j !<
    206        INTEGER(iwp) ::  k !<
    207 
    208        REAL(wp) ::  ukomp !<
    209        REAL(wp) ::  vkomp !<
    210        REAL(wp) ::  wkomp !<
     207       INTEGER(iwp) ::  i !< grid index along x-direction
     208       INTEGER(iwp) ::  j !< grid index along y-direction
     209       INTEGER(iwp) ::  k !< grid index along z-direction
     210
     211       REAL(wp) ::  ukomp !< advection velocity along x-direction
     212       REAL(wp) ::  vkomp !< advection velocity along y-direction
     213       REAL(wp) ::  wkomp !< advection velocity along z-direction
    211214       
    212215#if defined( __nopointer )
    213        REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  sk !<
     216       REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  sk !< treated scalar
    214217#else
    215        REAL(wp), DIMENSION(:,:,:), POINTER ::  sk
     218       REAL(wp), DIMENSION(:,:,:), POINTER ::  sk !< treated scalar
    216219#endif
    217220
  • palm/trunk/SOURCE/advec_u_pw.f90

    r3538 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3538 2018-11-20 10:55:41Z suehring
    2730! Remove unnecessary double-masking of topography
    2831!
     
    115118       IMPLICIT NONE
    116119
    117        INTEGER(iwp) ::  i !<
    118        INTEGER(iwp) ::  j !<
    119        INTEGER(iwp) ::  k !<
     120       INTEGER(iwp) ::  i !< grid index along x-direction
     121       INTEGER(iwp) ::  j !< grid index along y-direction
     122       INTEGER(iwp) ::  k !< grid index along z-direction
    120123       
    121        REAL(wp)    ::  gu !<
    122        REAL(wp)    ::  gv !<
     124       REAL(wp)    ::  gu !< Galilei-transformation velocity along x
     125       REAL(wp)    ::  gv !< Galilei-transformation velocity along y
    123126 
    124127       gu = 2.0_wp * u_gtrans
     
    167170       IMPLICIT NONE
    168171
    169        INTEGER(iwp) ::  i !<
    170        INTEGER(iwp) ::  j !<
    171        INTEGER(iwp) ::  k !<
     172       INTEGER(iwp) ::  i !< grid index along x-direction
     173       INTEGER(iwp) ::  j !< grid index along y-direction
     174       INTEGER(iwp) ::  k !< grid index along z-direction
    172175       
    173        REAL(wp)    ::  gu !<
    174        REAL(wp)    ::  gv !<
     176       REAL(wp)    ::  gu !< Galilei-transformation velocity along x
     177       REAL(wp)    ::  gv !< Galilei-transformation velocity along y
    175178
    176179       gu = 2.0_wp * u_gtrans
  • palm/trunk/SOURCE/advec_u_up.f90

    r3538 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3538 2018-11-20 10:55:41Z suehring
    2730! Remove unnecessary double-masking of topography
    2831!
     
    114117       IMPLICIT NONE
    115118
    116        INTEGER(iwp) ::  i !<
    117        INTEGER(iwp) ::  j !<
    118        INTEGER(iwp) ::  k !<
    119 
    120        REAL(wp) ::  ukomp !<
    121        REAL(wp) ::  vkomp !<
    122        REAL(wp) ::  wkomp !<
     119       INTEGER(iwp) ::  i !< grid index along x-direction
     120       INTEGER(iwp) ::  j !< grid index along y-direction
     121       INTEGER(iwp) ::  k !< grid index along z-direction
     122
     123       REAL(wp) ::  ukomp !< advection velocity along x-direction
     124       REAL(wp) ::  vkomp !< advection velocity along y-direction
     125       REAL(wp) ::  wkomp !< advection velocity along z-direction
    123126
    124127       
     
    190193       IMPLICIT NONE
    191194
    192        INTEGER(iwp) ::  i !<
    193        INTEGER(iwp) ::  j !<
    194        INTEGER(iwp) ::  k !<
    195 
    196        REAL(wp) ::  ukomp !<
    197        REAL(wp) ::  vkomp !<
    198        REAL(wp) ::  wkomp !<
     195       INTEGER(iwp) ::  i !< grid index along x-direction
     196       INTEGER(iwp) ::  j !< grid index along y-direction
     197       INTEGER(iwp) ::  k !< grid index along z-direction
     198
     199       REAL(wp) ::  ukomp !< advection velocity along x-direction
     200       REAL(wp) ::  vkomp !< advection velocity along y-direction
     201       REAL(wp) ::  wkomp !< advection velocity along z-direction
    199202
    200203
  • palm/trunk/SOURCE/advec_v_pw.f90

    r3538 r3547  
    2626! -----------------
    2727! $Id$
     28! variables documented
     29!
     30! 3538 2018-11-20 10:55:41Z suehring
    2831! Remove unnecessary double-masking of topography
    2932!
     
    116119       IMPLICIT NONE
    117120
    118        INTEGER(iwp) ::  i !<
    119        INTEGER(iwp) ::  j !<
    120        INTEGER(iwp) ::  k !<
     121       INTEGER(iwp) ::  i !< grid index along x-direction
     122       INTEGER(iwp) ::  j !< grid index along y-direction
     123       INTEGER(iwp) ::  k !< grid index along z-direction
    121124       
    122        REAL(wp)    ::  gu !<
    123        REAL(wp)    ::  gv !<
     125       REAL(wp)    ::  gu !< Galilei-transformation velocity along x
     126       REAL(wp)    ::  gv !< Galilei-transformation velocity along y
    124127 
    125128
     
    169172       IMPLICIT NONE
    170173
    171        INTEGER(iwp) ::  i !<
    172        INTEGER(iwp) ::  j !<
    173        INTEGER(iwp) ::  k !<
     174       INTEGER(iwp) ::  i !< grid index along x-direction
     175       INTEGER(iwp) ::  j !< grid index along y-direction
     176       INTEGER(iwp) ::  k !< grid index along z-direction
    174177       
    175        REAL(wp)    ::  gu !<
    176        REAL(wp)    ::  gv !<
     178       REAL(wp)    ::  gu !< Galilei-transformation velocity along x
     179       REAL(wp)    ::  gv !< Galilei-transformation velocity along y
    177180
    178181
  • palm/trunk/SOURCE/advec_v_up.f90

    r3538 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3538 2018-11-20 10:55:41Z suehring
    2730! Remove unnecessary double-masking of topography
    2831!
     
    114117       IMPLICIT NONE
    115118
    116        INTEGER(iwp) ::  i !<
    117        INTEGER(iwp) ::  j !<
    118        INTEGER(iwp) ::  k !<
    119 
    120        REAL(wp) ::  ukomp !<
    121        REAL(wp) ::  vkomp !<
    122        REAL(wp) ::  wkomp !<      
     119       INTEGER(iwp) ::  i !< grid index along x-direction
     120       INTEGER(iwp) ::  j !< grid index along y-direction
     121       INTEGER(iwp) ::  k !< grid index along z-direction
     122
     123       REAL(wp) ::  ukomp !< advection velocity along x-direction
     124       REAL(wp) ::  vkomp !< advection velocity along y-direction
     125       REAL(wp) ::  wkomp !< advection velocity along z-direction      
    123126
    124127
     
    190193       IMPLICIT NONE
    191194
    192        INTEGER(iwp) ::  i !<
    193        INTEGER(iwp) ::  j !<
    194        INTEGER(iwp) ::  k !<
    195 
    196        REAL(wp) ::  ukomp !<
    197        REAL(wp) ::  vkomp !<
    198        REAL(wp) ::  wkomp !<
     195       INTEGER(iwp) ::  i !< grid index along x-direction
     196       INTEGER(iwp) ::  j !< grid index along y-direction
     197       INTEGER(iwp) ::  k !< grid index along z-direction
     198
     199       REAL(wp) ::  ukomp !< advection velocity along x-direction
     200       REAL(wp) ::  vkomp !< advection velocity along y-direction
     201       REAL(wp) ::  wkomp !< advection velocity along z-direction
    199202
    200203
  • palm/trunk/SOURCE/advec_w_pw.f90

    r3538 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3538 2018-11-20 10:55:41Z suehring
    2730! Remove unnecessary double-masking of topography
    2831!
     
    115118       IMPLICIT NONE
    116119
    117        INTEGER(iwp) ::  i !<
    118        INTEGER(iwp) ::  j !<
    119        INTEGER(iwp) ::  k !<
     120       INTEGER(iwp) ::  i !< grid index along x-direction
     121       INTEGER(iwp) ::  j !< grid index along y-direction
     122       INTEGER(iwp) ::  k !< grid index along z-direction
    120123       
    121        REAL(wp)    ::  gu !<
    122        REAL(wp)    ::  gv !<
     124       REAL(wp)    ::  gu !< Galilei-transformation velocity along x
     125       REAL(wp)    ::  gv !< Galilei-transformation velocity along y
    123126
    124127 
     
    168171       IMPLICIT NONE
    169172
    170        INTEGER(iwp) ::  i !<
    171        INTEGER(iwp) ::  j !<
    172        INTEGER(iwp) ::  k !<
     173       INTEGER(iwp) ::  i !< grid index along x-direction
     174       INTEGER(iwp) ::  j !< grid index along y-direction
     175       INTEGER(iwp) ::  k !< grid index along z-direction
    173176       
    174        REAL(wp)    ::  gu !<
    175        REAL(wp)    ::  gv !<
     177       REAL(wp)    ::  gu !< Galilei-transformation velocity along x
     178       REAL(wp)    ::  gv !< Galilei-transformation velocity along y
    176179
    177180       gu = 2.0_wp * u_gtrans
  • palm/trunk/SOURCE/advec_w_up.f90

    r3538 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3538 2018-11-20 10:55:41Z suehring
    2730! Remove unnecessary double-masking of topography
    2831!
     
    117120       IMPLICIT NONE
    118121
    119        INTEGER(iwp) ::  i !<
    120        INTEGER(iwp) ::  j !<
    121        INTEGER(iwp) ::  k !<
    122 
    123        REAL(wp) ::  ukomp !<
    124        REAL(wp) ::  vkomp !<
     122       INTEGER(iwp) ::  i !< grid index along x-direction
     123       INTEGER(iwp) ::  j !< grid index along y-direction
     124       INTEGER(iwp) ::  k !< grid index along z-direction
     125
     126       REAL(wp) ::  ukomp !< advection velocity along x-direction
     127       REAL(wp) ::  vkomp !< advection velocity along y-direction
    125128
    126129
     
    191194       IMPLICIT NONE
    192195
    193        INTEGER(iwp) ::  i !<
    194        INTEGER(iwp) ::  j !<
    195        INTEGER(iwp) ::  k !<
    196 
    197        REAL(wp) ::  ukomp !<
    198        REAL(wp) ::  vkomp !<
     196       INTEGER(iwp) ::  i !< grid index along x-direction
     197       INTEGER(iwp) ::  j !< grid index along y-direction
     198       INTEGER(iwp) ::  k !< grid index along z-direction
     199
     200       REAL(wp) ::  ukomp !< advection velocity along x-direction
     201       REAL(wp) ::  vkomp !< advection velocity along y-direction
    199202
    200203
  • palm/trunk/SOURCE/advec_ws.f90

    r3467 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3467 2018-10-30 19:05:21Z suehring
    2730! Implementation of a new aerosol module salsa.
    2831!
     
    11941197       IMPLICIT NONE
    11951198
    1196        CHARACTER (LEN = *), INTENT(IN) ::  sk_char !<
     1199       CHARACTER (LEN = *), INTENT(IN) ::  sk_char !< string identifier, used for assign fluxes to the correct dimension in the analysis array
    11971200       
    1198        INTEGER(iwp) ::  i     !<
    1199        INTEGER(iwp) ::  ibit0 !<
    1200        INTEGER(iwp) ::  ibit1 !<
    1201        INTEGER(iwp) ::  ibit2 !<
    1202        INTEGER(iwp) ::  ibit3 !<
    1203        INTEGER(iwp) ::  ibit4 !<
    1204        INTEGER(iwp) ::  ibit5 !<
    1205        INTEGER(iwp) ::  ibit6 !<
    1206        INTEGER(iwp) ::  ibit7 !<
    1207        INTEGER(iwp) ::  ibit8 !<
    1208        INTEGER(iwp) ::  i_omp !<
    1209        INTEGER(iwp) ::  j     !<
    1210        INTEGER(iwp) ::  k     !<
    1211        INTEGER(iwp) ::  k_mm  !<
    1212        INTEGER(iwp) ::  k_pp  !<
    1213        INTEGER(iwp) ::  k_ppp !<
    1214        INTEGER(iwp) ::  tn    !<
     1201       INTEGER(iwp) ::  i     !< grid index along x-direction
     1202       INTEGER(iwp) ::  ibit0 !< flag indicating 1st-order scheme along x-direction
     1203       INTEGER(iwp) ::  ibit1 !< flag indicating 3rd-order scheme along x-direction
     1204       INTEGER(iwp) ::  ibit2 !< flag indicating 5th-order scheme along x-direction
     1205       INTEGER(iwp) ::  ibit3 !< flag indicating 1st-order scheme along y-direction
     1206       INTEGER(iwp) ::  ibit4 !< flag indicating 3rd-order scheme along y-direction
     1207       INTEGER(iwp) ::  ibit5 !< flag indicating 5th-order scheme along y-direction
     1208       INTEGER(iwp) ::  ibit6 !< flag indicating 1st-order scheme along z-direction
     1209       INTEGER(iwp) ::  ibit7 !< flag indicating 3rd-order scheme along z-direction
     1210       INTEGER(iwp) ::  ibit8 !< flag indicating 5th-order scheme along z-direction
     1211       INTEGER(iwp) ::  i_omp !< leftmost index on subdomain, or in case of OpenMP, on thread
     1212       INTEGER(iwp) ::  j     !< grid index along y-direction
     1213       INTEGER(iwp) ::  k     !< grid index along z-direction
     1214       INTEGER(iwp) ::  k_mm  !< k-2 index in disretization, can be modified to avoid segmentation faults
     1215       INTEGER(iwp) ::  k_pp  !< k+2 index in disretization, can be modified to avoid segmentation faults
     1216       INTEGER(iwp) ::  k_ppp !< k+3 index in disretization, can be modified to avoid segmentation faults
     1217       INTEGER(iwp) ::  tn    !< number of OpenMP thread
    12151218       
    1216        REAL(wp)     ::  diss_d !<
    1217        REAL(wp)     ::  div    !<
    1218        REAL(wp)     ::  flux_d !<
    1219        REAL(wp)     ::  u_comp !<
    1220        REAL(wp)     ::  v_comp !<
     1219       REAL(wp)     ::  diss_d !< artificial dissipation term at grid box bottom
     1220       REAL(wp)     ::  div    !< diverence on scalar grid
     1221       REAL(wp)     ::  flux_d !< 6th-order flux at grid box bottom
     1222       REAL(wp)     ::  u_comp !< advection velocity along x-direction
     1223       REAL(wp)     ::  v_comp !< advection velocity along y-direction
    12211224       
    12221225#if defined( __nopointer )
    1223        REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  sk !<
     1226       REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  sk !< advected scalar
    12241227#else
    1225        REAL(wp), DIMENSION(:,:,:), POINTER    ::  sk     !<
     1228       REAL(wp), DIMENSION(:,:,:), POINTER    ::  sk     !< advected scalar
    12261229#endif
    1227        REAL(wp), DIMENSION(nzb:nzt+1)         ::  diss_n !<
    1228        REAL(wp), DIMENSION(nzb:nzt+1)         ::  diss_r !<
    1229        REAL(wp), DIMENSION(nzb:nzt+1)         ::  diss_t !<
    1230        REAL(wp), DIMENSION(nzb:nzt+1)         ::  flux_n !<
    1231        REAL(wp), DIMENSION(nzb:nzt+1)         ::  flux_r !<
    1232        REAL(wp), DIMENSION(nzb:nzt+1)         ::  flux_t !<
     1230       REAL(wp), DIMENSION(nzb:nzt+1)         ::  diss_n !< discretized artificial dissipation at northward-side of the grid box
     1231       REAL(wp), DIMENSION(nzb:nzt+1)         ::  diss_r !< discretized artificial dissipation at rightward-side of the grid box
     1232       REAL(wp), DIMENSION(nzb:nzt+1)         ::  diss_t !< discretized artificial dissipation at top of the grid box
     1233       REAL(wp), DIMENSION(nzb:nzt+1)         ::  flux_n !< discretized 6th-order flux at northward-side of the grid box
     1234       REAL(wp), DIMENSION(nzb:nzt+1)         ::  flux_r !< discretized 6th-order flux at rightward-side of the grid box
     1235       REAL(wp), DIMENSION(nzb:nzt+1)         ::  flux_t !< discretized 6th-order flux at top of the grid box
    12331236       
    1234        REAL(wp), DIMENSION(nzb+1:nzt,0:threads_per_task-1) ::  swap_diss_y_local !<
    1235        REAL(wp), DIMENSION(nzb+1:nzt,0:threads_per_task-1) ::  swap_flux_y_local !<
     1237       REAL(wp), DIMENSION(nzb+1:nzt,0:threads_per_task-1) ::  swap_diss_y_local !< discretized artificial dissipation at southward-side of the grid box
     1238       REAL(wp), DIMENSION(nzb+1:nzt,0:threads_per_task-1) ::  swap_flux_y_local !< discretized 6th-order flux at northward-side of the grid box
    12361239       
    1237        REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn,0:threads_per_task-1) ::  swap_diss_x_local !<
    1238        REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn,0:threads_per_task-1) ::  swap_flux_x_local !<
     1240       REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn,0:threads_per_task-1) ::  swap_diss_x_local !< discretized artificial dissipation at leftward-side of the grid box
     1241       REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn,0:threads_per_task-1) ::  swap_flux_x_local !< discretized 6th-order flux at leftward-side of the grid box
    12391242       
    12401243
     
    17731776       IMPLICIT NONE
    17741777
    1775        INTEGER(iwp) ::  i      !<
    1776        INTEGER(iwp) ::  ibit9  !<
    1777        INTEGER(iwp) ::  ibit10 !<
    1778        INTEGER(iwp) ::  ibit11 !<
    1779        INTEGER(iwp) ::  ibit12 !<
    1780        INTEGER(iwp) ::  ibit13 !<
    1781        INTEGER(iwp) ::  ibit14 !<
    1782        INTEGER(iwp) ::  ibit15 !<
    1783        INTEGER(iwp) ::  ibit16 !<
    1784        INTEGER(iwp) ::  ibit17 !<
    1785        INTEGER(iwp) ::  i_omp  !<
    1786        INTEGER(iwp) ::  j      !<
    1787        INTEGER(iwp) ::  k      !<
    1788        INTEGER(iwp) ::  k_mm   !<
    1789        INTEGER(iwp) ::  k_pp   !<
    1790        INTEGER(iwp) ::  k_ppp  !<
    1791        INTEGER(iwp) ::  tn     !<
     1778       INTEGER(iwp) ::  i      !< grid index along x-direction
     1779       INTEGER(iwp) ::  ibit9  !< flag indicating 1st-order scheme along x-direction
     1780       INTEGER(iwp) ::  ibit10 !< flag indicating 3rd-order scheme along x-direction
     1781       INTEGER(iwp) ::  ibit11 !< flag indicating 5th-order scheme along x-direction
     1782       INTEGER(iwp) ::  ibit12 !< flag indicating 1st-order scheme along y-direction
     1783       INTEGER(iwp) ::  ibit13 !< flag indicating 3rd-order scheme along y-direction
     1784       INTEGER(iwp) ::  ibit14 !< flag indicating 5th-order scheme along y-direction
     1785       INTEGER(iwp) ::  ibit15 !< flag indicating 1st-order scheme along z-direction
     1786       INTEGER(iwp) ::  ibit16 !< flag indicating 3rd-order scheme along z-direction
     1787       INTEGER(iwp) ::  ibit17 !< flag indicating 5th-order scheme along z-direction
     1788       INTEGER(iwp) ::  i_omp  !< leftmost index on subdomain, or in case of OpenMP, on thread
     1789       INTEGER(iwp) ::  j      !< grid index along y-direction
     1790       INTEGER(iwp) ::  k      !< grid index along z-direction
     1791       INTEGER(iwp) ::  k_mm   !< k-2 index in disretization, can be modified to avoid segmentation faults
     1792       INTEGER(iwp) ::  k_pp   !< k+2 index in disretization, can be modified to avoid segmentation faults
     1793       INTEGER(iwp) ::  k_ppp  !< k+3 index in disretization, can be modified to avoid segmentation faults
     1794       INTEGER(iwp) ::  tn     !< number of OpenMP thread
    17921795       
    1793        REAL(wp)    ::  diss_d   !<
    1794        REAL(wp)    ::  div      !<
    1795        REAL(wp)    ::  flux_d   !<
    1796        REAL(wp)    ::  gu       !<
    1797        REAL(wp)    ::  gv       !<
    1798        REAL(wp)    ::  u_comp_l !<
     1796       REAL(wp)    ::  diss_d   !< artificial dissipation term at grid box bottom
     1797       REAL(wp)    ::  div      !< diverence on u-grid
     1798       REAL(wp)    ::  flux_d   !< 6th-order flux at grid box bottom
     1799       REAL(wp)    ::  gu       !< Galilei-transformation velocity along x
     1800       REAL(wp)    ::  gv       !< Galilei-transformation velocity along y
     1801       REAL(wp)    ::  u_comp_l !< advection velocity along x at leftmost grid point on subdomain
    17991802       
    1800        REAL(wp), DIMENSION(nzb:nzt+1) ::  diss_n !<
    1801        REAL(wp), DIMENSION(nzb:nzt+1) ::  diss_r !<
    1802        REAL(wp), DIMENSION(nzb:nzt+1) ::  diss_t !<
    1803        REAL(wp), DIMENSION(nzb:nzt+1) ::  flux_n !<
    1804        REAL(wp), DIMENSION(nzb:nzt+1) ::  flux_r !<
    1805        REAL(wp), DIMENSION(nzb:nzt+1) ::  flux_t !<
    1806        REAL(wp), DIMENSION(nzb:nzt+1) ::  u_comp !<
    1807        REAL(wp), DIMENSION(nzb:nzt+1) ::  v_comp !<
    1808        REAL(wp), DIMENSION(nzb:nzt+1) ::  w_comp !<
     1803       REAL(wp), DIMENSION(nzb:nzt+1) ::  diss_n !< discretized artificial dissipation at northward-side of the grid box
     1804       REAL(wp), DIMENSION(nzb:nzt+1) ::  diss_r !< discretized artificial dissipation at rightward-side of the grid box
     1805       REAL(wp), DIMENSION(nzb:nzt+1) ::  diss_t !< discretized artificial dissipation at top of the grid box
     1806       REAL(wp), DIMENSION(nzb:nzt+1) ::  flux_n !< discretized 6th-order flux at northward-side of the grid box
     1807       REAL(wp), DIMENSION(nzb:nzt+1) ::  flux_r !< discretized 6th-order flux at rightward-side of the grid box
     1808       REAL(wp), DIMENSION(nzb:nzt+1) ::  flux_t !< discretized 6th-order flux at top of the grid box
     1809       REAL(wp), DIMENSION(nzb:nzt+1) ::  u_comp !< advection velocity along x
     1810       REAL(wp), DIMENSION(nzb:nzt+1) ::  v_comp !< advection velocity along y
     1811       REAL(wp), DIMENSION(nzb:nzt+1) ::  w_comp !< advection velocity along z
    18091812
    18101813       gu = 2.0_wp * u_gtrans
     
    22682271       IMPLICIT NONE
    22692272
    2270        INTEGER(iwp)  ::  i      !<
    2271        INTEGER(iwp)  ::  ibit18 !<
    2272        INTEGER(iwp)  ::  ibit19 !<
    2273        INTEGER(iwp)  ::  ibit20 !<
    2274        INTEGER(iwp)  ::  ibit21 !<
    2275        INTEGER(iwp)  ::  ibit22 !<
    2276        INTEGER(iwp)  ::  ibit23 !<
    2277        INTEGER(iwp)  ::  ibit24 !<
    2278        INTEGER(iwp)  ::  ibit25 !<
    2279        INTEGER(iwp)  ::  ibit26 !<
    2280        INTEGER(iwp)  ::  i_omp  !<
    2281        INTEGER(iwp)  ::  j      !<
    2282        INTEGER(iwp)  ::  k      !<
    2283        INTEGER(iwp)  ::  k_mm   !<
    2284        INTEGER(iwp)  ::  k_pp   !<
    2285        INTEGER(iwp)  ::  k_ppp  !<
    2286        INTEGER(iwp)  ::  tn     !<
     2273       INTEGER(iwp)  ::  i      !< grid index along x-direction
     2274       INTEGER(iwp)  ::  ibit18 !< flag indicating 1st-order scheme along x-direction
     2275       INTEGER(iwp)  ::  ibit19 !< flag indicating 3rd-order scheme along x-direction
     2276       INTEGER(iwp)  ::  ibit20 !< flag indicating 5th-order scheme along x-direction
     2277       INTEGER(iwp)  ::  ibit21 !< flag indicating 1st-order scheme along y-direction
     2278       INTEGER(iwp)  ::  ibit22 !< flag indicating 3rd-order scheme along y-direction
     2279       INTEGER(iwp)  ::  ibit23 !< flag indicating 3rd-order scheme along y-direction
     2280       INTEGER(iwp)  ::  ibit24 !< flag indicating 1st-order scheme along z-direction
     2281       INTEGER(iwp)  ::  ibit25 !< flag indicating 3rd-order scheme along z-direction
     2282       INTEGER(iwp)  ::  ibit26 !< flag indicating 3rd-order scheme along z-direction
     2283       INTEGER(iwp)  ::  i_omp  !< leftmost index on subdomain, or in case of OpenMP, on thread
     2284       INTEGER(iwp)  ::  j      !< grid index along y-direction
     2285       INTEGER(iwp)  ::  k      !< grid index along z-direction
     2286       INTEGER(iwp)  ::  k_mm   !< k-2 index in disretization, can be modified to avoid segmentation faults
     2287       INTEGER(iwp)  ::  k_pp   !< k+2 index in disretization, can be modified to avoid segmentation faults
     2288       INTEGER(iwp)  ::  k_ppp  !< k+3 index in disretization, can be modified to avoid segmentation faults
     2289       INTEGER(iwp)  ::  tn     !< number of OpenMP thread
    22872290       
    2288        REAL(wp)     ::  diss_d   !<
    2289        REAL(wp)     ::  div      !<
    2290        REAL(wp)     ::  flux_d   !<
    2291        REAL(wp)     ::  gu       !<
    2292        REAL(wp)     ::  gv       !<
    2293        REAL(wp)     ::  v_comp_l !<
     2291       REAL(wp)     ::  diss_d   !< artificial dissipation term at grid box bottom
     2292       REAL(wp)     ::  div      !< divergence on v-grid
     2293       REAL(wp)     ::  flux_d   !< 6th-order flux at grid box bottom
     2294       REAL(wp)     ::  gu       !< Galilei-transformation velocity along x
     2295       REAL(wp)     ::  gv       !< Galilei-transformation velocity along y
     2296       REAL(wp)     ::  v_comp_l !< advection velocity along y on leftmost grid point on subdomain
    22942297       
    2295        REAL(wp), DIMENSION(nzb:nzt+1)  ::  diss_n !<
    2296        REAL(wp), DIMENSION(nzb:nzt+1)  ::  diss_r !<
    2297        REAL(wp), DIMENSION(nzb:nzt+1)  ::  diss_t !<
    2298        REAL(wp), DIMENSION(nzb:nzt+1)  ::  flux_n !<
    2299        REAL(wp), DIMENSION(nzb:nzt+1)  ::  flux_r !<
    2300        REAL(wp), DIMENSION(nzb:nzt+1)  ::  flux_t !<
    2301        REAL(wp), DIMENSION(nzb:nzt+1)  ::  u_comp !<
    2302        REAL(wp), DIMENSION(nzb:nzt+1)  ::  v_comp !<
    2303        REAL(wp), DIMENSION(nzb:nzt+1)  ::  w_comp !<
     2298       REAL(wp), DIMENSION(nzb:nzt+1)  ::  diss_n !< discretized artificial dissipation at northward-side of the grid box
     2299       REAL(wp), DIMENSION(nzb:nzt+1)  ::  diss_r !< discretized artificial dissipation at rightward-side of the grid box
     2300       REAL(wp), DIMENSION(nzb:nzt+1)  ::  diss_t !< discretized artificial dissipation at top of the grid box
     2301       REAL(wp), DIMENSION(nzb:nzt+1)  ::  flux_n !< discretized 6th-order flux at northward-side of the grid box
     2302       REAL(wp), DIMENSION(nzb:nzt+1)  ::  flux_r !< discretized 6th-order flux at rightward-side of the grid box
     2303       REAL(wp), DIMENSION(nzb:nzt+1)  ::  flux_t !< discretized 6th-order flux at top of the grid box
     2304       REAL(wp), DIMENSION(nzb:nzt+1)  ::  u_comp !< advection velocity along x
     2305       REAL(wp), DIMENSION(nzb:nzt+1)  ::  v_comp !< advection velocity along y
     2306       REAL(wp), DIMENSION(nzb:nzt+1)  ::  w_comp !< advection velocity along z
    23042307
    23052308       gu = 2.0_wp * u_gtrans
     
    27702773       IMPLICIT NONE
    27712774
    2772        INTEGER(iwp) ::  i      !<
    2773        INTEGER(iwp) ::  ibit27 !<
    2774        INTEGER(iwp) ::  ibit28 !<
    2775        INTEGER(iwp) ::  ibit29 !<
    2776        INTEGER(iwp) ::  ibit30 !<
    2777        INTEGER(iwp) ::  ibit31 !<
    2778        INTEGER(iwp) ::  ibit32 !<
    2779        INTEGER(iwp) ::  ibit33 !<
    2780        INTEGER(iwp) ::  ibit34 !<
    2781        INTEGER(iwp) ::  ibit35 !<
    2782        INTEGER(iwp) ::  i_omp  !<
    2783        INTEGER(iwp) ::  j      !<
    2784        INTEGER(iwp) ::  k      !<
    2785        INTEGER(iwp) ::  k_mm   !<
    2786        INTEGER(iwp) ::  k_pp   !<
    2787        INTEGER(iwp) ::  k_ppp  !<
    2788        INTEGER(iwp) ::  tn     !<
     2775       INTEGER(iwp) ::  i      !< grid index along x-direction
     2776       INTEGER(iwp) ::  ibit27 !< flag indicating 1st-order scheme along x-direction
     2777       INTEGER(iwp) ::  ibit28 !< flag indicating 3rd-order scheme along x-direction
     2778       INTEGER(iwp) ::  ibit29 !< flag indicating 5th-order scheme along x-direction
     2779       INTEGER(iwp) ::  ibit30 !< flag indicating 1st-order scheme along y-direction
     2780       INTEGER(iwp) ::  ibit31 !< flag indicating 3rd-order scheme along y-direction
     2781       INTEGER(iwp) ::  ibit32 !< flag indicating 5th-order scheme along y-direction
     2782       INTEGER(iwp) ::  ibit33 !< flag indicating 1st-order scheme along z-direction
     2783       INTEGER(iwp) ::  ibit34 !< flag indicating 3rd-order scheme along z-direction
     2784       INTEGER(iwp) ::  ibit35 !< flag indicating 5th-order scheme along z-direction
     2785       INTEGER(iwp) ::  i_omp  !< leftmost index on subdomain, or in case of OpenMP, on thread
     2786       INTEGER(iwp) ::  j      !< grid index along y-direction
     2787       INTEGER(iwp) ::  k      !< grid index along z-direction
     2788       INTEGER(iwp) ::  k_mm   !< k-2 index in disretization, can be modified to avoid segmentation faults
     2789       INTEGER(iwp) ::  k_pp   !< k+2 index in disretization, can be modified to avoid segmentation faults
     2790       INTEGER(iwp) ::  k_ppp  !< k+3 index in disretization, can be modified to avoid segmentation faults
     2791       INTEGER(iwp) ::  tn     !< number of OpenMP thread
    27892792       
    2790        REAL(wp)    ::  diss_d  !<
    2791        REAL(wp)    ::  div     !<
    2792        REAL(wp)    ::  flux_d  !<
    2793        REAL(wp)    ::  gu      !<
    2794        REAL(wp)    ::  gv      !<
     2793       REAL(wp)    ::  diss_d  !< discretized artificial dissipation at top of the grid box
     2794       REAL(wp)    ::  div     !< divergence on w-grid
     2795       REAL(wp)    ::  flux_d  !< discretized 6th-order flux at top of the grid box
     2796       REAL(wp)    ::  gu      !< Galilei-transformation velocity along x
     2797       REAL(wp)    ::  gv      !< Galilei-transformation velocity along y
    27952798       
    2796        REAL(wp), DIMENSION(nzb:nzt+1)  ::  diss_n !<
    2797        REAL(wp), DIMENSION(nzb:nzt+1)  ::  diss_r !<
    2798        REAL(wp), DIMENSION(nzb:nzt+1)  ::  diss_t !<
    2799        REAL(wp), DIMENSION(nzb:nzt+1)  ::  flux_n !<
    2800        REAL(wp), DIMENSION(nzb:nzt+1)  ::  flux_r !<
    2801        REAL(wp), DIMENSION(nzb:nzt+1)  ::  flux_t !<
    2802        REAL(wp), DIMENSION(nzb:nzt+1)  ::  u_comp !<
    2803        REAL(wp), DIMENSION(nzb:nzt+1)  ::  v_comp !<
    2804        REAL(wp), DIMENSION(nzb:nzt+1)  ::  w_comp !<
     2799       REAL(wp), DIMENSION(nzb:nzt+1)  ::  diss_n !< discretized artificial dissipation at northward-side of the grid box
     2800       REAL(wp), DIMENSION(nzb:nzt+1)  ::  diss_r !< discretized artificial dissipation at rightward-side of the grid box
     2801       REAL(wp), DIMENSION(nzb:nzt+1)  ::  diss_t !< discretized artificial dissipation at top of the grid box
     2802       REAL(wp), DIMENSION(nzb:nzt+1)  ::  flux_n !< discretized 6th-order flux at northward-side of the grid box
     2803       REAL(wp), DIMENSION(nzb:nzt+1)  ::  flux_r !< discretized 6th-order flux at rightward-side of the grid box
     2804       REAL(wp), DIMENSION(nzb:nzt+1)  ::  flux_t !< discretized 6th-order flux at top of the grid box
     2805       REAL(wp), DIMENSION(nzb:nzt+1)  ::  u_comp !< advection velocity along x
     2806       REAL(wp), DIMENSION(nzb:nzt+1)  ::  v_comp !< advection velocity along y
     2807       REAL(wp), DIMENSION(nzb:nzt+1)  ::  w_comp !< advection velocity along z
    28052808
    28062809       gu = 2.0_wp * u_gtrans
     
    32543257       IMPLICIT NONE
    32553258
    3256        CHARACTER (LEN = *), INTENT(IN)    ::  sk_char !<
     3259       CHARACTER (LEN = *), INTENT(IN)    ::  sk_char !< string identifier, used for assign fluxes to the correct dimension in the analysis array
    32573260       
    3258        INTEGER(iwp) ::  i      !<
    3259        INTEGER(iwp) ::  ibit0  !<
    3260        INTEGER(iwp) ::  ibit1  !<
    3261        INTEGER(iwp) ::  ibit2  !<
    3262        INTEGER(iwp) ::  ibit3  !<
    3263        INTEGER(iwp) ::  ibit4  !<
    3264        INTEGER(iwp) ::  ibit5  !<
    3265        INTEGER(iwp) ::  ibit6  !<
    3266        INTEGER(iwp) ::  ibit7  !<
    3267        INTEGER(iwp) ::  ibit8  !<
    3268        INTEGER(iwp) ::  j      !<
    3269        INTEGER(iwp) ::  k      !<
    3270        INTEGER(iwp) ::  k_mm   !<
    3271        INTEGER(iwp) ::  k_pp   !<
    3272        INTEGER(iwp) ::  k_ppp  !<
    3273        INTEGER(iwp) ::  tn = 0 !<
     3261       INTEGER(iwp) ::  i      !< grid index along x-direction
     3262       INTEGER(iwp) ::  ibit0  !< flag indicating 1st-order scheme along x-direction
     3263       INTEGER(iwp) ::  ibit1  !< flag indicating 3rd-order scheme along x-direction
     3264       INTEGER(iwp) ::  ibit2  !< flag indicating 5th-order scheme along x-direction
     3265       INTEGER(iwp) ::  ibit3  !< flag indicating 1st-order scheme along y-direction
     3266       INTEGER(iwp) ::  ibit4  !< flag indicating 3rd-order scheme along y-direction
     3267       INTEGER(iwp) ::  ibit5  !< flag indicating 5th-order scheme along y-direction
     3268       INTEGER(iwp) ::  ibit6  !< flag indicating 1st-order scheme along z-direction
     3269       INTEGER(iwp) ::  ibit7  !< flag indicating 3rd-order scheme along z-direction
     3270       INTEGER(iwp) ::  ibit8  !< flag indicating 5th-order scheme along z-direction
     3271       INTEGER(iwp) ::  j      !< grid index along y-direction
     3272       INTEGER(iwp) ::  k      !< grid index along z-direction
     3273       INTEGER(iwp) ::  k_mm   !< k-2 index in disretization, can be modified to avoid segmentation faults
     3274       INTEGER(iwp) ::  k_pp   !< k+2 index in disretization, can be modified to avoid segmentation faults
     3275       INTEGER(iwp) ::  k_ppp  !< k+3 index in disretization, can be modified to avoid segmentation faults
     3276       INTEGER(iwp) ::  tn = 0 !< number of OpenMP thread
    32743277       
    32753278#if defined( __nopointer )
    3276        REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  sk !<
     3279       REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  sk !< advected scalar
    32773280#else
    3278        REAL(wp), DIMENSION(:,:,:), POINTER ::  sk !<
     3281       REAL(wp), DIMENSION(:,:,:), POINTER ::  sk !< advected scalar
    32793282#endif
    32803283
    3281        REAL(wp) ::  diss_d !<
    3282        REAL(wp) ::  div    !<
    3283        REAL(wp) ::  flux_d !<
    3284        REAL(wp) ::  u_comp !<
    3285        REAL(wp) ::  v_comp !<
     3284       REAL(wp) ::  diss_d !< artificial dissipation term at grid box bottom
     3285       REAL(wp) ::  div    !< diverence on scalar grid
     3286       REAL(wp) ::  flux_d !< 6th-order flux at grid box bottom
     3287       REAL(wp) ::  u_comp !< advection velocity along x-direction
     3288       REAL(wp) ::  v_comp !< advection velocity along y-direction
    32863289       
    3287        REAL(wp), DIMENSION(nzb:nzt)   ::  diss_n !<
    3288        REAL(wp), DIMENSION(nzb:nzt)   ::  diss_r !<
    3289        REAL(wp), DIMENSION(nzb:nzt)   ::  diss_t !<
    3290        REAL(wp), DIMENSION(nzb:nzt)   ::  flux_n !<
    3291        REAL(wp), DIMENSION(nzb:nzt)   ::  flux_r !<
    3292        REAL(wp), DIMENSION(nzb:nzt)   ::  flux_t !<
     3290       REAL(wp), DIMENSION(nzb:nzt)   ::  diss_n !< discretized artificial dissipation at northward-side of the grid box
     3291       REAL(wp), DIMENSION(nzb:nzt)   ::  diss_r !< discretized artificial dissipation at rightward-side of the grid box
     3292       REAL(wp), DIMENSION(nzb:nzt)   ::  diss_t !< discretized artificial dissipation at rightward-side of the grid box
     3293       REAL(wp), DIMENSION(nzb:nzt)   ::  flux_n !< discretized 6th-order flux at northward-side of the grid box
     3294       REAL(wp), DIMENSION(nzb:nzt)   ::  flux_r !< discretized 6th-order flux at rightward-side of the grid box
     3295       REAL(wp), DIMENSION(nzb:nzt)   ::  flux_t !< discretized 6th-order flux at rightward-side of the grid box
    32933296       
    3294        REAL(wp), DIMENSION(nzb+1:nzt) ::  swap_diss_y_local !<
    3295        REAL(wp), DIMENSION(nzb+1:nzt) ::  swap_flux_y_local !<
     3297       REAL(wp), DIMENSION(nzb+1:nzt) ::  swap_diss_y_local !< discretized artificial dissipation term at southward-side of the grid box
     3298       REAL(wp), DIMENSION(nzb+1:nzt) ::  swap_flux_y_local !< discretized 6th-order flux at northward-side of the grid box
    32963299       
    3297        REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_diss_x_local !<
    3298        REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_flux_x_local !<
     3300       REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_diss_x_local !< discretized artificial dissipation term at leftward-side of the grid box
     3301       REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_flux_x_local !< discretized 6th-order flux at leftward-side of the grid box
    32993302       
    33003303
     
    38133816       IMPLICIT NONE
    38143817
    3815        INTEGER(iwp) ::  i      !<
    3816        INTEGER(iwp) ::  ibit9  !<
    3817        INTEGER(iwp) ::  ibit10 !<
    3818        INTEGER(iwp) ::  ibit11 !<
    3819        INTEGER(iwp) ::  ibit12 !<
    3820        INTEGER(iwp) ::  ibit13 !<
    3821        INTEGER(iwp) ::  ibit14 !<
    3822        INTEGER(iwp) ::  ibit15 !<
    3823        INTEGER(iwp) ::  ibit16 !<
    3824        INTEGER(iwp) ::  ibit17 !<
    3825        INTEGER(iwp) ::  j      !<
    3826        INTEGER(iwp) ::  k      !<
    3827        INTEGER(iwp) ::  k_mm   !<
    3828        INTEGER(iwp) ::  k_pp   !<
    3829        INTEGER(iwp) ::  k_ppp  !<
    3830        INTEGER(iwp) ::  tn = 0 !<
     3818       INTEGER(iwp) ::  i      !< grid index along x-direction
     3819       INTEGER(iwp) ::  ibit9  !< flag indicating 1st-order scheme along x-direction
     3820       INTEGER(iwp) ::  ibit10 !< flag indicating 3rd-order scheme along x-direction
     3821       INTEGER(iwp) ::  ibit11 !< flag indicating 5th-order scheme along x-direction
     3822       INTEGER(iwp) ::  ibit12 !< flag indicating 1st-order scheme along y-direction
     3823       INTEGER(iwp) ::  ibit13 !< flag indicating 3rd-order scheme along y-direction
     3824       INTEGER(iwp) ::  ibit14 !< flag indicating 5th-order scheme along y-direction
     3825       INTEGER(iwp) ::  ibit15 !< flag indicating 1st-order scheme along z-direction
     3826       INTEGER(iwp) ::  ibit16 !< flag indicating 3rd-order scheme along z-direction
     3827       INTEGER(iwp) ::  ibit17 !< flag indicating 5th-order scheme along z-direction
     3828       INTEGER(iwp) ::  j      !< grid index along y-direction
     3829       INTEGER(iwp) ::  k      !< grid index along z-direction
     3830       INTEGER(iwp) ::  k_mm   !< k-2 index in disretization, can be modified to avoid segmentation faults
     3831       INTEGER(iwp) ::  k_pp   !< k+2 index in disretization, can be modified to avoid segmentation faults
     3832       INTEGER(iwp) ::  k_ppp  !< k+3 index in disretization, can be modified to avoid segmentation faults
     3833       INTEGER(iwp) ::  tn = 0 !< number of OpenMP thread
    38313834       
    3832        REAL(wp)    ::  diss_d !<
    3833        REAL(wp)    ::  div    !<
    3834        REAL(wp)    ::  flux_d !<
    3835        REAL(wp)    ::  gu     !<
    3836        REAL(wp)    ::  gv     !<
    3837        REAL(wp)    ::  v_comp !<
    3838        REAL(wp)    ::  w_comp !<
     3835       REAL(wp)    ::  diss_d !< artificial dissipation term at grid box bottom
     3836       REAL(wp)    ::  div    !< diverence on u-grid
     3837       REAL(wp)    ::  flux_d !< 6th-order flux at grid box bottom
     3838       REAL(wp)    ::  gu     !< Galilei-transformation velocity along x
     3839       REAL(wp)    ::  gv     !< Galilei-transformation velocity along y
     3840       REAL(wp)    ::  v_comp !< advection velocity along y
     3841       REAL(wp)    ::  w_comp !< advection velocity along z
    38393842       
    3840        REAL(wp), DIMENSION(nzb+1:nzt) ::  swap_diss_y_local_u !<
    3841        REAL(wp), DIMENSION(nzb+1:nzt) ::  swap_flux_y_local_u !<
     3843       REAL(wp), DIMENSION(nzb+1:nzt) ::  swap_diss_y_local_u !< discretized artificial dissipation at southward-side of the grid box
     3844       REAL(wp), DIMENSION(nzb+1:nzt) ::  swap_flux_y_local_u !< discretized 6th-order flux at southward-side of the grid box
    38423845       
    3843        REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_diss_x_local_u !<
    3844        REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_flux_x_local_u !<
     3846       REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_diss_x_local_u !< discretized artificial dissipation at leftward-side of the grid box
     3847       REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_flux_x_local_u !< discretized 6th-order flux at leftward-side of the grid box
    38453848       
    3846        REAL(wp), DIMENSION(nzb:nzt) ::  diss_n !<
    3847        REAL(wp), DIMENSION(nzb:nzt) ::  diss_r !<
    3848        REAL(wp), DIMENSION(nzb:nzt) ::  diss_t !<
    3849        REAL(wp), DIMENSION(nzb:nzt) ::  flux_n !<
    3850        REAL(wp), DIMENSION(nzb:nzt) ::  flux_r !<
    3851        REAL(wp), DIMENSION(nzb:nzt) ::  flux_t !<
    3852        REAL(wp), DIMENSION(nzb:nzt) ::  u_comp !<
     3849       REAL(wp), DIMENSION(nzb:nzt) ::  diss_n !< discretized artificial dissipation at northward-side of the grid box
     3850       REAL(wp), DIMENSION(nzb:nzt) ::  diss_r !< discretized artificial dissipation at leftward-side of the grid box
     3851       REAL(wp), DIMENSION(nzb:nzt) ::  diss_t !< discretized artificial dissipation at top of the grid box
     3852       REAL(wp), DIMENSION(nzb:nzt) ::  flux_n !< discretized 6th-order flux at northward-side of the grid box
     3853       REAL(wp), DIMENSION(nzb:nzt) ::  flux_r !< discretized 6th-order flux at rightward-side of the grid box
     3854       REAL(wp), DIMENSION(nzb:nzt) ::  flux_t !< discretized 6th-order flux at top of the grid box
     3855       REAL(wp), DIMENSION(nzb:nzt) ::  u_comp !< advection velocity along x
    38533856 
    38543857       gu = 2.0_wp * u_gtrans
     
    43054308
    43064309
    4307        INTEGER(iwp) ::  i      !<
    4308        INTEGER(iwp) ::  ibit18 !<
    4309        INTEGER(iwp) ::  ibit19 !<
    4310        INTEGER(iwp) ::  ibit20 !<
    4311        INTEGER(iwp) ::  ibit21 !<
    4312        INTEGER(iwp) ::  ibit22 !<
    4313        INTEGER(iwp) ::  ibit23 !<
    4314        INTEGER(iwp) ::  ibit24 !<
    4315        INTEGER(iwp) ::  ibit25 !<
    4316        INTEGER(iwp) ::  ibit26 !<
    4317        INTEGER(iwp) ::  j      !<
    4318        INTEGER(iwp) ::  k      !<
    4319        INTEGER(iwp) ::  k_mm   !<
    4320        INTEGER(iwp) ::  k_pp   !<
    4321        INTEGER(iwp) ::  k_ppp  !<
    4322        INTEGER(iwp) ::  tn = 0 !<
     4310       INTEGER(iwp) ::  i      !< grid index along x-direction
     4311       INTEGER(iwp) ::  ibit18 !< flag indicating 1st-order scheme along x-direction
     4312       INTEGER(iwp) ::  ibit19 !< flag indicating 3rd-order scheme along x-direction
     4313       INTEGER(iwp) ::  ibit20 !< flag indicating 5th-order scheme along x-direction
     4314       INTEGER(iwp) ::  ibit21 !< flag indicating 1st-order scheme along y-direction
     4315       INTEGER(iwp) ::  ibit22 !< flag indicating 3rd-order scheme along y-direction
     4316       INTEGER(iwp) ::  ibit23 !< flag indicating 5th-order scheme along y-direction
     4317       INTEGER(iwp) ::  ibit24 !< flag indicating 1st-order scheme along z-direction
     4318       INTEGER(iwp) ::  ibit25 !< flag indicating 3rd-order scheme along z-direction
     4319       INTEGER(iwp) ::  ibit26 !< flag indicating 5th-order scheme along z-direction
     4320       INTEGER(iwp) ::  j      !< grid index along y-direction
     4321       INTEGER(iwp) ::  k      !< grid index along z-direction
     4322       INTEGER(iwp) ::  k_mm   !< k-2 index in disretization, can be modified to avoid segmentation faults
     4323       INTEGER(iwp) ::  k_pp   !< k+2 index in disretization, can be modified to avoid segmentation faults
     4324       INTEGER(iwp) ::  k_ppp  !< k+3 index in disretization, can be modified to avoid segmentation faults
     4325       INTEGER(iwp) ::  tn = 0 !< number of OpenMP thread
    43234326       
    4324        REAL(wp)    ::  diss_d !<
    4325        REAL(wp)    ::  div    !<
    4326        REAL(wp)    ::  flux_d !<
    4327        REAL(wp)    ::  gu     !<
    4328        REAL(wp)    ::  gv     !<
    4329        REAL(wp)    ::  u_comp !<
    4330        REAL(wp)    ::  w_comp !<
     4327       REAL(wp)    ::  diss_d !< artificial dissipation term at grid box bottom
     4328       REAL(wp)    ::  div    !< diverence on v-grid
     4329       REAL(wp)    ::  flux_d !< artificial 6th-order flux at grid box bottom
     4330       REAL(wp)    ::  gu     !< Galilei-transformation velocity along x
     4331       REAL(wp)    ::  gv     !< Galilei-transformation velocity along y
     4332       REAL(wp)    ::  u_comp !< advection velocity along x
     4333       REAL(wp)    ::  w_comp !< advection velocity along z
    43314334       
    4332        REAL(wp), DIMENSION(nzb+1:nzt) ::  swap_diss_y_local_v !<
    4333        REAL(wp), DIMENSION(nzb+1:nzt) ::  swap_flux_y_local_v !<
     4335       REAL(wp), DIMENSION(nzb+1:nzt) ::  swap_diss_y_local_v !< discretized artificial dissipation at southward-side of the grid box
     4336       REAL(wp), DIMENSION(nzb+1:nzt) ::  swap_flux_y_local_v !< discretized 6th-order flux at southward-side of the grid box
    43344337       
    4335        REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_diss_x_local_v !<
    4336        REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_flux_x_local_v !<
     4338       REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_diss_x_local_v !< discretized artificial dissipation at leftward-side of the grid box
     4339       REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_flux_x_local_v !< discretized 6th-order flux at leftward-side of the grid box
    43374340       
    4338        REAL(wp), DIMENSION(nzb:nzt) ::  diss_n !<
    4339        REAL(wp), DIMENSION(nzb:nzt) ::  diss_r !<
    4340        REAL(wp), DIMENSION(nzb:nzt) ::  diss_t !<
    4341        REAL(wp), DIMENSION(nzb:nzt) ::  flux_n !<
    4342        REAL(wp), DIMENSION(nzb:nzt) ::  flux_r !<
    4343        REAL(wp), DIMENSION(nzb:nzt) ::  flux_t !<
    4344        REAL(wp), DIMENSION(nzb:nzt) ::  v_comp !<
     4341       REAL(wp), DIMENSION(nzb:nzt) ::  diss_n !< discretized artificial dissipation at northward-side of the grid box
     4342       REAL(wp), DIMENSION(nzb:nzt) ::  diss_r !< discretized artificial dissipation at rightward-side of the grid box
     4343       REAL(wp), DIMENSION(nzb:nzt) ::  diss_t !< discretized artificial dissipation at top of the grid box
     4344       REAL(wp), DIMENSION(nzb:nzt) ::  flux_n !< discretized 6th-order flux at northward-side of the grid box
     4345       REAL(wp), DIMENSION(nzb:nzt) ::  flux_r !< discretized 6th-order flux at rightward-side of the grid box
     4346       REAL(wp), DIMENSION(nzb:nzt) ::  flux_t !< discretized 6th-order flux at top of the grid box
     4347       REAL(wp), DIMENSION(nzb:nzt) ::  v_comp !< advection velocity along y
    43454348
    43464349       gu = 2.0_wp * u_gtrans
     
    48074810       IMPLICIT NONE
    48084811
    4809        INTEGER(iwp) ::  i      !<
    4810        INTEGER(iwp) ::  ibit27 !<
    4811        INTEGER(iwp) ::  ibit28 !<
    4812        INTEGER(iwp) ::  ibit29 !<
    4813        INTEGER(iwp) ::  ibit30 !<
    4814        INTEGER(iwp) ::  ibit31 !<
    4815        INTEGER(iwp) ::  ibit32 !<
    4816        INTEGER(iwp) ::  ibit33 !<
    4817        INTEGER(iwp) ::  ibit34 !<
    4818        INTEGER(iwp) ::  ibit35 !<
    4819        INTEGER(iwp) ::  j      !<
    4820        INTEGER(iwp) ::  k      !<
    4821        INTEGER(iwp) ::  k_mm   !<
    4822        INTEGER(iwp) ::  k_pp   !<
    4823        INTEGER(iwp) ::  k_ppp  !<
    4824        INTEGER(iwp) ::  tn = 0 !<
     4812       INTEGER(iwp) ::  i      !< grid index along x-direction
     4813       INTEGER(iwp) ::  ibit27 !< flag indicating 1st-order scheme along x-direction
     4814       INTEGER(iwp) ::  ibit28 !< flag indicating 3rd-order scheme along x-direction
     4815       INTEGER(iwp) ::  ibit29 !< flag indicating 5th-order scheme along x-direction
     4816       INTEGER(iwp) ::  ibit30 !< flag indicating 1st-order scheme along y-direction
     4817       INTEGER(iwp) ::  ibit31 !< flag indicating 3rd-order scheme along y-direction
     4818       INTEGER(iwp) ::  ibit32 !< flag indicating 5th-order scheme along y-direction
     4819       INTEGER(iwp) ::  ibit33 !< flag indicating 1st-order scheme along z-direction
     4820       INTEGER(iwp) ::  ibit34 !< flag indicating 3rd-order scheme along z-direction
     4821       INTEGER(iwp) ::  ibit35 !< flag indicating 5th-order scheme along z-direction
     4822       INTEGER(iwp) ::  j      !< grid index along y-direction
     4823       INTEGER(iwp) ::  k      !< grid index along z-direction
     4824       INTEGER(iwp) ::  k_mm   !< k-2 index in disretization, can be modified to avoid segmentation faults
     4825       INTEGER(iwp) ::  k_pp   !< k+2 index in disretization, can be modified to avoid segmentation faults
     4826       INTEGER(iwp) ::  k_ppp  !< k+3 index in disretization, can be modified to avoid segmentation faults
     4827       INTEGER(iwp) ::  tn = 0 !< number of OpenMP thread
    48254828       
    4826        REAL(wp)    ::  diss_d !<
    4827        REAL(wp)    ::  div    !<
    4828        REAL(wp)    ::  flux_d !<
    4829        REAL(wp)    ::  gu     !<
    4830        REAL(wp)    ::  gv     !<
    4831        REAL(wp)    ::  u_comp !<
    4832        REAL(wp)    ::  v_comp !<
    4833        REAL(wp)    ::  w_comp !<
     4829       REAL(wp)    ::  diss_d !< artificial dissipation term at grid box bottom
     4830       REAL(wp)    ::  div    !< divergence on w-grid
     4831       REAL(wp)    ::  flux_d !< 6th-order flux at grid box bottom
     4832       REAL(wp)    ::  gu     !< Galilei-transformation velocity along x
     4833       REAL(wp)    ::  gv     !< Galilei-transformation velocity along y
     4834       REAL(wp)    ::  u_comp !< advection velocity along x
     4835       REAL(wp)    ::  v_comp !< advection velocity along y
     4836       REAL(wp)    ::  w_comp !< advection velocity along z
    48344837       
    4835        REAL(wp), DIMENSION(nzb:nzt)    ::  diss_t !<
    4836        REAL(wp), DIMENSION(nzb:nzt)    ::  flux_t !<
     4838       REAL(wp), DIMENSION(nzb:nzt)    ::  diss_t !< discretized artificial dissipation at top of the grid box
     4839       REAL(wp), DIMENSION(nzb:nzt)    ::  flux_t !< discretized 6th-order flux at top of the grid box
    48374840       
    4838        REAL(wp), DIMENSION(nzb+1:nzt)  ::  diss_n !<
    4839        REAL(wp), DIMENSION(nzb+1:nzt)  ::  diss_r !<
    4840        REAL(wp), DIMENSION(nzb+1:nzt)  ::  flux_n !<
    4841        REAL(wp), DIMENSION(nzb+1:nzt)  ::  flux_r !<
    4842        REAL(wp), DIMENSION(nzb+1:nzt)  ::  swap_diss_y_local_w !<
    4843        REAL(wp), DIMENSION(nzb+1:nzt)  ::  swap_flux_y_local_w !<
     4841       REAL(wp), DIMENSION(nzb+1:nzt)  ::  diss_n !< discretized artificial dissipation at northward-side of the grid box
     4842       REAL(wp), DIMENSION(nzb+1:nzt)  ::  diss_r !< discretized artificial dissipation at rightward-side of the grid box
     4843       REAL(wp), DIMENSION(nzb+1:nzt)  ::  flux_n !< discretized 6th-order flux at northward-side of the grid box
     4844       REAL(wp), DIMENSION(nzb+1:nzt)  ::  flux_r !< discretized 6th-order flux at rightward-side of the grid box
     4845       REAL(wp), DIMENSION(nzb+1:nzt)  ::  swap_diss_y_local_w !< discretized artificial dissipation at southward-side of the grid box
     4846       REAL(wp), DIMENSION(nzb+1:nzt)  ::  swap_flux_y_local_w !< discretized 6th-order flux at southward-side of the grid box
    48444847       
    4845        REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_diss_x_local_w !<
    4846        REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_flux_x_local_w !<
     4848       REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_diss_x_local_w !< discretized artificial dissipation at leftward-side of the grid box
     4849       REAL(wp), DIMENSION(nzb+1:nzt,nys:nyn) ::  swap_flux_x_local_w !< discretized 6th-order flux at leftward-side of the grid box
    48474850 
    48484851       gu = 2.0_wp * u_gtrans
  • palm/trunk/SOURCE/data_output_spectra.f90

    r3421 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3421 2018-10-24 18:39:32Z gronemeier
    2730! Renamed output variables
    2831!
     
    126129    IMPLICIT NONE
    127130
    128     INTEGER(iwp) ::  m       !<
    129     INTEGER(iwp) ::  pr      !<
     131    INTEGER(iwp) ::  m       !< running index over spectra output
     132    INTEGER(iwp) ::  pr      !< index used to assign default quantities to data output
    130133   
    131134    CALL cpu_log( log_point(31), 'data_output_spectra', 'start' )
     
    267270    IMPLICIT NONE
    268271
    269     CHARACTER (LEN=1), INTENT(IN) ::  direction     !<
    270 
    271     INTEGER(iwp), INTENT(IN)      ::  nsp           !<
    272 
    273     INTEGER(iwp)                  ::  i             !<
    274     INTEGER(iwp)                  ::  k             !<
    275 
    276     REAL(wp)                      ::  frequency     !<
    277 
    278     REAL(wp), DIMENSION(nx/2)     ::  netcdf_data_x !<
    279     REAL(wp), DIMENSION(ny/2)     ::  netcdf_data_y !<
     272    CHARACTER (LEN=1), INTENT(IN) ::  direction     !< directio of spectra evaluation
     273
     274    INTEGER(iwp), INTENT(IN)      ::  nsp           !< number of spectrum
     275
     276    INTEGER(iwp)                  ::  i             !< running index in frequency space
     277    INTEGER(iwp)                  ::  k             !< running index over number of spectrum
     278
     279    REAL(wp)                      ::  frequency     !< wavenumber
     280
     281    REAL(wp), DIMENSION(nx/2)     ::  netcdf_data_x !< normalized wavenumber along x written into NetCDF file
     282    REAL(wp), DIMENSION(ny/2)     ::  netcdf_data_y !< normalized wavenumber along y written into NetCDF file
    280283
    281284
  • palm/trunk/SOURCE/diffusion_s.f90

    r2759 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 2759 2018-01-17 16:24:59Z suehring
    2730! Major bugfix, horizontal diffusion at vertical surfaces corrected.
    2831!
     
    186189       REAL(wp), DIMENSION(1:surf_def_h(2)%ns) ::  s_flux_t           !< flux at model top
    187190#if defined( __nopointer )
    188        REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  s  !<
     191       REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  s  !< treated scalar
    189192#else
    190        REAL(wp), DIMENSION(:,:,:), POINTER ::  s  !<
     193       REAL(wp), DIMENSION(:,:,:), POINTER ::  s  !< treated scalar
    191194#endif
    192195
     
    497500       REAL(wp), DIMENSION(1:surf_def_h(2)%ns) ::  s_flux_t           !< flux at model top
    498501#if defined( __nopointer )
    499        REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  s !<
     502       REAL(wp), DIMENSION(nzb:nzt+1,nysg:nyng,nxlg:nxrg) ::  s !< treated scalar
    500503#else
    501        REAL(wp), DIMENSION(:,:,:), POINTER ::  s  !<
     504       REAL(wp), DIMENSION(:,:,:), POINTER ::  s  !< treated scalar
    502505#endif
    503506
  • palm/trunk/SOURCE/diffusion_u.f90

    r3241 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3241 2018-09-12 15:02:00Z raasch
    2730! unused variables removed
    2831!
     
    154157       INTEGER(iwp) ::  l             !< running index of surface type, south- or north-facing wall
    155158       INTEGER(iwp) ::  m             !< running index surface elements
    156        INTEGER(iwp) ::  surf_e        !< End index of surface elements at (j,i)-gridpoint
    157        INTEGER(iwp) ::  surf_s        !< Start index of surface elements at (j,i)-gridpoint
     159       INTEGER(iwp) ::  surf_e        !< end index of surface elements at (j,i)-gridpoint
     160       INTEGER(iwp) ::  surf_s        !< start index of surface elements at (j,i)-gridpoint
    158161
    159162       REAL(wp)     ::  flag          !< flag to mask topography grid points
    160        REAL(wp)     ::  kmym          !<
    161        REAL(wp)     ::  kmyp          !<
    162        REAL(wp)     ::  kmzm          !<
    163        REAL(wp)     ::  kmzp          !<
     163       REAL(wp)     ::  kmym          !< diffusion coefficient on southward side of the u-gridbox - interpolated onto xu-yv grid
     164       REAL(wp)     ::  kmyp          !< diffusion coefficient on northward side of the u-gridbox - interpolated onto xu-yv grid
     165       REAL(wp)     ::  kmzm          !< diffusion coefficient on bottom of the gridbox - interpolated onto xu-zw grid
     166       REAL(wp)     ::  kmzp          !< diffusion coefficient on top of the gridbox - interpolated onto xu-zw grid
    164167       REAL(wp)     ::  mask_bottom   !< flag to mask vertical upward-facing surface       
    165168       REAL(wp)     ::  mask_north    !< flag to mask vertical surface north of the grid point
     
    395398
    396399       REAL(wp)     ::  flag          !< flag to mask topography grid points
    397        REAL(wp)     ::  kmym          !<
    398        REAL(wp)     ::  kmyp          !<
    399        REAL(wp)     ::  kmzm          !<
    400        REAL(wp)     ::  kmzp          !<
     400       REAL(wp)     ::  kmym          !< diffusion coefficient on southward side of the u-gridbox - interpolated onto xu-yv grid
     401       REAL(wp)     ::  kmyp          !<diffusion coefficient on northward side of the u-gridbox - interpolated onto xu-yv grid
     402       REAL(wp)     ::  kmzm          !< diffusion coefficient on bottom of the gridbox - interpolated onto xu-zw grid
     403       REAL(wp)     ::  kmzp          !< diffusion coefficient on top of the gridbox - interpolated onto xu-zw grid
    401404       REAL(wp)     ::  mask_bottom   !< flag to mask vertical upward-facing surface       
    402405       REAL(wp)     ::  mask_north    !< flag to mask vertical surface north of the grid point
  • palm/trunk/SOURCE/diffusion_v.f90

    r3241 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3241 2018-09-12 15:02:00Z raasch
    2730! unused variables removed
    2831!
     
    153156
    154157       REAL(wp)     ::  flag          !< flag to mask topography grid points
    155        REAL(wp)     ::  kmxm          !<
    156        REAL(wp)     ::  kmxp          !<
    157        REAL(wp)     ::  kmzm          !<
    158        REAL(wp)     ::  kmzp          !<
     158       REAL(wp)     ::  kmxm          !< diffusion coefficient on leftward side of the v-gridbox - interpolated onto xu-yv grid
     159       REAL(wp)     ::  kmxp          !< diffusion coefficient on rightward side of the v-gridbox - interpolated onto xu-yv grid
     160       REAL(wp)     ::  kmzm          !< diffusion coefficient on bottom of the gridbox - interpolated onto yv-zw grid
     161       REAL(wp)     ::  kmzp          !< diffusion coefficient on top of the gridbox - interpolated onto yv-zw grid
    159162       REAL(wp)     ::  mask_bottom   !< flag to mask vertical upward-facing surface 
    160163       REAL(wp)     ::  mask_east     !< flag to mask vertical surface south of the grid point
     
    384387
    385388       REAL(wp)     ::  flag          !< flag to mask topography grid points
    386        REAL(wp)     ::  kmxm          !<
    387        REAL(wp)     ::  kmxp          !<
    388        REAL(wp)     ::  kmzm          !<
    389        REAL(wp)     ::  kmzp          !<
     389       REAL(wp)     ::  kmxm          !< diffusion coefficient on leftward side of the v-gridbox - interpolated onto xu-yv grid
     390       REAL(wp)     ::  kmxp          !< diffusion coefficient on rightward side of the v-gridbox - interpolated onto xu-yv grid
     391       REAL(wp)     ::  kmzm          !< diffusion coefficient on bottom of the gridbox - interpolated onto xu-zw grid
     392       REAL(wp)     ::  kmzp          !< diffusion coefficient on top of the gridbox - interpolated onto xu-zw grid
    390393       REAL(wp)     ::  mask_bottom   !< flag to mask vertical upward-facing surface 
    391394       REAL(wp)     ::  mask_east     !< flag to mask vertical surface south of the grid point
  • palm/trunk/SOURCE/diffusion_w.f90

    r3241 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3241 2018-09-12 15:02:00Z raasch
    2730! unused variables removed
    2831!
     
    152155       
    153156       REAL(wp) ::  flag              !< flag to mask topography grid points
    154        REAL(wp) ::  kmxm              !<
    155        REAL(wp) ::  kmxp              !<
    156        REAL(wp) ::  kmym              !<
    157        REAL(wp) ::  kmyp              !<
     157       REAL(wp) ::  kmxm              !< diffusion coefficient on leftward side of the w-gridbox - interpolated onto xu-y grid
     158       REAL(wp) ::  kmxp              !<diffusion coefficient on rightward side of the w-gridbox - interpolated onto xu-y grid
     159       REAL(wp) ::  kmym              !< diffusion coefficient on southward side of the w-gridbox - interpolated onto x-yv grid
     160       REAL(wp) ::  kmyp              !< diffusion coefficient on northward side of the w-gridbox - interpolated onto x-yv grid
    158161       REAL(wp) ::  mask_west         !< flag to mask vertical wall west of the grid point
    159162       REAL(wp) ::  mask_east         !< flag to mask vertical wall east of the grid point
     
    180183!
    181184!--             Interpolate eddy diffusivities on staggered gridpoints
    182                 kmxp = 0.25_wp * ( km(k,j,i)   +   km(k,j,i+1) +               &
     185                kmxp = 0.25_wp * ( km(k,j,i)   + km(k,j,i+1)  +               &
    183186                                   km(k+1,j,i) + km(k+1,j,i+1) )
    184187                kmxm = 0.25_wp * ( km(k,j,i)   + km(k,j,i-1)   +               &
     
    328331       
    329332       REAL(wp) ::  flag              !< flag to mask topography grid points
    330        REAL(wp) ::  kmxm              !<
    331        REAL(wp) ::  kmxp              !<
    332        REAL(wp) ::  kmym              !<
    333        REAL(wp) ::  kmyp              !<
     333       REAL(wp) ::  kmxm              !< diffusion coefficient on leftward side of the w-gridbox - interpolated onto xu-y grid
     334       REAL(wp) ::  kmxp              !< diffusion coefficient on rightward side of the w-gridbox - interpolated onto xu-y grid
     335       REAL(wp) ::  kmym              !< diffusion coefficient on southward side of the w-gridbox - interpolated onto x-yv grid
     336       REAL(wp) ::  kmyp              !< diffusion coefficient on northward side of the w-gridbox - interpolated onto x-yv grid
    334337       REAL(wp) ::  mask_west         !< flag to mask vertical wall west of the grid point
    335338       REAL(wp) ::  mask_east         !< flag to mask vertical wall east of the grid point
  • palm/trunk/SOURCE/init_3d_model.f90

    r3525 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3525 2018-11-14 16:06:14Z kanani
    2730! Changes related to clean-up of biometeorology (dom_dwd_user)
    2831!
     
    685688    IMPLICIT NONE
    686689
    687     INTEGER(iwp) ::  i             !<
    688     INTEGER(iwp) ::  ind_array(1)  !<
    689     INTEGER(iwp) ::  j             !<
    690     INTEGER(iwp) ::  k             !<
     690    INTEGER(iwp) ::  i             !< grid index in x direction
     691    INTEGER(iwp) ::  ind_array(1)  !< dummy used to determine start index for external pressure forcing
     692    INTEGER(iwp) ::  j             !< grid index in y direction
     693    INTEGER(iwp) ::  k             !< grid index in z direction
    691694    INTEGER(iwp) ::  k_surf        !< surface level index
    692695    INTEGER(iwp) ::  m             !< index of surface element in surface data type
    693696    INTEGER(iwp) ::  sr            !< index of statistic region
    694697
    695     INTEGER(iwp), DIMENSION(:), ALLOCATABLE   ::  ngp_2dh_l  !<
    696 
    697     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  ngp_2dh_outer_l    !<
    698     INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  ngp_2dh_s_inner_l  !<
     698    INTEGER(iwp), DIMENSION(:), ALLOCATABLE   ::  ngp_2dh_l  !< toal number of horizontal grid points in statistical region on subdomain
     699
     700    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  ngp_2dh_outer_l    !< number of horizontal non-wall bounded grid points on subdomain
     701    INTEGER(iwp), DIMENSION(:,:), ALLOCATABLE ::  ngp_2dh_s_inner_l  !< number of horizontal non-topography grid points on subdomain
    699702
    700703    REAL(wp)     ::  t_surface !< air temperature at the surface
     
    705708    INTEGER(iwp) ::  l       !< loop variable
    706709    INTEGER(iwp) ::  nzt_l   !< index of top PE boundary for multigrid level
     710   
    707711    REAL(wp) ::  dx_l !< grid spacing along x on different multigrid level
    708712    REAL(wp) ::  dy_l !< grid spacing along y on different multigrid level
    709713
    710     REAL(wp), DIMENSION(1:3) ::  volume_flow_area_l     !<
    711     REAL(wp), DIMENSION(1:3) ::  volume_flow_initial_l  !<
    712 
    713     REAL(wp), DIMENSION(:), ALLOCATABLE ::  mean_surface_level_height_l    !<
    714     REAL(wp), DIMENSION(:), ALLOCATABLE ::  ngp_3d_inner_l    !<
    715     REAL(wp), DIMENSION(:), ALLOCATABLE ::  ngp_3d_inner_tmp  !<
    716 
    717     INTEGER(iwp) ::  nz_u_shift   !<
    718     INTEGER(iwp) ::  nz_v_shift   !<
    719     INTEGER(iwp) ::  nz_w_shift   !<
    720     INTEGER(iwp) ::  nz_s_shift   !<
    721     INTEGER(iwp) ::  nz_u_shift_l !<
    722     INTEGER(iwp) ::  nz_v_shift_l !<
    723     INTEGER(iwp) ::  nz_w_shift_l !<
    724     INTEGER(iwp) ::  nz_s_shift_l !<
     714    REAL(wp), DIMENSION(1:3) ::  volume_flow_area_l     !< area of lateral and top model domain surface on local subdomain
     715    REAL(wp), DIMENSION(1:3) ::  volume_flow_initial_l  !< initial volume flow into model domain
     716
     717    REAL(wp), DIMENSION(:), ALLOCATABLE ::  mean_surface_level_height_l !< mean surface level height on subdomain
     718    REAL(wp), DIMENSION(:), ALLOCATABLE ::  ngp_3d_inner_l    !< total number of non-topography grid points on subdomain
     719    REAL(wp), DIMENSION(:), ALLOCATABLE ::  ngp_3d_inner_tmp  !< total number of non-topography grid points
     720
     721    INTEGER(iwp) ::  nz_u_shift   !< topography-top index on u-grid, used to vertically shift initial profiles
     722    INTEGER(iwp) ::  nz_v_shift   !< topography-top index on v-grid, used to vertically shift initial profiles
     723    INTEGER(iwp) ::  nz_w_shift   !< topography-top index on w-grid, used to vertically shift initial profiles
     724    INTEGER(iwp) ::  nz_s_shift   !< topography-top index on scalar-grid, used to vertically shift initial profiles
     725    INTEGER(iwp) ::  nz_u_shift_l !< topography-top index on u-grid, used to vertically shift initial profiles
     726    INTEGER(iwp) ::  nz_v_shift_l !< topography-top index on v-grid, used to vertically shift initial profiles
     727    INTEGER(iwp) ::  nz_w_shift_l !< topography-top index on w-grid, used to vertically shift initial profiles
     728    INTEGER(iwp) ::  nz_s_shift_l !< topography-top index on scalar-grid, used to vertically shift initial profiles
    725729
    726730    CALL location_message( 'allocating arrays', .FALSE. )
  • palm/trunk/SOURCE/surface_layer_fluxes_mod.f90

    r3361 r3547  
    2626! -----------------
    2727! $Id$
     28! variables documented
     29!
     30! 3361 2018-10-16 20:39:37Z knoop
    2831! Modularization of all bulk cloud physics code components
    2932!
     
    337340       IMPLICIT NONE
    338341
    339        surf_vertical = .FALSE.
    340        downward      = .FALSE.
     342       surf_vertical = .FALSE. !< flag indicating vertically orientated surface elements
     343       downward      = .FALSE. !< flag indicating downward-facing surface elements
    341344!
    342345!--    Derive potential temperature and specific humidity at first grid level
     
    10361039       IMPLICIT NONE
    10371040
    1038        INTEGER(iwp) ::  i             !< running index x direction
    1039        INTEGER(iwp) ::  j             !< running index y direction
    1040        INTEGER(iwp) ::  k             !< running index z direction
    1041        INTEGER(iwp) ::  m             !< running index surface elements
    1042 
    1043        REAL(wp)     ::  u_i
    1044        REAL(wp)     ::  w_i
     1041       INTEGER(iwp) ::  i   !< running index x direction
     1042       INTEGER(iwp) ::  j   !< running index y direction
     1043       INTEGER(iwp) ::  k   !< running index z direction
     1044       INTEGER(iwp) ::  m   !< running index surface elements
     1045
     1046       REAL(wp)     ::  u_i !< u-component on xu-grid
     1047       REAL(wp)     ::  w_i !< w-component on xu-grid
    10451048
    10461049
     
    10711074       IMPLICIT NONE
    10721075
    1073        INTEGER(iwp) ::  i             !< running index x direction
    1074        INTEGER(iwp) ::  j             !< running index y direction
    1075        INTEGER(iwp) ::  k             !< running index z direction
    1076        INTEGER(iwp) ::  m             !< running index surface elements
    1077 
    1078        REAL(wp)     ::  v_i
    1079        REAL(wp)     ::  w_i
     1076       INTEGER(iwp) ::  i   !< running index x direction
     1077       INTEGER(iwp) ::  j   !< running index y direction
     1078       INTEGER(iwp) ::  k   !< running index z direction
     1079       INTEGER(iwp) ::  m   !< running index surface elements
     1080
     1081       REAL(wp)     ::  v_i !< v-component on yv-grid
     1082       REAL(wp)     ::  w_i !< w-component on yv-grid
    10801083
    10811084
     
    11051108       IMPLICIT NONE
    11061109
    1107        INTEGER(iwp) ::  i             !< running index x direction
    1108        INTEGER(iwp) ::  j             !< running index y direction
    1109        INTEGER(iwp) ::  k             !< running index z direction
    1110        INTEGER(iwp) ::  m             !< running index surface elements
    1111 
    1112        REAL(wp)     ::  u_i
    1113        REAL(wp)     ::  v_i
    1114        REAL(wp)     ::  w_i
     1110       INTEGER(iwp) ::  i   !< running index x direction
     1111       INTEGER(iwp) ::  j   !< running index y direction
     1112       INTEGER(iwp) ::  k   !< running index z direction
     1113       INTEGER(iwp) ::  m   !< running index surface elements
     1114
     1115       REAL(wp)     ::  u_i !< u-component on x-zw-grid
     1116       REAL(wp)     ::  v_i !< v-component on y-zw-grid
     1117       REAL(wp)     ::  w_i !< w-component on zw-grid
    11151118!
    11161119!--    North- (l=0) and south-facing (l=1) surfaces
     
    11571160       IMPLICIT NONE
    11581161
    1159        INTEGER(iwp) ::  i             !< running index x direction
    1160        INTEGER(iwp) ::  j             !< running index y direction
    1161        INTEGER(iwp) ::  k             !< running index z direction
    1162        INTEGER(iwp) ::  m             !< running index surface elements
    1163 
    1164        REAL(wp)     ::  u_i
    1165        REAL(wp)     ::  v_i
    1166        REAL(wp)     ::  w_i
     1162       INTEGER(iwp) ::  i   !< running index x direction
     1163       INTEGER(iwp) ::  j   !< running index y direction
     1164       INTEGER(iwp) ::  k   !< running index z direction
     1165       INTEGER(iwp) ::  m   !< running index surface elements
     1166
     1167       REAL(wp)     ::  u_i !< u-component on scalar grid
     1168       REAL(wp)     ::  v_i !< v-component on scalar grid
     1169       REAL(wp)     ::  w_i !< w-component on scalar grid
    11671170
    11681171!
  • palm/trunk/SOURCE/surface_mod.f90

    r3467 r3547  
    2626! -----------------
    2727! $Id$
     28! variables documented
     29!
     30! 3467 2018-10-30 19:05:21Z suehring
    2831! Bugfix in surface_wrd_local for cloud and rain water flux
    2932!
     
    473476       REAL(wp), DIMENSION(:), ALLOCATABLE ::  surfinlw            !< longwave radiation falling to local surface including radiation from reflections
    474477       REAL(wp), DIMENSION(:), ALLOCATABLE ::  surfoutlw           !< total longwave radiation outgoing from nonvirtual surfaces surfaces after all reflection
     478       
     479       REAL(wp), DIMENSION(:), ALLOCATABLE ::  n_vg_green          !< vangenuchten parameters
     480       REAL(wp), DIMENSION(:), ALLOCATABLE ::  alpha_vg_green      !< vangenuchten parameters
     481       REAL(wp), DIMENSION(:), ALLOCATABLE ::  l_vg_green          !< vangenuchten parameters
    475482
    476483
     
    497504       REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  tt_green_m        !< t_green prognostic array
    498505       REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  zw_green          !< green layer depths (m)
    499        REAL(wp), DIMENSION(:), ALLOCATABLE ::  n_vg_green        !< vangenuchten parameters
    500        REAL(wp), DIMENSION(:), ALLOCATABLE ::  alpha_vg_green    !< vangenuchten parameters
    501        REAL(wp), DIMENSION(:), ALLOCATABLE ::  l_vg_green    !< vangenuchten parameters
    502        REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  gamma_w_green_sat    !< hydraulic conductivity
     506
     507       REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  gamma_w_green_sat !< hydraulic conductivity
    503508       REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  lambda_w_green
    504        REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  gamma_w_green    !< hydraulic conductivity
    505        REAL(wp), DIMENSION(:,:), ALLOCATABLE :: tswc_h_m
     509       REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  gamma_w_green     !< hydraulic conductivity
     510       REAL(wp), DIMENSION(:,:), ALLOCATABLE ::  tswc_h_m
    506511
    507512
     
    625630       IMPLICIT NONE
    626631
    627        INTEGER(iwp) ::  i         !<
    628        INTEGER(iwp) ::  j         !<
    629        INTEGER(iwp) ::  k         !<
    630 
    631        INTEGER(iwp), DIMENSION(0:1) ::  num_h         !<
    632        INTEGER(iwp), DIMENSION(0:1) ::  num_h_kji     !<
    633        INTEGER(iwp), DIMENSION(0:1) ::  start_index_h !<
     632       INTEGER(iwp) ::  i         !< loop index along x-direction
     633       INTEGER(iwp) ::  j         !< loop index along y-direction
     634       INTEGER(iwp) ::  k         !< loop index along y-direction
     635
     636       INTEGER(iwp), DIMENSION(0:1) ::  num_h         !< number of surfaces
     637       INTEGER(iwp), DIMENSION(0:1) ::  num_h_kji     !< number of surfaces for (j,i)-grid point
     638       INTEGER(iwp), DIMENSION(0:1) ::  start_index_h !< local start index
    634639
    635640!
  • palm/trunk/SOURCE/virtual_flight_mod.f90

    r3421 r3547  
    2525! -----------------
    2626! $Id$
     27! variables documented
     28!
     29! 3421 2018-10-24 18:39:32Z gronemeier
    2730! Renamed output variables
    2831!
     
    963966       IMPLICIT NONE
    964967       
    965        LOGICAL, INTENT(OUT)  ::  found
     968       LOGICAL, INTENT(OUT)  ::  found !< flag indicating if a variable string is found
    966969
    967970
Note: See TracChangeset for help on using the changeset viewer.