Changeset 759 for palm/trunk/SOURCE/init_grid.f90

Timestamp:

Sep 15, 2011 1:58:31 PM (14 years ago)

Author:

raasch

Message:

New:
---

The number of parallel I/O operations can be limited with new mrun-option -w.
(advec_particles, data_output_2d, data_output_3d, header, init_grid, init_pegrid, init_3d_model, modules, palm, parin, write_3d_binary)

Changed:

---

mrun option -T is obligatory

Errors:

---

Bugfix: No zero assignments to volume_flow_initial and volume_flow_area in
case of normal restart runs. (init_3d_model)

initialization of u_0, v_0. This is just to avoid access of uninitialized
memory in exchange_horiz_2d, which causes respective error messages
when the Intel thread checker (inspector) is used. (production_e)

Bugfix for ts limitation (prandtl_fluxes)

File:

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

palm/trunk/SOURCE/init_grid.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! Splitting of parallel I/O in blocks of PEs
 !
 ! Former revisions:
@@ -80 +80 @@
 
     INTEGER ::  bh, blx, bly, bxl, bxr, byn, bys, ch, cwx, cwy, cxl, cxr, cyn, &
-                cys, gls, i, inc, i_center, j, j_center, k, l, nxl_l, nxr_l, &
-                nyn_l, nys_l, nzb_si, nzt_l, vi
+                cys, gls, i, ii, inc, i_center, j, j_center, k, l, nxl_l,      &
+                nxr_l, nyn_l, nys_l, nzb_si, nzt_l, vi
 
     INTEGER, DIMENSION(:), ALLOCATABLE   ::  vertical_influence
@@ -507 +507 @@
 
        CASE ( 'read_from_file' )
-!
-!--       Arbitrary irregular topography data in PALM format (exactly matching
-!--       the grid size and total domain size)
-          OPEN( 90, FILE='TOPOGRAPHY_DATA', STATUS='OLD', FORM='FORMATTED',  &
-               ERR=10 )
-          DO  j = ny, 0, -1
-             READ( 90, *, ERR=11, END=11 )  ( topo_height(j,i), i = 0, nx )
-          ENDDO
+
+          DO  ii = 0, io_blocks-1
+             IF ( ii == io_group )  THEN
+
+!
+!--             Arbitrary irregular topography data in PALM format (exactly
+!--             matching the grid size and total domain size)
+                OPEN( 90, FILE='TOPOGRAPHY_DATA', STATUS='OLD', &
+                      FORM='FORMATTED', ERR=10 )
+                DO  j = ny, 0, -1
+                   READ( 90, *, ERR=11, END=11 )  ( topo_height(j,i), i = 0,nx )
+                ENDDO
+
+                GOTO 12
+          
+ 10             message_string = 'file TOPOGRAPHY_DATA does not exist'
+                CALL message( 'init_grid', 'PA0208', 1, 2, 0, 6, 0 )
+
+ 11             message_string = 'errors in file TOPOGRAPHY_DATA'
+                CALL message( 'init_grid', 'PA0209', 1, 2, 0, 6, 0 )
+
+ 12             CLOSE( 90 )
+
+             ENDIF
+#if defined( __parallel )
+             CALL MPI_BARRIER( comm2d, ierr )
+#endif
+          ENDDO
+
 !
 !--       Calculate the index height of the topography
@@ -523 +544 @@
           ENDDO
 !
-!--       Add cyclic boundaries (additional layers are for calculating flag
-!--       arrays needed for the multigrid sover)
+!--       Add cyclic boundaries (additional layers are for calculating
+!--       flag arrays needed for the multigrid sover)
           nzb_local(-gls:-1,0:nx)     = nzb_local(ny-gls+1:ny,0:nx)
           nzb_local(ny+1:ny+gls,0:nx) = nzb_local(0:gls-1,0:nx)
           nzb_local(:,-gls:-1)        = nzb_local(:,nx-gls+1:nx)
           nzb_local(:,nx+1:nx+gls)    = nzb_local(:,0:gls-1)
-
-
-      
-          GOTO 12
-          
- 10       message_string = 'file TOPOGRAPHY_DATA does not exist'
-          CALL message( 'init_grid', 'PA0208', 1, 2, 0, 6, 0 )
-
- 11       message_string = 'errors in file TOPOGRAPHY_DATA'
-          CALL message( 'init_grid', 'PA0209', 1, 2, 0, 6, 0 )
-
- 12       CLOSE( 90 )
 
        CASE DEFAULT