Changeset 1600 for trunk/LMDZ.GENERIC

Timestamp:

Sep 15, 2016, 4:33:30 PM (8 years ago)

Author:

emillour

Message:

Generic model:

• bug fix in tpindex.F : when computing weights for abundance of species one should not return the maximum index (leads to out-of-bounds index use in the caller).

TB

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• libf/phystd/tpindex.F (modified) (4 diffs)

trunk/LMDZ.GENERIC/README

@@ -1257 +1257 @@
 - Further cleanup to harmonize with LMDZ.COMMON turn "idissip" into
   "dissip_period".
+
+== 15/09/2016 == TB
+- bug fix in tpindex.F : when computing weights for abundance of species one
+  should not return the maximum index (leads to out-of-bounds index use in 
+  the caller).

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

@@ -69 +69 @@
 !     requested temperature. If the requested temperature is outside
 !     the T-grid, set up to extrapolate from the appropriate end.
-
+!     TW : temperature to be interpolated
+!     TREF : grid array
+!     MT : index of TREF for bounding new temperature
+!     U : new index (real) for temperature interpolated
 
       IF(TW.LE.TREF(1)) THEN
@@ -85 +88 @@
    10   continue
       END IF
+
+      !TB15 : case low temp : MT=1: fixed TW right above tref(1)
+      IF (MT.eq.1) THEN
+         TW=tref(1)*1.00
+         write(*,*) 'tpindex: Caution! Pressure of upper levels lower 
+     $than ref pressure for k-coef: k-coeff fixed for upper levels'
+      ENDIF
 
       U = (TW-TREF(MT))/(TREF(MT+1)-TREF(MT))
@@ -104 +114 @@
 
    20 continue
+      
+      !TB15 : case low pressure : n=2 : fixed pwl, right above pref(1)
+      IF (MP.eq.1) THEN
+         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'
+      ENDIF
 
+!     interpolated pressure
       T = (PWL-PREF(MP))/(PREF(MP+1)-PREF(MP))
 
@@ -118 +136 @@
       IF(QVAR.le.WREFVAR(1)) then
         NVAR   = 1
-        WRATIO = 0.0D0
+        WRATIO = 0.0D0      ! put all the weight on the first point
       ELSEIF(QVAR.ge.WREFVAR(L_REFVAR)) then
-        NVAR   = L_REFVAR
-        WRATIO = 0.0D0
+        NVAR   = L_REFVAR-1 ! TB16 in order to not oversize NVAr when doing
+                                  !NVAR+1
+        WRATIO = 1.00D0     ! put all the weight on the last point
       ELSE
         DO N=2,L_REFVAR