Changeset 2821 for trunk/LMDZ.GENERIC

Timestamp:

Nov 17, 2022, 3:42:23 PM (2 years ago)

Author:

aslmd

Message:

Fixed the extrapolation in tpindex.F

if above the Tmax or Pmax of the grid, we go back to the maximum value of the grid. We do not extrapolate anymore. LT

File:

• trunk/LMDZ.GENERIC/libf/phystd/tpindex.F (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/tpindex.F

@@ -83 +83 @@
 !     Get the upper and lower temperature grid indicies that bound the
 !     requested temperature. If the sought temperature is below
-!     the T-grid minimum, assume that lower value. If the sought teperature
-!     is above the T-grid maximum, extrapolate.
+!     the T-grid minimum, assume that lower value. If the sought temperature
+!     is above the T-grid maximum, assume that maximum value.
 !     TW : temperature to interpolate to
 !     TREF : reference temperature array
@@ -96 +96 @@
         MT = 1
         TWL=TREF(1)*1.00
-        write(*,*) 'tpindex: Caution! Temperature of upper levels lower 
-     $ than ref temperature for k-coef: k-coeff fixed for upper levels'
-        write(*,*) "         TW=",TWL
-        write(*,*) "         TREF(1)=",TREF(1)
       ELSEIF(TWL.GE.TREF(L_NTREF)) THEN
        ! Above highest tabulated temperature
-       MT = L_NTREF-1
+        MT=L_NTREF-1
+        TWL = TREF(L_NTREF)*1.00
       ELSE 
         ! Regular case, find encompassing tabulated temperatures
@@ -121 +118 @@
 !     requested pressure. If the sought pressure is below
 !     the P-grid minimum, assume that lower value. If the sought pressure
-!     is above the P-grid maximum, extrapolate.
+!     is above the P-grid maximum, assume that maximum value.
 
       PWL = log10(PW)
@@ -129 +126 @@
         MP = 1
         PWL=Pref(1)*1.00
-        write(*,*) 'tpindex: Caution! Pressure of upper levels lower 
-     $than ref pressure for k-coef: k-coeff fixed for upper levels'
-        write(*,*) "         PWL=",PWL
-        write(*,*) "         PREF(1)=",PREF(1)
       ELSEIF(PWL.GE.PREF(L_PINT)) THEN
         ! Above highest tabulated pressure
         MP = L_PINT-1
+        PWL = PREF(L_PINT)
       ELSE 
         ! Regular case, find encompassing tabulated pressures