Changeset 4099 for palm/trunk/SOURCE

Timestamp:

Jul 15, 2019 3:29:37 PM (5 years ago)

Author:

suehring

Message:

Bugfix in opening the parameter file (unit 11) in case of ocean precursor runs

File:

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

palm/trunk/SOURCE/check_open.f90

@@ -25 +25 @@
 ! -----------------
 ! $Id$
+! Bugfix in opening the parameter file (unit 11) in case of ocean precursor
+! runs. 
+! 
+! 4069 2019-07-01 14:05:51Z Giersch
 ! Masked output running index mid has been introduced as a local variable to 
 ! avoid runtime error (Loop variable has been modified) in time_integration
@@ -233 +237 @@
     INTEGER(iwp) ::  mid         !< masked output running index
     
+    LOGICAL ::  file_exist       !< file check
     LOGICAL ::  netcdf_extend    !<
 
@@ -311 +316 @@
 
        CASE ( 11 )
-
-          OPEN ( 11, FILE='PARIN'//TRIM( coupling_char ), FORM='FORMATTED',    &
-                     STATUS='OLD' )
+!
+!--       Read the parameter file. Therefore, inquire whether the file exist or 
+!--       not. This is required for the ocean-atmoshere coupling. For an ocean
+!--       precursor run palmrun provides a PARIN_O file instead of a PARIN 
+!--       file. Actually this should be considered in coupling_char, however,
+!--       in pmc_init the parameter file is already opened to read the 
+!--       nesting parameters and decide whether it is a nested run or not, 
+!--       but coupling_char is still not set at that moment (must be set after
+!-        the nesting setup is read).
+!--       This, however, leads to the situation that for ocean
+!--       precursor runs PARIN is not available and the run crashes. Thus, 
+!--       if the file is not there, PARIN_O will be read. An ocean precursor 
+!--       run will be the only situation where this can happen. 
+          INQUIRE( FILE = 'PARIN' // TRIM( coupling_char ),                    &
+                   EXIST = file_exist )
+                   
+          IF ( file_exist )  THEN
+             filename = 'PARIN' // TRIM( coupling_char )
+          ELSE
+             filename = 'PARIN_O'
+          ENDIF
+
+          OPEN ( 11, FILE= TRIM( filename ), FORM='FORMATTED', STATUS='OLD' )
 
        CASE ( 13 )