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

Timestamp:

Apr 22, 2014 3:03:56 PM (10 years ago)

Author:

boeske

Message:

large scale forcing enabled

File:

• palm/trunk/SOURCE/check_parameters.f90 (modified) (3 diffs)

palm/trunk/SOURCE/check_parameters.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Usage of large scale forcing for pt and q enabled:
+! output for profiles of large scale advection (td_lsa_lpt, td_lsa_q), 
+! large scale subsidence (td_sub_lpt, td_sub_q) 
+! and nudging tendencies (td_nud_lpt, td_nud_q, td_nud_u and td_nud_v) added,
+! check if use_subsidence_tendencies is used correctly
 ! 
 ! Former revisions:
@@ -1271 +1275 @@
        ENDIF
          
+    ENDIF
+
+!
+!-- Check if the control parameter use_subsidence_tendencies is used correctly
+    IF ( use_subsidence_tendencies  .AND.  .NOT. large_scale_subsidence )  THEN
+       message_string = 'The usage of use_subsidence_tendencies ' // &
+                            'requires large_scale_subsidence = .T..'
+       CALL message( 'check_parameters', 'PA0396', 1, 2, 0, 6, 0 )
+    ELSEIF ( use_subsidence_tendencies  .AND.  .NOT. large_scale_forcing )  THEN
+       message_string = 'The usage of use_subsidence_tendencies ' // &
+                            'requires large_scale_forcing = .T..'
+       CALL message( 'check_parameters', 'PA0397', 1, 2, 0, 6, 0 )
     ENDIF
 
@@ -2661 +2677 @@
              ENDIF
 
+          CASE ( 'td_lsa_lpt' )
+             IF ( .NOT. large_scale_forcing )  THEN
+                message_string = 'data_output_pr = ' // &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for large_scale_forcing = .FALSE.'
+                CALL message( 'check_parameters', 'PA0393', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 81
+                dopr_unit(i)  = 'K/s'
+                hom(:,2,81,:) = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
+          CASE ( 'td_lsa_q' )
+             IF ( .NOT. large_scale_forcing )  THEN
+                message_string = 'data_output_pr = ' // &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for large_scale_forcing = .FALSE.'
+                CALL message( 'check_parameters', 'PA0393', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 82
+                dopr_unit(i)  = 'kg/kgs'
+                hom(:,2,82,:) = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
+          CASE ( 'td_sub_lpt' )
+             IF ( .NOT. large_scale_forcing )  THEN
+                message_string = 'data_output_pr = ' // &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for large_scale_forcing = .FALSE.'
+                CALL message( 'check_parameters', 'PA0393', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 83
+                dopr_unit(i)  = 'K/s'
+                hom(:,2,83,:) = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
+          CASE ( 'td_sub_q' )
+             IF ( .NOT. large_scale_forcing )  THEN
+                message_string = 'data_output_pr = ' // &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for large_scale_forcing = .FALSE.'
+                CALL message( 'check_parameters', 'PA0393', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 84
+                dopr_unit(i)  = 'kg/kgs'
+                hom(:,2,84,:) = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
+          CASE ( 'td_nud_lpt' )
+             IF ( .NOT. nudging )  THEN
+                message_string = 'data_output_pr = ' // &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for nudging = .FALSE.'
+                CALL message( 'check_parameters', 'PA0394', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 85
+                dopr_unit(i)  = 'K/s'
+                hom(:,2,85,:) = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
+          CASE ( 'td_nud_q' )
+             IF ( .NOT. nudging )  THEN
+                message_string = 'data_output_pr = ' // &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for nudging = .FALSE.'
+                CALL message( 'check_parameters', 'PA0394', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 86
+                dopr_unit(i)  = 'kg/kgs'
+                hom(:,2,86,:) = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
+          CASE ( 'td_nud_u' )
+             IF ( .NOT. nudging )  THEN
+                message_string = 'data_output_pr = ' // &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for nudging = .FALSE.'
+                CALL message( 'check_parameters', 'PA0394', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 87
+                dopr_unit(i)  = 'm/s2'
+                hom(:,2,87,:) = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
+          CASE ( 'td_nud_v' )
+             IF ( .NOT. nudging )  THEN
+                message_string = 'data_output_pr = ' // &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for nudging = .FALSE.'
+                CALL message( 'check_parameters', 'PA0394', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 88
+                dopr_unit(i)  = 'm/s2'
+                hom(:,2,88,:) = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
+
           CASE DEFAULT