1subroutine tm2(day,xlat4,xlon4,xl4,b4)
2
3  implicit real*8 (a-h,o-z)
4  parameter (RADS=0.0174532925199433d0)
5
6  real*4 xlat4,xlon4,xl4,b4
7
8  glat=xlat4*RADS
9  glong=xlon4*RADS
10  call tmoonsub(day,glat,glong,el,rv,xl,b,pax)
11  xl4=xl
12  b4=b
13
14end subroutine tm2
15