!
! $Id: hgardfou.F 1795 2013-07-18 08:20:28Z emillour $
      SUBROUTINE hgardfou (t,tsol,text)
      USE dimphy
      USE phys_state_var_mod
      USE indice_sol_mod
      IMPLICIT none
!c======================================================================
!c Verifier la temperature
!c======================================================================
#include "dimensions.h"
#include "YOMCST.h"
#include "iniprint.h"
      REAL t(klon,klev), tsol(klon,nbsrf)
      CHARACTER*(*) text
      character (len=20) :: modname = 'hgardfou'
      character (len=80) :: abort_message
!C
      INTEGER i, k, nsrf
      REAL zt(klon)
      INTEGER jadrs(klon), jbad
      LOGICAL ok
!c
      LOGICAL firstcall
      SAVE firstcall
      DATA firstcall /.TRUE./
!$OMP THREADPRIVATE(firstcall)

! Lluis
      INTEGER                                            :: lp
      lp = 885

      IF (firstcall) THEN
         WRITE(lunout,*)                                                             &
       &  'hgardfou garantit la temperature dans [100,370] K'
         firstcall = .FALSE.
!c        DO i = 1, klon
!c         WRITE(lunout,*)'i=',i,'rlon=',rlon(i),'rlat=',rlat(i)
!c        ENDDO
!c
      ENDIF
!c
      ok = .TRUE.
      DO k = 1, klev
         DO i = 1, klon
            zt(i) = t(i,k)
         ENDDO
#ifdef CRAY
         CALL WHENFGT(klon, zt, 1, 370.0, jadrs, jbad)
#else
         jbad = 0
         DO i = 1, klon
         IF (zt(i) > 370.) THEN
            jbad = jbad + 1
            jadrs(jbad) = i
         ENDIF
         ENDDO
#endif
         IF (jbad .GT. 0) THEN
           ok = .FALSE.
           DO i = 1, jbad
             WRITE(lunout,*)                                                         &
       &       'i,k,temperature,lon,lat,pourc ter,lic,oce,sic =',                    &
       &       jadrs(i),k,zt(jadrs(i)),rlon(jadrs(i)),rlat(jadrs(i)),                &
       &       (pctsrf(jadrs(i),nsrf),nsrf=1,nbsrf)
           ENDDO
         ENDIF
#ifdef CRAY
         CALL WHENFLT(klon, zt, 1, 100.0, jadrs, jbad)
#else
         jbad = 0
         DO i = 1, klon
!         IF (zt(i).LT.100.0) THEN
         IF (zt(i).LT.50.0) THEN
            jbad = jbad + 1
            jadrs(jbad) = i
         ENDIF
         ENDDO
#endif
         IF (jbad .GT. 0) THEN
           ok = .FALSE.
           DO i = 1, jbad
             WRITE(lunout,*)                                                         &
       &       'i,k,temperature,lon,lat,pourc ter,lic,oce,sic =',                    &
       &       jadrs(i),k,zt(jadrs(i)),rlon(jadrs(i)),rlat(jadrs(i)),                &
       &       (pctsrf(jadrs(i),nsrf),nsrf=1,nbsrf)
           ENDDO
         ENDIF
      ENDDO
!c
      DO nsrf = 1, nbsrf
         PRINT *,'  Lluis nsrf: ',nsrf
         DO i = 1, klon
            IF (pctsrf(i,nsrf) > 0.) THEN
              zt(i) = tsol(i,nsrf)
            ELSE
              zt(i) = 300.0
            END IF
         ENDDO
#ifdef CRAY
         CALL WHENFGT(klon, zt, 1, 370.0, jadrs, jbad)
#else
         jbad = 0
         DO i = 1, klon
         IF (zt(i).GT.370.0) THEN
            jbad = jbad + 1
            jadrs(jbad) = i
         ENDIF
         ENDDO
#endif
         IF (jbad .GT. 0) THEN
           ok = .FALSE.
           DO i = 1, jbad
            WRITE(lunout,*)                                                          &
       &      'i,nsrf,temperature,lon,lat,pourc ter,lic,oce,sic ='                   &
       &      ,jadrs(i),nsrf,zt(jadrs(i)),rlon(jadrs(i)),rlat(jadrs(i))              &
       &      ,pctsrf(jadrs(i),nsrf)
           ENDDO
         ENDIF
#ifdef CRAY
         CALL WHENFLT(klon, zt, 1, 100.0, jadrs, jbad)
#else
         jbad = 0
         DO i = 1, klon
!         IF (zt(i).LT.100.0) THEN
         IF (zt(i).LT.50.0) THEN
            jbad = jbad + 1
            jadrs(jbad) = i
         ENDIF
         ENDDO
#endif
         IF (jbad .GT. 0) THEN
           ok = .FALSE.
           DO i = 1, jbad
            WRITE(lunout,*)                                                          &
       &      'i,nsrf,temperature,lon,lat,pourc ter,lic,oce,sic ='                   &
       &      ,jadrs(i),nsrf,zt(jadrs(i)),rlon(jadrs(i)),rlat(jadrs(i))              &
       &      ,pctsrf(jadrs(i),nsrf)
             PRINT *,'    Lluis: nsrf k isfc tsol zt pctsrf_____'
             DO k=1,nbsrf
               PRINT *,'    ',nsrf,k,jadrs(i),tsol(jadrs(i),k),zt(jadrs(i)),pctsrf(jadrs(i),k)
             END DO
           ENDDO
         ENDIF
      ENDDO
!c
      IF (.NOT. ok) THEN
         abort_message= 'hgardfou s arrete '//text
         CALL abort_gcm (modname,abort_message,1)
      ENDIF

      RETURN
      END