首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 296 毫秒
1.
The semi-Lagrangian advection scheme is implemented on a new quasi-uniform overset (Yin-Yang) grid on the sphere. The Yin-Yang grid is a newly developed grid system in spherical geometry with two perpendicularly-oriented latitude-longitude grid components (called Yin and Yang respectively) that overlapp each other, and this effectively avoids the coordinate singularity and the grid convergence near the poles. In this overset grid, the way of transferring data between the Yin and Yang components is the key to maintaining the accuracy and robustness in numerical solutions. A numerical interpolation for boundary data exchange, which maintains the accuracy of the original advection scheme and is computationally efficient, is given in this paper. A standard test of the solid-body advection proposed by Williamson is carried out on the Yin-Yang grid. Numerical results show that the quasi-uniform Yin-Yang grid can get around the problems near the poles, and the numerical accuracy in the original semi-Lagrangian scheme is effectively maintained in the Yin-Yang grid.  相似文献   

2.
Designed for grid point systems, the traditional semi-Lagrangian semi-implicit scheme is not mass-conserving and can lead to significant solution errors. In the present study, a finite-volume semi-Lagrangian semi-implicit scheme (hereafter “FVSLSI”) is designed for the Yin-Yang mesh and tested in a barotropic shallow water model in the spherical coordinate system. Three test cases, i.e. the advection of a solid body, a steady state nonlinear zonal geostrophic flow and the deformation flow, are simulated to compare the performance of the FVSLSI with that of the traditional semi-Lagrangian scheme (hereafter “SL”) from perspectives of shape preservation, mass conservation, normalized bias, and convergence rate. Results indicate that the FVSLSI performs better than the SL in mass conservation and shape preservation. The bias by the FVSLSI is smaller than that by the SL, while the rate of convergence by the FVSLSI is larger than that by the SL. The FVSLSI also allows large time step. Therefore, the FVSLSI is suggested to be distributed to communities that are developing atmospheric/oceanic models.  相似文献   

3.
ABSTRACT The Global/Regional Assimilation and PrEdiction System (GRAPES) is the newgeneration numerical weather predic- tion (NWP) system developed by the China Meteorological Administration. It is a fully compressible non-hydrostatical global/regional unified model that uses a traditional semi-Lagrangian advection scheme with cubic Lagrangian interpola tion (referred to as the SL_CL scheme). The SL_CL scheme has been used in many operational NWP models, but there are still some deficiencies, such as the damping effects due to the interpolation and the relatively low accuracy. Based on Reich's semi-Lagrangian advection scheme (referred to as the R2007 scheme), the Re_R2007 scheme that uses the low- and high-order B-spline function for interpolation at the departure point, is developed in this paper. One- and two-dimensional idealized tests in the rectangular coordinate system with uniform grid cells were conducted to compare the Re..R2007 scheme and the SL_CL scheme. The numerical results showed that: (1) the damping effects were remarkably reduced with the Re_R2007 scheme; and (2) the normalized errors of the Re_R2007 scheme were about 7.5 and 3 times smaller than those of the SL_CL scheme in one- and two-dimensional tests, respectively, indicating the higher accuracy of the Re..R2007 scheme. Furthermore, two solid-body rotation tests were conducted in the latitude-longitude spherical coordinate system with non uniform grid cells, which also verified the Re_R2007 scheme's advantages. Finally, in comparison with other global advection schemes, the Re_R2007 scheme was competitive in terms of accuracy and flow independence. An encouraging possibility for the application of the Re_R2007 scheme to the GRAPES model is provided.  相似文献   

4.
This paper further explores the estimating and expressing of dynamic balance constraints using statistical methods in GRAPES-3DVAR (Version GM). Unlike the single-level scheme which only considers the coupling between mass and wind at one level, the multi-level scheme considers the coupling between their vertical profiles and calculates the balanced mass field at each layer using the rotational wind at all model levels. A reformed ridge regression method is used in the new scheme to avoid the multicollinearity problem and reduce the noises caused by unbalanced mesoscale disturbances. The results of numerical experiments show that the new scheme can get more reasonable vertical mass field, reduce the magnitude of the adjustment by the initialization, and improve the potential temperature analysis performance. Furthermore, the results of forecast verification in January (winter) and July (summer) both confirm that the new scheme can significantly improve the temperature forecast accuracy and bring slight positive effects to the pressure and wind forecast.  相似文献   

5.
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.  相似文献   

