Changeset 4169 for trunk/LMDZ.MARS/util/zrecast.F90

Timestamp:

Apr 2, 2026, 6:19:13 PM (30 hours ago)

Author:

jmauxion

Message:

Mars PCM:

• Switching output variable name from "phisfi" to "phisinit" in physiq_mod to ensure compatibility with utils.
• Updating xml files accordingly.
• Switching localtime and zrecast back to previous version (undoing changes in r4161).

JM

File:

• trunk/LMDZ.MARS/util/zrecast.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/util/zrecast.F90

@@ -654 +654 @@
 if (ierr.ne.NF_NOERR) then
   write(*,*) "Warning: Failed to get phisinit ID from file ",trim(infile)
-  write(*,*) "Trying name phisfi (xios name)..."
-  ierr=NF_INQ_VARID(infid,"phisfi",tmpvarid)
+  infile2="diagfi.nc"
+  write(*,*) "         Trying file ",trim(infile2)
+  ierr=NF_OPEN(infile2,NF_NOWRITE,infid2)
   if (ierr.ne.NF_NOERR) then
-    write(*,*) "Failed to find phisfi..."
-    infile2="diagfi.nc"
+    write(*,*) "Problem: Could not find/open that file"
+    infile2="diagfi1.nc"
     write(*,*) "         Trying file ",trim(infile2)
     ierr=NF_OPEN(infile2,NF_NOWRITE,infid2)
     if (ierr.ne.NF_NOERR) then
       write(*,*) "Problem: Could not find/open that file"
-      infile2="diagfi1.nc"
+      infile2="phisinit.nc"
       write(*,*) "         Trying file ",trim(infile2)
       ierr=NF_OPEN(infile2,NF_NOWRITE,infid2)
       if (ierr.ne.NF_NOERR) then
-        write(*,*) "Problem: Could not find/open that file"
-        infile2="phisinit.nc"
-        write(*,*) "         Trying file ",trim(infile2)
-        ierr=NF_OPEN(infile2,NF_NOWRITE,infid2)
-        if (ierr.ne.NF_NOERR) then
-          write(*,*) "Error: Could not open that file either"
-          write(*,*) "Might as well stop here"
-          stop
-        endif
+        write(*,*) "Error: Could not open that file either"
+        write(*,*) "Might as well stop here"
+        stop
       endif
     endif
-
-    ! Get ID for phisinit
-    ierr=NF_INQ_VARID(infid2,"phisinit",tmpvarid)
-    if (ierr.ne.NF_NOERR) then
-      write(*,*) "Error: Failed to get phisinit ID"
-      stop
-    endif
-  
-    ! Get physinit
-    ierr=NF_GET_VAR_REAL(infid2,tmpvarid,phisinit)
-    if (ierr.ne.NF_NOERR) then
-      write(*,*) "Error: Failed reading phisinit"
-      stop
-    endif
-    ! Close file
-    write(*,*) 'OK, got phisinit'
-    ierr=NF_CLOSE(infid2)
-  else
-    ierr=NF_GET_VAR_REAL(infid,tmpvarid,phisinit)
-    if (ierr.ne.NF_NOERR) then
-      write(*,*) "Error: Failed reading phisinit"
-      stop
-    endif
-  endif
+  endif
+
+  ! Get ID for phisinit
+  ierr=NF_INQ_VARID(infid2,"phisinit",tmpvarid)
+  if (ierr.ne.NF_NOERR) then
+    write(*,*) "Error: Failed to get phisinit ID"
+    stop
+  endif
+  ! Get physinit
+  ierr=NF_GET_VAR_REAL(infid2,tmpvarid,phisinit)
+  if (ierr.ne.NF_NOERR) then
+    write(*,*) "Error: Failed reading phisinit"
+    stop
+  endif
+  ! Close file
+  write(*,*) 'OK, got phisinit'
+  ierr=NF_CLOSE(infid2)
 else
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,phisinit)
@@ -709 +697 @@
   endif
 endif
+
 !===============================================================================
 ! 1.4 Choose and build the new vertical coordinate