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

Timestamp:

Sep 15, 2011 1:58:31 PM (10 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/palm.f90 (modified) (3 diffs)

palm/trunk/SOURCE/palm.f90

@@ -4 +4 @@
 ! Current revisions:
 ! -----------------
-! 
+! Splitting of parallel I/O, cpu measurement for write_3d_binary and opening
+! of unit 14 moved to here
 !
 ! Former revisions:
@@ -175 +176 @@
 !-- If required, write binary data for restart runs
     IF ( write_binary(1:4) == 'true' )  THEN
-!
-!--    Write flow field data
-       CALL write_3d_binary
+
+       CALL cpu_log( log_point(22), 'write_3d_binary', 'start' )
+
+       CALL check_open( 14 )
+
+       DO  i = 0, io_blocks-1
+          IF ( i == io_group )  THEN
+!
+!--          Write flow field data
+             CALL write_3d_binary
+          ENDIF
+#if defined( __parallel )
+          CALL MPI_BARRIER( comm2d, ierr )
+#endif
+       ENDDO
+
+       CALL cpu_log( log_point(22), 'write_3d_binary', 'stop' )
+
 !
 !--    If required, write particle data
@@ -192 +208 @@
 !-- unit in routine user_last_actions.
     CALL cpu_log( log_point(4), 'last actions', 'start' )
-    CALL user_last_actions
-    IF ( write_binary(1:4) == 'true' )  CALL close_file( 14 )
+    DO  i = 0, io_blocks-1
+       IF ( i == io_group )  THEN
+          CALL user_last_actions
+          IF ( write_binary(1:4) == 'true' )  CALL close_file( 14 )
+       ENDIF
+#if defined( __parallel )
+       CALL MPI_BARRIER( comm2d, ierr )
+#endif
+    ENDDO
     CALL close_file( 0 )
     CALL close_dvrp