Changeset 264 for palm/trunk/SOURCE/init_dvrp.f90

Timestamp:

Mar 21, 2009 8:14:44 AM (16 years ago)

Author:

raasch

Message:

new dvrp features added

File:

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

palm/trunk/SOURCE/init_dvrp.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
+! Definition of a colortable to be used for particles.
 ! Output names are changed: surface=groundplate, buildings=topography
 ! Output of messages replaced by message handling routine.
@@ -112 +113 @@
        dvrp_overlap       = .FALSE.
        dvrp_total_overlap = .FALSE.
-       nxl_dvrp = nx_dvrp_l
-       nxr_dvrp = nxl_dvrp + 4
-       nys_dvrp = ny_dvrp_s
-       nyn_dvrp = nys_dvrp + 4
+!       nxl_dvrp = nx_dvrp_l
+!       nxr_dvrp = nxl_dvrp + 4
+!       nys_dvrp = ny_dvrp_s
+!       nyn_dvrp = nys_dvrp + 4
+       nxl_dvrp = nxl
+       nxr_dvrp = MIN( nxl+4, nxr )
+       nys_dvrp = nys
+       nyn_dvrp = MIN( nys+4, nyn )
 
     ENDIF
@@ -250 +255 @@
           tmp_3 = clip_dvrp_s * superelevation_y
           tmp_4 = 0.0
-          tmp_5 = clip_dvrp_r * superelevation_x
-          tmp_6 = clip_dvrp_n * superelevation_y
+          tmp_5 = (clip_dvrp_r+dx) * superelevation_x
+          tmp_6 = (clip_dvrp_n+dy) * superelevation_y
           tmp_7 = zu(nz_do3d) * superelevation
           CALL DVRP_BOUNDINGBOX( m-1, 1, tmp_1, tmp_2, tmp_3, tmp_4, tmp_5, &
@@ -644 +649 @@
        ENDIF
 
+       IF ( mode_dvrp(m) == 'particles' )  THEN
+!
+!--       Define a default colourtable for particles
+          DO  i = 1, 11
+             interval_values_dvrp_prt(1,i) = i - 1.0
+             interval_values_dvrp_prt(2,i) = REAL( i )
+             interval_h_dvrp_prt(:,i) = 270.0 - ( i - 1.0 ) * 9.0
+          ENDDO
+
+          DO  j = 12, 22
+             interval_values_dvrp_prt(1,i) = i - 1.0
+             interval_values_dvrp_prt(2,i) = REAL( i )
+             interval_h_dvrp_prt(:,i) = 70.0 - ( i - 12.0 ) * 9.5
+          ENDDO
+
+          dvrp_colortable_entries_prt = 22
+
+       ENDIF
+
        m = m + 1