Changeset 97 for palm/trunk/SOURCE/flow_statistics.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/flow_statistics.f90 (modified) (6 diffs)

palm/trunk/SOURCE/flow_statistics.f90

@@ -5 +5 @@
 ! -----------------
 ! Statistics for ocean version (salinity, density) added,
-! calculation of Deardorff velocity scale adjusted to be used with the ocean
-! version (HAS STILL TO BE COMPLETED!!!)
+! calculation of z_i and Deardorff velocity scale adjusted to be used with
+! the ocean version
 !
 ! Former revisions:
@@ -558 +558 @@
 !
 !--             Salinity flux and density (density does not belong to here,
-!--             but so far there is no suitable place to calculate)
+!--             but so far there is no other suitable place to calculate)
                 IF ( ocean )  THEN
                    pts = 0.5 * ( sa(k,j,i)   - hom(k,1,23,sr) + &
@@ -611 +611 @@
 
 !
+!--    Density at top follows Neumann condition
+       IF ( ocean )  sums_l(nzt+1,64,tn) = sums_l(nzt,64,tn)
+
+!
 !--    Divergence of vertical flux of resolved scale energy and pressure
 !--    fluctuations. First calculate the products, then the divergence.
@@ -709 +713 @@
 !
 !--       Fluxes at the surface must be zero (e.g. due to the Prandtl-layer)
-          sums(nzb,58) = 0.0
-          sums(nzb,59) = 0.0
-          sums(nzb,60) = 0.0
-          sums(nzb,61) = 0.0
-          sums(nzb,62) = 0.0
-          sums(nzb,63) = 0.0
+          sums_l(nzb,58,tn) = 0.0
+          sums_l(nzb,59,tn) = 0.0
+          sums_l(nzb,60,tn) = 0.0
+          sums_l(nzb,61,tn) = 0.0
+          sums_l(nzb,62,tn) = 0.0
+          sums_l(nzb,63,tn) = 0.0
 
        ENDIF
@@ -857 +861 @@
        z_i(1) = 0.0
        first = .TRUE.
-!       IF ( .NOT. ocean )  THEN
+       IF ( ocean )  THEN
+          DO  k = nzt, nzb+1, -1
+             IF ( first .AND. hom(k,1,18,sr) < 0.0 )  THEN
+                first = .FALSE.
+                height = zw(k)
+             ENDIF
+             IF ( hom(k,1,18,sr) < 0.0  .AND. &
+                  hom(k-1,1,18,sr) > hom(k,1,18,sr) )  THEN
+                IF ( zw(k) < 1.5 * height )  THEN
+                   z_i(1) = zw(k)
+                ELSE
+                   z_i(1) = height
+                ENDIF
+                EXIT
+             ENDIF
+          ENDDO
+       ELSE
           DO  k = nzb, nzt-1
              IF ( first .AND. hom(k,1,18,sr) < 0.0 )  THEN
@@ -873 +893 @@
              ENDIF
           ENDDO
-!       ELSE
-!       ENDIF
-
-!
-!--    Second scheme: Starting from the top model boundary, look for the first 
-!--    characteristic kink in the temperature profile, where the originally 
-!--    stable stratification notably weakens.
+       ENDIF
+
+!
+!--    Second scheme: Starting from the top/bottom model boundary, look for
+!--    the first characteristic kink in the temperature profile, where the
+!--    originally stable stratification notably weakens.
        z_i(2) = 0.0
-       DO  k = nzt-1, nzb+1, -1
-          IF ( ( hom(k+1,1,4,sr) - hom(k,1,4,sr) ) > &
-               2.0 * ( hom(k,1,4,sr) - hom(k-1,1,4,sr) ) )  THEN
-             z_i(2) = zu(k)
-             EXIT
-          ENDIF
-       ENDDO
+       IF ( ocean )  THEN
+          DO  k = nzb+1, nzt-1
+             IF ( ( hom(k,1,4,sr) - hom(k-1,1,4,sr) ) > &
+                  2.0 * ( hom(k+1,1,4,sr) - hom(k,1,4,sr) ) )  THEN
+                z_i(2) = zu(k)
+                EXIT
+             ENDIF
+          ENDDO
+       ELSE
+          DO  k = nzt-1, nzb+1, -1
+             IF ( ( hom(k+1,1,4,sr) - hom(k,1,4,sr) ) > &
+                  2.0 * ( hom(k,1,4,sr) - hom(k-1,1,4,sr) ) )  THEN
+                z_i(2) = zu(k)
+                EXIT
+             ENDIF
+          ENDDO
+       ENDIF
 
        hom(nzb+6,1,pr_palm,sr) = z_i(1)