Changeset 2334 for palm/trunk

Timestamp:

Aug 4, 2017 11:57:04 AM (7 years ago)

Author:

gronemeier

Message:

set c_m to 0.4 within 1D model

File:

• palm/trunk/SOURCE/init_1d_model.f90 (modified) (7 diffs)

palm/trunk/SOURCE/init_1d_model.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! set c_m = 0.4 according to Detering and Etling (1985)
+! 
+! 2299 2017-06-29 10:14:38Z maronga
 ! Removed german text
 !
@@ -273 +276 @@
     REAL(wp) ::  a            !<
     REAL(wp) ::  b            !<
+    REAL(wp) ::  c_m = 0.4_wp !< model constant, 0.4 according to Detering and Etling (1985)
     REAL(wp) ::  dissipation  !<
     REAL(wp) ::  dpt_dz       !<
@@ -473 +477 @@
 !--       Boundary conditions for the prognostic variables.
 !--       At the top boundary (nzt+1) u,v and e keep their initial values 
-!--       (ug(nzt+1), vg(nzt+1), 0), at the bottom boundary the mirror 
-!--       boundary condition applies to u and v.
-!--       The boundary condition for e is set further below ( (u*/cm)**2 ).
-         ! u1d_p(nzb) = -u1d_p(nzb+1)
-         ! v1d_p(nzb) = -v1d_p(nzb+1)
-
+!--       (ug(nzt+1), vg(nzt+1), 0).
+!--       At the bottom boundary, Dirichlet condition is used for u and v (0)
+!--       and Neumann condition for e (e(nzb)=e(nzb+1)).
           u1d_p(nzb) = 0.0_wp
           v1d_p(nzb) = 0.0_wp
@@ -499 +500 @@
 !
 !--             Compute theta* using Rif numbers of the previous time step
-                IF ( rif1d(1) >= 0.0_wp )  THEN
+                IF ( rif1d(nzb+1) >= 0.0_wp )  THEN
 !
 !--                Stable stratification
@@ -628 +629 @@
 !--             compatibility with the 3D model.
                 IF ( ibc_e_b == 2 )  THEN
-                   e1d(nzb+1) = ( us1d / 0.1_wp )**2
-!                  e1d(nzb+1) = ( us1d / 0.4_wp )**2  !not used so far, see also
-                                                      !prandtl_fluxes
+                   e1d(nzb+1) = ( us1d / c_m )**2
                 ENDIF
                 e1d(nzb) = e1d(nzb+1)
@@ -637 +636 @@
 !
 !--                Compute q* 
-                   IF ( rif1d(1) >= 0.0_wp )  THEN
+                   IF ( rif1d(nzb+1) >= 0.0_wp )  THEN
 !
 !--                   Stable stratification
@@ -714 +713 @@
              ENDIF
              DO  k = nzb_diff, nzt
-!                km1d(k) = 0.4 * SQRT( e1d(k) ) !changed: adjustment to 3D-model
-                km1d(k) = 0.1_wp * SQRT( e1d(k) )
+                km1d(k) = c_m * SQRT( e1d(k) )
+
                 IF ( rif1d(k) >= 0.0_wp )  THEN
                    km1d(k) = km1d(k) * l1d(k)