Changeset 308 for palm/trunk

Timestamp:

Apr 29, 2009 9:26:10 AM (15 years ago)

Author:

letzel

Message:

• Check for illegal entries in section_xy|xz|yz that exceed nz+1|ny+1|nx+1 (check_parameters)

Location:

palm/trunk

Files:

• DOC/tec/message_identifiers (modified) (1 diff)
• SOURCE/CURRENT_MODIFICATIONS (modified) (1 diff)
• SOURCE/check_parameters.f90 (modified) (2 diffs)

palm/trunk/DOC/tec/message_identifiers

@@ -583 +583 @@
                       because &the inversion height calculated by the prerun
                       is zero
-    
+    PA0319            section_xy must be <= nz + 1 = ...    
+    PA0320            section_xz must be <= ny + 1 = ...    
+    PA0321            section_yz must be <= nx + 1 = ...    
 
     UI0001            unknown location "..."

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -1 +1 @@
 New:
 ---
+Check for illegal entries in section_xy|xz|yz that exceed nz+1|ny+1|nx+1 
+(check_parameters)
+
 Clipping of dvrp output implemented. Default colourtable for particles
 implemented, particle attributes (color, dvrp_size) can be set with new

palm/trunk/SOURCE/check_parameters.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! Check for illegal entries in section_xy|xz|yz that exceed nz+1|ny+1|nx+1
 ! Coupling with independent precursor runs.
 ! Check particle_color, particle_dvrpsize, color_interval, dvrpsize_interval
@@ -2459 +2460 @@
 
 !
-!-- Store sectional planes in one shared array
+!-- Check sectional planes and store them in one shared array
+    IF ( ANY( section_xy > nz + 1 ) )  THEN
+       WRITE( message_string, * )  'section_xy must be <= nz + 1 = ', nz + 1
+       CALL message( 'check_parameters', 'PA0319', 1, 2, 0, 6, 0 )
+    ENDIF
+    IF ( ANY( section_xz > ny + 1 ) )  THEN
+       WRITE( message_string, * )  'section_xz must be <= ny + 1 = ', ny + 1
+       CALL message( 'check_parameters', 'PA0320', 1, 2, 0, 6, 0 )
+    ENDIF
+    IF ( ANY( section_yz > nx + 1 ) )  THEN
+       WRITE( message_string, * )  'section_yz must be <= nx + 1 = ', nx + 1
+       CALL message( 'check_parameters', 'PA0321', 1, 2, 0, 6, 0 )
+    ENDIF
     section(:,1) = section_xy
     section(:,2) = section_xz