6.
A 3D dynamic core of the non-hydrostatic model GRAPES(Global/Regional Assimilation and Prediction System) is developed on the Yin-Yang grid to address the polar problem and to enhance the computational efficiency. Three-dimensional Coriolis forcing is introduced to the new core, and full representation of the Coriolis forcing makes it straightforward to share code between the Yin and Yang subdomains. Similar to that in the original GRAPES model, a semi-implicit semi-Lagrangian scheme is adopted for temporal integration and advection with additional arrangement for cross-boundary transport. Under a non-centered second-order temporal and spatial discretization, the dry nonhydrostatic frame is summarized as the solution of an elliptical problem. The resulting Helmholtz equation is solved with the Generalized Conjugate Residual solver in cooperation with the classic Schwarz method. Even though the coefficients of the equation are quite different from those in the original model, the computational procedure of the new core is just the same. The bi-cubic Lagrangian interpolation serves to provide Dirichlet-type boundary conditions with data transfer between the subdomains. The dry core is evaluated with several benchmark test cases, and all the tests display reasonable numerical stability and computing performance. Persistency of the balanced flow and development of both the mountain-induced Rossby wave and Rossby–Haurwitz wave confirms the appropriate installation of the 3D Coriolis terms in the semi-implicit semi-Lagrangian dynamic core on the Yin-Yang grid.  相似文献   

7.
A series of data assimilation and forecast test have been carried out with a hemispheric spectralmodel(T42L9H).It is found that the numerical scheme for determining hemispheric initial wind isimportant to data assimilation and forecast.An inappropriate scheme may cause computationalsources of divergence near the equator,which are responsible for the spurious strong precipitationand corresponding latent heat release.Obviously,this problem differs from either thehemispheric/global domain effect or the tropical data effect pointed by Dalley et al.(1981).Basedon the previous studies,the new scheme of divergence and vorticity correction is presented,andthe difference with other schemes and its effects on the data assimilation are discussed against thecontrol test.Preliminary tests have shown that the new divergence correction scheme proposed in thispaper may be a preferable choice to overcome the initial computational errors in the hemisphericdata assimilation.  相似文献   

8.
Axisymmetric bogus vortexes at sea level are usually used in the traditional bogus data assimilation (BDA) scheme. In the traditional scheme, the vortex could not accurately describe the specific characteristics of a typhoon, and the evolving real typhoon is forced to unreasonably adapt to this changeless vortex. For this reason, an asymmetrical typhoon bogus method with information blended from the analysis and the observation is put forward in this paper, in which the impact of the Subtropical High is also taken into consideration. With the fifth-generation Penn State/NCAR Mesoscale Model (MM5) and its adjoint model, a four-dimensional variational data assimilation (4D-Var) technique is employed to build a dynamic asymmetrical BDA scheme to assimilate different asymmetrical bogus vortexes at different time. The track and intensity of six summer typhoons much influenced by the Subtropical High are simulated and the results are compared. It is shown that the improvement in track simulation in the new scheme is more significant than that in the traditional scheme. Moreover, the periods for which the track cannot be simulated well by the traditional scheme can be improved with the new scheme. The results also reveal that although the simulated typhoon intensity in the new scheme is generally weaker than that in the traditional scheme, this trend enables the new scheme to simulate, in the later period, closer-to-observation intensity than the traditional scheme. However, despite the fact that the observed intensity has been largely weakened, the simulated intensity at later periods of the BDA schemes is still very intensive, resulting in overly development of the typhoon during the simulation. The limitation to the simulation effect of the BDA scheme due to this condition needs to be further studied.  相似文献   

9.
1. IntroductionRecently, the improvement of accuracy in the out-puts of a numerical mesoscale model by the physi-cal dissipative technique is reached (Liu et al., 2002;Liu and Liu, 2003). The effect of improvement ofthis technique differs not only from model to model,but from scheme to scheme of parameterization em-ployed in the same model. The rapid developmentof the computer technology makes possible the com-plicated numerical experiments by a model with highresolution and multiple domains …  相似文献   

10.
Canopy interception of incident precipitation, as a critical component of a forest's water budget, can affect the amount of water available to the soil, and ultimately vegetation distribution and function. In this paper, a statistical-dynamic approach based on leaf area index and statistical canopy interception is used to parameterize the canopy interception process. The statistical-dynamic canopy interception scheme is implemented into the Community Land Model with dynamic global vegetation model (CLM-DGVM) to improve its dynamic vegetation simulation. The simulation for continental China by the land surface model with the new canopy interception scheme shows that the new one reasonably represents the precipitation intercepted by the canopy. Moreover, the new scheme enhances the water availability in the root zone for vegetation growth, especially in the densely vegetated and semi-arid areas, and improves the model's performance of potential vegetation simulation.  相似文献   

