Changeset 97 for palm/trunk/SOURCE/eqn_state_seawater.f90

Timestamp:

Jun 21, 2007 8:23:15 AM (17 years ago)

Author:

raasch

Message:

New:
---
ocean version including prognostic equation for salinity and equation of state for seawater. Routine buoyancy can be used with both temperature and density.
+ inipar-parameters bc_sa_t, bottom_salinityflux, ocean, sa_surface, sa_vertical_gradient, sa_vertical_gradient_level, top_salinityflux

advec_s_bc, average_3d_data, boundary_conds, buoyancy, check_parameters, data_output_2d, data_output_3d, diffusion_e, flow_statistics, header, init_grid, init_3d_model, modules, netcdf, parin, production_e, prognostic_equations, read_var_list, sum_up_3d_data, swap_timelevel, time_integration, user_interface, write_var_list, write_3d_binary

New:
eqn_state_seawater, init_ocean

Changed:

---

inipar-parameter use_pt_reference renamed use_reference

hydro_press renamed hyp, routine calc_mean_pt_profile renamed calc_mean_profile

format adjustments for the ocean version (run_control)

advec_particles, buoyancy, calc_liquid_water_content, check_parameters, diffusion_e, diffusivities, header, init_cloud_physics, modules, production_e, prognostic_equations, run_control

Errors:

---

Bugfix: height above topography instead of height above level k=0 is used for calculating the mixing length (diffusion_e and diffusivities).

Bugfix: error in boundary condition for TKE removed (advec_s_bc)

advec_s_bc, diffusion_e, prognostic_equations

File:

• palm/trunk/SOURCE/eqn_state_seawater.f90 (modified) (6 diffs)

palm/trunk/SOURCE/eqn_state_seawater.f90

@@ -8 +8 @@
 ! Former revisions:
 ! -----------------
-! $Id:$
+! $Id$
 !
 ! Initial revision
@@ -18 +18 @@
 ! salinity, and pressure.
 ! For coefficients see Jackett et al., 2006: J. Atm. Ocean Tech.
+! eqn_state_seawater calculates the potential density referred at hyp(0).
+! eqn_state_seawater_func calculates density.
 !------------------------------------------------------------------------------!
 
@@ -70 +72 @@
        DO  i = nxl, nxr
           DO  j = nys, nyn
-             DO  k = nzb_u_inner(j,i)+1, nzt
+             DO  k = nzb_s_inner(j,i)+1, nzt
 !
 !--             Pressure is needed in dbar
-                p1 = hyp(k) * 1E-4
+                p1 = hyp(0) * 1E-4
                 p2 = p1 * p1
                 p3 = p2 * p1
@@ -102 +104 @@
 
              ENDDO
+!
+!--          Neumann conditions are assumed at bottom and top boundary
+             rho(nzt+1,j,i)            = rho(nzt,j,i)
+             rho(nzb_s_inner(j,i),j,i) = rho(nzb_s_inner(j,i)+1,j,i)
           ENDDO
        ENDDO
@@ -122 +128 @@
        REAL ::  p1, p2, p3, pt1, pt2, pt3, pt4, sa1, sa15, sa2
 
-       DO  k = nzb_u_inner(j,i)+1, nzt
+       DO  k = nzb_s_inner(j,i)+1, nzt
 !
 !--       Pressure is needed in dbar
-          p1 = hyp(k) * 1E-4
+          p1 = hyp(0) * 1E-4
           p2 = p1 * p1
           p3 = p2 * p1
@@ -152 +158 @@
                        )
        ENDDO
+!
+!--    Neumann conditions are assumed at bottom and top boundary
+       rho(nzt+1,j,i)            = rho(nzt,j,i)
+       rho(nzb_s_inner(j,i),j,i) = rho(nzb_s_inner(j,i)+1,j,i)
 
     END SUBROUTINE eqn_state_seawater_ij