SUBROUTINE thresh_supersat(t,RH_crit)

        IMPLICIT NONE

        REAL, intent(in) :: t   ! temperature (K)
        REAL, intent(out) :: RH_crit ! seuil de sursaturation en %
                                     

! Selon Kaercher and Lohmann (JGR 2002)

        RH_crit = 2.583-t/207.83

        END SUBROUTINE thresh_supersat