Ignore:
Timestamp:
Jun 9, 2016 4:25:25 PM (5 years ago)
Author:
suehring
Message:

several bugfixes in particle model and serial mode

File:
1 edited

Legend:

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

    r1823 r1929  
    1919! Current revisions:
    2020! ------------------
    21 !
     21! sgs_wfu_par, sgs_wfv_par and sgs_wfw_par are replaced by sgs_wf_par
    2222!
    2323! Former revisions:
     
    7272
    7373    USE particle_attributes,                                                   &
    74         ONLY:  sgs_wfu_part, sgs_wfv_part, sgs_wfw_part
     74        ONLY:  sgs_wf_part
    7575
    7676    USE pegrid
     
    9595                  k  > nzb_s_inner(j,i+1) )                                    &
    9696             THEN
    97                 de_dx(k,j,i) = 2.0_wp * sgs_wfu_part *                         &
     97                de_dx(k,j,i) = 2.0_wp * sgs_wf_part *                          &
    9898                               ( e(k,j,i+1) - e(k,j,i) ) * ddx
    9999             ELSEIF ( k  > nzb_s_inner(j,i-1)  .AND.  k > nzb_s_inner(j,i)     &
    100100                      .AND.  k <= nzb_s_inner(j,i+1) )                         &
    101101             THEN
    102                 de_dx(k,j,i) = 2.0_wp * sgs_wfu_part *                         &
     102                de_dx(k,j,i) = 2.0_wp * sgs_wf_part *                          &
    103103                               ( e(k,j,i) - e(k,j,i-1) ) * ddx
    104104             ELSEIF ( k < nzb_s_inner(j,i)  .AND.  k < nzb_s_inner(j,i+1) )    &
     
    109109                de_dx(k,j,i) = 0.0_wp
    110110             ELSE
    111                 de_dx(k,j,i) = sgs_wfu_part * ( e(k,j,i+1) - e(k,j,i-1) ) * ddx
     111                de_dx(k,j,i) = sgs_wf_part * ( e(k,j,i+1) - e(k,j,i-1) ) * ddx
    112112             ENDIF
    113113
     
    115115                  k  > nzb_s_inner(j+1,i) )                                    &
    116116             THEN
    117                 de_dy(k,j,i) = 2.0_wp * sgs_wfv_part *                         &
     117                de_dy(k,j,i) = 2.0_wp * sgs_wf_part *                          &
    118118                               ( e(k,j+1,i) - e(k,j,i) ) * ddy
    119119             ELSEIF ( k  > nzb_s_inner(j-1,i)  .AND.  k  > nzb_s_inner(j,i)    &
    120120                      .AND.  k <= nzb_s_inner(j+1,i) )                         &
    121121             THEN
    122                 de_dy(k,j,i) = 2.0_wp * sgs_wfv_part *                         &
     122                de_dy(k,j,i) = 2.0_wp * sgs_wf_part *                          &
    123123                               ( e(k,j,i) - e(k,j-1,i) ) * ddy
    124124             ELSEIF ( k < nzb_s_inner(j,i)  .AND.  k < nzb_s_inner(j+1,i) )    &
     
    129129                de_dy(k,j,i) = 0.0_wp
    130130             ELSE
    131                 de_dy(k,j,i) = sgs_wfv_part * ( e(k,j+1,i) - e(k,j-1,i) ) * ddy
     131                de_dy(k,j,i) = sgs_wf_part * ( e(k,j+1,i) - e(k,j-1,i) ) * ddy
    132132             ENDIF
    133133
     
    142142
    143143          DO  k = nzb_s_inner(j,i)+2, nzt-1
    144              de_dz(k,j,i)  = 2.0_wp * sgs_wfw_part *                           &
     144             de_dz(k,j,i)  = 2.0_wp * sgs_wf_part *                            &
    145145                             ( e(k+1,j,i) - e(k-1,j,i) ) / ( zu(k+1)-zu(k-1) )
    146146          ENDDO
     
    148148          k = nzb_s_inner(j,i)
    149149          de_dz(nzb:k,j,i) = 0.0_wp
    150           de_dz(k+1,j,i)   = 2.0_wp * sgs_wfw_part *                           &
     150          de_dz(k+1,j,i)   = 2.0_wp * sgs_wf_part *                            &
    151151                             ( e(k+2,j,i) - e(k+1,j,i) ) / ( zu(k+2) - zu(k+1) )
    152152          de_dz(nzt,j,i)   = 0.0_wp
Note: See TracChangeset for help on using the changeset viewer.