首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
We report on our implementation of EULAG as a dynamical core in the Community Atmospheric Model (CAM). EULAG is a non-hydrostatic, parallel computational model for all-scale geophysical flows. EULAG’s name derives from its two computational options: EULerian (flux form) or semi-LAGrangian (advective form). The model combines nonoscillatory forward-in-time (NFT) numerical algorithms with a robust elliptic Krylov solver. A signature feature of EULAG is that it is formulated in generalized coordinates. In particular, this enables grid adaptivity. In total, these features give EULAG novel advantages over the existing dynamical cores in CAM. This paper uses a series of aqua-planet simulations to demonstrate that CAM-EULAG results compare favorably with those from CAM simulations at standard CAM resolution that use current finite volume or Eulerian-spectral dynamical core options. We also show that the grid adaptivity implemented in CAM3-EULAG allows higher resolution in selected regions without causing anomalous behavior such as spurious wave reflection.  相似文献   

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

3.
The Korea Institute of Atmospheric Prediction Systems (KIAPS) began a national project to develop a new global atmospheric model system in 2011. The ultimate goal of this 9-year project is to replace the current operational model at the Korea Meteorological Administration (KMA), which was adopted from the United Kingdom’s Meteorological Office’s unified model (UM) in 2010. The 12-km Korean Integrated Model (KIM) system, consisting of a spectral-element non-hydrostatic dynamical core on a cubed sphere grid and a state-of-the-art physics parameterization package, has been launched in a real-time forecast framework, with initial conditions obtained via the advanced hybrid four-dimensional ensemble variational data assimilation (4DEnVar) over its native grid. A development strategy for KIM and the evolution of its performance in medium-range forecasts toward a world-class global forecast system are described. Outstanding issues in KIM 3.1 as of February 2018 are discussed, along with a future plan for operational deployment in 2020.  相似文献   

4.
GRAPES模式动力框架的长期积分特征   总被引:4,自引:2,他引:4       下载免费PDF全文
通过考虑动量表面拖曳并利用牛顿松弛方法将温度松弛到纬向对称的温度场, 对GRAPES模式的动力框架进行了长期时间积分试验。通过统计分析其积分结果表明:GRAPES模式的动力框架可以模拟出大气环流的基本特征; 随着分辨率的提高, GRAPES模式动力框架的模拟结果显示出收敛的特性; 虽然GRAPES模式动力框架采用的是能量不守衡的半隐式半拉格朗日时间积分方案, 但长期时间积分试验表明其能量基本保持稳定。这些结果显示使用GRAPES模式动力框架作大气环流和气候研究的大气动力框架是可行的, 同时也为进一步改进GRAPES模式动力框架提供了线索和依据。  相似文献   

5.
A global non-hydrostatic atmospheric model, i.e., GRAPES_YY (Global/Regional Assimilation and Prediction System on the Yin–Yang grid), with a semi-implicit semi-Lagrangian (SISL) dynamical core developed on the Yin–Yang grid was coupled with the physical parameterization package of the operational version of GRAPES. A 3.5-yr integration was carried out on an aqua planet to assess the numerical performance of this non-hydrostatic model relative to other models. Specific aspects of precipitation and general circulation under two different sea surface temperature (SST) conditions (CONTROL and FLAT) were analyzed. The CONTROL SST peaked at the equator. The FLAT SST had its maximum gradient at about 20° latitude, giving a broad equatorial SST maximum in the tropics and flat profile approaching the equator. The tropical precipitation showed different propagation features in the CONTROL and FLAT simulations. The CONTROL showed tropical precipitation bands moving eastward with some envelopes of westward convective-scale disturbance. Less organized westward-propagating rainfall cells and bands were seen in the FLAT and the propagation of the tropical wave varied with the SST gradient. The Inter Tropical Convergence Zone (ITCZ), Hadley cell, and westerly jet core were weaker and more poleward as the SST profile flattened from the CONTROL to FLAT. The climatological structures simulated by GRAPES_YY, such as the distribution of precipitation and the large-scale circulation, fell within the bounds from other models. The stronger ITCZ precipitation, accompanied with stronger Hadley cells and convective heating in the CONTROL simulation, may be summed up as a result of stronger parameterized convection and the non-hydrostatic effects in GRAPES_YY. In addition, mechanism of the zonal mean circulation maintaining is analyzed for the different SST patterns referring the transient eddy flux.  相似文献   

