Ignore:
Timestamp:
Sep 7, 2012 10:41:47 AM (9 years ago)
Author:
raasch
Message:

parameter use_prior_plot1d_parameters removed; little reformatting

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/init_masks.f90

    r979 r996  
    44! Current revisions:
    55! -----------------
    6 !
     6! little reformatting
    77!
    88! Former revisions:
     
    8080!-- Parallel mask output not yet tested
    8181    IF ( netcdf_data_format > 2 )  THEN
    82        message_string = 'NetCDF file formats '//                         &
    83                         '3 (NetCDF 4) and 4 (NetCDF 4 Classic model)'//  &
     82       message_string = 'netCDF file formats '//                         &
     83                        '3 (netCDF 4) and 4 (netCDF 4 Classic model)'//  &
    8484                        '&are currently not supported (not yet tested).'
    8585       CALL message( 'init_masks', 'PA0328', 1, 2, 0, 6, 0 )
     
    8888!
    8989!-- Store data output parameters for masked data output in few shared arrays
    90     DO mid = 1, masks
     90    DO  mid = 1, masks
    9191   
    9292       do_mask     (mid,:) = data_output_masks(mid,:)
     
    9696       mask      (mid,3,:) = mask_z(mid,:)
    9797       
    98        IF ( mask_x_loop(mid,1) == -1.0 .AND. mask_x_loop(mid,2) == -1.0  &
    99             .AND. mask_x_loop(mid,3) == -1.0 ) THEN
    100           mask_loop(mid,1,1:2) = -1.0
     98       IF ( mask_x_loop(mid,1) == -1.0  .AND. mask_x_loop(mid,2) == -1.0  &
     99            .AND.  mask_x_loop(mid,3) == -1.0 ) THEN
     100          mask_loop(mid,1,1:2) = -1.0
    101101          mask_loop(mid,1,3) = 0.0
    102102       ELSE
    103103          mask_loop(mid,1,:) = mask_x_loop(mid,:)
    104104       ENDIF
    105        IF ( mask_y_loop(mid,1) == -1.0 .AND. mask_y_loop(mid,2) == -1.0  &
    106             .AND. mask_y_loop(mid,3) == -1.0 ) THEN
    107           mask_loop(mid,2,1:2) = -1.0
     105       IF ( mask_y_loop(mid,1) == -1.0  .AND. mask_y_loop(mid,2) == -1.0  &
     106            .AND.  mask_y_loop(mid,3) == -1.0 ) THEN
     107          mask_loop(mid,2,1:2) = -1.0
    108108          mask_loop(mid,2,3) = 0.0
    109109       ELSE
    110110          mask_loop(mid,2,:) = mask_y_loop(mid,:)
    111111       ENDIF
    112        IF ( mask_z_loop(mid,1) == -1.0 .AND. mask_z_loop(mid,2) == -1.0  &
    113             .AND. mask_z_loop(mid,3) == -1.0 ) THEN
    114           mask_loop(mid,3,1:2) = -1.0
     112       IF ( mask_z_loop(mid,1) == -1.0  .AND. mask_z_loop(mid,2) == -1.0  &
     113            .AND.  mask_z_loop(mid,3) == -1.0 ) THEN
     114          mask_loop(mid,3,1:2) = -1.0
    115115          mask_loop(mid,3,3) = 0.0
    116116       ELSE
     
    158158!--    Check and set steering parameters for mask data output and averaging
    159159       i   = 1
    160        DO  WHILE ( do_mask(mid,i) /= ' '  .AND.  i <= 100 )
     160       DO WHILE ( do_mask(mid,i) /= ' '  .AND.  i <= 100 )
    161161!
    162162!--       Check for data averaging
     
    265265                     'for horizontal cross section'
    266266                CALL message( 'init_masks', 'PA0111', 1, 2, 0, 6, 0 )
    267 !                IF ( TRIM( var ) == 'lwp*'  .AND.  .NOT. cloud_physics )  THEN
    268 !                   WRITE ( message_string, * ) 'output of "', TRIM( var ), &
    269 !                         '" requires cloud_physics = .TRUE.'
    270 !                   CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
    271 !                ENDIF
    272 !                IF ( TRIM( var ) == 'pra*'  .AND.  .NOT. precipitation )  THEN
    273 !                   WRITE ( message_string, * ) 'output of "', TRIM( var ), &
    274 !                         '" requires precipitation = .TRUE.'
    275 !                   CALL message( 'init_masks', 'PA0112', 1, 2, 0, 6, 0 )
    276 !                ENDIF
    277 !                IF ( TRIM( var ) == 'pra*'  .AND.  j == 1 )  THEN
    278 !                   WRITE ( message_string, * ) 'temporal averaging of ', &
    279 !                         ' precipitation amount "', TRIM( var ),         &
    280 !                         '" not possible'
    281 !                   CALL message( 'init_masks', 'PA0113', 1, 2, 0, 6, 0 )
    282 !                ENDIF
    283 !                IF ( TRIM( var ) == 'prr*'  .AND.  .NOT. precipitation )  THEN
    284 !                   WRITE ( message_string, * ) 'output of "', TRIM( var ), &
    285 !                         '" requires precipitation = .TRUE.'
    286 !                   CALL message( 'init_masks', 'PA0112', 1, 2, 0, 6, 0 )
    287 !                ENDIF
    288 !
    289 !                IF ( TRIM( var ) == 'u*'   )  unit = 'm/s'
    290 !                IF ( TRIM( var ) == 't*'   )  unit = 'K'
    291 !                IF ( TRIM( var ) == 'lwp*' )  unit = 'kg/kg*m'
    292 !                IF ( TRIM( var ) == 'pra*' )  unit = 'mm'
    293 !                IF ( TRIM( var ) == 'prr*' )  unit = 'mm/s'
    294 !                IF ( TRIM( var ) == 'z0*'  )  unit = 'm'
    295267
    296268             CASE ( 'p', 'pt', 'u', 'v', 'w' )
     
    404376!--       If at least part of the mask resides on the PE, send the index limits
    405377!--       for the target array, otherwise send -9999 to PE0.
    406           IF ( mask_size_l(mid,1) > 0 .AND.  mask_size_l(mid,2) > 0 .AND. &
    407                mask_size_l(mid,3) > 0  ) &
    408                THEN
     378          IF ( mask_size_l(mid,1) > 0  .AND.  mask_size_l(mid,2) > 0  .AND.  &
     379               mask_size_l(mid,3) > 0  )  THEN
    409380             ind(1) = mask_start_l(mid,1)
    410381             ind(2) = mask_start_l(mid,1) + mask_size_l(mid,1) - 1
     
    599570          DO  m = loop_begin, loop_end, loop_stride
    600571             count = count + 1
    601              IF ( m >= lb .AND. m <= ub )  THEN
     572             IF ( m >= lb  .AND. m <= ub )  THEN
    602573                IF ( count_l == 0 )  mask_start_l(mid,dim) = count
    603574                count_l = count_l + 1
Note: See TracChangeset for help on using the changeset viewer.