首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
定常情况下,本文利用球坐标系（λ,φ,r）的动力、热力学方程导出三维速度场（vλ,vφ,vr）的动力系统.这种包括摩擦力和热传导的不可压缩大气运动的动力系统,无论从定性上还是从定量上都能描述由赤道和极地间的加热不均匀造成的三圈环流.定性上表明,在北半球经向速度vφ和纬向速度vλ同符号,地表刮北风（vφ<0）和刮东风（vλ<0）相对应,刮南风（vφ>0）和刮西风（vλ>0）相对应.在南半球,经向速度vφ和纬向速度vλ符号相反,刮北风（vφ<0）和刮西风（vλ>0）相对应,刮南风（vφ>0）和刮东风（vλ<0）相对应.定量分析表明球面上的压力场p可以用球调和函数plm（sinφ）cosmλ来表达.当取l=6,m=0时即可导出三圈环流.在经圈剖面（φ,r）上,地表的φ1=±56°和φ2=±28°左右,以及赤道是速度场的奇点,它们都是鞍点,说明在副热带是下沉运动,在中纬度是上升运动,这正是三圈环流中的Ferrel环流的特征.这样经向速度vφ和纬向速度vλ相联系,经向速度vφ又和垂直速度相联系,那么三圈环流的三维速度场就构成了一个整体.  相似文献   

2.
Nonhydrostatic Atmospheric Normal Modes on Beta-Planes   总被引：1，自引：0，他引：1  
--To facilitate the understanding of nonhydrostatic effect in global and regional nonhydrostatic models, the normal modes of a nonhydrostatic, stratified, and compressible atmosphere are studied using Cartesian coordinates on midlatitude and equatorial #-planes. The dynamical equations without forcing and dissipation are linearized around the basic state at rest, and solved by using the method of separation of variables. An eigenvalue-eigenfunction problem is formulated, consisting of the horizontal and vertical structure equations with suitable boundary conditions. The wave frequency and the separation parameter, referred to as "equivalent height," appear in both the horizontal and vertical characteristic equations as a coupled problem, unlike the hydrostatic case. Therefore, the nonhydrostatic equivalent height depends not only on the vertical modal scale, as in the hydrostatic case, but also on the zonal and meridional modal scales. Numerical resu lts on the dispersion relations are presented for an isothermal atmosphere. Three kinds of normal modes, namely acoustic, gravity, and Rossby modes, are solved and compared with the corresponding global solutions. Nonhydrostatic effects are studied in terms of normal modes in a wide range of wavelengths from small to planetary scales. It is demonstrated that Rossby modes are hardly affected by nonhydrostatic effects regardless of wavelengths. However, nonhydrostatic effects on gravity modes become significant for smaller horizontal and deeper vertical scales of motion. The equivalent height plays a particularly important role in evaluating nonhydrostatic effects of normal modes on the equatorial #-plane, because the equivalent height appears in the scaling of meridional distance variable of the eigenfunctions. The implementation of nonhydrostatic normal mode analysis on high-resolution numerical modeling is also discussed.  相似文献   

3.
1 Introduction Antarctic Oscillation is a major mode of Southern Hemispheric (SH) extratropical atmospheric circula- tion. The SH Annual Mode represents a zonally sym- metric exchange of mass between polarward of 60°S and 40°S[1], therefore AAO indicates both the intensity of circumpolar low and zonal mean west wind at mid-high latitudes in SH. Positive phase of AAO tends to deepen circumpolar low and enhance west wind at mid-high latitudes in SH. Many studies show that AAO has bar…  相似文献   

4.
— The mean zonal velocity in the atmosphere is taken as being created continually by the global scale Hadley circulation produced by the differential solar heating through the balance between the Coriolis effect and vertical diffusion, and not by conservation of absolute momentum. Hence a proper determination of the diffusion coefficient becomes the key to the solution of the zonal flow problem. In this study we take the flow field as composed of a primary global scale Hadley circulation, and a secondary flow created by the convergences of the eddy transports of heat and momentum and surface friction, which give rise to the classical three cell structure of the meridional circulation but which only modifies the zonal velocity distribution slightly.¶Finally, we use the equilibrium solution of the perturbation potential vorticity equation to obtain the eddy transports of momentum and heat, with the zonal velocity given by the primary Hadley flow as the basic flow, and we found that they are close to the statistically observed values, demonstrating that the system can maintain itself.  相似文献   