6.
Summary The dynamical effect of land surface heterogeneity on heat fluxes in the atmospheric boundary layer (ABL) is investigated using numerical simulations with a non-hydrostatic model over a wide range of grid resolutions. It is commonly assumed that mesoscale or dynamical fluxes associated with mesoscale and convective circulations simulated by a high-resolution model (subgrid (SG) model) on the subgrid scale of a climate model (large-scale (LS) model) represent additional processes in the ABL, which are not considered by the turbulence scheme of the LS-model, and which can be parameterized using the SG-model. The present study investigates the usefulness of this methodology for small-scale and large-scale idealized heterogeneities using a SG-model resolving mesoscale or even microscale circulations to compute the mesoscale fluxes on the scale of the LS-model. It is shown that the dynamical transports as derived from the SG-model should not be used to correct the parameterized turbulent fluxes of the LS-model. The reason is that the subgrid circulations simulated by the SG-model interact with the fields of wind and scalars in the ABL, which results in reduced turbulent fluxes in the ABL. Thus the methodology of previous studies to use mesoscale/dynamical fluxes for the correction of flux profiles simulated by climate models seems to be questionable.  相似文献   

7.
程锐  宇如聪  徐幼平  王斌 《大气科学》2018,42(6):1286-1296
AREM(Advanced Regional Eta-coordinate Model)对中国暴雨、台风等中尺度天气系统的模拟、预报能力突出。但是伴随模式分辨率的提高,制约该模式发展的一个问题日渐突出,即"静力平衡近似"的约束。本文通过对原静力平衡系统进行修正,引入高阶订正参数定义第三运动方程来构建该模式的非静力动力框架。我们基于Euler原始方程组,有效结合原静力平衡模式的标准层结扣除及IAP(Institute of Atmospheric Physics)变换方法,推导出了球面余纬坐标下的非静力框架,并在E网格和η坐标下进行了时空离散。采用时间两部分离技术进行积分运算以提高计算效率,并通过"追赶法"结合迭代法计算声波。此框架可方便地继承静力平衡框架的特点,最大限度地保留静力平衡框架的优势。理论推导和数值试验表明,当非静力框架退化为静力平衡框架后,方程形式及其模拟结果一致。在文章第二部分将通过理想和实例试验检验非静力模式性能。  相似文献   

8.
This paper deals with a non-hydrostatic mesoscale model that achieves full vectorization on computers like the CYBER 205. The model formulation ensures the conservation of all fluxes and takes into account the terrain inhomogeneities by the aid of suitable transformations. The diagnostic equation for the pressure change is solved using a very efficient vectorized elliptic solver. By imposing appropriate boundary conditions no additional precautions at the boundaries are necessary to achieve meaningful results. As an application, the steady-state inviscid flow over a single mountain is simulated.  相似文献   

9.
多尺度非静力通用模式框架的设计策略   总被引:20,自引:8,他引:12       下载免费PDF全文
针对目前国际上数值预报模式的最新发展趋势,综述了新形势下多尺度非静力通用数值预报模式框架的设计策略。目前世界上开发通用模式主要有两个途经,一是建立一个离散化方案和源程序代码共享的、全球和有限区通用的模式,二是建立一个单一的全球可变分辨率的模式。还从通用模式方程组的选取策略、模式网格属性的构造,时间积分方案、空间离散方案,垂直坐标的选取等方面进行了分析。  相似文献   

10.
JFNK方法概述及其在大气全隐式非静力模式中的应用方案   总被引:3,自引:3,他引:0  
首先介绍了近年来新发展的非线性方程全隐式数值求解的JFNK方法,及其在地球流体力学方面应用计算实例.可看到,无论在计算精度还是计算效率方面,全隐式数值求解远远超过常规的半隐式计算格式.其次,还讨论了JFNK方法在气象非静力模式中应用方案,并提出了用静力假定和半隐式差分格式来构造预条件处理器,变三维求解为二维求解,简化了方程组求解难度.该方案不仅可用于差分模式,也为用譜方法求解非静力模式提供可能.  相似文献   

11.
非静力中尺度高分辨率模式模拟中的垂直坐标影响研究   总被引:1,自引:1,他引:1  
李兴良  陈德辉 《气象学报》2005,63(2):161-172
近几年来,随着高性能计算机技术的高速发展,甚高分辨率的中尺度数值预报模式业务应用已成为可能,与之相关的一系列模式技术的新问题也随之提了出来,垂直坐标系的影响就是其中之一。文中借助于美国新一代数值预报模式WRF(WeatherReseachandForecast),比较了非静力中尺度模式高分辨率模拟应用的垂直坐标影响问题。研究表明,当选用几何高度(z)和气压(p)来构造地形追随坐标时,低层两坐标引起的误差基本一样,中高层高度地形追随坐标引起的误差小于气压地形追随坐标;而且分辨率越高差异愈大。高分辨率模拟结果也表明,这种差异趋势是存在的;此外,本文对天气过程的预报要素场进行了相关的分析。  相似文献   

