Changeset 996 for palm/trunk/SOURCE/init_masks.f90

Timestamp:

Sep 7, 2012 10:41:47 AM (12 years ago)

Author:

raasch

Message:

parameter use_prior_plot1d_parameters removed; little reformatting

File:

• palm/trunk/SOURCE/init_masks.f90 (modified) (8 diffs)

palm/trunk/SOURCE/init_masks.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-!
+! little reformatting
 !
 ! Former revisions:
@@ -80 +80 @@
 !-- Parallel mask output not yet tested
     IF ( netcdf_data_format > 2 )  THEN
-       message_string = 'NetCDF file formats '//                         &
-                        '3 (NetCDF 4) and 4 (NetCDF 4 Classic model)'//  &
+       message_string = 'netCDF file formats '//                         &
+                        '3 (netCDF 4) and 4 (netCDF 4 Classic model)'//  &
                         '&are currently not supported (not yet tested).'
        CALL message( 'init_masks', 'PA0328', 1, 2, 0, 6, 0 )
@@ -88 +88 @@
 !
 !-- Store data output parameters for masked data output in few shared arrays
-    DO mid = 1, masks
+    DO  mid = 1, masks
     
        do_mask     (mid,:) = data_output_masks(mid,:)
@@ -96 +96 @@
        mask      (mid,3,:) = mask_z(mid,:) 
        
-       IF ( mask_x_loop(mid,1) == -1.0 .AND. mask_x_loop(mid,2) == -1.0  &
-            .AND. mask_x_loop(mid,3) == -1.0 ) THEN
-          mask_loop(mid,1,1:2) = -1.0
+       IF ( mask_x_loop(mid,1) == -1.0  .AND.  mask_x_loop(mid,2) == -1.0  &
+            .AND.  mask_x_loop(mid,3) == -1.0 )  THEN
+          mask_loop(mid,1,1:2) = -1.0
           mask_loop(mid,1,3) = 0.0
        ELSE
           mask_loop(mid,1,:) = mask_x_loop(mid,:)
        ENDIF
-       IF ( mask_y_loop(mid,1) == -1.0 .AND. mask_y_loop(mid,2) == -1.0  &
-            .AND. mask_y_loop(mid,3) == -1.0 ) THEN
-          mask_loop(mid,2,1:2) = -1.0
+       IF ( mask_y_loop(mid,1) == -1.0  .AND.  mask_y_loop(mid,2) == -1.0  &
+            .AND.  mask_y_loop(mid,3) == -1.0 )  THEN
+          mask_loop(mid,2,1:2) = -1.0
           mask_loop(mid,2,3) = 0.0
        ELSE
           mask_loop(mid,2,:) = mask_y_loop(mid,:)
        ENDIF
-       IF ( mask_z_loop(mid,1) == -1.0 .AND. mask_z_loop(mid,2) == -1.0  &
-            .AND. mask_z_loop(mid,3) == -1.0 ) THEN
-          mask_loop(mid,3,1:2) = -1.0
+       IF ( mask_z_loop(mid,1) == -1.0  .AND.  mask_z_loop(mid,2) == -1.0  &
+            .AND.  mask_z_loop(mid,3) == -1.0 )  THEN
+          mask_loop(mid,3,1:2) = -1.0
           mask_loop(mid,3,3) = 0.0
        ELSE
@@ -158 +158 @@
 !--    Check and set steering parameters for mask data output and averaging
        i   = 1
-       DO  WHILE ( do_mask(mid,i) /= ' '  .AND.  i <= 100 )
+       DO WHILE ( do_mask(mid,i) /= ' '  .AND.  i <= 100 )
 !
 !--       Check for data averaging
@@ -265 +265 @@
                      'for horizontal cross section'
                 CALL message( 'init_masks', 'PA0111', 1, 2, 0, 6, 0 )
-!                IF ( TRIM( var ) == 'lwp*'  .AND.  .NOT. cloud_physics )  THEN
-!                   WRITE ( message_string, * ) 'output of "', TRIM( var ), &
-!                         '" requires cloud_physics = .TRUE.'
-!                   CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
-!                ENDIF
-!                IF ( TRIM( var ) == 'pra*'  .AND.  .NOT. precipitation )  THEN
-!                   WRITE ( message_string, * ) 'output of "', TRIM( var ), &
-!                         '" requires precipitation = .TRUE.'
-!                   CALL message( 'init_masks', 'PA0112', 1, 2, 0, 6, 0 )
-!                ENDIF
-!                IF ( TRIM( var ) == 'pra*'  .AND.  j == 1 )  THEN
-!                   WRITE ( message_string, * ) 'temporal averaging of ', &
-!                         ' precipitation amount "', TRIM( var ),         &
-!                         '" not possible' 
-!                   CALL message( 'init_masks', 'PA0113', 1, 2, 0, 6, 0 )
-!                ENDIF
-!                IF ( TRIM( var ) == 'prr*'  .AND.  .NOT. precipitation )  THEN
-!                   WRITE ( message_string, * ) 'output of "', TRIM( var ), &
-!                         '" requires precipitation = .TRUE.'
-!                   CALL message( 'init_masks', 'PA0112', 1, 2, 0, 6, 0 )
-!                ENDIF
-!
-!                IF ( TRIM( var ) == 'u*'   )  unit = 'm/s'
-!                IF ( TRIM( var ) == 't*'   )  unit = 'K'
-!                IF ( TRIM( var ) == 'lwp*' )  unit = 'kg/kg*m'
-!                IF ( TRIM( var ) == 'pra*' )  unit = 'mm'
-!                IF ( TRIM( var ) == 'prr*' )  unit = 'mm/s'
-!                IF ( TRIM( var ) == 'z0*'  )  unit = 'm'
 
              CASE ( 'p', 'pt', 'u', 'v', 'w' )
@@ -404 +376 @@
 !--       If at least part of the mask resides on the PE, send the index limits
 !--       for the target array, otherwise send -9999 to PE0.
-          IF ( mask_size_l(mid,1) > 0 .AND.  mask_size_l(mid,2) > 0 .AND. &
-               mask_size_l(mid,3) > 0  ) &
-               THEN
+          IF ( mask_size_l(mid,1) > 0  .AND.  mask_size_l(mid,2) > 0  .AND.  &
+               mask_size_l(mid,3) > 0  )  THEN
              ind(1) = mask_start_l(mid,1)
              ind(2) = mask_start_l(mid,1) + mask_size_l(mid,1) - 1
@@ -599 +570 @@
           DO  m = loop_begin, loop_end, loop_stride
              count = count + 1
-             IF ( m >= lb .AND. m <= ub )  THEN
+             IF ( m >= lb  .AND.  m <= ub )  THEN
                 IF ( count_l == 0 )  mask_start_l(mid,dim) = count
                 count_l = count_l + 1