首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 797 毫秒
1.
By using IAP 9L AGCM, two sets of long-term climatological integration have been per-formed with the two different interpolation procedures for generating the daily surface boundary conditions. One interpolation procedure is the so-called “traditional” scheme, for which the daily surface boundary conditions are obtained by linearly interpolating between the observed monthly mean values, however the observed monthly means cannot be preserved after interpolation. The other one is the “new” scheme, for which the daily surface boundary conditions are obtained by linearly interpolating between the "artificial" monthly mean values which are based on, but are dif-ferent from the observed ones, after interpolating with this new scheme, not only the observed monthly mean values are preserved, the time series of the new generated daily values is also more consistent with the observation. Comparison of the model results shows that the differences of the globally or zonally averaged fields between these two integrations are quite small, and this is due to the compensating effect between the different regions. However, the differences of the two patterns (the global or regional geographical distributions), are quite significant, for example, the magni-tude of the difference in the JJA mean rainfall between these two integrations can exceed 2 mm/day over Asian monsoon regions, and the difference in DJF mean surface air temperature can also exceed 2oC over this region. The fact that the model climatology depends quite strongly on the method of prescribing the daily surface boundary conditions suggests that in order to validate the climate model or to predict the short-term climate anomalies, either the " new* interpolation scheme or the high frequency surface boundary conditions (e.g., daily or weekly data instead of the monthly data) should be introduced. Meanwhile, as for the coupled model, the daily coupling scheme between the different component cli?mate models (e.g., atmospheric and oceanic general circulation models) is preferred in order to partly eliminate the “climate drift” problem which may appear during the course of direct coupling.  相似文献   

