Changeset 3159 for palm/trunk/SOURCE/check_open.f90

Timestamp:

Jul 20, 2018 11:20:01 AM (6 years ago)

Author:

sward

Message:

Added multi agent system

File:

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

palm/trunk/SOURCE/check_open.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Added multi agent system
+! 
+! 3045 2018-05-28 07:55:41Z Giersch
 ! Output case 108 disabled
 ! 
@@ -190 +193 @@
 
     USE netcdf_interface,                                                      &
-        ONLY:  id_set_fl, id_set_mask, id_set_pr, id_set_prt, id_set_pts,      &
-               id_set_sp, id_set_ts, id_set_xy, id_set_xz, id_set_yz,          &
-               id_set_3d, nc_stat, netcdf_create_file, netcdf_data_format,     &
-               netcdf_define_header, netcdf_handle_error, netcdf_open_write_file
+        ONLY:  id_set_agt, id_set_fl, id_set_mask, id_set_pr, id_set_prt,      &
+               id_set_pts, id_set_sp, id_set_ts, id_set_xy, id_set_xz,         &
+               id_set_yz, id_set_3d, nc_stat, netcdf_create_file,              &
+               netcdf_data_format, netcdf_define_header, netcdf_handle_error,  &
+               netcdf_open_write_file
 
     USE particle_attributes,                                                   &
@@ -1077 +1081 @@
           ENDIF
 
+       CASE ( 118 )
+
+          IF ( myid == 0 )  THEN
+             filename = 'DATA_AGT_NETCDF'
+!
+!--       Inquire, if there is a netCDF file from a previuos run. This should
+!--       be opened for extension, if its variables match the actual run.
+          INQUIRE( FILE=filename, EXIST=netcdf_extend )
+
+!
+!--          Create a new netCDF output file with requested netCDF format
+             CALL netcdf_create_file( filename, id_set_agt, .FALSE., 43 )
+
+!
+!--          Define the header
+             CALL netcdf_define_header( 'ag', netcdf_extend, 0 )
+          ENDIF
+
+!           IF ( netcdf_extend )  THEN
+! !
+! !--          Open an existing netCDF file for output
+!              CALL netcdf_open_write_file( filename, id_set_agt, .FALSE., 41 )
+! !
+! !--          Read header information and set all ids. If there is a mismatch
+! !--          between the previuos and the actual run, netcdf_extend is returned
+! !--          as .FALSE.
+!              CALL netcdf_define_header( 'ag', netcdf_extend, 0 )
+! 
+! !
+! !--          Remove the local file, if it can not be extended
+!              IF ( .NOT. netcdf_extend )  THEN
+!                 nc_stat = NF90_CLOSE( id_set_agt )
+!                 CALL netcdf_handle_error( 'check_open', 42 )
+!                 CALL local_system( 'rm ' // TRIM( filename ) )
+!              ENDIF
+! 
+!           ENDIF
+
+          IF ( .NOT. netcdf_extend )  THEN
+
+!
+! !--          For runs on multiple processors create the subdirectory
+!              IF ( myid_char /= '' )  THEN
+!                 IF ( myid == 0  .AND. .NOT. openfile(file_id)%opened_before )  &
+!                 THEN    ! needs modification in case of non-extendable sets
+!                    CALL local_system( 'mkdir  DATA_PRT_NETCDF' //              &
+!                                        TRIM( coupling_char ) // '/' )
+!                 ENDIF
+! #if defined( __parallel )
+! ! 
+! !--             Set a barrier in order to allow that all other processors in the
+! !--             directory created by PE0 can open their file
+!                 CALL MPI_BARRIER( comm2d, ierr )
+! #endif
+!              ENDIF
+
+          ENDIF
+
 !
 !--    nc-file for virtual flight measurements