Changeset 3834

Timestamp:

Jul 7, 2025, 1:39:14 PM (4 days ago)

Author:

afalco

Message:

Pluto: fix for zrecast, reading phisfi when phisinit unavailable.
AF

Location:

trunk/LMDZ.PLUTO/util

Files:

• compile (modified) (1 diff)
• zrecast.F90 (modified) (3 diffs)

trunk/LMDZ.PLUTO/util/compile

@@ -11 +11 @@
 
 FLAGS=
-# FLAGS=-g -O0
+# FLAGS="-g -O0"
 
 #pgf90 $SOURCE.F90 \

trunk/LMDZ.PLUTO/util/zrecast.F90

@@ -596 +596 @@
     write(*,*) 'OK, got phisinit or phsifi'
     ierr=NF_CLOSE(infid2)
+  else
+    ierr=NF_GET_VAR_REAL(infid,tmpvarid,phisinit)
+    if (ierr.ne.NF_NOERR) then
+        stop "Error: Failed reading phisinit/phisfi"
+    endif
   endif
 else
@@ -603 +608 @@
   endif
 endif
+print*,"phisinit range:", minval(phisinit), maxval(phisinit)
 
 !===============================================================================
@@ -953 +959 @@
     call crude_gcm_za(lonlength,latlength,altlength,timelength, &
                       phisinit,ps,press,temp,za_gcm)
+    print*, "GCM altitude range:", minval(za_gcm),maxval(za_gcm)
   else if (ztype.eq.3) then ! above local surface altitudes
     allocate(zs_gcm(lonlength,latlength,altlength,timelength),stat=ierr)