a 3-d finite-volume nonhydrostatic icosahedral model … earth system research laboratory - boulder,...

A 3-D Finite-Volume Nonhydrostatic Icosahedral Model (NIM)

Jin Lee

Earth System Research Laboratory (ESRL)

Aeronomy Lab. Climate Diagnostic center

Climate Monitoring and Diagnostic Lab Environmental Technology Lab

Forecast Systems Lab

Chemical Sciences Div Global Monitoring Div Physical Sciences Div Global Systems Div

GFDL,NSSL,ARL,AOML,GLERL,PMEL

Earth System Research Laboratory(ESRL)

Aeronomy Lab. Climate Diagnostic center

Climate Monitoring and Diagnostic Lab Environmental Technology Lab

Forecast Systems Lab

Chemical Sciences Div Global Monitoring Div Physical Sciences Div Global Systems Div

Modeling goal: to develop a non-hydrostatic icosahedral global model for weather and climate applications

GFDL,NSSL,ARL,AOML,GLERL,PMEL

Director, Dr. A.E. (Sandy) MacDonald

ESRL finite-volume Icos- models (FIM/NIM)

FIM (flow-following finite- volume Icosahedral model):

Target resolution ≥ 10 km A hydrostatic model consists of 2-D

finite-volume SWM coupled with hybrid σ-θ vertical solver.

A hydrostatic model for NCEP global model ensemble.

NIM (Nonhydrostatic Icosahedral model):

Target resolution : O(1 km) and beyond Extension of 2-D finite-volume integration

into 3-D integration on control volume defined on the height coordinate.

Use the latest GPU technology to speed up high-resolution model calculations.

ESMF

19 Apr 2011 4 Global Modeling Program Review

Presenter

Presentation Notes

Preliminary tests of FIM done on NCEP CCS

NOAA Earth System Research Laboratory - Boulder, Colorado

• Finite-volume Integrations on Local Coordinate

Novel features of FIM/NIM:

x

y

z

Ο

S

P

••

•

•

•

• ikl ,

ikn ,

2-D f.-v. operator carried out on straight lines, rather than along the 3-D curved lines on the sphere

Lee and MacDonald (MWR, 2009): A Finite-Volume Icosahedral Shallow Water Model on Local Coordinate.


• Finite-volume Integrations on Local Coordinate • Conservative and Monotonic Adams-Bashforth 3rd-order FCT Scheme - Lee, Bleck, and MacDonald (2010, JCP): A Multistep Flux-Corrected Transport Scheme.

.



• Finite-volume Integrations on Local Coordinate • Conservative and Monotonic Adams-Bashforth 3rd-order FCT Scheme • FIM: Hybrid σ-θ Coordinate w/ GFS Physics - Bleck, Benjamin, Lee and MacDonald (2010, MWR): On the Use of an Arbitrary

Lagrangian-Eulerian Vertical Coordinate in Global Atmospheric Modeling. .


FIM forecasts initialized with GSI initial condition at 00 and 12 UTC

(comparisons of FIM and GFS 500 mb ACC)

• Interpolate GFS initial data to FIM model grid.

• Perform hydrostatic initialization.

• Perform 7-day fcst with dx~30km.

• Same initial condition, terrain & sfc parameters, physics package runing at similar model resolutions (dx~30 km)

Main dropouts from GFS in NH Improvement over GFS from FIM using same initial conditions.

120h NH 500 hPa fcsts -- Nov 2009 – Sept 2010

GFS better

FIM better


• Finite-volume Integrations on Local Coordinate • Conservative and Monotonic Adams-Bashforth 3rd-order FCT Scheme • FIM: Hybrid σ-θ Coordinate w/ GFS Physics • Efficient Indirect Addressing Scheme on Irregular Grid - MacDonald, Middlecoff, Henderson, and Lee (2010, IJHPC) : A General Method for

Modeling on Irregular Grids. .



• Finite-volume Integrations on Local Coordinate • Conservative and Monotonic Adams-Bashforth 3rd-order FCT Scheme • FIM: Hybrid σ-θ Coordinate w/ GFS Physics • Efficient Indirect Addressing Scheme on Irregular Grid • Grid Optimization for Efficiency and Accuracy - Wang and Lee (2011, SIAM): Geometric Properties of Icosahedral-Hexagonal Grid on Sphere.


SBiR

SPDN

MGCL

