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

Timestamp:

Mar 6, 2014 1:15:21 PM (10 years ago)

Author:

heinze

Message:

enable usage of large_scale subsidence in combination with large_scale_forcing

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! enable usage of large_scale subsidence in combination with large_scale_forcing
+! output for profile of large scale vertical velocity w_subs added
 !
 ! Former revisions:
@@ -1354 +1355 @@
 !
 !-- Initialize large scale subsidence if required
-    IF ( subs_vertical_gradient_level(1) /= -9999999.9 )  THEN
-       large_scale_subsidence = .TRUE.
-       CALL init_w_subsidence
-    ENDIF
+    If ( large_scale_subsidence )  THEN
+       IF ( subs_vertical_gradient_level(1) /= -9999999.9 .AND. &
+                                     .NOT. large_scale_forcing )  THEN
+          CALL init_w_subsidence
+       ENDIF
+!
+!--    In case large_scale_forcing is used, profiles for subsidence velocity
+!--    are read in from file LSF_DATA
+
+       IF ( subs_vertical_gradient_level(1) == -9999999.9 .AND. &
+                                     .NOT. large_scale_forcing )  THEN
+          message_string = 'There is no default large scale vertical ' // &
+                           'velocity profile set. Specify the subsidence ' // &
+                           'velocity profile via subs_vertical_gradient and ' // &
+                           'subs_vertical_gradient_level.'
+          CALL message( 'check_parameters', 'PA0380', 1, 2, 0, 6, 0 )
+       ENDIF
+    ELSE
+        IF ( subs_vertical_gradient_level(1) /= -9999999.9 )  THEN
+           message_string = 'Enable usage of large scale subsidence by ' // &
+                            'setting large_scale_subsidence = .T..'
+          CALL message( 'check_parameters', 'PA0381', 1, 2, 0, 6, 0 )
+        ENDIF
+    ENDIF    
 
 !
@@ -2723 +2744 @@
              hom(:,2,79,:) = SPREAD( zu, 2, statistic_regions+1 )
 
+          CASE ( 'w_subs' )
+             IF ( .NOT. large_scale_subsidence )  THEN
+                message_string = 'data_output_pr = ' // &
+                                 TRIM( data_output_pr(i) ) // ' is not imp' // &
+                                 'lemented for large_scale_subsidence = .FALSE.'
+                CALL message( 'check_parameters', 'PA0382', 1, 2, 0, 6, 0 )
+             ELSE
+                dopr_index(i) = 80
+                dopr_unit(i)  = 'm/s'
+                hom(:,2,80,:) = SPREAD( zu, 2, statistic_regions+1 )
+             ENDIF
+
           CASE DEFAULT