#include "dims.h"
      SUBROUTINE DFACT(T,J,LEN)
      use cons_module,only: hor
      implicit none
#include "params.h"
!
! Args:
      integer,intent(in) :: j,len
      real,intent(out) :: T(len)
!
! Local:
      integer :: i
      real :: rlat,fack
C
C     ****     T3 = A = .5*(1.+SIN(PI*(ABS(RLATM)-PI/6.)/(PI/3.)))
C     ****               FOR  ABS(RLATM).LT.PI/3.
C     ****          A = 1.    FOR  ABS(RLATM).GE.PI/3
C
      RLAT =(-87.5+(J-1)*5.)/57.295
!
!     FACK = merge(1.,.5*(1.+SIN(pi*(ABS(RLAT)-pi/9.)
!    1          /(pi/4.5))),ABS(RLAT)-pi/4.5>=0.)
!     fack = 1.
!     if (abs(rlat) < pi/4.5) fack = .5*(1.+sin(pi*(abs(rlat)-
!    |  pi/9.)/(pi/4.5)))
!
      DO 1 I=1,LEN
! 1/4/01: do not use fack, as per ~roble/timegcm/tgcm22mt/modsrc.kibo12:
!         (HOR is set in runmds.f)
      T(I)=HOR(J)
C     T(I)=HOR(J)+FACK
    1 CONTINUE
      RETURN
      END
C