5.
A mean climatology is studied to examine atmospheric circulation characteristics to assess the wintertime (December, January, February and March - DJFM) synoptic weather system affecting northern India. The main objective is to study the mean circulation and mean energetics distribution pertaining to the winter season, which are embedded with an eastward moving synoptic weather system in westerlies, called Western Disturbances (WDs). Forty years (1958–1997) of uninitialized daily re-analysis data of the National Center for Environmental Prediction - National Center for Atmospheric Research (NCEP- NCAR, henceafter NCEP), U.S. has been considered for this study. Winter circulations are considered over the domain 15°S–45°N and 30°E–120°E. This domain is considered particularly to illustrate the impact of wintertime synoptic weather system Western Disturbances (WDs), which travel towards the east over the western Himalayas during winter and yield an enormous amount of precipitation in the form of snow. Large-scale balances of kinetic energy, vorticity, angular momentum, heat and moisture budget terms are analyzed. The main findings of the study show that strong rising motion in the extratropical region brings a significant amount of precipitation over the region of study. Also, horizontal flux of kinetic energy converges in the tropical region and diverges over the extratropical region. It is seen that both the zonal and meridional component of kinetic energy contributes to the production of kinetic energy in the upper troposphere. Vorticity budget shows that wintertime circulation over the western Himalayas is characterized by a negative generation of vorticity. The relative and planetary vorticity advection contributes to the horizontal transport of vorticity. The moisture flux transported into the region shows that in the middle tropospheric levels moisture undergoes phase transformation due to turbulent exchange and hence releases latent heat.  相似文献   

6.
Meridional circulation in stellar convection zones is not generally well observed, but may be critical for the workings of MHD dynamos operating in these domains. Coriolis forces from differential rotation play a large role in determining what the meridional circulation is. Here, we consider the question of whether a stellar differential rotation that is constant on cylinders concentric with the rotation axis can drive a meridional circulation. Conventional wisdom says that it can not. Using two related forms of the governing equations that respectively estimate the longitudinal components of the curl of the meridional mass flux and the vorticity, we show that such differential rotation will drive a meridional flow. This is because to satisfy anelastic mass conservation, non-spherically symmetric pressure contours must be present for all differential rotations, not just ones that depart from constancy on cylinders concentric with the rotation axis. Therefore, the fluid is always baroclinic if differential rotation is present. This is because, in anelastic systems, the perturbation pressure must satisfy a Poisson type equation, as well as an equation of state and a thermodynamic equation. We support our qualitative reasoning with numerical examples, and show that meridional circulation is sensitive to the magnitude and form of departures from rotation constant on cylinders. The effect should be present in 3D global anelastic convection simulations, particularly those for which the differential rotation driven by global convection is nearly cylindrical in profile. For solar-like differential rotation, Coriolis forces generally drive a two-celled circulation in each hemisphere, with a second, reversed flow at high latitudes. For solar like turbulent viscosities, the meridional circulation produced by Coriolis forces is much larger than observed on the Sun. Therefore, there must be at least one additional force, probably a buoyancy force, which opposes the meridional flow to bring its amplitude down to observed values.  相似文献   

7.
Based on numerical simulation, this study explored the characteristics and interactions of surface sensible heating and atmospheric latent heating over the main part of the Tibetan Plateau, i.e., terrain at elevations >2 km in summer. The impacts of these two types of heating on local vertical motion and monsoonal meridional circulation were compared. Theoretical analysis and numerical experimentation demonstrated that by changing the configuration of the upper-tropospheric air temperature and circulation, the two types of heating could generate both minimum absolute vorticity and abnormal potential vorticity forcing near the tropopause, enhance the meridional circulation of the Asian summer monsoon, and produce an eastward- propagating Rossby wave train within the mid-latitude westerly flow. Consequently, the manifestations of these features were shown to influence the circulation of the Northern Hemisphere.  相似文献   

