Ignore:
Timestamp:
Jul 22, 2014 2:14:06 PM (7 years ago)
Author:
heinze
Message:

+nr, qc, qr for mask output

File:
1 edited

Legend:

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

    r1415 r1438  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! +nr, qc, qr
    2323!
    2424! Former revisions:
     
    8585               data_output_masks, data_output_masks_user,                      &
    8686               doav, doav_n, domask, domask_no, dz, dz_stretch_level, humidity,&
    87                mask, masks, mask_scale, mask_i,                                &
     87               icloud_scheme, mask, masks, mask_scale, mask_i,                 &
    8888               mask_i_global, mask_j, mask_j_global, mask_k, mask_k_global,    &
    8989               mask_loop, mask_size, mask_size_l, mask_start_l, mask_x,        &
    9090               mask_x_loop, mask_xyz_dimension, mask_y, mask_y_loop, mask_z,   &
    9191               mask_z_loop, max_masks,  message_string, mid,                   &
    92                netcdf_data_format, passive_scalar, ocean
     92               netcdf_data_format, passive_scalar, precipitation, ocean
    9393
    9494    USE grid_variables,                                                        &
     
    252252                unit = 'K'
    253253
     254             CASE ( 'nr' )
     255                IF ( .NOT. cloud_physics )  THEN
     256                   WRITE ( message_string, * ) 'output of "', TRIM( var ), &
     257                        '" requires cloud_physics = .TRUE.'
     258                   CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
     259                 ELSEIF ( icloud_scheme /= 0 ) THEN
     260                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
     261                         'res cloud_scheme = seifert_beheng'
     262                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
     263                 ELSEIF ( .NOT. precipitation )  THEN
     264                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
     265                                 'res precipitation = .TRUE.'
     266                   CALL message( 'check_parameters', 'PA0112', 1, 2, 0, 6, 0 )
     267                ENDIF
     268                unit = '1/m3'
     269
    254270             CASE ( 'pc', 'pr' )
    255271                IF ( .NOT. particle_advection )  THEN
     
    271287                IF ( TRIM( var ) == 'vpt' )  unit = 'K'
    272288
     289             CASE ( 'qc' )
     290                IF ( .NOT. cloud_physics )  THEN
     291                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
     292                            'res cloud_physics = .TRUE.'
     293                   CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
     294                ELSEIF ( icloud_scheme /= 0 ) THEN
     295                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
     296                            'res cloud_scheme = seifert_beheng'
     297                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
     298                ENDIF
     299                unit = 'kg/kg'
     300
    273301             CASE ( 'ql' )
    274302                IF ( .NOT. ( cloud_physics  .OR.  cloud_droplets ) )  THEN
     
    295323                        '" requires cloud_physics = .TRUE.'
    296324                   CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
     325                ENDIF
     326                unit = 'kg/kg'
     327
     328             CASE ( 'qr' )
     329                IF ( .NOT. cloud_physics )  THEN
     330                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
     331                            'res cloud_physics = .TRUE.'
     332                   CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
     333                ELSEIF ( icloud_scheme /= 0 ) THEN
     334                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
     335                            'res cloud_scheme = seifert_beheng'
     336                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
     337                ELSEIF ( .NOT. precipitation )  THEN
     338                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
     339                                    'res precipitation = .TRUE.'
     340                   CALL message( 'check_parameters', 'PA0112', 1, 2, 0, 6, 0 )
    297341                ENDIF
    298342                unit = 'kg/kg'
Note: See TracChangeset for help on using the changeset viewer.