11.
The second author studied the nonlinear stability of N-layer quasi-geostrophic flow subject to perturbations of parameters and initial data, and established the stability criteria for the flow in question, which involve finding out the lowest eigenvalue of an elliptic boundary value problem.In this paper when the domain is a periodic zonal channel, a formula of the lowest eigenvalue is established, which is useful for further studies and practical applications.  相似文献   

12.
由中国气象科学研究院,中国科学院地球环境研究所,中国科学院大气物理研究所,北京大学等单位参加,周秀骥院士主持的国家自然科学基金委员会重大项目“中国地区树轮和千年气候变化研究”于2008年批准立项,执行时间为2009年1月-2012年12月。本重大项目拟从建立中国典型气候区域的树轮气候网络入手,重建中国区域千年气候变化时间序列;通过非线性分析诊断,  相似文献   

13.
SAMIL大气环流模式海面湍流通量参数化方案研究   总被引:1,自引:1,他引:0  
王自强  缪启龙  高志球 《大气科学》2010,34(6):1155-1167
将中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室发展的大气环流谱模式SAMIL_R42L26 2.08中的海气通量参数化方案——Louis方案, 与新发展的一种新型近海层湍流通量参数化方案——LGLC方案进行比较和分析。离线测试结果表明, LGLC方案计算的通量结果与观测数据具有更好的一致性, 且由于其区分热力粗糙度和动力粗糙度, 使得对热量通量的计算更加准确。在线测试则证明, 引入LGLC方案的SAMIL模式对洋面风应力、 感热通量、 潜热通量和降水率的模拟能力有了进一步的提高, 尤其对北半球夏季印度季风和南海季风区的降水改善明显。  相似文献   

14.
大气环流模式和耦合模式模拟的降水-海温关系之比较   总被引:9,自引:6,他引:3  
李博  周天军  吴春强 《大气科学》2009,33(5):1071-1086
本文讨论了气候系统耦合模式FGOALS_s模拟的局地降水和海温的关系, 并通过与单独大气模式SAMIL的AMIP试验结果进行对比分析, 考察了海气耦合过程对局地降水和海温关系模拟的影响。结果表明, 耦合模式FGOALS_s和单独大气模式SAMIL在模拟局地降水和海温关系上各有优势。在赤道中东太平洋地区, 观测中局地降水和海温的关系是海洋强迫为主, FGOALS_s模拟的海洋对大气的强迫比观测偏弱, 因此, SAMIL相对于FGOALS_s更有优势。在西北太平洋东部地区, 观测中夏秋季节降水和海温的关系是大气强迫为主, 由于考虑了海气相互作用过程, FGOALS_s对降水和海温关系的模拟能力要优于单独的大气模式。此外, 由于大气模式SAMIL的云参数化方案导致的模拟偏差, 在赤道中东太平洋地区 (9~11月), 降水增加时入射的短波辐射通量也是增加的, 并且这种模拟的偏差在耦合后仍然保留了下来, 导致了与观测不符的云[CD*2]辐射反馈过程。因此, 改进大气模式的云参数化方案是未来工作重点之一。此外, 分析发现耦合模式模拟的潜热通量的变化过分依赖海气湿度差。  相似文献   

15.
SAMIL模式中Tiedtke 积云对流方案对热带降水模拟的影响   总被引:2,自引:1,他引:1  
刘琨  刘屹岷  吴国雄 《大气科学》2010,34(1):163-174
目前, 大多数全球耦合模式及大气环流模式在降水模拟中普遍存在不同程度的“热带偏差”问题, 中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室所发展的全球大气环流谱模式SAMIL-R42L26也存在这一现象, 主要表现在SPCZ (南太平洋辐合区) 降水过强且过分东伸、 赤道附近降水偏少等方面。本文通过修改SAMIL中的积云对流方案有效地削弱了这一偏差, 并进一步探讨其原因, 发现对流方案修改后, 改变了对流层低层至地面的温度分布状态, 进而影响了风速及散度场的模拟, 最终通过垂直速度的调整反作用于对流过程。比较修改前后对流过程云底质量通量, 发现修改后的方案主要通过削弱浅对流来提高热带降水的模拟性能。  相似文献   

