Changeset 2546 for trunk/LMDZ.MARS/util/zrecast.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/zrecast.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/util/zrecast.F90

@@ -400 +400 @@
   ierr=NF_INQ_DIMID(infid,"time",tmpdimid)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to get lon dimension ID"
+    write(*,*) "Could not get time dimension ID"
+    write(*,*) "  looking for time_counter dimension instead... "
+    ierr=NF_INQ_DIMID(infid,"time_counter",tmpdimid)
+    if (ierr.ne.NF_NOERR) then
+      stop "Error: Failed to get time_counter dimension ID"
+    else
+      ierr=NF_INQ_VARID(infid,"time_counter",tmpvarid)
+      if (ierr.ne.NF_NOERR) then
+        stop "Error: Failed to get time_counter ID"
+      else
+        ierr=NF_INQ_DIMLEN(infid,tmpdimid,timelength)
+        if (ierr.ne.NF_NOERR) then
+          stop "Error: Failed to get time_counter length"
+        else
+          allocate(time(timelength),stat=ierr)
+          if (ierr.ne.0) then
+            write(*,*) "Error: Failed to allocate time(timelength)"
+            write(*,*) "     timelength=",timelength
+            stop
+          endif
+          ierr=NF_GET_VAR_REAL(infid,tmpvarid,time)
+          if (ierr.ne.NF_NOERR) then
+            stop "Error: Failed reading time"
+          endif
+        endif
+      endif
+    endif
   else
     ierr=NF_INQ_VARID(infid,"time",tmpvarid)
@@ -471 +497 @@
 if (ierr.ne.NF_NOERR) then
   write(*,*) "Could not get controle dimension ID"
-  ctllength = 0
+  write(*,*) "  looking for controle_axe dimension instead... "
+  ierr=NF_INQ_DIMID(infid,"controle_axe",tmpdimid)
+  if (ierr.ne.NF_NOERR) then
+    ctllength = 0
+  else
+    ierr=NF_INQ_VARID(infid,"controle",tmpvarid)
+    if (ierr.ne.NF_NOERR) then
+      stop "Error: Failed to get controle ID"
+    else
+      ierr=NF_INQ_DIMLEN(infid,tmpdimid,ctllength)
+      if (ierr.ne.NF_NOERR) then
+        stop "Error: Failed to get controle length"
+      else
+        allocate(ctl(ctllength),stat=ierr)
+        if (ierr.ne.0) then
+          write(*,*) "Error: Failed to allocate ctl(ctllength)"
+          write(*,*) "     ctllength=",ctllength
+          stop
+        endif
+        ierr=NF_GET_VAR_REAL(infid,tmpvarid,ctl)
+        if (ierr.ne.NF_NOERR) then
+          stop "Error: Failed reading controle"
+        endif
+      endif
+    endif
+  endif
 else
   ierr=NF_INQ_VARID(infid,"controle",tmpvarid)
@@ -943 +994 @@
      ierr=NF_INQ_VARID(infid,"teta",tmpvarid)
      if (ierr.ne.NF_NOERR) then
-       write(*,*) "Error: Failed to get temp, t or teta ID"
-       have_temp=.false.
-       if (ztype.ne.1) stop
+       ierr=NF_INQ_VARID(infid,"temperature",tmpvarid)
+       if (ierr.ne.NF_NOERR) then
+          write(*,*) "Error: Failed to get temp, t, teta or temperature ID"
+          have_temp=.false.
+          if (ztype.ne.1) stop
+       endif
      endif
        ierr=NF_GET_VAR_REAL(infid,tmpvarid,teta)