Ignore:
Timestamp:
Aug 19, 2019 5:12:31 PM (2 years ago)
Author:
gronemeier
Message:

changes in turbulence_closure_mod:

  • add performance optimizations according to K. Ketelsen to diffusion_e and tcm_diffusivities_default
  • bugfix in calculating l_wall for vertical walls
  • bugfix in using l_wall in initialization (consider wall_adjustment_factor)
  • always initialize diss and save the dissipation to that array

related changes in time_integration:

  • copy diss, diss_p, tdiss_m to GPU
File:
1 edited

Legend:

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

    r4144 r4170  
    2525! -----------------
    2626! $Id$
     27! copy diss, diss_p, tdiss_m to GPU
     28!
     29! 4144 2019-08-06 09:11:47Z raasch
    2730! relational operators .EQ., .NE., etc. replaced by ==, /=, etc.
    2831!
     
    718721        ONLY:  d, dd2zu, ddzu, ddzw, drho_air, drho_air_zw, dzw, heatflux_output_conversion, kh,   &
    719722               km, momentumflux_output_conversion, p, ptdf_x, ptdf_y, rdf, rdf_sc, rho_air,        &
    720                rho_air_zw, te_m, tpt_m, tu_m, tv_m, tw_m, ug, u_init, u_stokes_zu, vg, v_init,     &
    721                v_stokes_zu, zu
     723               rho_air_zw, tdiss_m, te_m, tpt_m, tu_m, tv_m, tw_m, ug, u_init, u_stokes_zu, vg,    &
     724               v_init, v_stokes_zu, zu
    722725
    723726    USE control_parameters,                                                                        &
     
    761764!$ACC DATA &
    762765!$ACC COPY(d(nzb+1:nzt,nys:nyn,nxl:nxr)) &
     766!$ACC COPY(diss(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
    763767!$ACC COPY(e(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
    764768!$ACC COPY(u(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
     
    771775
    772776!$ACC DATA &
     777!$ACC COPY(diss_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
    773778!$ACC COPY(e_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
    774779!$ACC COPY(u_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
     
    777782!$ACC COPY(pt_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
    778783!$ACC COPY(tend(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
     784!$ACC COPY(tdiss_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
    779785!$ACC COPY(te_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
    780786!$ACC COPY(tu_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
Note: See TracChangeset for help on using the changeset viewer.