Changeset 107 for palm/trunk/SOURCE

Timestamp:

Aug 17, 2007 1:54:45 PM (17 years ago)

Author:

raasch

Message:

further bugfix for non-cyclic BCs

Location:

palm/trunk/SOURCE

Files:

• CURRENT_MODIFICATIONS (modified) (1 diff)
• boundary_conds.f90 (modified) (2 diffs)
• init_rankine.f90 (modified) (2 diffs)

palm/trunk/SOURCE/CURRENT_MODIFICATIONS

@@ -17 +17 @@
 Quantities w*p* and w"e can be output as vertical profiles.
 
-boundary_conds, check_open, check_parameters, diffusion_u, diffusion_v, flow_statistics, header, init_pegrid, init_3d_model, modules, palm, parin, pres, prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, time_integration, write_var_list, write_3d_binary
+Initial profiles are reset to constant profiles in case that initializing_actions /= 'set_constant_profiles'. (init_rankine)
+
+boundary_conds, check_open, check_parameters, diffusion_u, diffusion_v, flow_statistics, header, init_pegrid, init_rankine, init_3d_model, modules, palm, parin, pres, prognostic_equations, read_var_list, read_3d_binary, swap_timelevel, time_integration, write_var_list, write_3d_binary
 
 New:

palm/trunk/SOURCE/boundary_conds.f90

@@ -285 +285 @@
                                        ( u(k,-1,i) - u(k,0,i) ) * ddy
 
-             v_p(k,-1,i) = v(k,-1,i) - dt_3d * tsc(2) * c_v(k,i) * &
+             v_p(k,0,i)  = v(k,0,i)  - dt_3d * tsc(2) * c_v(k,i) * &
                                        ( v(k,0,i) - v(k,1,i) ) * ddy
 
@@ -433 +433 @@
                 IF ( denom /= 0.0 )  THEN
                    c_u(k,j) = -c_max * ( u(k,j,1) - u_m_l(k,j,1) ) / denom
-                   IF ( c_u(k,j) > 0.0 )  THEN
+                   IF ( c_u(k,j) < 0.0 )  THEN
                       c_u(k,j) = 0.0
-                   ELSEIF ( c_u(k,j) < -c_max )  THEN
-                      c_u(k,j) = -c_max
-                   ENDIF
-                ELSE
-                   c_u(k,j) = -c_max
+                   ELSEIF ( c_u(k,j) > c_max )  THEN
+                      c_u(k,j) = c_max
+                   ENDIF
+                ELSE
+                   c_u(k,j) = c_max
                 ENDIF
 

palm/trunk/SOURCE/init_rankine.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
-! 
+! Initial profiles are reset to constant profiles
 !
 ! Former revisions:
@@ -48 +48 @@
     kc1 = nzb
     kc2 = nzt+1
+
+!
+!-- Reset initial profiles to constant profiles
+    IF ( INDEX(initializing_actions, 'set_constant_profiles') /= 0 )  THEN
+       DO  i = nxl-1, nxr+1
+          DO  j = nys-1, nyn+1
+             pt(:,j,i) = pt_init
+             u(:,j,i)  = u_init
+             v(:,j,i)  = v_init
+          ENDDO
+       ENDDO
+    ENDIF
 
 !