Changeset 1179 for palm/trunk/SOURCE/time_integration.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/time_integration.f90 (modified) (3 diffs)

palm/trunk/SOURCE/time_integration.f90

@@ -20 +20 @@
 ! Current revisions:
 ! ------------------
-!
+! mean profiles for reference state are only calculated if required,
+! small bugfix for background communication
 !
 ! Former revisions:
@@ -230 +231 @@
 !--       Calculate those variables needed in the tendency terms which need
 !--       global communication
-          IF ( .NOT. neutral )  CALL calc_mean_profile( pt, 4 )
-          IF ( ocean         )  CALL calc_mean_profile( rho, 64 )
-          IF ( humidity      )  CALL calc_mean_profile( vpt, 44 )
+          IF ( .NOT. use_single_reference_value  .AND. &
+               .NOT. use_initial_profile_as_reference )  THEN
+!
+!--          Horizontally averaged profiles to be used as reference state in
+!--          buoyancy terms (WARNING: only the respective last call of
+!--          calc_mean_profile defines the reference state!)
+             IF ( .NOT. neutral )  CALL calc_mean_profile( pt, 4 )
+             IF ( ocean         )  CALL calc_mean_profile( rho, 64 )
+             IF ( humidity      )  CALL calc_mean_profile( vpt, 44 )
+          ENDIF
+
           IF ( .NOT. constant_diffusion )  CALL production_e_init
           IF ( ( ws_scheme_mom .OR. ws_scheme_sca )  .AND.  &
@@ -438 +447 @@
              IF (humidity  .OR.  passive_scalar)  THEN
                 CALL exchange_horiz( q_p, nbgp )
-                IF ( cloud_physics .AND. icloud_scheme == 0 )  THEN
+                IF ( cloud_physics .AND. icloud_scheme == 0  .AND.  &
+                     precipitation )  THEN
                    CALL exchange_horiz( qr_p, nbgp )
                    CALL exchange_horiz( nr_p, nbgp )