Changeset 4819

Timestamp:

Dec 10, 2020 12:57:52 PM (4 years ago)

Author:

gronemeier

Message:

Check if netCDF file actually exists before opening it for reading (netcdf_data_input_mod).

File:

• palm/trunk/SOURCE/netcdf_data_input_mod.f90 (modified) (2 diffs)

palm/trunk/SOURCE/netcdf_data_input_mod.f90

@@ -24 +24 @@
 ! -----------------
 ! $Id$
+! Check if netCDF file actually exists before opening it for reading.
+!
+! 4807 2020-12-02 21:02:28Z gronemeier
 ! Deactivated reading of building_obstruction_full due to improper implementation (conflicts with
 ! building_obstruction_f and it is also not used anywhere else in the code).
@@ -3589 +3592 @@
  SUBROUTINE open_read_file( filename, id )
 
+    USE control_parameters,                                                                        &
+        ONLY: message_string
+
     USE pegrid
 
     IMPLICIT NONE
 
-    CHARACTER (LEN=*), INTENT(IN) ::  filename  !< filename
-    INTEGER(iwp), INTENT(INOUT)   ::  id        !< file id
+    CHARACTER (LEN=*), INTENT(IN) ::  filename    !< filename
+    INTEGER(iwp), INTENT(INOUT)   ::  id          !< file id
+    LOGICAL                       ::  file_exists !< true if file exists
 
 #if defined( __netcdf )
+!
+!-- Check if requested file exists
+    INQUIRE( FILE=filename, EXIST=file_exists )
+
+    IF ( .NOT. file_exists )  THEN
+       WRITE( message_string, * ) 'Required input file "' // filename // '" does not exist!'
+       CALL message( 'open_read_file', 'PA0516', 2, 2, 0, 6, 1 )
+    ENDIF
 
 #if defined( __netcdf4_parallel )