Changeset 410 for palm/trunk/SOURCE/check_parameters.f90

Timestamp:

Dec 4, 2009 5:05:40 PM (15 years ago)

Author:

letzel

Message:

• reintegrate branch letzel/masked_output into trunk; new funtionality: masked data output (not yet documented)

Location:

palm/trunk

Files:

• . (modified) (1 prop)
• SOURCE/check_parameters.f90 (modified) (4 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -5 +5 @@
 ! -----------------
 ! Enabled passive scalar/humidity wall fluxes for non-flat topography
+!
+! Branch revisions:
+! -----------------
+! Masked data output
 !
 ! Former revisions:
@@ -1477 +1481 @@
        IF ( dt_do3d           == 9999999.9 )  dt_do3d           = dt_data_output
        IF ( dt_data_output_av == 9999999.9 )  dt_data_output_av = dt_data_output
+       DO  mid = 1, 20
+          IF ( dt_domask(mid) == 9999999.9 )  dt_domask(mid)    = dt_data_output
+       ENDDO
     ENDIF
 
@@ -1495 +1502 @@
     IF ( skip_time_data_output_av == 9999999.9 ) &
                                 skip_time_data_output_av = skip_time_data_output
+    DO  mid = 1, 20
+       IF ( skip_time_domask(mid) == 9999999.9 ) &
+                                skip_time_domask(mid)    = skip_time_data_output
+    ENDDO
 
 !
@@ -2647 +2658 @@
 
 !
+!-- Check mask conditions
+    IF ( masks < 0 .OR. masks > max_masks )  THEN
+       WRITE( message_string, * )  'illegal value: masks must be >= 0 and ', &
+            '<= ', max_masks, ' (=max_masks)'
+       CALL message( 'check_parameters', 'PA9998', 1, 2, 0, 6, 0 )
+    ENDIF
+    IF ( masks > 0 )  THEN
+       mask_scale(1) = mask_scale_x
+       mask_scale(2) = mask_scale_y
+       mask_scale(3) = mask_scale_z
+       IF ( ANY( mask_scale <= 0.0 ) )  THEN
+          WRITE( message_string, * )  &
+               'illegal value: mask_scale_x, mask_scale_y and mask_scale_z', &
+               'must be > 0.0'
+          CALL message( 'check_parameters', 'PA9998', 1, 2, 0, 6, 0 )
+       ENDIF
+!
+!--    Generate masks for masked data output
+       CALL init_masks
+    ENDIF
+
+!
 !-- Check netcdf precison
     ldum = .FALSE.