Comparisons of Icosa-grids

Uniformity Regularity

SBiR 1.195 1.476

MBiR 1.175 1.405

SGCL 1.476 1.194

MGCL 1.446 1.135

STD

SPN MGC

Williamson etal.(1992) Case V: Zonal flow over Mountain

SBiR (G8/dt=45 sec)

MGCL (G8/dt=36 sec, blow up with dt=45)

Geostrophic adjustment of finite-volume A-grid model in inviscid and viscid flows

Geostrophic adjustment Experiment : Finite-volume A-Grid SWM

Inviscide model simulation on hexagonal A-grid (dx~60 km)

19



20

Presenter

Presentation Notes

12 hrs?



21

Geostrophic adjustment Experiment : Finite-volume A-Grid Model

Inviscide model simulation on hexagonal A-grid (dx~60 km) ( repeat same experiment)

22

typical dissipation operator (small coefficient)

2∇


+ typical dissipation operator (small coefficient)

2∇


23


+ typical dissipation operator (small coefficient)

2∇


24

25

The 2nd-order C-grid scheme over-amplifies mountain waves. Reinecke and Durran (MWR,2009)


• Finite-volume Integrations on Local Coordinate • Conservative and Monotonic Adams-Bashforth 3rd-order FCT Scheme • FIM: Hybrid σ-θ Coordinate w/ GFS Physics • Efficient Indirect Addressing Scheme on Irregular Grid • Grid Optimization for Efficiency and Accuracy • Novel Features of NIM: -Three-dimensional finite-volume integration.


ikl ,

ikn ,

3-D control volume box


• Finite-volume Integrations on Local Coordinate • Conservative and Monotonic Adams-Bashforth 3rd-order FCT Scheme • FIM: Hybrid σ-θ Coordinate w/ GFS Physics • Efficient Indirect Addressing Scheme on Irregular Grid • Grid Optimization for Efficiency and Accuracy • Novel Features of NIM: -Three-dimensional finite-volume integration. - Conservative flux formulation on height coordinate.



3-D control volume box ( ) ( ) ( )

( ) ( ) ( )

( ) ( ) ( )

( ) ( ) ( )

( ) ( ) ( )

( ) ( ) ( ) ( )

( ) ( ) ( )( )

κ

π

γ

πγρρρ

ρρθρρρρ

ρ

θθθ

ρ

ρππρπγ

πγ

πγ

=∂

∂+

∂∂

+∂

∂+

∂∂

Θ=

∂∂

+∂

∂+

∂∂

+∂Θ∂

=∂

∂+

∂∂

+∂∂

+∂∂

=

′+

′−

∂Θ′∂

+∂

∂+

∂∂

+∂

∂+

∂∂

=∂Θ′∂

+∂

∂+

∂∂

+∂

∂+

∂∂

=∂Θ′∂

+∂

∂+

∂∂

+∂

∂+

∂∂

=Θ

=

Θ∇=∇′+=

Θ′+Θ=Θ=Θ

0;

