Changeset 1328 for trunk/LMDZ.GENERIC

Timestamp:

Aug 21, 2014, 5:32:24 PM (10 years ago)

Author:

sglmd

Message:

minor change to have a more consistent pressure interpolation for uneven pressure grids

File:

• trunk/LMDZ.GENERIC/libf/phystd/sugas_corrk.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/sugas_corrk.F90

@@ -220 +220 @@
       pgasmax = 10.0**pgasref(L_NPREF)
 
-      ! interpolate to finer grid
+      ! interpolate to finer grid, adapted to uneven grids
       do n=1,L_NPREF-1
          do m=1,5
-            pfgasref((n-1)*5+m) = pgasref(n)+(m-1)*0.2
+            pfgasref((n-1)*5+m) = pgasref(n)+(m-1)*(pgasref(n+1) - pgasref(n))/5.
          end do
       end do
       pfgasref(L_PINT) = pgasref(L_NPREF)
-      ! Warning, pfgasref needs generalisation to uneven grids!
 
 !-----------------------------------------------------------------------