12.
为了深入理解非静力近似下的波-湍相互作用问题,本研究在σ坐标的海洋环境研究和预报模型(MERF)中引入常用的Mellor-Yamada两方程湍混合参数化方案(MY2.5),评估垂向湍混合对小尺度背风波传播过程的影响.瞬时状态场的模拟结果表明,无论是否为非静力近似条件,上述湍参数化方案的引入都会减弱背风波传播的模拟效果.从时间平均场的试验结果来看,垂向湍混合过程会显著减小非静力近似和静力近似之间的差异.此外,能量收支分析的诊断结果表明,MY2.5方案会显著抑制陆坡地形下的背风波传播过程,进而将更多的潮能转化到不可逆的湍混合过程中.  相似文献   

13.
设计了适用于四维变分同化系统的扰动预报模式GRAPES_PF。根据GRAPES的地形追随坐标非静力原始方程组,采用小扰动分离方法推导微分形式的线性扰动预报方程组,并利用与GRAPES非线性模式相似的数值求解方案求解线性扰动微分方程组。在设计扰动预报模式时采用了两个时间层半隐式半拉格朗日方案对动量方程、热力学方程、水汽方程和连续方程进行时间差分,空间差分方案的变量分布水平方向采用Arakawa C跳点网格,垂直方向采用Charney/Phillips跳层。利用代数消元法进一步推导得到只包含未来时刻扰动Exner气压的亥姆霍兹方程,进而通过广义共轭余差法(GCR)求解,在此基础上得到未来时刻扰动量的预报值。基于所开发的扰动模式开展了数值试验。首先在非线性模式中施加一个中尺度初始扰动高压,得到初始扰动在非线性模式中的后续演变,然后将相同的初始扰动作为扰动模式的初值进行时间积分,将扰动模式预报的结果与非线性模式的结果做了对比。结果表明,所开发的扰动模式GRAPES_PF较好地模拟了惯性重力内波的传播过程:初始高压扰动激发了一个迅速向外传播的惯性重力内波,在气压场向风场适应的过程中,水平风场、垂直运动、位温和湿度等变量均出现了扰动增量,与非线性模式得到的结果相当接近。GRAPES_PF作为GRAPES非线性模式的合理线性模式为建立基于线性扰动预报的区域四维变分同化系统奠定了科学基础。   相似文献   

14.
Abstract

An important objective of middle atmosphere global climate modelling is the development of the capability of predicting the response of the middle atmosphere to natural or anthropogenic perturbations. To achieve this, a comprehensive chemistry package interactively coupled with radiative and dynamical modules is required. This paper presents preliminary results obtained with a photochemistry module which has been incorporated in the Canadian Middle Atmosphere Model (CMAM). The module contains 42 species including necessary oxygen, hydrogen, nitrogen, chlorine, bromine and methane oxidation cycle species. Photochemical balance equations are solved on‐line throughout the middle atmosphere at every dynamical time step. A full diurnal cycle is simulated with photolysis rates provided by a look‐up table. The chemistry solver is a mass conserving, fully implicit, backward difference scheme which currently uses less than 10% of the GCM run time. We present the results obtained from short integrations and compare them with UARS measurements. The model ozone distribution appears in quantitative agreement with observations showing peak values near 10 ppmv and confined to the 35‐km region. The abundance of nitrogen, chlorine, bromine oxides and their respective contributions to the overall ozone budget is realistic. The study illustrates the capability of the model to simulate middle atmosphere photochemistry for the disparate conditions occurring throughout the region.  相似文献   

15.
程锐等(2018)中,我们完成了非静力AREM(Advanced Regional Eta-coordinate Model)模式动力框架设计。本文将通过理想和实例试验检验其模拟能力。设计理想试验并通过与国际成熟的中尺度非静力框架比较,直接检验非静力AREM三维动力框架在细致分辨率(约1 km)下的模拟性能。可以看出,非静力AREM与ARPS(Advanced Regional Prediction System)、WRF(Weather Research and Forecasting Model)模拟出类似的积云对流结构及演变特征,从而基本确证了本文发展的非静力框架的正确性。结合原静力平衡模式的初始化和物理参数化过程,形成非静力AREM模式系统。台风实例模拟表明,粗分辨率下静力、非静力AREM模式性能接近;但在高分辨率下,非静力明显优于静力模式。我们还开展了批量降水试验检验,对非静力AREM模式性能进行了进一步的验证。  相似文献   

