Ignore:
Timestamp:
Apr 11, 2014 5:15:14 PM (7 years ago)
Author:
hoffmann
Message:

new Lagrangian particle structure integrated

File:
1 edited

Legend:

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

    r1321 r1359  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! New particle structure integrated.
     23! Kind definition added to all floating point numbers.
    2324!
    2425! Former revisions:
     
    8485                  k  > nzb_s_inner(j,i+1) )                                    &
    8586             THEN
    86                 de_dx(k,j,i) = 2.0 * sgs_wfu_part * ( e(k,j,i+1) - e(k,j,i) )  &
    87                                * ddx
     87                de_dx(k,j,i) = 2.0_wp * sgs_wfu_part *                         &
     88                               ( e(k,j,i+1) - e(k,j,i) ) * ddx
    8889             ELSEIF ( k  > nzb_s_inner(j,i-1)  .AND.  k > nzb_s_inner(j,i)     &
    8990                      .AND.  k <= nzb_s_inner(j,i+1) )                         &
    9091             THEN
    91                 de_dx(k,j,i) = 2.0 * sgs_wfu_part * ( e(k,j,i) - e(k,j,i-1) )  &
    92                                * ddx
     92                de_dx(k,j,i) = 2.0_wp * sgs_wfu_part *                         &
     93                               ( e(k,j,i) - e(k,j,i-1) ) * ddx
    9394             ELSEIF ( k < nzb_s_inner(j,i)  .AND.  k < nzb_s_inner(j,i+1) )    &
    9495             THEN
    95                 de_dx(k,j,i) = 0.0
     96                de_dx(k,j,i) = 0.0_wp
    9697             ELSEIF ( k < nzb_s_inner(j,i-1)  .AND.  k < nzb_s_inner(j,i) )    &
    9798             THEN
    98                 de_dx(k,j,i) = 0.0
     99                de_dx(k,j,i) = 0.0_wp
    99100             ELSE
    100101                de_dx(k,j,i) = sgs_wfu_part * ( e(k,j,i+1) - e(k,j,i-1) ) * ddx
     
    104105                  k  > nzb_s_inner(j+1,i) )                                    &
    105106             THEN
    106                 de_dy(k,j,i) = 2.0 * sgs_wfv_part * ( e(k,j+1,i) - e(k,j,i) )  &
    107                                * ddy
     107                de_dy(k,j,i) = 2.0_wp * sgs_wfv_part *                         &
     108                               ( e(k,j+1,i) - e(k,j,i) ) * ddy
    108109             ELSEIF ( k  > nzb_s_inner(j-1,i)  .AND.  k  > nzb_s_inner(j,i)    &
    109110                      .AND.  k <= nzb_s_inner(j+1,i) )                         &
    110111             THEN
    111                 de_dy(k,j,i) = 2.0 * sgs_wfv_part * ( e(k,j,i) - e(k,j-1,i) )  &
    112                                * ddy
     112                de_dy(k,j,i) = 2.0_wp * sgs_wfv_part *                         &
     113                               ( e(k,j,i) - e(k,j-1,i) ) * ddy
    113114             ELSEIF ( k < nzb_s_inner(j,i)  .AND.  k < nzb_s_inner(j+1,i) )    &
    114115             THEN
    115                 de_dy(k,j,i) = 0.0
     116                de_dy(k,j,i) = 0.0_wp
    116117             ELSEIF ( k < nzb_s_inner(j-1,i)  .AND.  k < nzb_s_inner(j,i) )    &
    117118             THEN
    118                 de_dy(k,j,i) = 0.0
     119                de_dy(k,j,i) = 0.0_wp
    119120             ELSE
    120121                de_dy(k,j,i) = sgs_wfv_part * ( e(k,j+1,i) - e(k,j-1,i) ) * ddy
     
    131132
    132133          DO  k = nzb_s_inner(j,i)+2, nzt-1
    133              de_dz(k,j,i)  = 2.0 * sgs_wfw_part * &
     134             de_dz(k,j,i)  = 2.0_wp * sgs_wfw_part *                          &
    134135                             ( e(k+1,j,i) - e(k-1,j,i) ) / ( zu(k+1)-zu(k-1) )
    135136          ENDDO
    136137
    137138          k = nzb_s_inner(j,i)
    138           de_dz(nzb:k,j,i) = 0.0
    139           de_dz(k+1,j,i)   = 2.0 * sgs_wfw_part * ( e(k+2,j,i) - e(k+1,j,i) ) &
    140                                                / ( zu(k+2) - zu(k+1) )
    141           de_dz(nzt,j,i)   = 0.0
    142           de_dz(nzt+1,j,i) = 0.0
     139          de_dz(nzb:k,j,i) = 0.0_wp
     140          de_dz(k+1,j,i)   = 2.0_wp * sgs_wfw_part *                          &
     141                             ( e(k+2,j,i) - e(k+1,j,i) ) / ( zu(k+2) - zu(k+1) )
     142          de_dz(nzt,j,i)   = 0.0_wp
     143          de_dz(nzt+1,j,i) = 0.0_wp
    143144       ENDDO
    144145    ENDDO
     
    162163!--    First calculate horizontally averaged profiles of the horizontal
    163164!--    velocities.
    164        sums_l(:,1,0) = 0.0
    165        sums_l(:,2,0) = 0.0
     165       sums_l(:,1,0) = 0.0_wp
     166       sums_l(:,2,0) = 0.0_wp
    166167
    167168       DO  i = nxl, nxr
     
    197198!--    Now calculate the profiles of SGS TKE and the resolved-scale
    198199!--    velocity variances
    199        sums_l(:,8,0)  = 0.0
    200        sums_l(:,30,0) = 0.0
    201        sums_l(:,31,0) = 0.0
    202        sums_l(:,32,0) = 0.0
     200       sums_l(:,8,0)  = 0.0_wp
     201       sums_l(:,30,0) = 0.0_wp
     202       sums_l(:,31,0) = 0.0_wp
     203       sums_l(:,32,0) = 0.0_wp
    203204       DO  i = nxl, nxr
    204205          DO  j = nys, nyn
Note: See TracChangeset for help on using the changeset viewer.