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