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

Timestamp:

Jul 7, 2016 10:43:48 AM (7 years ago)

Author:

suehring

Message:

flight module added

File:

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

palm/trunk/SOURCE/check_open.f90

@@ -19 +19 @@
 ! Current revisions:
 ! -----------------
-! 
+! flight module added
 ! 
 ! Former revisions:
@@ -128 +128 @@
 
     USE netcdf_interface,                                                      &
-        ONLY:  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_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,                                                   &
@@ -1108 +1108 @@
           OPEN ( 117, FILE='PROGRESS'//TRIM( coupling_char ),                  &
                       STATUS='REPLACE', FORM='FORMATTED' )
+!
+!--    nc-file for virtual flight measurements
+       CASE ( 199 )
+!
+!--       Set filename
+          filename = 'DATA_1D_FL_NETCDF' // TRIM( coupling_char )
+
+!
+!--       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 )
+
+          IF ( netcdf_extend )  THEN
+!
+!--          Open an existing netCDF file for output
+             CALL netcdf_open_write_file( filename, id_set_fl, .FALSE., 532 )
+!
+!--          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( 'fl', netcdf_extend, 0 )
+
+!
+!--          Remove the local file, if it can not be extended
+             IF ( .NOT. netcdf_extend )  THEN
+                nc_stat = NF90_CLOSE( id_set_fl )
+                CALL netcdf_handle_error( 'check_open', 533 )
+                CALL local_system( 'rm ' // TRIM( filename ) )
+             ENDIF
+
+          ENDIF          
+
+          IF ( .NOT. netcdf_extend )  THEN
+!
+!--          Create a new netCDF output file with requested netCDF format
+             CALL netcdf_create_file( filename, id_set_fl, .FALSE., 534 )
+!
+!--          Define the header
+             CALL netcdf_define_header( 'fl', netcdf_extend, 0 )
+
+          ENDIF