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

Timestamp:

Mar 7, 2007 8:38:00 AM (17 years ago)

Author:

raasch

Message:

preliminary version, several changes to be explained later

File:

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

palm/trunk/SOURCE/init_3d_model.f90

@@ -8 +8 @@
 ! -----------------
 ! New initializing action "by_user" calls user_init_3d_model,
-! ts_value is allocated
+! ts_value is allocated, +module netcdf_control,
+! initial velocities at nzb+1 are regarded for volume
+! flow control in case they have been set zero before (to avoid small timesteps)
 !
 ! Former revisions:
@@ -44 +46 @@
     USE interfaces
     USE model_1d
+    USE netcdf_control
     USE particle_attributes
     USE pegrid
@@ -183 +186 @@
 
 !
+!-- 4D-array for storing the Rif-values at vertical walls
+    IF ( topography /= 'flat' )  THEN
+       ALLOCATE( rif_wall(nzb:nzt+1,nys-1:nyn+1,nxl-1:nxr+1,1:4) )
+       rif_wall = 0.0
+    ENDIF
+
+!
+!-- Velocities at nzb+1 needed for volume flow control
+    IF ( conserve_volume_flow )  THEN
+       ALLOCATE( u_nzb_p1_for_vfc(nys:nyn), v_nzb_p1_for_vfc(nxl:nxr) )
+       u_nzb_p1_for_vfc = 0.0
+       v_nzb_p1_for_vfc = 0.0
+    ENDIF
+
+!
 !-- Initial assignment of the pointers
     IF ( timestep_scheme(1:5) /= 'runge' )  THEN
@@ -382 +400 @@
           ENDDO
 !
-!--       Set initial horizontal velocities at the lowest grid levels to zero
-!--       in order to avoid too small time steps caused by the diffusion
+!--       Set initial horizontal velocities at the lowest computational grid levels
+!--       to zero in order to avoid too small time steps caused by the diffusion
 !--       limit in the initial phase of a run (at k=1, dz/2 occurs in the
-!--       limiting formula!)
+!--       limiting formula!). The original values are stored to be later used for
+!--       volume flow control.
           DO  i = nxl-1, nxr+1
              DO  j = nys-1, nyn+1
@@ -392 +411 @@
              ENDDO
           ENDDO
+          IF ( conserve_volume_flow )  THEN
+             IF ( nxr == nx )  THEN
+                DO  j = nys, nyn
+                   k = nzb_u_inner(j,nx) + 1
+                   u_nzb_p1_for_vfc(j) = u_init(k) * dzu(k)
+                ENDDO
+             ENDIF
+             IF ( nyn == ny )  THEN
+                DO  i = nxl, nxr
+                   k = nzb_v_inner(ny,i) + 1
+                   v_nzb_p1_for_vfc(i) = v_init(k) * dzu(k)
+                ENDDO
+             ENDIF
+          ENDIF
 
           IF ( moisture  .OR.  passive_scalar )  THEN
@@ -568 +601 @@
                    volume_flow_area_l(1)    = volume_flow_area_l(1) + dzu(k)
                 ENDDO
+!
+!--             Correction if velocity at nzb+1 has been set zero further above
+                volume_flow_initial_l(1) = volume_flow_initial_l(1) + &
+                                           u_nzb_p1_for_vfc(j)
              ENDDO
           ENDIF
@@ -578 +615 @@
                    volume_flow_area_l(2)    = volume_flow_area_l(2) + dzu(k)
                 ENDDO
+!
+!--             Correction if velocity at nzb+1 has been set zero further above
+                volume_flow_initial_l(2) = volume_flow_initial_l(2) + &
+                                           v_nzb_p1_for_vfc(i)
              ENDDO
           ENDIF
@@ -704 +745 @@
 !
 !-- If required, initialize dvrp-software
-!    WRITE ( 9, * ) '*** myid=', myid, ' vor init_dvrp'
-!    CALL FLUSH_( 9 )
     IF ( dt_dvrp /= 9999999.9 )  CALL init_dvrp
-!    WRITE ( 9, * ) '*** myid=', myid, ' nach init_dvrp'
-!    CALL FLUSH_( 9 )
 
 !
@@ -719 +756 @@
 !
 !-- If required, initialize particles
-!    WRITE ( 9, * ) '*** myid=', myid, ' vor init_particles'
-!    CALL FLUSH_( 9 )
     CALL init_particles
-!    WRITE ( 9, * ) '*** myid=', myid, ' nach init_particles'
-!    CALL FLUSH_( 9 )
 
 !
@@ -825 +858 @@
 
 !
-!-- User-defined initializing actions
+!-- User-defined initializing actions. Check afterwards, if maximum number
+!-- of allowed timeseries is not exceeded
     CALL user_init
+
+    IF ( dots_num > dots_max )  THEN
+       IF ( myid == 0 )  THEN
+          PRINT*, '+++ user_init: number of time series quantities exceeds', &
+                  ' its maximum of dots_max = ', dots_max
+          PRINT*, '    Please increase dots_max in modules.f90.'
+       ENDIF
+       CALL local_stop
+    ENDIF
 
 !