Changeset 2122 for palm/trunk/SOURCE/lpm_set_attributes.f90

Timestamp:

Jan 18, 2017 12:22:54 PM (8 years ago)

Author:

hoffmann

Message:

introduction of a particle ID, improvement of equilibrium radius calculation, and reformatting

File:

• palm/trunk/SOURCE/lpm_set_attributes.f90 (modified) (2 diffs)

palm/trunk/SOURCE/lpm_set_attributes.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! DVRP routine removed
 ! 
 ! Former revisions:
@@ -397 +397 @@
     ENDIF
 
-!
-!-- Set particle size for dvrp graphics
-    IF ( particle_dvrpsize == 'absw' )  THEN
-
-       DO  ip = nxl, nxr
-          DO  jp = nys, nyn
-             DO  kp = nzb+1, nzt
-
-                number_of_particles = prt_count(kp,jp,ip)
-                particles => grid_particles(kp,jp,ip)%particles(1:number_of_particles)
-                IF ( number_of_particles <= 0 )  CYCLE
-                start_index = grid_particles(kp,jp,ip)%start_index
-                end_index   = grid_particles(kp,jp,ip)%end_index
-
-                ALLOCATE( xv(1:number_of_particles), &
-                          yv(1:number_of_particles) )
-
-                xv = particles(1:number_of_particles)%x
-                yv = particles(1:number_of_particles)%y
-                zv = particles(1:number_of_particles)%z
-
-                DO  nb = 0,7
-
-                   i = ip + block_offset(nb)%i_off
-                   j = jp + block_offset(nb)%j_off
-                   k = kp-1
-
-                   DO  n = start_index(nb), end_index(nb)
-!
-!--                   Interpolate w-component to the current particle position
-                      x  = xv(n) - i * dx
-                      y  = yv(n) - j * dy
-                      aa = x**2          + y**2
-                      bb = ( dx - x )**2 + y**2
-                      cc = x**2          + ( dy - y )**2
-                      dd = ( dx - x )**2 + ( dy - y )**2
-                      gg = aa + bb + cc + dd
-
-                      w_int_l = ( ( gg - aa ) * w(k,j,i)   + ( gg - bb ) *     &
-                                  w(k,j,i+1) + ( gg - cc ) * w(k,j+1,i) +      &
-                                  ( gg - dd ) * w(k,j+1,i+1)                   &
-                                ) / ( 3.0_wp * gg )
-
-                      IF ( k+1 == nzt+1 )  THEN
-                         w_int = w_int_l
-                      ELSE
-                         w_int_u = ( ( gg - aa ) * w(k+1,j,i)   + ( gg - bb ) *  &
-                                     w(k+1,j,i+1) + ( gg - cc ) * w(k+1,j+1,i) + &
-                                     ( gg - dd ) * w(k+1,j+1,i+1)                &
-                                   ) / ( 3.0_wp * gg )
-                         w_int = w_int_l + ( zv(n) - zw(k) ) / dz *     &
-                                           ( w_int_u - w_int_l )
-                      ENDIF
-
-!
-!--                   Limit values by the given interval and normalize to
-!--                   interval [0,1]
-                      w_int = ABS( w_int )
-                      w_int = MIN( w_int, dvrpsize_interval(2) )
-                      w_int = MAX( w_int, dvrpsize_interval(1) )
-
-                      w_int = ( w_int - dvrpsize_interval(1) ) / &
-                              ( dvrpsize_interval(2) - dvrpsize_interval(1) )
-
-                      particles(n)%dvrp_psize = ( 0.25_wp + w_int * 0.6_wp ) * &
-                                                dx
-
-                   ENDDO
-                ENDDO
-
-                DEALLOCATE( xv, yv, zv )
-
-             ENDDO
-          ENDDO
-       ENDDO
-
-    ENDIF
-
     CALL cpu_log( log_point_s(49), 'lpm_set_attributes', 'stop' )