Changeset 4170 for palm/trunk/SOURCE/time_integration.f90

Timestamp:

Aug 19, 2019 5:12:31 PM (5 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:

• palm/trunk/SOURCE/time_integration.f90 (modified) (5 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! copy diss, diss_p, tdiss_m to GPU
+! 
+! 4144 2019-08-06 09:11:47Z raasch
 ! relational operators .EQ., .NE., etc. replaced by ==, /=, etc.
 ! 
@@ -718 +721 @@
         ONLY:  d, dd2zu, ddzu, ddzw, drho_air, drho_air_zw, dzw, heatflux_output_conversion, kh,   &
                km, momentumflux_output_conversion, p, ptdf_x, ptdf_y, rdf, rdf_sc, rho_air,        &
-               rho_air_zw, te_m, tpt_m, tu_m, tv_m, tw_m, ug, u_init, u_stokes_zu, vg, v_init,     &
-               v_stokes_zu, zu
+               rho_air_zw, tdiss_m, te_m, tpt_m, tu_m, tv_m, tw_m, ug, u_init, u_stokes_zu, vg,    &
+               v_init, v_stokes_zu, zu
 
     USE control_parameters,                                                                        &
@@ -761 +764 @@
 !$ACC DATA & 
 !$ACC COPY(d(nzb+1:nzt,nys:nyn,nxl:nxr)) &
+!$ACC COPY(diss(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
 !$ACC COPY(e(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
 !$ACC COPY(u(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
@@ -771 +775 @@
 
 !$ACC DATA &
+!$ACC COPY(diss_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
 !$ACC COPY(e_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
 !$ACC COPY(u_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
@@ -777 +782 @@
 !$ACC COPY(pt_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
 !$ACC COPY(tend(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
+!$ACC COPY(tdiss_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
 !$ACC COPY(te_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &
 !$ACC COPY(tu_m(nzb:nzt+1,nysg:nyng,nxlg:nxrg)) &