Changeset 1575 for palm/trunk/SOURCE/lpm_init.f90

Timestamp:

Mar 27, 2015 9:56:27 AM (9 years ago)

Author:

raasch

Message:

optimized multigrid method installed, new parameter seed_follows_topography for particle release, small adjustment in subjob for HLRN

File:

• palm/trunk/SOURCE/lpm_init.f90 (modified) (6 diffs)

palm/trunk/SOURCE/lpm_init.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! initial vertical particle position is allowed to follow the topography
 ! 
 ! Former revisions:
@@ -101 +101 @@
 
     USE indices,                                                               &
-        ONLY:  nx, nxl, nxlg, nxrg, nxr, ny, nyn, nys, nyng, nysg, nz, nzb, nzt
+        ONLY:  nx, nxl, nxlg, nxrg, nxr, ny, nyn, nys, nyng, nysg, nz, nzb,    &
+               nzb_w_inner, nzt
 
     USE kinds
@@ -126 +127 @@
                 prt_count, psb, psl, psn, psr, pss, pst,                       &
                 radius, random_start_position, read_particles_from_restartfile,&
-                skip_particles_for_tail, sort_count,                           &
+                seed_follows_topography, skip_particles_for_tail, sort_count,  &
                 tail_mask, total_number_of_particles, total_number_of_tails,   &
                 use_particle_tails, use_sgs_for_particles,                     &
@@ -649 +650 @@
                    pos_x = psl(i)
 
-                   DO WHILE ( pos_x <= psr(i) )
+            xloop: DO WHILE ( pos_x <= psr(i) )
 
                       IF ( pos_x >= ( nxl - 0.5_wp ) * dx  .AND.  &
@@ -738 +739 @@
                                tmp_particle%tail_id    = 0
                             ENDIF
+!
+!--                         Determine the grid indices of the particle position
                             ip = ( tmp_particle%x + 0.5_wp * dx ) * ddx
                             jp = ( tmp_particle%y + 0.5_wp * dy ) * ddy
                             kp = tmp_particle%z / dz + 1 + offset_ocean_nzt_m1
+
+                            IF ( seed_follows_topography )  THEN
+!
+!--                            Particle height is given relative to topography
+                               kp = kp + nzb_w_inner(jp,ip)
+                               tmp_particle%z = tmp_particle%z + zw(kp)
+                               IF ( kp > nzt )  THEN
+                                  pos_x = pos_x + pdx(i)
+                                  CYCLE xloop
+                               ENDIF
+                            ENDIF
 
                             local_count(kp,jp,ip) = local_count(kp,jp,ip) + 1
@@ -758 +772 @@
                       pos_x = pos_x + pdx(i)
 
-                   ENDDO
+                   ENDDO xloop
 
                 ENDIF