2.
The synoptic climatology of monthly mean surface air temperature anomaly (MMATA) at 47 surface stations in Asia is assessed by two types of maps. The first consists of simple linear correlation coefficients between the MMATA at a station and the monthly mean 700 hPa height anomaly (MMHA) at 134 grid points in the Northern Hemisphere. In winter, these fields generally contain two equally important centers: a center of maximum positive correlation approximately 1300 km from the reference station and a negative center about 2800 km to its northwest. In summer, the positive center is located an average of 800 km from the reference station and strongly dominates the MMATA, while the negative center is no longer as important as in winter. The second kind of map is a special type of anomaly composite, constructed by multiplying the correlation coefficient at each grid point by the standard deviation of MMHA at the same point. This map shows the optimum distribution of height anomalies for abnormally warm (or c  相似文献   

3.
The skill of probability density function (PDF) prediction of summer rainfall over East China using optimal ensemble schemes is evaluated based on the precipitation data from ˉve coupled atmosphere-ocean general circulation models that participate in the ENSEMBLES project. The optimal ensemble scheme in each region is the scheme with the highest skill among the four commonly-used ones: the equally-weighted ensemble (EE), EE for calibrated model-simulations (Cali-EE), the ensemble scheme based on multiple linear regression analysis (MLR), and the Bayesian ensemble scheme (Bayes). The results show that the optimal ensemble scheme is the Bayes in the southern part of East China; the Cali-EE in the Yangtze River valley, the Yangtze-Huaihe River basin, and the central part of northern China; and the MLR in the eastern part of northern China. Their PDF predictions are well calibrated, and are sharper than or have approximately equal interval-width to the climatology prediction. In all regions, these optimal ensemble schemes outperform the climatology prediction, indicating that current commonly-used multi-model ensemble schemes are able to produce skillful PDF prediction of summer rainfall over East China, even though more information for other model variables is not derived.  相似文献   

4.
In order to meet the needs of work in numerical weather forecast and in numerical simulations for climate change and ocean current, a kind of difference scheme in high precision in the time direction developed from the completely square-conservative difference scheme in explicit way is built by means of the Taylor expansion. A numerical test with 4-wave Rossby-Haurwitz waves on them and an application of them on the monthly mean current the of South China Sea are carried out, from which, it is found that not only do the new schemes have high harmony and approximate precision but also can the time step of the schemes be lengthened and can much computational time be saved. Therefore, they are worth generalizing and applying.  相似文献   

5.
To improve the simulation of the surface radiation budget and related thermal processes in arid regions, three sophisticated surface albedo schemes designed for such regions were incorporated into the Biosphere- Atmosphere Transfer Scheme (BATS). Two of these schemes are functions of the solar zenith angle (SZA), where the first one has one adjustable parameter defined as SZA1 scheme, and the second one has two empirical parameters defined as SZA2 scheme. The third albedo scheme is a function of solar angle and soil water that were developed based on arid-region observations from the Dunhuang field experiment (DHEX) (defined as DH scheme). We evaluated the performance of the original and newly-incorporated albedo schemes within BATS using the in-situ data from the Oasis System Energy and Water Cycle Field Experiment that was carried out in JinTa, Gansu arid area (JTEX). The results indicate that a control run by the original version of the BATS generates a constant albedo, while the SZA1 and SZA2 schemes basically can reproduce the observed diurnal cycle of surface albedo, although these two schemes still underestimate the albedo when SZA is high in the early morning and late afternoon, and overestimate it when SZA is low during noontime. The SZA2 scheme has a better overall performance than the SZA1 scheme. In addition, BATS with the DH scheme slightly improves the albedo simulation in magnitude as compared to that from the control run, but a diurnal cycle of albedo is not produced by this scheme. The SZA1 and SZA2 schemes significantly increase the surface absorbed solar radiation by nearly 70 W m^-2, which further raises the ground temperature by 6 K and the sensible heat flux by 35 W m^-2. The increased solar radiation, heat flux, and temperature are more consistent with the observations that those from the control run. However, a significant improvement in these three variables is not found in BATS with the DH scheme due to the neglect of the diurnal cycle of albedo. Further analysis indi  相似文献   

6.
Nepal lies on the southern slope of Himalaya in Asia. In a width ranging between 150 and 250 km, the altitude varies greatly from about 100 m at its southern border to a maximum of 8848 min the northern part. Like the variation in altitude, climatic condition varies quite a lot. Long-term monthly mean erythemal UV daily dose values for Nepal are evaluated using Total Ozone Mapping Spectrometer (TOMS) estimation from the time of its overpass between 1996 and 2003. The results are presented as summer and winter maps of mean UV levels in each satellite grid. The mean winter erythemal UV daily dose ranges between 2.1 and 3.6 kJ m-2 whereas summer values are found to lie between 4.6 and 9.7 kJ m-2. The altitude variation increases the UV levels by about 0.2 kJ km-1 in winter months, and 0.9 kJ km-1 in summer. A multiyear monthly average erythemal daily dose in most of the areas shows that the summer value is about three times higher than that in winter. Although year-to-year variation is not pronounced in high- and mid-elevation regions, UV levels seemed to decrease from 1997 to 2002 in the southern part of the country in the low elevation region by about 5.35%. Due to the combined effects of the altitude, low ozone concentration in the troposphere, and thin air, surface UV radiation at higher altitudes is found to be higher than in the surrounding regions.  相似文献   

7.
In this paper, a special three-step difference scheme is applied to the solution of nonlinear time-evolution equations, whose coefficients are determined according to accuracy constraints, necessary conditions of square conservation, and historical observation information under the linear supposition. As in the linear case, the schemes also have obvious superiority in overall performance in the nonlinear case compared with traditional finite difference schemes, e.g., the leapfrog(LF) scheme and the complete square conservation difference(CSCD) scheme that do not use historical observations in determining their coefficients, and the retrospective time integration(RTI) scheme that does not consider compatibility and square conservation. Ideal numerical experiments using the one-dimensional nonlinear advection equation with an exact solution show that this three-step scheme minimizes its root mean square error(RMSE) during the first 2500 integration steps when no shock waves occur in the exact solution, while the RTI scheme outperforms the LF scheme and CSCD scheme only in the first 1000 steps and then becomes the worst in terms of RMSE up to the 2500th step. It is concluded that reasonable consideration of accuracy, square conservation, and historical observations is also critical for good performance of a finite difference scheme for solving nonlinear equations.  相似文献   

8.
A single-column model (SCM) is developed in the regional climate model RegCM4. The evolution of a dry convection boundary layer (DCBL) is used to evaluate this SCM. Moreover, four planetary boundary layer (PBL) schemes, namely the Holtslag-Boville scheme (HB), Yonsei University scheme (YSU), and two University of Washington schemes (UW01, Grenier-Bretherton-McCaa scheme and UW09, Bretherton-Park scheme), are compared by using the SCM approach. A large-eddy simulation (LES) of the DCBL is performed as a benchmark to examine how well a PBL parameterization scheme reproduces the LES results, and several diagnostic outputs are compared to evaluate the schemes. The results show that the SCM is proper constructed. In general, with the DCBL case, the YSU scheme performs best for reproducing the LES results, which include well-mixed features and vertical sensible heat fluxes; the simulated wind speed, turbulent kinetic energy, entrainment flux, and height of the entrainment zone are all underestimated in the UW09; the UW01 has all those biases of the UW09 but larger, and the simulated potential temperature is not well mixed; the HB is the least skillful scheme, by which the PBL height, entrainment flux, height of the entrainment zone, and the vertical gradients within the mixed layer are all overestimated, and a inversion layer near the top of the surface layer is wrongly simulated.Although more cases and further testing are required, these simulations show encouraging results towards the use of this SCM framework for evaluating the simulated physical processes by the RegCM4.  相似文献   

9.
In this paper, the East Asia summer monsoon onset date lines in East China are calculated by the definition similar to the traditional one, with the ECMWF reanalyzed 850 hPa daily wind and observed, reana-lyzed and combined daily rainfall during 1980~1993. To make the onset date line as close as possible to the previous work, the earliest onset date limits have to be applied for the regions with different latitude and the daily mean datasets have to be smoothed by space before calculation, therefore their space-resolution is reduced to about 3 longitude ×1°latitude. The results show that the multiyear mean summer monsoon onset date lines are quite similar to each other. Compared with the one from the reanalysis, the 14-year average onset date line form combination is obviously improved in the southern Sichuan Basin and the correlation between observed and combined onset date is also slightly higher over the Huaihe valley and Northeast China. Since daily rainfall combination also improved the long term daily mean and standard deviation through the pentad CMAP, if no better daily dataset is available, such a kind of daily rainfall combination can be used to get reasonable result in the Indian monsoon region without sufficient observatories or over the North Pacific without any ground observation at all in future study.  相似文献   

10.
In this paper, the dynamic effect of oceanic upwelling on the intensity of El Nio-Southern Oscillation (ENSO) is studied using a simple coupled model (Zebiak-Cane Model). The term balance analysis in the temperature variability equation shows that the anomalous upwelling of the mean vertical temperature gradient and the mean advection of the anomalous meridional temperature gradient are the two of most important factors that determine the intensity of ENSO events, in which the "vertical oceanic heat flux" in the eastern equatorial Pacific (EEP) is the primary influencing factor. The lag correlation between "vertical heat flux (VHF)" and ENSO intensity shows that the highest correlation occurs when the former leads the latter by one to two weeks. The VHF is positively correlated with the background thermocline strength in the EEP, and an increase of both could result in strong ENSO variability. Comparison of the forced and coupled experiments suggests that the coupled process can affect both the intensity and frequency of ENSO.  相似文献   

11.
This paper analyzed the variations of latent heat flux (LHF) over the tropical Pacific in the period 1978-1988 by using COADS (Comprehensive Ocean and Atmospheric Data Set). It has been founded that the interannual variabili ty of LHF exhibits strong ENSO signal, with the significant increasing LHF during the recent two warm events, i.e., 1982 / 83 and 1986 / 87 and decreasing LHF in the cold episodes. However the longitudinal distribution of the LHF departures varies from event to event. In the eastern Pacific, the specific humidity difference at air-sea interface (qs -qa) makes a dominant contribution to the interannual variability of LHF ( r = 0.73 ), while in the western Pacific the surface wind speed, W and the qs - qa make nearly equal contribution to that of LHF.  相似文献   

12.
The regional climate model (RegCM4) is customized for 10-year climate simulation over Indian region through sensitivity studies on cumulus convection and land surface parameterization schemes. The model is configured over 30° E–120° E and 15° S–45° N at 30-km horizontal resolution with 23 vertical levels. Six 10-year (1991–2000) simulations are conducted with the combinations of two land surface schemes (BATS, CLM3.5) and three cumulus convection schemes (Kuo, Grell, MIT). The simulated annual and seasonal climatology of surface temperature and precipitation are compared with CRU observations. The interannual variability of these two parameters is also analyzed. The results indicate that the model simulated climatology is sensitive to the convection as well as land surface parameterization. The analysis of surface temperature (precipitation) climatology indicates that the model with CLM produces warmer (dryer) climatology, particularly over India. The warmer (dryer) climatology is due to the higher sensible heat flux (lower evapotranspiration) in CLM. The model with MIT convection scheme simulated wetter and warmer climatology (higher precipitation and temperature) with smaller Bowen ratio over southern India compared to that with the Grell and Kuo schemes. This indicates that a land surface scheme produces warmer but drier climatology with sensible heating contributing to warming where as a convection scheme warmer but wetter climatology with latent heat contributing to warming. The climatology of surface temperature over India is better simulated by the model with BATS land surface model in combination with MIT convection scheme while the precipitation climatology is better simulated with BATS land surface model in combination with Grell convection scheme. Overall, the modeling system with the combination of Grell convection and BATS land surface scheme provides better climate simulation over the Indian region.  相似文献   

13.
Comprehensive diagnostic comparisons and evaluations have been carried out with the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) and European Centre for Medium Range Weather Forecasts (ECMWF) reanalyses of the vertically integrated atmospheric energy budgets. For 1979 to 1993 the focus is on the monthly means of the divergence of the atmospheric energy transports. For February 1985 to April 1989, when there are reliable top-of-the-atmosphere (TOA) radiation data from the Earth Radiation Budget Experiment (ERBE), the implied monthly mean surface fluxes are derived and compared with those from the assimilating models and from the Comprehensive Ocean Atmosphere Data Set (COADS), both locally and zonally integrated, to deduce the implied ocean meridional heat transports. While broadscale aspects and some details of both the divergence of atmospheric energy and the surface flux climatological means are reproducible, especially in the zonal means, differences are also readily apparent. Systematic differences are typically ∼20 W m−2. The evaluation highlights the poor results over land. Land imbalances indicate local errors in the divergence of the atmospheric energy transports for monthly means on scales of 500 km (T31) of 30 W m−2 in both reanalyses and ∼50 W m−2 in areas of high topography and over Antarctica for NCEP/NCAR. Over the oceans in the extratropics, the monthly mean anomaly time series of the vertically integrated total energy divergence from the two reanalyses correspond reasonably well, with correlations exceeding 0.7. A common monthly mean climate signal of about 40 W m−2 is inferred along with local errors of 25 to 30 W m−2 in most extratropical regions. Except for large scales, there is no useful common signal in the tropics, and reproducibility is especially poor in regions of active convection and where stratocumulus prevails. Although time series of monthly anomalies of surface bulk fluxes from the two models and COADS agree very well over the northern extratropical oceans, the total fields all contain large systematic biases which make them unsuitable for determining ocean heat transports. TOA biases in absorbed shortwave, outgoing longwave and net radiation from both reanalysis models are substantial (>20 W m−2 in the tropics) and indicate that clouds are a primary source of problems in the model fluxes, both at the surface and the TOA. Time series of monthly COADS surface fluxes are shown to be unreliable south of about 20N where there are fewer than 25 observations per 5 square per month. Only the derived surface fluxes give reasonable implied meridional ocean heat transports. Received: 21 March 2000 / Accepted: 21 June 2000  相似文献   

