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

palm/trunk/SOURCE/init_ocean.f90

@@ -20 +20 @@
 ! Current revisions:
 ! -----------------
-! 
+! Initial density profile is stored in array hom
 !
 ! Former revisions:
@@ -47 +47 @@
     USE control_parameters
     USE eqn_state_seawater_mod
-    USE pegrid
     USE grid_variables
     USE indices
+    USE pegrid
+    USE statistics
 
     IMPLICIT NONE
@@ -126 +127 @@
     CALL eqn_state_seawater
 
+!
+!-- Store initial density profile
+    hom(:,1,77,:)  = SPREAD( rho_init(:), 2, statistic_regions+1 )
+
+!
+!-- Set the reference state to be used in the buoyancy terms
+    IF ( use_single_reference_value )  THEN
+       ref_state(:) = prho_reference
+    ELSE
+       ref_state(:) = rho_init(:)
+    ENDIF
+
 
  END SUBROUTINE init_ocean