Ignore:
Timestamp:
May 7, 2014 11:15:00 AM (10 years ago)
Author:
heinze
Message:

adjustments in case of nudging + bugfix for KIND in CMPLX function

File:
1 edited

Legend:

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

    r1381 r1398  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! Rayleigh-damping for horizontal velocity components changed: instead of damping
     23! against ug and vg, damping against u_init and v_init is used to allow for a
     24! homogenized treatment in case of nudging
    2325!
    2426! Former revisions:
     
    149151               ref_state, rho, sa, sa_init, sa_p, saswsb, saswst, shf, tend,   &
    150152               te_m, tnr_m, tpt_m, tq_m, tqr_m, tsa_m, tswst, tu_m, tv_m,      &
    151                tw_m, u, ug, u_p, v, vg, vpt, v_p, w, w_p
     153               tw_m, u, ug, u_init, u_p, v, vg, vpt, v_init, v_p, w, w_p
    152154       
    153155    USE control_parameters,                                                    &
     
    383385                u_p(k,j,i) = u(k,j,i) + dt_3d * ( tsc(2) * tend(k,j,i) +       &
    384386                                                  tsc(3) * tu_m(k,j,i) )       &
    385                                       - tsc(5) * rdf(k) * ( u(k,j,i) - ug(k) )
     387                                      - tsc(5) * rdf(k) * ( u(k,j,i) - u_init(k) )
    386388             ENDDO
    387389
     
    442444                v_p(k,j,i) = v(k,j,i) + dt_3d * ( tsc(2) * tend(k,j,i) +       &
    443445                                                  tsc(3) * tv_m(k,j,i) )       &
    444                                       - tsc(5) * rdf(k) * ( v(k,j,i) - vg(k) )
     446                                      - tsc(5) * rdf(k) * ( v(k,j,i) - v_init(k) )
    445447             ENDDO
    446448
     
    975977             u_p(k,j,i) = u(k,j,i) + dt_3d * ( tsc(2) * tend(k,j,i) +          &
    976978                                               tsc(3) * tu_m(k,j,i) )          &
    977                                    - tsc(5) * rdf(k) * ( u(k,j,i) - ug(k) )
     979                                   - tsc(5) * rdf(k) * ( u(k,j,i) - u_init(k) )
    978980          ENDDO
    979981       ENDDO
     
    10511053             v_p(k,j,i) = v(k,j,i) + dt_3d * ( tsc(2) * tend(k,j,i) +          &
    10521054                                               tsc(3) * tv_m(k,j,i) )          &
    1053                                    - tsc(5) * rdf(k) * ( v(k,j,i) - vg(k) )
     1055                                   - tsc(5) * rdf(k) * ( v(k,j,i) - v_init(k) )
    10541056          ENDDO
    10551057       ENDDO
     
    18111813!
    18121814!-- Prognostic equation for u-velocity component
    1813     !$acc kernels present( nzb_u_inner, rdf, tend, tu_m, u, ug, u_p )
     1815    !$acc kernels present( nzb_u_inner, rdf, tend, tu_m, u, u_init, u_p )
    18141816    !$acc loop independent
    18151817    DO  i = i_left, i_right
     
    18211823                u_p(k,j,i) = u(k,j,i) + dt_3d * ( tsc(2) * tend(k,j,i) +       &
    18221824                                                  tsc(3) * tu_m(k,j,i) )       &
    1823                                       - tsc(5) * rdf(k) * ( u(k,j,i) - ug(k) )
     1825                                      - tsc(5) * rdf(k) * ( u(k,j,i) - u_init(k) )
    18241826!
    18251827!--             Tendencies for the next Runge-Kutta step
     
    18781880!
    18791881!-- Prognostic equation for v-velocity component
    1880     !$acc kernels present( nzb_v_inner, rdf, tend, tv_m, v, vg, v_p )
     1882    !$acc kernels present( nzb_v_inner, rdf, tend, tv_m, v, v_init, v_p )
    18811883    !$acc loop independent
    18821884    DO  i = i_left, i_right
     
    18881890                v_p(k,j,i) = v(k,j,i) + dt_3d * ( tsc(2) * tend(k,j,i) +       &
    18891891                                                  tsc(3) * tv_m(k,j,i) )       &
    1890                                       - tsc(5) * rdf(k) * ( v(k,j,i) - vg(k) )
     1892                                      - tsc(5) * rdf(k) * ( v(k,j,i) - v_init(k) )
    18911893!
    18921894!--             Tendencies for the next Runge-Kutta step
Note: See TracChangeset for help on using the changeset viewer.