16.
Summary Quantitatively comparative experiments of moist convection using hydrostatic and non-hydrostatic models are reviewed and a further study is made of the suitability of the hydrostatic approximation for a high-resolution model when the grid size falls below 20 km. Idealized moist convection is treated, and then the torrential rain that occurred on 6 August 1993 in Kagoshima, southern Kyushu, Japan is simulated by each model. An explicit warm-rain process predicting cloud water and rainwater and the scheme of moist convective adjustment are individually or conjunctively employed in the model. The effect of hydrostatic water loading is also examined in detall.For the simulation of idealized convection, the hydrostatic simulation tends to overestimate and overexpand precipitation in comparison with the non-hydrostatic counterpart, and the drag effect of hydrostatic water loading is more significant for convective development than the non-hydrostatic effect. In the 20-km simulations, however, the hydrostatic simulation with hydrostatic water loading produces results that are comparable to the nonhydrostatic counterpart. For the simulation with real data, the comparative results well correspond to those of idealized convection. Furthermore, the 5 km hydrostatic simulation overestimates total precipitation more than that of dealized convection. On the basis of these results, when developing 1020 km numerical weather prediction (NWP) models, hydrostatic water loading should be evaluated in preference to adopting non-hydrostatic models, and a non-hydrostatic model with hydrostatic water loading is thought to be recommendable for a high-resolution NWP model.With 7 Figures  相似文献   

17.
Summary ?A code for the simulation of atmospheric flows in 3D is presented. The underlying mathematical model is fully compressible, it takes gravity into account but Coriolis forces, turbulence and viscosity are neglected. The general numerical code consists of a finite volume discretization on unstructured hexahedral grids in 3D. The code is presently being investigated on applications to the calculation of atmospheric gravity waves on a mesh which has a structured type and is locally refined near the orography. We develop two schemes, the main difference between them lies in the different discretizations for the mass fluxes. We show that both schemes resolve typical structures of gravity waves in potential flow, linear hydrostatic motion and nonlinear non-hydrostatic regime. We compare advantages and disadvantages of the developed schemes. Received April 20, 2001; revised September 10, 2001  相似文献   

18.
华北回流暴雪发展机理个例研究   总被引:23,自引:5,他引:18  
周雪松  谈哲敏 《气象》2008,34(1):18-26
利用新一代中尺度数值模式WRF对2004年11月24-25日发生在华北的一次回流暴雪过程进行了数值模拟,并着重对回流暴雪形成发展的机制进行了分析研究.在回流暴雪形成与发展的过程中,动力锋生机制发挥了重要作用;其产生的锋面次级环流是造成这次暴雪的主要原因;降水过程中,雨雪相态的转变与温度廓线有直接关系,随着降水由雨转为雪,雪水混合比增大迅速且其高度有所降低.地形在此回流暴雪过程中,通过其强迫的、分别位于迎风坡和背风坡的正负垂直速度中心,对降雪起明显的增幅作用,同时对雨水和云水分布产生影响.  相似文献   

19.
Summary The full-physics adjoint of the FSU Global Spectral Model of version T42L12 is applied to carry out sensitivity analysis of the localized model forecast error to the initial conditions for a case test occurring on June 8, 1988 during the Indian summer monsoon. The results show that adjoint sensitivity based on ECMWF analysis can be used to identify regions with large analysis uncertainties. The conclusion is that more observations are required over the northern Bay of Bengal to improve the quality of analyses so as to ameliorate the model forecast skill.With 17 Figures  相似文献   

20.
In order to simulate the plume produced by large natural draft cooling towers, a semi-spectral warm cloud parameterization has been implemented in an anelastic and non-hydrostatic 3D micro-scale meteorological code. The model results are compared to observations from a detailed field experiment carried out in 1980 at Bugey (location of an electrical nuclear power plant in the Rhône valley in East Central France) including airborne dynamical and microphysical measurements. Although we observe a slight overestimation of the liquid-water content, the results are satisfactory for all the 15 different cases simulated, which include different meteorological conditions ranging from low wind speed and convective conditions in clear sky to high wind and very cloudy. Such parameterization, which includes semi-spectral determination for droplet spectra, seems to be promising to describe plume interaction with atmosphere especially for aerosols and cloud droplets.  相似文献   

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

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