14.
The surface heat flux feedback is estimated in the Atlantic and the North Pacific, using monthly heat flux and sea surface temperature (SST) anomaly data from COADS and the NCEP reanalysis. In the Atlantic, the heat flux feedback is dominated by the turbulent flux. It is negative, mostly ranging between 10 and 35 W m-2 K-1, but larger in the western subtropical gyre and part of the subpolar gyre, and weaker in the tropics. The radiative feedback is generally weak. In the North Pacific, the heat flux feedback is also dominated by the turbulent flux and is negative nearly everywhere, peaking in the subtropics. In both oceans, the turbulent heat flux feedback remains primarily negative in each season, and is stronger in fall and winter; patches of positive feedback can be seen, but often with limited correspondence between COADS and NCEP. The radiative feedback remains weak, and is positive in spring and summer at mid-latitudes. It is also shown that the turbulent heat flux feedback is weaker over large-scale areas, that no positive heat flux feedback sustains the SST anomaly "dipole" in the tropical Atlantic, and that the main SST anomaly mode in the North Pacific is damped by a negative heat flux feedback. The energy exchange with the atmosphere that results from the heat flux feedback can be substantial at mid-latitudes, but does not exceed 7 W m-2 at basin scale.  相似文献   

15.
Summary The seasonal and intraseasonal variation of tropical climate in National Center for Atmospheric Research (NCAR) Community Climate Model Version 2 (CCM2) General Circulation Model (GCM) has been examined using two different cumulus parameterization schemes, the moist convective adjustment scheme of Manabe et al. (1965) and the mass-flux scheme of Hack (1994). Ten-year simulations have been undertaken with each of these schemes with SST prescribed according to the monthly mean climatology. The seasonal mean rainfall in the tropics simulated by the moist convective adjustment scheme (MCA) scheme was found to be more realistic than the mass-flux (Hack) scheme. The more realistic simulation by the MCA scheme was found to be on account of the fact that the mean moist static energy of the lower troposphere in the MCA scheme was closer to the observations than in the Hack scheme. In both the schemes, the precipitation in the tropics increases montonically with precipitable water vapour when the precipitable water vapour is above 40 mm. This is consistent with relationship between precipitation and precipitable water in the observations. The Hack scheme tends to simulate lower precipitation (for a given amount of precipitable water) when compared to observations. The MCA scheme simulates the eastward migration of convective systems along the equator quite well, although the speed of propagation is somewhat low. The poleward migration of convective systems in the Indian region is more realistically simulated by the MCA scheme than the Hack scheme. This is because the latitudinal gradient of the mean moist static energy in the MCA scheme is more realistic than in the Hack scheme. Over most of the tropics, simulation by the MCA scheme is more realistic on both seasonal and intraseasonal timescales. Received November 1, 2000 Revised June 20, 2001  相似文献   