8.
Summary The time-dependent primitive equations for a shallow homogeneous ocean with a free surface are solved for a bounded basin on the sphere, driven by a steady zonal wind stress and subject to lateral viscous dissipation. These are the vertically integrated equations for a free-surface model, and are integrated to 60 days from an initial state of rest by an explicit centered-difference method with zero-slip lateral boundary conditions. In a series of comparative numerical solutions it is shown that at least a 2-deg resolution is needed to resolve the western boundary currents adequately and to avoid undue distortion of the transient (Rossby waves. The -plane formulation is shown to be an adequate approximation for the mean circulation in the lower and middle latitudes, but noticeably intensifies the transports poleward of about 50 deg and both slows and distorts the transients in the central basin. The influence of the (southern) zonal boundary on the transport solutions is confined to the southernmost gyre, except in the region of the western boundary currents where its influence spreads to the northern edge of the basin by 30 days. The total boundary current transport is shown to be approximately proportional to the zonal width of the basin and independent of the basin's (uniform) depth, while the elevation of the free water surface is inversely proportional to the basin depth, in accordance with linear theory. The introduction of bottom friction has a marked damping effect on the transient Rossby waves, and also reduces the maximum boundary-current transport. The solutions throughout are approximately geostrophic and are only slightly nonlinear.The root-mean-square (rms) transport variability during the period 30 to 60 days is concentrated in the southwest portion of the basin through the reflection of the transient Rossby waves from the western shore and has a maximum corresponding to an rms current variability of about 3 cm sec–1. The transport variabilities are about 10 percent of the mean zonal transport and more than 100 percent of the mean meridional transport over a considerable region of the western basin (outside the western boundary current regime). Some 99 percent of the total kinetic energy is associated with the zonal mean and standing zonal waves, which are also responsible for the bulk of the meridional transport of zonal angular momentum. Although the transient Rossby waves systematically produce a momentum flux convergence at the latitude of the maximum eastward current, much in the manner of their atmospheric counterparts, this is only a relatively small contribution to the zonal oceanic momentum balance; the bulk of the mean zonal stress is here balanced by a nearly stationary net pressure torque exerted against the meridional boundaries by the wind-raised water. In an ocean without such boundaries the role of the transient circulations may be somewhat more important.  相似文献   

9.
two-dimensional time dependent model of the stratosphere incorporating the major interactions between radiative-photochemical and dynamical processes is described. The main prognostic equations considered are the thermodynamic equation and the general conservation equation for the minor chemical constituents representing the odd oxygen (O x =O+('D)+O3), odd hydrogen (HO x =HO+HO2), N2O, odd nitrogen (NO x =NO+NO2+HNO3), CF2Cl2, CFCl3 and odd chlorine (Cl x =Cl+ClO+HCl). The zonal wind and mean meridional circulations are determined diagnostically by the integration of the thermal wind equation and the stream function equation in the meridional plane espectively. The large scale eddy processes are parameterized in terms of zonal mean quantities using the generalized diffusion formulation on a sloping surface. The radiative heating and cooling and the hotochemical sources and sinks are incorporated in a form which allows for the major interactions among the minor trace constituents, temperature and mean circulation.Two integrations consisting of natural stratosphere and a stratosphere contaminated by the chlorofluoromethanes through lower boundary fluxes are carried out for 23 model years by changing the declination of the sun every day and using 6-hour time step. The model simulations of temperature, mean circulation, ozone, HO x , N2O and NO x in the meridional plane for the normal stratosphere, show satisfactory agreement with the available observations. Based on the results of second integration it is found that the injection of chlorofluoromethanes in the atmosphere at the estimated current production rates can lead to significant changes in the meridional distribution of ozone, temperature and NO x in the middle and upper stratosphere. The results also indicate that the percentage total ozone depletion increases from tropics to high latitudes and from summer to winter high latitudes. Also discussed are the results of additional experiments incorporating the reaction of HO2 with NO and the reactions involving ClNO3.  相似文献   

10.
Three-dimensional structures of the ionospheric dynamo currents are examined using the neutral winds in a general circulation model of the middle atmosphere at Kyushu University. A quasi-three-dimensional ionospheric dynamo model is constructed assuming an infinite parallel conductivity in the ionosphere. This model is able to simulate both the equatorial electrojet and the global Sq current system successfully. The simulated results reveal that the equatorial electrojet is confined in quite narrow latitudes around the equator accompanied with meridional current circulations and satisfies a non-divergent structure mainly within the E region. A vertically stratified double layered structure is seen in the east–west current density near the focus latitude of the global Sq current system. It is shown that the stratified structure mainly consists of the east–west Hall current associated with the eastward wind of zonal wavenumbers 1 and 2 in the lower altitudes and the westward wind of zonal wavenumber 2 in the upper altitudes. The day-to-day variation of the neutral winds can significantly vary the induced ionospheric dynamo current system, which is recognized as changes of the focus latitude and/or the maximum value of the equatorial electrojet.  相似文献   

11.
Summary The mean zonal and meridional wind components and the mean mass circulation in different latitudes are discussed in relation to previous studies. Divergence and vertical motion are calculated for various latitude belts. There is evidence for a strong tropical Hadley cell with a temperate latitude indirect circulation during the winter season. During summer, the northern Hadley cell is weaker and displaced poleward; a circulation in the opposite sense appears in equatorial latitudes.The regional patterns of divergence and vertical motion appear related particularly to the position of the subtropical high pressure cells. Subsidence and lower-layer divergence are characteristic of the eastern flanks of the subtropical anticyclones, while the opposite pattern prevails on its western flanks. These longitudinal contrasts appear particularly pronounced during the summer season. The three-dimensional flow pattern in the tropics is illustrated by selected trajectories for the winter and summer seasons.  相似文献   

12.
近60年全球大气环流经向模态的气候变化   总被引：1，自引：0，他引：1       下载免费PDF全文
本文根据1948～2004年NCEP/NCAR 1000 hPa、500 hPa、100 hPa高度场逐月再分析资料，分析了近60年全球大气环流经向模态的气候变化. 结果表明：近60年来第一模态从低层到高层都表现出高纬与低纬地区之间明显的反向变化关系，且随时间有明显的增强趋势. 第一模态位相发生了相反的改变，低纬地区由负距平演变为正距平，高纬地区由正距平演变为负距平. 1000 hPa和500 hPa高度场上的南半球比北半球变化激烈，而100 hPa高度场上的北半球比南半球变化激烈. 第二模态在1000 hPa高度场上，主要表现为南极涛动（AAO）和北极涛动(AO)，且两涛动在年际、年代际尺度上表现出明显的负相关关系；在100 hPa高度场上，主要表现为南北半球高纬度地区之间的反向变化；500 hPa高度场是1000 hPa和100 hPa的一个过渡层次，主要表现出明显的南极涛动(AAO). 第二模态可能是南北半球中高纬环流相互作用的桥梁.  相似文献   

13.
By adopting the idea of three-dimensional Walker, Hadley and Rossby stream functions, the global atmospheric circulation can be considered as the sum of three stream functions from a global per- spective. Therefore, a mathematical model of three-dimensional decomposition of global atmospheric circulation is proposed and the existence and uniqueness of the model are proved. Besides, the model includes a numerical method leading to no truncation error in the discrete three-dimensional grid points. Results also show that the three-dimensional stream functions exist and are unique for a given velocity field. The mathematical model shows the generalized form of three-dimensional stream func- tions equal to the velocity field in representing the features of atmospheric motion. Besides, the vertical velocity calculated through the model can represent the main characteristics of the vertical motion. In sum, the three-dimensional decomposition of atmospheric circulation is convenient for the further in- vestigation of the features of global atmospheric motions.  相似文献   

14.
We present results from the Numerical Spectral Model (NSM), which focus on the temperature environment of the mesopause region where polar mesospheric clouds (PMC) form. The PMC occur in summer and are observed varying on time scales from months to years, and the NSM describes the dynamical processes that can generate the temperature variations involved. The NSM simulates the quasi-biennial oscillation (QBO), which dominates the zonal circulation of the lower stratosphere at equatorial latitudes. The modeled QBO extends into the upper mesosphere, due to gravity wave (GW) filtering, consistent with UARS zonal wind and TIMED temperature measurements. While the QBO zonal winds are confined to equatorial latitudes, the associated temperature variations extend to high latitudes. The meridional circulation redistributes the QBO energy—and the resulting temperature oscillations away from the equator produce inter-annual variations that can exceed 5 K in the polar mesopause region, with considerable differences between the two hemispheres. The NSM shows that the 30-month QBO produces a 5-year or semi-decadal (SD) oscillation, and stratospheric NCEP data provide observational evidence for that. This SD oscillation extends in the temperature to the upper mesosphere, where it could contribute to the long-term variations of the region.  相似文献   

15.
The NCEP/NCAR R1 reanalysis data are employed to investigate the impact of forced and inertial instability in the lower troposphere over the Arabian Sea on the onset process of Indian summer monsoon(ISM),and to reveal the important role of zonal advection of zonal geostrophic momentum played in the forced unstable convection.Results show that during the ISM onset the zero absolute vorticity contour(??=0)shifts northward due to the strong cross-equatorial pressure gradient in the lower troposphere over southern Arabian Sea.Thus a region with negative absolute vorticity is generated near the equator in the northern hemisphere,manifesting the evident free inertial instability.When a southerly passes through this region,under the influence of friction a lower convergence that facilitates the convection flourishing at the lower latitudes appears to the north of zero absolute vorticity contour.However,owing to such a traditional inertial instability,the convection is confined near the equator which does not have direct influence on the ISM onset.On the contrary in the region to the north of the zero absolute vorticity contour and to the south of the low pressure center near the surface,although the atmosphere there is inertially stable,the lower westerly jet can develop and bring on the apparent zonal advection of zonal geostrophic momentum.Both theoretical study and diagnosing analysis present that such a zonal advection of geostrophic momentum is closely associated with the zonal asymmetric distribution of meridional land-sea thermal contrast,which induces a convergence center near and further north of the westerly jet in the lower troposphere over the southwestern coast of the Indian Peninsula,providing a favorable lower circulation for the ISM onset.It illustrates that the development of convection over the Arabian Sea in late spring and early summer is not only due to the frictional inertial instability but also strongly affected by the zonal asymmetric distribution of land-sea thermal contrast.Moreover,before the ISM onset due to the eastward development of the South Asian High(SAH)in the upper troposphere,high potential vorticity is transported to the region over the Arabian Sea.Then a local trumpet-shaped stream field is generated to cause the evident upper divergence-pumping effect which favors the ISM onset.When the upper divergence is vertically coupled with the lower convergence resulted from the aforementioned forced unstable convection development near the southwestern coast of Indian Peninsula,the atmospheric baroclinic unstable development is stimulated and the ISM onset is triggered.  相似文献   

16.
选用每天12∶00UTC时次的逐日ERA-Interim再分析资料,根据transformed Eulerian-mean(TEM)方程通过积分剩余速度珔v*,研究了1979—2011年间Brewer-Dobson(BD)环流的时空演变规律.并将其与downward control(DC)原理研究的结果进行比较,同时还探讨了平流层温度与BD环流之间的相互联系.结果表明,由TEM方程通过积分剩余速度珔v*估算的BD环流与利用DC原理估算的环流相比较,在热带地区的形势更加明显.环流在热带对流层中上层上升至平流层中下层,最高可达1hPa等压面附近.然后在热带外向极向下运动,最后在中高纬度下沉回到对流层.BD环流的上升中心及质量通量均随季节的变化产生变动,环流在冬半球的形势显著地强于夏半球.在春季和秋季期间,环流呈现出南北两半球的对称形势.从全球尺度物质输送的角度来看,在过去的33a间平流层BD环流的长期变化趋势是减弱的,且在平流层中下层减弱是明显的.环流的减弱趋势与纬向平均温度的长期变化趋势相匹配.  相似文献   

17.
This study is diagnostic in character and entails a systematic investigation of the kinematic parameters of the mean monthly 300 mb horizontal-motion field for the global tropical region during January, April, and October. The kinematic parameters analysed include the horizontal velocity divergence, relative vorticity, velocity potential, and stream function. The magnitude of relative voritcity has been observed to be nearly equal to that of the horizontal velocity divergence in the equatorial latitudes and an order of magnitude larger in the subtropics. The magnitudes of?u/?x, ?v/?y, and horizontal velocity divergence are comparable and vary very little with latitude, while that of?u/?y is an order of magnitude larger than that of?v/?x, confirming its dominance in the computation of relative vorticity. The latitudinal position of zero absolute vorticity isopleth fluctuates within 5 degrees of the equator and is confined more to the winter hemisphere than to the summer hemisphere. The rotational component of wind is large compared with the divergent wind component. The magnitude of the rotational vector wind is four times that of the divergent vector wind in the equatorial latitudes and an order of magnitude larger in the subtropical latitudes.  相似文献   

18.
A continuing goal in the diagnostic studies of the atmospheric general circulation is to estimate various quantities that cannot be directly observed. Evaluation of all the dynamical terms in the budget equations for kinetic energy, vorticity, heat and moisture provide estimates of kinetic energy and vorticity generation, diabatic heating and source/sinks of moisture. All these are important forcing factors to the climate system. In this paper, diagnostic aspects of the dynamics and energetics of the Asian summer monsoon and its spatial variability in terms of contrasting features of surplus and deficient summer monsoon seasons over India are studied with reanalysis data sets. The daily reanalysis data sets from the National Centre for Environmental Prediction/National Centre for Atmospheric Research (NCEP/NCAR) are used for a fifty-two year (1948–1999) period to investigate the large-scale budget of kinetic energy, vorticity, heat and moisture. The primary objectives of the study are to comprehend the climate diagnostics of the Asian summer monsoon and the role of equatorial convection of the summer monsoon activity over India.It is observed that the entrance/exit regions of the Tropical Easterly Jet (TEJ) are characterized by the production/destruction of the kinetic energy, which is essential to maintain outflow/inflow prevailing at the respective location of the TEJ. Both zonal and meridional components contribute to the production of kinetic energy over the monsoon domain, though the significant contribution to the adiabatic generation of kinetic energy originates from the meridional component over the Bay of Bengal in the upper level and over the Somali Coast in the low level. The results indicate that the entire Indian peninsula including the Bay of Bengal is quite unstable during the summer monsoon associated with the production of vorticity within the domain itself and maintain the circulation. The summer monsoon evinces strong convergence of heat and moisture over the monsoon domain. Also, considerable heat energy is generated through the action of the adiabatic process. The combined effect of these processes leads to the formation of a strong diabatic heat source in the region to maintain the monsoon circulation. The interesting aspect noted in this study is that the large-scale budgets of heat and moisture indicate excess magnitudes over the Arabian Sea and the western equatorial Indian Ocean during surplus monsoon. On the other hand, the east equatorial Indian Ocean and the Bay of Bengal region show stronger activity during deficient monsoon. This is reflected in various budget terms considered in this study.  相似文献   

19.
Early in the1970s,Madden and Julian[1,2]first re-vealed the existence of eastward propagation of the intraseasonal(or30-60-day)oscillation(ISO)in the zonal wind fields over the tropics.Later on,the northward propagation of the ISO activities over the Indian Ocean was discovered by Yasunari[3,4],which accounts for the seasonal variations of the Indian summer monsoon.Similar studies were carried out on the low-frequency oscillation activities related to the ENSO and the East Asian summer …  相似文献   

20.
Mechanisms of the meridional heat transport in the Southern Ocean   总被引：1，自引：0，他引：1  
The Southern Ocean (SO) transports heat towards Antarctica and plays an important role in determining the heat budget of the Antarctic climate system. A global ocean data synthesis product at eddy-permitting resolution from the Estimating the Circulation and Climate of the Ocean, Phase II (ECCO2) project is used to estimate the meridional heat transport (MHT) in the SO and to analyze its mechanisms. Despite the intense eddy activity, we demonstrate that most of the poleward MHT in the SO is due to the time-mean fields of the meridional velocity, V, and potential temperature, θ. This is because the mean circulation in the SO is not strictly zonal. The Antarctic Circumpolar Current carries warm waters from the region south of the Agulhas Retroflection to the lower latitudes of the Drake Passage and the Malvinas Current carries cold waters northward along the Argentinian shelf. Correlations between the time-varying fields of V and θ (defined as transient processes) significantly contribute to the horizontal-gyre heat transport, but not the overturning heat transport. In the highly energetic regions of the Agulhas Retroflection and the Brazil-Malvinas Confluence the contribution of the horizontal transient processes to the total MHT exceeds the contribution of the mean horizontal flow. We show that the southward total MHT is mainly maintained by the meridional excursion of the mean geostrophic horizontal shear flow (i.e., deviation from the zonal average) associated with the Antarctic Circumpolar Current that balances the equatorward MHT due to the Ekman transport and provides a net poleward MHT in the SO. The Indian sector of the SO serves as the main pathway for the poleward MHT.  相似文献   

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

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