Changeset 2215

Timestamp:

Jan 10, 2020, 10:59:33 AM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Improve localtime utility: handle the case when fields to be
interpolated may contain "undefined" (a.k.a "missing") values by
refraining from doing the interpolation in these cases.
EM

File:

• trunk/LMDZ.MARS/util/localtime.F90 (modified) (1 diff)

trunk/LMDZ.MARS/util/localtime.F90

@@ -1275 +1275 @@
      if (((x.ge.xd(i)).and.(x.le.xd(i+1))).or.&
           ((x.le.xd(i)).and.(x.ge.xd(i+1)))) then
-        y=yd(i)+(x-xd(i))*(yd(i+1)-yd(i))/(xd(i+1)-xd(i))
+        if ((yd(i).eq.missing).or.(yd(i+1).eq.missing)) then
+          ! cannot perform the interpolation if an encompasing value
+          ! is already set to 'missing'
+        else
+          !linear interpolation based on encompassing values
+          y=yd(i)+(x-xd(i))*(yd(i+1)-yd(i))/(xd(i+1)-xd(i))
+        endif
         exit
      endif