Changeset 6058 for IOIPSL

Timestamp:

02/10/22 13:44:36 (2 years ago)

Author:

jgipsl

Message:

Add error checking after read of attributes related to time axis. Do nothing if there are errors. This is needed for the file ECDYN.nc without time axis, file read by create_etat0_limit in LMDZ. Note that it is not new that this file is read by flincom but maybe the debug options was not tested for this case before.

Done by Adriana Sima, LMD

Location:

IOIPSL

Files:

• tags/v2_2_5/src/flincom.f90 (modified) (2 diffs)
• trunk/src/flincom.f90 (modified) (2 diffs)

IOIPSL/tags/v2_2_5/src/flincom.f90

@@ -377 +377 @@
   REAL, PARAMETER :: eps = 1.e-4
 !-
-  INTEGER :: iret, vid, fid, nbdim, i, iilast, jjlast
+  INTEGER :: iret, nberr, vid, fid, nbdim, i, iilast, jjlast
   INTEGER :: gdtt_id, old_id, iv, gdtmaf_id
   CHARACTER(LEN=250) :: name
@@ -653 +653 @@
     ELSE IF (old_id > 0) THEN
       iret = NF90_GET_ATT (fid, NF90_GLOBAL, 'delta_tstep_sec', dt)
+      nberr = 0 
       iret = NF90_GET_ATT (fid, NF90_GLOBAL, 'day0', r_day)
+      IF (iret /= NF90_NOERR) nberr=nberr+1
       iret = NF90_GET_ATT (fid, NF90_GLOBAL, 'sec0', sec)
+      IF (iret /= NF90_NOERR) nberr=nberr+1 
       iret = NF90_GET_ATT (fid, NF90_GLOBAL, 'year0', r_year)
+      IF (iret /= NF90_NOERR) nberr=nberr+1
       iret = NF90_GET_ATT (fid, NF90_GLOBAL, 'month0', r_month)
-!-----
-      day = INT(r_day)
-      month = INT(r_month)
-      year = INT(r_year)
-!-----
-      CALL ymds2ju (year, month, day, sec, date0)
+      IF (iret /= NF90_NOERR) nberr=nberr+1
+
+      IF (nberr == 0) THEN
+         day = INT(r_day)
+         month = INT(r_month)
+         year = INT(r_year)
+
+         CALL ymds2ju (year, month, day, sec, date0)
+      ELSE
+          WRITE(*,*) 'Subroutine : flinopen_work ; time axis missing or not recognized in file', filename
+      ENDIF
+
     ENDIF
   ENDIF

IOIPSL/trunk/src/flincom.f90

@@ -377 +377 @@
   REAL, PARAMETER :: eps = 1.e-4
 !-
-  INTEGER :: iret, vid, fid, nbdim, i, iilast, jjlast
+  INTEGER :: iret, nberr, vid, fid, nbdim, i, iilast, jjlast
   INTEGER :: gdtt_id, old_id, iv, gdtmaf_id
   CHARACTER(LEN=250) :: name
@@ -653 +653 @@
     ELSE IF (old_id > 0) THEN
       iret = NF90_GET_ATT (fid, NF90_GLOBAL, 'delta_tstep_sec', dt)
+      nberr = 0 
       iret = NF90_GET_ATT (fid, NF90_GLOBAL, 'day0', r_day)
+      IF (iret /= NF90_NOERR) nberr=nberr+1
       iret = NF90_GET_ATT (fid, NF90_GLOBAL, 'sec0', sec)
+      IF (iret /= NF90_NOERR) nberr=nberr+1 
       iret = NF90_GET_ATT (fid, NF90_GLOBAL, 'year0', r_year)
+      IF (iret /= NF90_NOERR) nberr=nberr+1
       iret = NF90_GET_ATT (fid, NF90_GLOBAL, 'month0', r_month)
-!-----
-      day = INT(r_day)
-      month = INT(r_month)
-      year = INT(r_year)
-!-----
-      CALL ymds2ju (year, month, day, sec, date0)
+      IF (iret /= NF90_NOERR) nberr=nberr+1
+
+      IF (nberr == 0) THEN
+         day = INT(r_day)
+         month = INT(r_month)
+         year = INT(r_year)
+
+         CALL ymds2ju (year, month, day, sec, date0)
+      ELSE
+          WRITE(*,*) 'Subroutine : flinopen_work ; time axis missing or not recognized in file', filename
+      ENDIF
+
     ENDIF
   ENDIF