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

Timestamp:

Jun 14, 2013 5:57:58 AM (11 years ago)

Author:

raasch

Message:

New:
---
Initial profiles can be used as reference state in the buoyancy term. New parameter
reference_state introduced. Calculation and handling of reference state in buoyancy term revised.
binary version for restart files changed from 3.9 to 3.9a (no downward compatibility!),
initial profile for rho added to hom (id=77)

Errors:

---

small bugfix for background communication (time_integration)

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! checks and settings of buoyancy parameters and switches revised,
+! initial profile for rho added to hom (id=77)
 !
 ! Former revisions:
@@ -1323 +1324 @@
        large_scale_subsidence = .TRUE.
        CALL init_w_subsidence
-    END IF
- 
-             
+    ENDIF
 
 !
@@ -1333 +1332 @@
 
 !
-!-- Ocean runs always use reference values in the buoyancy term. Therefore
-!-- set the reference temperature equal to the surface temperature.
-    IF ( ocean  .AND.  pt_reference == 9999999.9 )  pt_reference = pt_surface
-
-!
-!-- Reference value has to be used in buoyancy terms
-    IF ( pt_reference /= 9999999.9 )  use_reference = .TRUE.
+!-- Check and set buoyancy related parameters and switches
+    IF ( reference_state == 'horizontal_average' )  THEN
+       CONTINUE
+    ELSEIF ( reference_state == 'initial_profile' )  THEN
+       use_initial_profile_as_reference = .TRUE.
+    ELSEIF ( reference_state == 'single_value' )  THEN
+       use_single_reference_value = .TRUE.
+       IF ( pt_reference == 9999999.9 )  pt_reference = pt_surface
+       vpt_reference = pt_reference * ( 1.0 + 0.61 * q_surface )
+    ELSE
+       message_string = 'illegal value for reference_state: "' // &
+                        TRIM( reference_state ) // '"'
+       CALL message( 'check_parameters', 'PA0056', 1, 2, 0, 6, 0 )
+    ENDIF
+
+!
+!-- Ocean runs always use reference values in the buoyancy term
+    IF ( ocean )  THEN
+       reference_state = 'single_value'
+       use_single_reference_value = .TRUE.
+    ENDIF
 
 !
@@ -2463 +2476 @@
                 dopr_unit(i)  = 'kg/m3'
                 hom(:,2,64,:) = SPREAD( zu, 2, statistic_regions+1 )
+                IF ( data_output_pr(i)(1:1) == '#' )  THEN
+                   dopr_initial_index(i) = 77
+                   hom(:,2,77,:)         = SPREAD( zu, 2, statistic_regions+1 )
+                   hom(nzb,2,77,:)       = 0.0    ! because zu(nzb) is negative
+                   data_output_pr(i)     = data_output_pr(i)(2:)
+                ENDIF
              ENDIF