00

);,,,(),,,(

),,,(),,,(;,,,,,,,,

.0

0

0

0

p

p

p

Rpp

Rptzyxztzyx

tzyxztzyxwvuWVU

q

p

Sz

Wqy

Vqx

Uqtq

TCH

zW

yV

xU

t

zW

yV

xU

t

ggz

Rz

Wwy

Vwx

Uwt

Wy

Rz

Wvy

Vvx

UvtV

xR

zWu

yVu

xUu

tU

Flux form GEs on 3-D control volume on height coord.


• Finite-volume Integrations on Local Coordinate • Conservative and Monotonic Adams-Bashforth 3rd-order FCT Scheme • FIM: Hybrid σ-θ Coordinate w/ GFS Physics • Efficient Indirect Addressing Scheme on Irregular Grid • Grid Optimization for Efficiency and Accuracy • Novel Features of NIM: -Three-dimensional finite-volume integration. - Conservative flux formulation on z-coordiante. - 3-D volume Integration to calculate PGF.



• Finite-volume Integrations on Local Coordinate • Conservative and Monotonic Adams-Bashforth 3rd-order FCT Scheme • FIM: Hybrid σ-θ Coordinate w/ GFS Physics • Efficient Indirect Addressing Scheme on Irregular Grid • Grid Optimization for Efficiency and Accuracy • Novel Features of NIM: -Three-dimensional finite-volume integration. - Conservative flux formulation on z-coordiante. - 3-D volume Integration to calculate PGF. - Horizontal explicit, semi-implicit tri-diagonal solver for vertically propagating acoustic waves.



• Finite-volume Integrations on Local Coordinate • Conservative and Monotonic Adams-Bashforth 3rd-order FCT Scheme • FIM: Hybrid σ-θ Coordinate w/ GFS Physics • Efficient Indirect Addressing Scheme on Irregular Grid • Grid Optimization for Efficiency and Accuracy • Novel Features of NIM: -Three-dimensional finite-volume integration. - Conservative flux formulation on z-coordiante. - 3-D volume Integration to calculate PGF. - Horizontal explicit, semi-implicit tri-diag solver for vertically propagating acoustic waves. - Runge-Kutta (RK4) for time discretization



• Finite-volume Integrations on Local Coordinate • Conservative and Monotonic Adams-Bashforth 3rd-order FCT Scheme • FIM: Hybrid σ-θ Coordinate w/ GFS Physics • Efficient Indirect Addressing Scheme on Irregular Grid • Grid Optimization for Efficiency and Accuracy • Novel Features of NIM: -Three-dimensional finite-volume integration. - Conservative flux formulation on z-coordiante. - 3-D volume Integration to calculate PGF. - Horizontal explicit, semi-implicit tri-diag solver for vertically propagating acoustic waves. - Runge-Kutta (RK4) for time discretization - Fast GPUs to speed up calculation.


NIM 2-D test cases: (in 3-D model grid)

heat forced circulation,

warm bubble, density current (cold bubble),

mountain waves, Internal gravity waves (icos-grid)

multi-months aqua-planet simulations (icos-grid)

•

• •

• •

• •

•

• •

•

•

•

X

Y

.0)( =yf

35

mzx

mHCo

10ground theabove m 260at centered

,125/2'

=∆=∆

== αθ

Comparisons of vertical solvers with the warm bubble simulation. A rising thermal in an isentropic atmosphere.

t= 0.0 min t= 0.0 min θ ′ θ ′

Explicit .vs. Implicit tri-diag solvers

t= 14.0 min t= 14.0 min θ ′ θ ′

Explicit .vs. Implicit tri-diag solvers

38

NIM density current simulation. Negative buoyancy to initiate density current.

39

2-D Mountain Waves (DCMIP case 2-1 with h=100 m, dx=500m, a=7km)

(Schar etal, 2002, MWR, Klemp 2011, MWR)

Internal Gravity Waves (DCMIP case 3-1, 3-D simulation on icosahedral grid)

41

NIM

Same as that of FIM To be used for medium-

range weather forecast

Global/Regional Integrated Modeling system (YSU, Korea)

Reference : Hong et al. (APJAS 2012, in review)

Climate research high-resolution run

GFS physics GRIMs physics

* In Namelist, Physics= ‘gfs’ or ‘grims’

Physics packages

NIM 800-day aqua-planet simulation MODEL NIM/GRIMs

SST Zonally uniform, max. temp. on equator Resolution G5 (∆x ~ 240 km)

Vertical 32 Stretch layers Model top 25 km

∆t 20 min

42

Preliminary Results

Zonal Mean SST

Presenter

Presentation Notes

To verify the interface Physics verification

Hoskins et al. (1999), Tellus NIM

43

Mean Zonal Wind (3-month)

perturbation pressure

meridional wind

vertical velocity

MSLP + precip

NIM aqua-planet simulation global precipitation map NIM precipitation

56

Microphysics sensitivities on model resolutions

• Microphysics : Reduction of precipitation in WSM3 • Resolution : G7 has more precipitation • Pick : Large on the equator : Small on mid latitudes

G7

G6

G5

Final remarks and future outlook •A 3-D f.-v. Nonhydrostatic Icosahedral Model (NIM) dycore has been developed and tested w/ mesoscale benchmarks

•Incorporated GFS and GRIMs physics packages into NIM dycore.

•Implemented NIM on GPU and CPU clusters. •Aqua-planet simulations to test dynamics/physics interfaces. •Test NIM for medium-range NWP at 30-km resolution initialized w/ GSI. •Future NIM applications for hurricane track/intensity forecasts and intra-seasonal predictions.

a 3-d finite-volume nonhydrostatic icosahedral model … earth system research laboratory - boulder,...

Documents