1 double precision function pinorm(x)
5 c this function takes an angle (in radians) and puts it in the range of
11 pinorm = x - n * dwapi
12 if ( pinorm .gt. pi ) then
13 pinorm = pinorm - dwapi
14 else if ( pinorm .lt. - pi ) then
15 pinorm = pinorm + dwapi