Changeset 1212 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

Aug 15, 2013 8:46:27 AM (11 years ago)

Author:

raasch

Message:

tridia-solver moved to seperate module; the tridiagonal matrix coefficients of array tri are calculated only once at the beginning

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (4 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -23 +23 @@
 ! Current revisions:
 ! ------------------
-!
+! array tri is allocated and included in data copy statement
 !
 ! Former revisions:
@@ -364 +364 @@
 !-- Array for storing constant coeffficients of the tridiagonal solver
     IF ( psolver == 'poisfft' )  THEN
+       ALLOCATE( tri(nxl_z:nxr_z,nys_z:nyn_z,0:nz-1,2) )
        ALLOCATE( tric(nxl_z:nxr_z,nys_z:nyn_z,0:nz-1) )
     ENDIF
@@ -518 +519 @@
     IF ( use_sgs_for_particles  .OR.  wang_kernel  .OR.  turbulence  .OR.  &
          num_acc_per_node > 0 )  THEN
-       print*, '*** allocating diss'
        ALLOCATE( diss(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     ENDIF
@@ -1496 +1496 @@
        CALL disturb_field( nzb_v_inner, tend, v )
        n_sor = nsor_ini
-       !$acc data copy( d, ddzu, ddzw, nzb_s_inner, nzb_u_inner, nzb_v_inner, nzb_w_inner, p, tric, u, v, w, weight_pres, weight_substep, tend )
+       !$acc data copy( d, ddzu, ddzw, nzb_s_inner, nzb_u_inner, nzb_v_inner, nzb_w_inner, p, tri, tric, u, v, w, weight_pres, weight_substep, tend )
        CALL pres
        !$acc end data