theta_constr0(i)=deg2rad*theta_constr0(i)
theta_drange(i)=deg2rad*theta_drange(i)
enddo
C if(me.eq.king.or..not.out1file)
C & write (iout,*) 'FTORS',ftors
theta_constr0(i)=deg2rad*theta_constr0(i)
theta_drange(i)=deg2rad*theta_drange(i)
enddo
C if(me.eq.king.or..not.out1file)
C & write (iout,*) 'FTORS',ftors