Changeset 623 for trunk/NEMO/OPA_SRC/SBC/flxfwb.F90

Timestamp:

2007-02-26T18:26:08+01:00 (17 years ago)

Author:

opalod

Message:

nemo_v2_bugfix_024:RB: Use of ctlopn for all files (except for dimg and coupled part)

File:

• trunk/NEMO/OPA_SRC/SBC/flxfwb.F90 (modified) (4 diffs)

trunk/NEMO/OPA_SRC/SBC/flxfwb.F90

@@ -80 +80 @@
 
       !! * Local declarations
+      CHARACTER (len=32) ::   &
+         clname = 'EMPav.dat'
       INTEGER  ::   ji, jj, jk        ! dummy loop indices
-      INTEGER  ::   inum = 11         ! temporary logical unit
+      INTEGER  ::   inum              ! temporary logical unit
       INTEGER  ::   ikty              ! 
       REAL(wp) ::   &
@@ -194 +196 @@
       IF( kt == nitend ) THEN
 
-         OPEN( inum, FILE='EMPave.dat' )
+         CALL ctlopn( inum, clname, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL',   &
+            &           1, numout, .FALSE., 0 )
          WRITE(inum, "(24X,I8,2ES24.16)" ) nyear, a_fwb_b, a_fwb
          WRITE(inum,*)
@@ -249 +252 @@
       !!----------------------------------------------------------------------
       !! * Local declarations
-      LOGICAL ::   llbon
       CHARACTER (len=32) ::   &
          clname = 'EMPave_old.dat'
-      INTEGER ::   inum = 11         ! temporary logical unit
-      INTEGER ::   iyear
+      INTEGER ::   inum, iyear
 
       NAMELIST/namfwb/ ln_fwb 
@@ -290 +291 @@
          !                                     ! ==============================
          ! Read the corrective factor on precipitations (empold)
-         INQUIRE( FILE=clname, EXIST=llbon )
-         IF( llbon ) THEN
-            OPEN ( inum, FILE=clname)
-            READ ( inum, "(24X,I8,2ES24.16)" ) iyear, a_fwb_b, a_fwb
-            CLOSE( inum )
-            empold = a_fwb                  ! current year freshwater budget correction
-            !                               ! estimate from the previous year budget
-            IF(lwp)WRITE(numout,*)
-            IF(lwp)WRITE(numout,*)'flx_fwb_init : year = ',iyear  , ' freshwater budget correction = ', empold
-            IF(lwp)WRITE(numout,*)'               year = ',iyear-1, ' freshwater budget read       = ', a_fwb
-            IF(lwp)WRITE(numout,*)'               year = ',iyear-2, ' freshwater budget read       = ', a_fwb_b
-         ELSE
-            WRITE(ctmp1,*)'flx_fwb_init : unable to read the file', clname
-            CALL ctl_stop( ctmp1 )
-         ENDIF
+         CALL ctlopn( inum, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL',   &
+            &           1, numout, .FALSE., 1 )
+         READ ( inum, "(24X,I8,2ES24.16)" ) iyear, a_fwb_b, a_fwb
+         CLOSE( inum )
+         empold = a_fwb                  ! current year freshwater budget correction
+         !                               ! estimate from the previous year budget
+         IF(lwp)WRITE(numout,*)
+         IF(lwp)WRITE(numout,*)'flx_fwb_init : year = ',iyear  , ' freshwater budget correction = ', empold
+         IF(lwp)WRITE(numout,*)'               year = ',iyear-1, ' freshwater budget read       = ', a_fwb
+         IF(lwp)WRITE(numout,*)'               year = ',iyear-2, ' freshwater budget read       = ', a_fwb_b
          !                                    ! ==============================
       ELSE                                    !      NO  budget correction