16.
海气耦合模式FGOALS_gl模拟的水汽和云辐射反馈过程   总被引:1,自引:0,他引:1  
本文分析了中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室(LASG)发展的快速耦合模式FGOALS_gl对低纬太平洋区域水汽温室效应(Ga)和云辐射强迫特征的模拟能力,讨论了模拟偏差的成因.结果表明,FGOALS_gl能合理再现Ga、云辐射强迫的气候态空间分布特征,但也存在明显的偏差.模式低估了冷...  相似文献   

17.
    
Arnol′d′s second nonlinear stability criterion for motions governed by a general multilayer quasi-geostrophic model is established. The model allows arbitrary density jumps and layer thickness, and at the top and the bottom of the fluid, the boundary condition is either free or rigid. The criterion is obtained by the establishment of the upper bounds of disturbance energy and potential enstrophy in terms of the initial disturbance field. Project supported by the State Education Commission of China Project supported by the National Natural Science Foundation of China and LASG, IAP, the Chinese Academy of Sciences  相似文献   

18.
积云对流和云物理过程调整对气候模拟的影响   总被引:1,自引:1,他引:0  
本文利用中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室 (LASG/IAP) 的大气环流谱模式SAMIL, 结合观测与政府间气候变化专门委员会第四次评估报告 (IPCC AR4) 大气模式集合平均结果, 以大气辐射通量为例, 诊断分析了物理过程调整前后模式对气候模拟的影响。旧版本SAMIL对大气辐射通量的模拟存在较大偏差, 经过大气辐射过程、 积云对流和诊断云等物理过程的调整后, 新版本SAMIL模拟的全球辐射通量的年平均结果与观测的偏差大幅减小, 其中大气顶能量收支的年变化及其平均值与观测更为接近。在积云对流方案调整基础上, 通过对诊断云物理方案的进一步调整, 新版本SAMIL对云物理量模拟更为合理, 在赤道辐合带等区域, 在很大程度上克服了单一积云对流物理过程调整引起的云宏观和微观属性不匹配问题, 能模拟出夏季气候平均辐射通量的全球分布特征, 尤其在东亚区域有较好的模拟能力。研究还表明, 在热带和副热带对流活跃区域, 当前SAMIL中积云对流过程偏差对辐射通量的模拟偏差有很大影响, 而模式中较为简单的诊断云方案也会将云宏观物理量模拟偏差带入云微观量模拟中, 也是主要偏差源之一。本文结果表明, 要继续提高SAMIL的模拟性能, 急需更新云物理参数化方案以改进云辐射过程的模拟, 同时也需要有针对性的研究积云对流和云物理过程之间相互作用, 并作进一步协同调整。  相似文献   

19.
正1.Introduction The Atmosphere Profiling Synthetic Observation System(APSOS)is the first ground-based facility for profiling atmospheric variables and multiple constituents in the whole(neutral)atmosphere from the surface up to the lower thermosphere.It enables simultaneous observations and extensive studies of the atmospheric vertical structure and constituent  相似文献   

20.
积云参数化和分辨率对MJO数值模拟的影响   总被引:3,自引:1,他引:2  
用中国科学院大气物理研究所发展的一个大气环流模式,使用不同的积云参数化方案和分辨率进行了6个模拟试验,考察了积云参数化方案和模式分辨率对热带大气季节内振荡(MJO)模拟的影响。结果显示:积云参数化方案和分辨率都会影响MJO的模拟。但积云参数化方案决定了模式对MJO模拟的基本能力,决定了模拟的MJO的基本特征。分辨率的变化并不能使模拟的MJO发生本质的改变,分辨率的作用更多的是对MJO的模拟起一定的调制作用,而这种调制作用又受到积云参数化方案的制约。在改进积云参数化方案的基础上提高模式的分辨率会在某些方面改善MJO的模拟。但是分辨率的提高需要同时提高水平分辨率和垂直分辨率,单独提高水平分辨率会降低模式模拟MJO的能力,引入更多的小尺度的高频扰动。非绝热加热垂直廓线对模式模拟MJO有重要的影响,而非绝热加热廓线很大程度上取决于所使用的积云参数化方案。模式水平分辨率的变化不会对加热廓线的结构产生明显的影响,而垂直分辨率的变化会对加热廓线的垂直结构产生一定的调制作用,进而对模拟的MJO起到调制作用。  相似文献   

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

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