Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/data_output_2d.f90

    r730 r759  
    44! Current revisions:
    55! -----------------
     6! Splitting of parallel I/O
    67!
    78! Former revisions:
     
    737738                         ENDIF
    738739#endif
    739                          WRITE ( 21 )  nxlg, nxrg, nysg, nyng
    740                          WRITE ( 21 )  local_2d
     740                         DO  i = 0, io_blocks-1
     741                            IF ( i == io_group )  THEN
     742                               WRITE ( 21 )  nxlg, nxrg, nysg, nyng
     743                               WRITE ( 21 )  local_2d
     744                            ENDIF
     745#if defined( __parallel )
     746                            CALL MPI_BARRIER( comm2d, ierr )
     747#endif
     748                         ENDDO
    741749
    742750                      ELSE
     
    10331041                         ENDIF
    10341042#endif
    1035                          IF ( ( section(is,s) >= nys  .AND.                  &
    1036                                 section(is,s) <= nyn )  .OR.                 &
    1037                               ( section(is,s) == -1  .AND.  nys-1 == -1 ) )  &
    1038                          THEN
    1039                             WRITE (22)  nxlg, nxrg, nzb, nzt+1
    1040                             WRITE (22)  local_2d
    1041                          ELSE
    1042                             WRITE (22)  -1, -1, -1, -1
    1043                          ENDIF
     1043                         DO  i = 0, io_blocks-1
     1044                            IF ( i == io_group )  THEN
     1045                               IF ( ( section(is,s) >= nys  .AND.   &
     1046                                      section(is,s) <= nyn )  .OR.  &
     1047                                    ( section(is,s) == -1  .AND.    &
     1048                                      nys-1 == -1 ) )               &
     1049                               THEN
     1050                                  WRITE (22)  nxlg, nxrg, nzb, nzt+1
     1051                                  WRITE (22)  local_2d
     1052                               ELSE
     1053                                  WRITE (22)  -1, -1, -1, -1
     1054                               ENDIF
     1055                            ENDIF
     1056#if defined( __parallel )
     1057                            CALL MPI_BARRIER( comm2d, ierr )
     1058#endif
     1059                         ENDDO
    10441060
    10451061                      ELSE
     
    13341350                         ENDIF
    13351351#endif
    1336                          IF ( ( section(is,s) >= nxl  .AND.                  &
    1337                                 section(is,s) <= nxr )  .OR.                 &
    1338                               ( section(is,s) == -1  .AND.  nxl-1 == -1 ) )  &
    1339                          THEN
    1340                             WRITE (23)  nysg, nyng, nzb, nzt+1
    1341                             WRITE (23)  local_2d
    1342                          ELSE
    1343                             WRITE (23)  -1, -1, -1, -1
    1344                          ENDIF
     1352                         DO  i = 0, io_blocks-1
     1353                            IF ( i == io_group )  THEN
     1354                               IF ( ( section(is,s) >= nxl  .AND.   &
     1355                                      section(is,s) <= nxr )  .OR.  &
     1356                                    ( section(is,s) == -1  .AND.    &
     1357                                      nxl-1 == -1 ) )               &
     1358                               THEN
     1359                                  WRITE (23)  nysg, nyng, nzb, nzt+1
     1360                                  WRITE (23)  local_2d
     1361                               ELSE
     1362                                  WRITE (23)  -1, -1, -1, -1
     1363                               ENDIF
     1364                            ENDIF
     1365#if defined( __parallel )
     1366                            CALL MPI_BARRIER( comm2d, ierr )
     1367#endif
     1368                         ENDDO
    13451369
    13461370                      ELSE
     
    15041528
    15051529    IF ( data_output_2d_on_each_pe )  THEN
    1506        CALL close_file( file_id )
     1530       DO  i = 0, io_blocks-1
     1531          IF ( i == io_group )  THEN
     1532             CALL close_file( file_id )
     1533          ENDIF
     1534#if defined( __parallel )
     1535          CALL MPI_BARRIER( comm2d, ierr )
     1536#endif
     1537       ENDDO
    15071538    ELSE
    15081539       IF ( myid == 0 )  CALL close_file( file_id )
Note: See TracChangeset for help on using the changeset viewer.