Changeset 1616

Timestamp:

Sep 27, 2016, 9:31:17 AM (8 years ago)

Author:

emillour

Message:

Generic physics: better control of warning messages in tpindex, only write these if strictly out of range.
EM

File:

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

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

@@ -76 +76 @@
       IF(TW.LE.TREF(1)) THEN
         MT = 1
+        IF (TW.LT.TREF(1)) THEN
+         write(*,*) 'tpindex: Caution! Temperature of upper levels lower 
+     $ than ref temperature for k-coef: k-coeff fixed for upper levels'
+         write(*,*) "         TW=",TW
+         write(*,*) "         TREF(1)=",TREF(1)
+        ENDIF
       ELSE
         do n=1,L_NTREF-1
@@ -92 +98 @@
       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'
+!         write(*,*) 'tpindex: Caution! Temperature of upper levels lower 
+!     $than ref temperature for k-coef: k-coeff fixed for upper levels'
+!         write(*,*) "         TW=",TW
+!         write(*,*) "         TREF(1)=",TREF(1)
       ENDIF
 
@@ -117 +125 @@
       !TB15 : case low pressure : n=2 : fixed pwl, right above pref(1)
       IF (MP.eq.1) THEN
-         PWL=Pref(1)*1.00
+        IF (PWL.LT.PREF(1)) THEN
          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)
+        ENDIF
+        PWL=Pref(1)*1.00
       ENDIF