Changeset 240 for palm/trunk/SOURCE/init_3d_model.f90

Timestamp:

Feb 18, 2009 5:50:38 PM (15 years ago)

Author:

letzel

Message:

• External pressure gradient (check_parameters, init_3d_model, header, modules, parin, prognostic_equations)
• New topography case 'single_street_canyon'

File:

• palm/trunk/SOURCE/init_3d_model.f90 (modified) (4 diffs)

palm/trunk/SOURCE/init_3d_model.f90

@@ -7 +7 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Set the starting level and the vertical smoothing factor used for 
+! the external pressure gradient
 !
 ! Former revisions:
@@ -91 +92 @@
     IMPLICIT NONE
 
-    INTEGER ::  i, j, k, sr
+    INTEGER ::  i, ind_array(1), j, k, sr
 
     INTEGER, DIMENSION(:), ALLOCATABLE ::  ngp_2dh_l, ngp_3d_inner_l
@@ -111 +112 @@
               sums_divnew_l(0:statistic_regions),                           &
               sums_divold_l(0:statistic_regions) )
-    ALLOCATE( rdf(nzb+1:nzt) )
+    ALLOCATE( dp_smooth_factor(nzb:nzt), rdf(nzb+1:nzt) )
     ALLOCATE( ngp_2dh_outer(nzb:nzt+1,0:statistic_regions),                 &
               ngp_2dh_outer_l(nzb:nzt+1,0:statistic_regions),               &
@@ -1264 +1265 @@
 
 !
+!-- Initialize the starting level and the vertical smoothing factor used for 
+!-- the external pressure gradient
+    dp_smooth_factor = 1.0
+    IF ( dp_external )  THEN
+!
+!--    Set the starting level dp_level_ind_b only if it has not been set before
+!--    (e.g. in init_grid).
+       IF ( dp_level_ind_b == 0 )  THEN
+          ind_array = MINLOC( ABS( dp_level_b - zu ) )
+          dp_level_ind_b = ind_array(1) - 1 + nzb 
+                                        ! MINLOC uses lower array bound 1
+       ENDIF
+       IF ( dp_smooth )  THEN
+          dp_smooth_factor(:dp_level_ind_b) = 0.0
+          DO  k = dp_level_ind_b+1, nzt
+             dp_smooth_factor(k) = 0.5 * ( 1.0 + SIN( pi * &
+                  ( REAL( k - dp_level_ind_b ) /  &
+                    REAL( nzt - dp_level_ind_b ) - 0.5 ) ) )
+          ENDDO
+       ENDIF
+    ENDIF
+
+!
 !-- Initialize diffusivities used within the outflow damping layer in case of
 !-- non-cyclic lateral boundaries. A linear increase is assumed over the first