Changeset 2546 for trunk/LMDZ.MARS/util/localtime.F90

Timestamp:

Jul 8, 2021, 4:29:44 PM (3 years ago)

Author:

romain.vande

Message:

MARS Dynamico:

Utilitaire concatnc, localtime, zrecast can process Xios output
RV

File:

• trunk/LMDZ.MARS/util/localtime.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/util/localtime.F90

@@ -247 +247 @@
    ierr=NF_INQ_DIMID(nid,"latitude",latdim)
    if (ierr.NE.NF_NOERR) then
-      write(*,*) 'ERROR: Dimension <latitude> is missing in file '//trim(file)
-      stop ""  
+      write(*,*) 'Dimension <latitude> is missing, looking for lat'
+      ierr=NF_INQ_DIMID(nid,"lat",latdim)
+         if (ierr.NE.NF_NOERR) then
+            write(*,*) 'ERROR: Dimension <latitude> and <lat> is missing in file '//trim(file)
+            stop ""  
+         endif
    endif
    ierr=NF_INQ_VARID(nid,"latitude",latvar)
    if (ierr.NE.NF_NOERR) then
-      write(*,*) 'ERROR: Field <latitude> is missing in file '//trim(file)
-      stop ""  
+      write(*,*) 'Field <latitude> is missing, looking for lat'
+      ierr=NF_INQ_VARID(nid,"lat",latvar)
+      if (ierr.NE.NF_NOERR) then
+        write(*,*) 'ERROR: Field <latitude> and <lat> is missing in file '//trim(file)
+        stop ""  
+      endif
    endif
    ierr=NF_INQ_DIMLEN(nid,latdim,latlen)
@@ -260 +268 @@
    ierr=NF_INQ_DIMID(nid,"longitude",londim)
    if (ierr.NE.NF_NOERR) then
-      write(*,*) 'ERROR: Dimension <longitude> is missing in file '//trim(file)
-      stop ""  
+      write(*,*) 'Field <longitude> is missing, looking for lon'
+      ierr=NF_INQ_DIMID(nid,"lon",londim)
+         if (ierr.NE.NF_NOERR) then
+            write(*,*) 'ERROR: Dimension <longitude> and <lon> is missing in file '//trim(file)
+            stop ""  
+         endif 
    endif
    ierr=NF_INQ_VARID(nid,"longitude",lonvar)
    if (ierr.NE.NF_NOERR) then
-      write(*,*) 'ERROR: Field <longitude> is missing in file '//trim(file)
-      stop "" 
+      write(*,*) 'Field <longitude> is missing, looking for lon'
+      ierr=NF_INQ_VARID(nid,"lon",lonvar)
+      if (ierr.NE.NF_NOERR) then
+        write(*,*) 'ERROR: Field <longitude> and <lon> is missing in file '//trim(file)
+        stop ""  
+      endif
    endif
    ierr=NF_INQ_DIMLEN(nid,londim,lonlen)
@@ -300 +316 @@
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'Dimension <index> is missing in file '//trim(file)
-      ctllen=0
+      write(*,*) 'Looking for controle_axe'
+      ierr=NF_INQ_DIMID(nid,"controle_axe",ctldim)
+      if (ierr.NE.NF_NOERR) then
+         ctllen=0
+      endif
       !stop ""  
    endif
+!   ierr=NF_INQ_DIMID(nid,"controle_axe",ctldim)
+!   if (ierr.NE.NF_NOERR) then
+!      !stop ""  
+!   endif
    ierr=NF_INQ_VARID(nid,"controle",ctlvar)
    if (ierr.NE.NF_NOERR) then
@@ -386 +410 @@
    ierr=NF_INQ_DIMID(nid,"Time",timedim)
    if (ierr.NE.NF_NOERR) then
-      write(*,*) 'ERROR: Dimension <Time> is missing in file'//trim(file)
-      stop ""
+     write(*,*) 'ERROR: Dimension <Time> is missing looking for time_counter'
+     ierr=NF_INQ_DIMID(nid,"time_counter",timedim)
+     if (ierr.NE.NF_NOERR) then
+        write(*,*) 'ERROR: Dimension <Time> or is missing in file'//trim(file)
+        stop ""
+     endif
    endif
    ierr=NF_INQ_VARID(nid,"Time",timevar)
    if (ierr.NE.NF_NOERR) then
-      write(*,*) 'ERROR: Field <Time> is missing in file'//trim(file)
-      stop ""
+     write(*,*) 'ERROR: Field <Time> is missing looking for time_counter'
+     ierr=NF_INQ_VARID(nid,"time_counter",timevar)
+     if (ierr.NE.NF_NOERR) then
+        write(*,*) 'ERROR: Field <Time> is missing in file '//trim(file)
+        stop ""
+     endif
    endif
    ierr=NF_INQ_DIMLEN(nid,timedim,timelen)