Ignore:
Timestamp:
Jun 14, 2013 5:57:58 AM (8 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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/init_ocean.f90

    r1037 r1179  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! Initial density profile is stored in array hom
    2323!
    2424! Former revisions:
     
    4747    USE control_parameters
    4848    USE eqn_state_seawater_mod
    49     USE pegrid
    5049    USE grid_variables
    5150    USE indices
     51    USE pegrid
     52    USE statistics
    5253
    5354    IMPLICIT NONE
     
    126127    CALL eqn_state_seawater
    127128
     129!
     130!-- Store initial density profile
     131    hom(:,1,77,:)  = SPREAD( rho_init(:), 2, statistic_regions+1 )
     132
     133!
     134!-- Set the reference state to be used in the buoyancy terms
     135    IF ( use_single_reference_value )  THEN
     136       ref_state(:) = prho_reference
     137    ELSE
     138       ref_state(:) = rho_init(:)
     139    ENDIF
     140
    128141
    129142 END SUBROUTINE init_ocean
Note: See TracChangeset for help on using the changeset viewer.