Changeset 2918 for palm/trunk/SOURCE/model_1d_mod.f90

Timestamp:

Mar 21, 2018 3:52:14 PM (6 years ago)

Author:

gronemeier

Message:

moved initialization of mixing length to turbulence_closure_mod; consider walls in calculation of ml in rans-mode

Location:

palm/trunk

Files:

• . (modified) (1 prop)
• SOURCE (modified) (1 prop)
• SOURCE/model_1d_mod.f90 (modified) (11 diffs)

palm/trunk/SOURCE/model_1d_mod.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! - rename l_black into l1d_init
+! - calculate l_grid within init_1d_model and save it as l1d_init
+! 
+! 2718 2018-01-02 08:49:38Z maronga
 ! Corrected "Former revisions" section
 ! 
@@ -128 +132 @@
 !> @todo harmonize code with new surface_layer_fluxes module
 !> @bug 1D model crashes when using small grid spacings in the order of 1 m
+!> @fixme option "as_in_3d_model" seems to eb an inappropriate option because
+!>        the 1D model uses different turbulence closure approaches at least if
+!>        the 3D model is set to LES-mode.
 !------------------------------------------------------------------------------!
  MODULE model_1d_mod
 
     USE arrays_3d,                                                             &
-        ONLY:  dd2zu, ddzu, ddzw, dzu, l_grid, pt_init, q_init, ug, u_init,    &
+        ONLY:  dd2zu, ddzu, ddzw, dzu, dzw, pt_init, q_init, ug, u_init,       &
                vg, v_init, zu
     
@@ -190 +197 @@
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  kh1d     !< turbulent diffusion coefficient for heat (1d-model)
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  km1d     !< turbulent diffusion coefficient for momentum (1d-model)
-    REAL(wp), DIMENSION(:), ALLOCATABLE ::  l_black  !< mixing length Blackadar (1d-model)
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  l1d      !< mixing length for turbulent diffusion coefficients (1d-model)
+    REAL(wp), DIMENSION(:), ALLOCATABLE ::  l1d_init !< initial mixing length (1d-model)
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  l1d_diss !< mixing length for dissipation (1d-model)
     REAL(wp), DIMENSION(:), ALLOCATABLE ::  rif1d    !< Richardson flux number (1d-model)
@@ -255 +262 @@
  SUBROUTINE init_1d_model
  
+    USE grid_variables,                                                        &
+        ONLY:  dx, dy
+
     IMPLICIT NONE
 
@@ -267 +277 @@
     ALLOCATE( diss1d(nzb:nzt+1), diss1d_p(nzb:nzt+1),                          &
               e1d(nzb:nzt+1), e1d_p(nzb:nzt+1), kh1d(nzb:nzt+1),               &
-              km1d(nzb:nzt+1), l_black(nzb:nzt+1), l1d(nzb:nzt+1),             &
+              km1d(nzb:nzt+1), l1d(nzb:nzt+1), l1d_init(nzb:nzt+1),            &
               l1d_diss(nzb:nzt+1), rif1d(nzb:nzt+1), te_diss(nzb:nzt+1),       &
               te_dissm(nzb:nzt+1), te_e(nzb:nzt+1),                            &
@@ -286 +296 @@
 !
 !--    Compute the mixing length
-       l_black(nzb) = 0.0_wp
+       l1d_init(nzb) = 0.0_wp
 
        IF ( TRIM( mixing_length_1d ) == 'blackadar' )  THEN
@@ -299 +309 @@
 
           DO  k = nzb+1, nzt+1
-             l_black(k) = kappa * zu(k) / ( 1.0_wp + kappa * zu(k) / lambda )
+             l1d_init(k) = kappa * zu(k) / ( 1.0_wp + kappa * zu(k) / lambda )
           ENDDO
 
        ELSEIF ( TRIM( mixing_length_1d ) == 'as_in_3d_model' )  THEN
 !
-!--       Use the same mixing length as in 3D model
-          l_black(1:nzt) = l_grid
-          l_black(nzt+1) = l_black(nzt)
+!--       Use the same mixing length as in 3D model (LES-mode)
+          !### TODO: rename (delete?) this option
+          !###  As the mixing length is different between RANS and LES mode, it
+          !###  must be distinguished here between these modes. For RANS mode,
+          !###  the mixing length is calculated accoding to Blackadar, which is
+          !###  the other option at this point.
+          !###  Maybe delete this option entirely (not appropriate in LES case)
+          !###  2018-03-20, gronemeier
+          DO  k = nzb+1, nzt
+             l1d_init(k)  = ( dx * dy * dzw(k) )**0.33333333333333_wp
+          ENDDO
+          l1d_init(nzt+1) = l1d_init(nzt)
 
        ENDIF
     ENDIF
-    l1d      = l_black
-    l1d_diss = l_black
+    l1d      = l1d_init
+    l1d_diss = l1d_init
     u1d      = u_init
     u1d_p    = u_init
@@ -457 +476 @@
                    diss1d(k) = c_m**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)
                 ELSEIF ( dissipation_1d == 'as_in_3d_model' )  THEN
-                   diss1d(k) = ( 0.19_wp + 0.74_wp * l1d_diss(k) /  l_grid(k)  &
+                   diss1d(k) = ( 0.19_wp + 0.74_wp * l1d_diss(k) / l1d_init(k) &
                                ) * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)
                 ENDIF
@@ -520 +539 @@
                 diss1d(k) = c_m**3 * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)
              ELSEIF ( dissipation_1d == 'as_in_3d_model' )  THEN
-                diss1d(k) = ( 0.19_wp + 0.74_wp * l1d_diss(k) / l_grid(k) )    &
+                diss1d(k) = ( 0.19_wp + 0.74_wp * l1d_diss(k) / l1d_init(k) )  &
                             * e1d(k) * SQRT( e1d(k) ) / l1d_diss(k)
              ENDIF
@@ -804 +823 @@
                 DO  k = nzb+1, nzt
                    IF ( rif1d(k) >= 0.0_wp )  THEN
-                      l1d(k) = l_black(k) / ( 1.0_wp + 5.0_wp * rif1d(k) )
+                      l1d(k) = l1d_init(k) / ( 1.0_wp + 5.0_wp * rif1d(k) )
                       l1d_diss(k) = l1d(k)
                    ELSE
-                      l1d(k) = l_black(k)
-                      l1d_diss(k) = l_black(k) *                               &
+                      l1d(k) = l1d_init(k)
+                      l1d_diss(k) = l1d_init(k) *                              &
                                     SQRT( 1.0_wp - 16.0_wp * rif1d(k) )
                    ENDIF
@@ -819 +838 @@
                                      SQRT( g / pt_init(k) * dpt_dz ) + 1E-5_wp
                    ELSE
-                      l_stable = l_grid(k)
+                      l_stable = l1d_init(k)
                    ENDIF
-                   l1d(k) = MIN( l_grid(k), l_stable )
+                   l1d(k) = MIN( l1d_init(k), l_stable )
                    l1d_diss(k) = l1d(k)
                 ENDDO