16.
为准确描述我国最大的固定/半固定沙漠-古尔班通古特沙漠区域的大气边界层结构,本文利用该沙漠腹地2017年的梯度铁塔和通量观测数据,基于中尺度气象模式WRF (Weather Research and Forecast v3.7.1),分析了5种边界层参数化方案在古尔班通古特沙漠的适用性。结果表明:1)采用WRF模拟沙漠腹地近地层内的边界层特征时,2m气温的模拟存在冷偏差,5种边界层参数化方案均能较好地模拟出四个季节2m气温的日变化特征,其中非局地方案ACM2(Asymmetric Convective Model version 2)对2m气温效果最好,局地方案BL方案的模拟偏差最大;2)5种边界层参数化方案均能够模拟出10m风速的日变化特征,其中局地方案BL(Bougeault-Lacarrere)对10m风速效果最佳;3)采用WRF模拟沙漠近地层内的地表通量特征时,感热通量存在高估现象,潜热通量存在低估现象,5种边界层参数化方案均能较好地模拟出四个季节模拟时间段内地表净辐射通量的日变化特征,其中局地方案MYJ(Mellor-Yamada-Janjie)的模拟精度最高。  相似文献   

17.
NASA/Goddard长波辐射方案在GRAPES_Meso模式中的应用研究   总被引:2,自引:0,他引:2  
张梦  王宏  黄兴友 《大气科学》2014,38(3):603-614
本文将NASA(National Aeronautics and Space Administration)/Goddard长波辐射方案引入到GRAPES_ Meso(Global/Regional Assimilation and PrEdiction System_Meso)模式中,对2006年4月中国地区进行了一个月的模拟试验,并与相应的NCEP(National Centers for Environmental Prediction)再分析资料进行了对比分析。试验结果表明:在模拟区域内,使用GRAPES_Meso模式进行24 h、48 h预报得到的晴空大气顶向外长波辐射通量(the clear sky outgoing longwave radiation flux,OLRC)、地面接收到向下长波辐射通量(the clear sky downward longwave radiation flux at ground,GLWC)分布形势与NCEP再分析资料具有较好的对应关系;模式预报24 h、48 h OLRC和NCEP再分析资料月平均误差百分比控制在-10%~+10%以内,GLWC月平均误差百分比比OLRC略大,但总体上两者误差都在合理和可接受范围之内。OLRC和GLWC 24 h、48 h的预报和NCEP再分析资料的逐日距平相关系数及标准误差的对比显示,模式24 h预报OLRC、GLWC的距平相关系数月平均值分别为0.96、0.98,标准误差月平均值分别为24.54 W m-2、27.23 W m-2;模式48 h预报OLRC、GLWC的距平相关系数月平均值分别为0.9521、0.9804,标准误差月平均值分别为22.43 W m-2、27.64 W m-2。总体上,模式24 h、48 h预报OLRC和GLWC的距平相关系数都在0.93以上,标准误差都在31 W m-2以内,且GLWC预报和NCEP再分析资料的相关性比OLRC略好,OLRC预报与NCEP再分析资料的的标准误差比GLWC略小。通过和RRTM长波辐射方案对比可知,两者的预报水平基本一致。本文研究结果表明,引入NASA/Goddard长波辐射方案后的GRAPES_Meso模式整体上能够较好地预报OLRC和GLWC,该辐射方案可以作为模式GRAPES_Meso的备选辐射方案之一。  相似文献   

