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

trunk/LMDZ.MARS/util/concatnc.F90

@@ -131 +131 @@
 ! valid_range(2): [netcdf] interval in which a value is considered valid
 real :: year_day = 669.
+integer :: XIOS =0
 
 !==============================================================================
@@ -328 +329 @@
    ierr=NF_INQ_VARID(nid,"latitude",latvar)
    if (ierr.NE.NF_NOERR) then
-      write(*,*) 'ERROR: Field <latitude> is missing in file '//file(i)
+     ierr=NF_INQ_DIMID(nid,"lat",latdim)
+     ierr=NF_INQ_VARID(nid,"lat",latvar)
+     if (ierr.EQ.NF_NOERR) then
+        write(*,*) 'XIOS concat spotted'
+        XIOS=1
+     endif
+   endif  
+   if (ierr.NE.NF_NOERR) then
+      write(*,*) 'ERROR: Field <latitude> or <lat> is missing in file '//file(i)
       stop ""  
    endif
@@ -334 +343 @@
 !  write(*,*) "latlen: ",latlen
 
-   ierr=NF_INQ_DIMID(nid,"longitude",londim)
-   ierr=NF_INQ_VARID(nid,"longitude",lonvar)
-   if (ierr.NE.NF_NOERR) then
-      write(*,*) 'ERROR: Field <longitude> is missing in file '//file(i)
-      stop "" 
+   if (XIOS.EQ.0) then
+       ierr=NF_INQ_DIMID(nid,"longitude",londim)
+       ierr=NF_INQ_VARID(nid,"longitude",lonvar)
+       if (ierr.NE.NF_NOERR) then
+          write(*,*) 'ERROR: Field <longitude> is missing in file '//file(i)
+          stop "" 
+       endif
+       ierr=NF_INQ_DIMLEN(nid,londim,lonlen)
+       !  write(*,*) "lonlen: ",lonlen
+   else
+       ierr=NF_INQ_DIMID(nid,"lon",londim)
+       ierr=NF_INQ_VARID(nid,"lon",lonvar)
+       if (ierr.NE.NF_NOERR) then
+          write(*,*) 'ERROR: Field <lon> is missing in file '//file(i)
+          stop "" 
+       endif
+       ierr=NF_INQ_DIMLEN(nid,londim,lonlen)
+       !  write(*,*) "lonlen: ",lonlen
    endif
-   ierr=NF_INQ_DIMLEN(nid,londim,lonlen)
-!  write(*,*) "lonlen: ",lonlen
 
    ierr=NF_INQ_DIMID(nid,"altitude",altdim)
@@ -368 +388 @@
 ! 2.2.1 Check presence and read "controle" dimension & variable from input file
 !==============================================================================
-   ierr=NF_INQ_DIMID(nid,"index",ctldim)
-   if (ierr.NE.NF_NOERR) then
-      write(*,*) 'Dimension <index> is missing in file '//file(i)
-      write(*,*) "The program continues..."
-      ctllen=0
-      !stop ""
+   if (XIOS.EQ.0) then
+        ierr=NF_INQ_DIMID(nid,"index",ctldim)
+        if (ierr.NE.NF_NOERR) then
+          write(*,*) 'Dimension <index> is missing in file '//file(i)
+          write(*,*) "The program continues..."
+          ctllen=0
+          !stop ""
+       else
+          ierr=NF_INQ_VARID(nid,"controle",ctlvar)
+          if (ierr.NE.NF_NOERR) then
+            write(*,*) 'Field <controle> is missing in file '//file(i)
+            write(*,*) "The program continues..."
+            ctllen=0
+            !stop ""
+          else
+            ierr=NF_INQ_DIMLEN(nid,ctldim,ctllen)
+          endif
+       endif
    else
-      ierr=NF_INQ_VARID(nid,"controle",ctlvar)
-      if (ierr.NE.NF_NOERR) then
-         write(*,*) 'Field <controle> is missing in file '//file(i)
-         write(*,*) "The program continues..."
-         ctllen=0
-         !stop ""
-      else
-         ierr=NF_INQ_DIMLEN(nid,ctldim,ctllen)
-      endif
+     ierr=NF_INQ_DIMID(nid,"controle_axe",ctldim)
+     if (ierr.NE.NF_NOERR) then
+       write(*,*) 'Dimension <controle_axe> is missing in file '//file(i)
+       write(*,*) "The program continues..."
+       ctllen=0
+       !stop ""
+     else
+       ierr=NF_INQ_VARID(nid,"controle",ctlvar)
+       if (ierr.NE.NF_NOERR) then
+         write(*,*) 'Field <controle> is missing in file '//file(i)
+         write(*,*) "The program continues..."
+         ctllen=0
+         !stop ""
+        else
+         ierr=NF_INQ_DIMLEN(nid,ctldim,ctllen)
+       endif
+     endif
    endif
 
@@ -395 +435 @@
 ! 2.3 Read "Time" dimension & variable from input file
 !==============================================================================
+   if (XIOS.EQ.0) then
    ierr=NF_INQ_DIMID(nid,"Time",timedim)
    if (ierr.NE.NF_NOERR) then
@@ -407 +448 @@
    ierr=NF_INQ_DIMLEN(nid,timedim,timelen)
 !  write(*,*) "timelen: ",timelen
+   else
+   ierr=NF_INQ_DIMID(nid,"time_counter",timedim)
+   if (ierr.NE.NF_NOERR) then
+      write(*,*) 'ERROR: Dimension <time_counter> is missing in file '//file(i)
+      stop ""
+   endif
+   ierr=NF_INQ_VARID(nid,"time_counter",timevar)
+   if (ierr.NE.NF_NOERR) then
+      write(*,*) 'ERROR: Field <time_counter> is missing in file '//file(i)
+      stop ""
+   endif
+   ierr=NF_INQ_DIMLEN(nid,timedim,timelen)
+!  write(*,*) "timelen: ",timelen
+   endif
 
    ! allocate time() array and fill it with values from input file
@@ -513 +568 @@
 !         correction if the file to concatenate seems "before" previous file
 !         (for instance to concatenate diagfi from the previous year at the enf of a year)
+          print *, "time(rep)", time(rep), rep
+          print *, "ctlsol", ctlsol
+          print *, "starttimeoffset", starttimeoffset
+          print *, "output_time", output_time
+          print *, "previous_last_output_time", previous_last_output_time
           do while (output_time.lt.previous_last_output_time)  
                    output_time = output_time + year_day