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

Timestamp:

May 27, 2016 2:35:57 PM (8 years ago)

Author:

raasch

Message:

bugfixes for calculating run control quantities, bugfix for calculating pressure with fft-method in case of Neumann conditions both at bottom and top, steering of pres modified, ocean mode now uses initial density profile as reference in the buoyancy term

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -19 +19 @@
 ! Current revisions:
 ! ------------------
-! 
+! intermediate_timestep_count is set 0 instead 1 for first call of pres,
+! bugfix: initialization of local sum arrays are moved to the beginning of the
+!         routine because otherwise results from pres are overwritten
 ! 
 ! Former revisions:
@@ -709 +711 @@
     weight_substep = 1.0_wp
     weight_pres    = 1.0_wp
-    intermediate_timestep_count = 1  ! needed when simulated_time = 0.0
+    intermediate_timestep_count = 0  ! needed when simulated_time = 0.0
        
     CALL location_message( 'finished', .TRUE. )
+
+!
+!-- Initialize local summation arrays for routine flow_statistics.
+!-- This is necessary because they may not yet have been initialized when they
+!-- are called from flow_statistics (or - depending on the chosen model run -
+!-- are never initialized)
+    sums_divnew_l      = 0.0_wp
+    sums_divold_l      = 0.0_wp
+    sums_l_l           = 0.0_wp
+    sums_up_fraction_l = 0.0_wp
+    sums_wsts_bc_l     = 0.0_wp
+
+
 !
 !-- Initialize model variables
@@ -1821 +1836 @@
        ENDDO
     ENDIF
-
-!
-!-- Initialize local summation arrays for routine flow_statistics.
-!-- This is necessary because they may not yet have been initialized when they
-!-- are called from flow_statistics (or - depending on the chosen model run -
-!-- are never initialized)
-    sums_divnew_l      = 0.0_wp
-    sums_divold_l      = 0.0_wp
-    sums_l_l           = 0.0_wp
-    sums_up_fraction_l = 0.0_wp
-    sums_wsts_bc_l     = 0.0_wp
 
 !