18.
刘昕晔  李煜斌 《气象科技》2024,52(3):434-445
稻麦轮作农田是我国典型农田类型,其模拟效果对我国农田气候模拟具有重要参考价值。气候中尺度模拟结果对近地层通量极为敏感,选择合适的通量计算方案对模拟效果至关重要。因此,对比分析稻麦轮作农田下不同的通量计算方案具有重要意义。本文选取了8种具有代表性的近地层湍流通量计算方案,采用寿县国家气候观象台实测资料对比分析了各方案在稻麦轮作农田的计算特征和差异。结果表明,在不同稳定度和不同风速情况下,各方案的误差特征各异。本文基于归一化标准差综合评价了各方案的准确度,总体而言,所有方案的动量通量总体平均归一化平均差为0.536,其中SS14(Sharan和Srivastava,2014年)方案最大为0.575,SS20(Sharan和Srivastava等,2020年)方案最小为0.517;所有方案的感热通量总体平均归一化标准差为0.638,其中GLGS20(Gryanik等,2020年)方案最大为0.871,SS14方案最小为0.476。此外,本研究还给出了稻麦轮作农田不同稳定层结和不同风速情况下各通量计算方案的误差特征。本文的研究结论,可为准确计算近地层湍流通量提供支撑。  相似文献   

