SUBROUTINE DFACT(T,J,LEN) implicit none include "params.h" include "cons.h" ! integer :: nerd real :: hor COMMON/RUNTIM/NERD(ZJMX+2),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 FACK = merge(1.,.5*(1.+SIN(C(110)*(ABS(RLAT)-C(110)/9.) 1 /(C(110)/4.5))),ABS(RLAT)-C(110)/4.5>=0.) DO 1 I=1,LEN C T(I)=HOR(J) T(I)=HOR(J)+FACK 1 CONTINUE RETURN END C