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

Timestamp:

Jul 22, 2014 2:14:06 PM (10 years ago)

Author:

heinze

Message:

+nr, qc, qr for mask output

File:

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

palm/trunk/SOURCE/init_masks.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! +nr, qc, qr
 ! 
 ! Former revisions:
@@ -85 +85 @@
                data_output_masks, data_output_masks_user,                      &
                doav, doav_n, domask, domask_no, dz, dz_stretch_level, humidity,&
-               mask, masks, mask_scale, mask_i,                                &
+               icloud_scheme, mask, masks, mask_scale, mask_i,                 &
                mask_i_global, mask_j, mask_j_global, mask_k, mask_k_global,    &
                mask_loop, mask_size, mask_size_l, mask_start_l, mask_x,        &
                mask_x_loop, mask_xyz_dimension, mask_y, mask_y_loop, mask_z,   &
                mask_z_loop, max_masks,  message_string, mid,                   &
-               netcdf_data_format, passive_scalar, ocean
+               netcdf_data_format, passive_scalar, precipitation, ocean
 
     USE grid_variables,                                                        &
@@ -252 +252 @@
                 unit = 'K'
 
+             CASE ( 'nr' )
+                IF ( .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 )
+                 ELSEIF ( icloud_scheme /= 0 ) THEN
+                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
+                         'res cloud_scheme = seifert_beheng'
+                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
+                 ELSEIF ( .NOT. precipitation )  THEN
+                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
+                                 'res precipitation = .TRUE.'
+                   CALL message( 'check_parameters', 'PA0112', 1, 2, 0, 6, 0 )
+                ENDIF
+                unit = '1/m3'
+
              CASE ( 'pc', 'pr' )
                 IF ( .NOT. particle_advection )  THEN
@@ -271 +287 @@
                 IF ( TRIM( var ) == 'vpt' )  unit = 'K'
 
+             CASE ( 'qc' )
+                IF ( .NOT. cloud_physics )  THEN
+                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
+                            'res cloud_physics = .TRUE.'
+                   CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
+                ELSEIF ( icloud_scheme /= 0 ) THEN
+                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
+                            'res cloud_scheme = seifert_beheng'
+                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
+                ENDIF
+                unit = 'kg/kg'
+
              CASE ( 'ql' )
                 IF ( .NOT. ( cloud_physics  .OR.  cloud_droplets ) )  THEN
@@ -295 +323 @@
                         '" requires cloud_physics = .TRUE.'
                    CALL message( 'init_masks', 'PA0108', 1, 2, 0, 6, 0 )
+                ENDIF
+                unit = 'kg/kg'
+
+             CASE ( 'qr' )
+                IF ( .NOT. cloud_physics )  THEN
+                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
+                            'res cloud_physics = .TRUE.'
+                   CALL message( 'check_parameters', 'PA0108', 1, 2, 0, 6, 0 )
+                ELSEIF ( icloud_scheme /= 0 ) THEN
+                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
+                            'res cloud_scheme = seifert_beheng'
+                   CALL message( 'check_parameters', 'PA0359', 1, 2, 0, 6, 0 )
+                ELSEIF ( .NOT. precipitation )  THEN
+                   message_string = 'output of "' // TRIM( var ) // '" requi' //  &
+                                    'res precipitation = .TRUE.'
+                   CALL message( 'check_parameters', 'PA0112', 1, 2, 0, 6, 0 )
                 ENDIF
                 unit = 'kg/kg'