& vpsipred(3,maxdih_constr),sdihed(2,maxdih_constr),wdihc
common /torcnstr/ phi0,drange,ftors,theta_constr0,theta_drange,
& for_thet_constr,vpsipred,sdihed,wdihc,
& vpsipred(3,maxdih_constr),sdihed(2,maxdih_constr),wdihc
common /torcnstr/ phi0,drange,ftors,theta_constr0,theta_drange,
& for_thet_constr,vpsipred,sdihed,wdihc,