19.
We demonstrate a method to improve the performance of commonly used parameterizations to calculate the cloud-free down-welling long-wave radiation at the surface. The method uses a monthly climatology of the effective radiating temperature of the atmosphere instead of the instantaneous screen-level temperature. The climatology of the effective radiating temperature can be derived from pyrgeometer measurements and was incorporated into two commonly used schemes. We compared the calculated cloud-free down-welling long-wave irradiances to high-quality pyrgeometer measurements from four Swiss sites. The discrepancies between observations and modified schemes can be reduced by up to 35 %, resulting in a model uncertainty close to 5 W m?2 which corresponds to the measurement uncertainty of pyrgeometers. Furthermore, we introduce a new long-wave model which is based on radiative transfer calculations in the 8–14- $\upmu $ m wavelength range. In the remaining long-wave spectrum, the radiation is calculated using the Planck function with the effective radiating temperature of the atmosphere. The performance of this new model is consistent with the modified parameterizations.  相似文献   

20.
南疆沙漠腹地大气边界层湍流通量特征的观测研究   总被引:4,自引:0,他引:4  
利用新疆塔中站2006年4月、8月的三维风速。温度和水汽脉动资料,运用涡旋相关法计算得到了春、夏季塔中10m高度的动量、感热和潜热通量。结果表明,塔中地区地表热量输送以感热输送为主。春季每天的最大感热通量变化范围为120—320W·m^-2,月平均值为220W·m^-2;夏季最大感热通量的变化范围为140—340W·m^-2,月平均值为230W·m^-2。感热通量值在夜间为负,白天为正,符号的改变出现在日出、日落前后。夏季潜热通量最大值一般为20—60W·m^-2,平均值为27W·m^-2,潜热通量比感热通量小一个量级。春季动量通量的平均值为-0.063W·m^-2,夏季动量通量的平均值为-0.091W·m^-2。日变化规律比较明显,日出后,动量向下传输增大,在09-10时(地方时)出现一个最大值,随后动量向下传输并开始减小。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号