In view of the difficulties in data processing and adjustment of coordinate calculation perplexed large-scale GNSS reference station network,this paper proposes a coordinate calculation method.The method is elaborated in four aspects including selection of global framework sites,division of reference station network subnets,single-day data processing of base stations,and global network adjustment calculation.Finally,the internal and external consistence accuracies of the calculation results are evaluated.The analysis results show that the method can solve reference coordinates for network with more than 3 000 stations,and the RMSE of the results is better than 3 mm in both north-south and east-west directions,and better than 6 mm in elevation.The accuracy of the calculation is higher than those of most IGS analysis centers.