Changeset 767 for palm/trunk/SOURCE/check_parameters.f90

Timestamp:

Oct 14, 2011 6:39:12 AM (13 years ago)

Author:

raasch

Message:

New:
---

Flow field initialization with given (e.g. measured) profiles. Profile data
for u-,v-velocity components + respective heights are given with new
inipar-parameters u_profile, v_profile, and uv_heights. Final profiles are
calculated from these given profiles by linear interpolation.
(check_parameters, header, init_3d_model, modules, parin)

Changed:

---

ug,vg replaced by u_init,v_init as the Dirichlet top boundary condition
(boundary_conds)

dirichlet_0 conditions moved from init_3d_model to
check_parameters (check_parameters, init_3d_model)

Errors:

---

bugfix: dirichlet_0 conditions moved from init_3d_model to
check_parameters (check_parameters, init_3d_model)

File:

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

palm/trunk/SOURCE/check_parameters.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! Calculating u,v-profiles from given profiles by linear interpolation.
+! bugfix: dirichlet_0 conditions for ug/vg moved from init_3d_model to here
 !
 ! Former revisions:
@@ -176 +177 @@
     CHARACTER (LEN=100) ::  action
 
-    INTEGER ::  i, ilen, intervals, iremote = 0, iter, j, k, nnxh, nnyh, &
-         position, prec
+    INTEGER ::  i, ilen, intervals, iremote = 0, iter, j, k, kk, nnxh, nnyh, &
+                position, prec
     LOGICAL ::  found, ldum
     REAL    ::  gradient, maxn, maxp, remote = 0.0, &
@@ -771 +772 @@
 
 !
-!--    Initial profiles for 1D and 3D model, respectively
-       u_init  = ug_surface
-       v_init  = vg_surface
+!--    Initial profiles for 1D and 3D model, respectively (u,v further below)
        pt_init = pt_surface
        IF ( humidity )        q_init  = q_surface
@@ -838 +837 @@
        ENDIF
 
-       u_init = ug
-
-!
-!--    In case of no given gradients for ug, choose a vanishing gradient
+!
+!--    In case of no given gradients for ug, choose a zero gradient
        IF ( ug_vertical_gradient_level(1) == -9999999.9 )  THEN
           ug_vertical_gradient_level(1) = 0.0
@@ -904 +901 @@
        ENDIF
 
-       v_init = vg
- 
-!
-!--    In case of no given gradients for vg, choose a vanishing gradient
+!
+!--    In case of no given gradients for vg, choose a zero gradient
        IF ( vg_vertical_gradient_level(1) == -9999999.9 )  THEN
           vg_vertical_gradient_level(1) = 0.0
+       ENDIF
+
+!
+!--    Let the initial wind profiles be the calculated ug/vg profiles or
+!--    interpolate them from wind profile data (if given)
+       IF ( u_profile(1) == 9999999.9  .AND.  v_profile(1) == 9999999.9 )  THEN
+
+          u_init = ug
+          v_init = vg
+
+       ELSEIF ( u_profile(1) == 0.0  .AND.  v_profile(1) == 0.0 )  THEN
+
+          IF ( uv_heights(1) /= 0.0 )  THEN
+             message_string = 'uv_heights(1) must be 0.0'
+             CALL message( 'check_parameters', 'PA0345', 1, 2, 0, 6, 0 )
+          ENDIF
+
+          use_prescribed_profile_data = .TRUE.
+
+          kk = 1
+          u_init(0) = 0.0
+          v_init(0) = 0.0
+
+          DO  k = 1, nz+1
+
+             IF ( kk < 100 )  THEN
+                DO WHILE ( uv_heights(kk+1) <= zu(k) )
+                   kk = kk + 1
+                   IF ( kk == 100 )  EXIT
+                ENDDO
+             ENDIF
+
+             IF ( kk < 100 )  THEN
+                u_init(k) = u_profile(kk) + ( zu(k) - uv_heights(kk) ) /       &
+                                       ( uv_heights(kk+1) - uv_heights(kk) ) * &
+                                       ( u_profile(kk+1) - u_profile(kk) )
+                v_init(k) = v_profile(kk) + ( zu(k) - uv_heights(kk) ) /       &
+                                       ( uv_heights(kk+1) - uv_heights(kk) ) * &
+                                       ( v_profile(kk+1) - v_profile(kk) )
+             ELSE
+                u_init(k) = u_profile(kk)
+                v_init(k) = v_profile(kk)
+             ENDIF
+
+          ENDDO
+
+       ELSE
+
+          message_string = 'u_profile(1) and v_profile(1) must be 0.0'
+          CALL message( 'check_parameters', 'PA0346', 1, 2, 0, 6, 0 )
+
        ENDIF
 
@@ -1539 +1585 @@
        IF ( bc_uv_t == 'dirichlet' .OR. bc_uv_t == 'dirichlet_0' )  THEN
           ibc_uv_t = 0
+          IF ( bc_uv_t == 'dirichlet_0' )  THEN
+!
+!--          Velocities for the initial u,v-profiles are set zero at the top
+!--          in case of dirichlet_0 conditions
+             u_init(nzt+1)    = 0.0
+             v_init(nzt+1)    = 0.0
+          ENDIF
        ELSEIF ( bc_uv_t == 'neumann' )  THEN
           ibc_uv_t = 1