#include "dims.h"
      SUBROUTINE DFACT(T,J,LEN)
      use cons_module,only: pi,hor
      implicit none
#include "params.h"
!
!     real :: ttbound,nob,avto,hor
!     COMMON/RUNTIM/TTBOUND,NOB(ZJMX),AVTO,HOR(ZJMX)
!
! Args:
      real,intent(out) :: T(1)
      integer,intent(in) :: j,len
!
! Local:
      real :: rlat,fack
      integer :: i
C
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
      if (ABS(RLAT)-pi/4.5>=0.) then
        fack = 1.
      else
        fack = .5*(1.+SIN(pi*(ABS(RLAT)-pi/9.)/(pi/4.5)))
      endif
      DO 1 I=1,LEN
C     T(I)=HOR(J)
      T(I)=HOR(J)+FACK
    1 CONTINUE
      RETURN
      END
C