a 3-d finite-volume nonhydrostatic icosahedral model (nim) jin lee
Post on 26-Jan-2016
16 Views
Preview:
DESCRIPTION
TRANSCRIPT
A 3-D Finite-VolumeNonhydrostatic Icosahedral Model (NIM)
Jin Lee
Earth System Research Laboratory (ESRL)
Aeronomy Lab.Climate Diagnostic center
Climate Monitoring and Diagnostic LabEnvironmental Technology Lab
Forecast Systems Lab
Chemical Sciences DivGlobal Monitoring DivPhysical Sciences DivGlobal Systems Div
GFDL,NSSL,ARL,AOML,GLERL,PMEL
Earth System Research Laboratory(ESRL)
Aeronomy Lab.Climate Diagnostic center
Climate Monitoring and Diagnostic LabEnvironmental Technology Lab
Forecast Systems Lab
Chemical Sciences DivGlobal Monitoring DivPhysical Sciences DivGlobal 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 4Global Modeling Program Review
NOAA Earth System Research Laboratory - Boulder, Colorado Page 5
• 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.
NOAA Earth System Research Laboratory - Boulder, Colorado Page 6
• 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.Lee, Bleck, and MacDonald (2010, JCP): A Multistep Flux-Corrected Transport Scheme.
.
Novel features of FIM/NIM:
NOAA Earth System Research Laboratory - Boulder, Colorado Page 7
• 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.
.
Novel features of FIM/NIM:
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 NHImprovement over GFS from FIM using same initial conditions.
120h NH 500 hPa fcsts -- Nov 2009 – Sept 2010
GFS better
FIM better
NOAA Earth System Research Laboratory - Boulder, Colorado Page 12
• 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 MacDonald, Middlecoff, Henderson, and Lee (2010, IJHPC) : A General Method for
Modeling on Irregular Grids.Modeling on Irregular Grids..
Novel features of FIM/NIM:
NOAA Earth System Research Laboratory - Boulder, Colorado Page 13
• 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.
Novel features of FIM/NIM:
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 Geostrophic adjustment of finite-volume A-grid model in inviscid and viscid flowsA-grid model in inviscid and viscid flows
Geostrophic adjustment ExperimentGeostrophic adjustment Experiment : Finite-volume A-Grid SWM: Finite-volume A-Grid SWM
Inviscide model simulation on hexagonal A-grid (dx~60 km)
19
Geostrophic adjustment ExperimentGeostrophic adjustment Experiment : Finite-volume A-Grid SWM: Finite-volume A-Grid SWM
Inviscide model simulation on hexagonal A-grid (dx~60 km)
20
Geostrophic adjustment ExperimentGeostrophic adjustment Experiment : Finite-volume A-Grid SWM: Finite-volume A-Grid SWM
Inviscide model simulation on hexagonal A-grid (dx~60 km)
21
Geostrophic adjustment ExperimentGeostrophic adjustment Experiment : Finite-volume : Finite-volume AA-Grid Model-Grid Model
Inviscide model simulation on hexagonal A-grid (dx~60 km)( repeat same experiment)
22
typical dissipation operator (small coefficient)
2
Geostrophic adjustment ExperimentGeostrophic adjustment Experiment : Finite-volume : Finite-volume AA-Grid Model-Grid Model
+ typical dissipation operator (small coefficient)
2
Inviscide model simulation on hexagonal A-grid (dx~60 km)
23
Geostrophic adjustment ExperimentGeostrophic adjustment Experiment : Finite-volume : Finite-volume AA-Grid Model-Grid Model
+ typical dissipation operator (small coefficient)
2
Inviscide model simulation on hexagonal A-grid (dx~60 km)
24
25
The 2nd-order C-grid scheme over-amplifies mountain waves. Reinecke and Durran (MWR,2009)
NOAA Earth System Research Laboratory - Boulder, Colorado Page 26
• 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.
Novel features of FIM/NIM:
ikl ,
ikn ,
3-D control volume box
NOAA Earth System Research Laboratory - Boulder, Colorado Page 27
• 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.
Novel features of FIM/NIM:
NOAA Earth System Research Laboratory - Boulder, Colorado Page 28
3-D control volume box
0
;
00
);,,,(),,,(
),,,(),,,(;,,,,,,,,
.0
0
0
0
p
p
p
Rpp
Rptzyxztzyx
tzyxztzyxwvuWVU
q
p
Sz
Wq
y
Vq
x
Uq
t
q
TC
H
z
W
y
V
x
U
t
z
W
y
V
x
U
t
ggz
Rz
Ww
y
Vw
x
Uw
t
W
yR
z
Wv
y
Vv
x
Uv
t
V
xR
z
Wu
y
Vu
x
Uu
t
U
Flux form GEs on 3-D control volume on height coord.
NOAA Earth System Research Laboratory - Boulder, Colorado Page 29
• 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.
Novel features of FIM/NIM:
NOAA Earth System Research Laboratory - Boulder, Colorado Page 30
• 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.
Novel features of FIM/NIM:
NOAA Earth System Research Laboratory - Boulder, Colorado Page 31
• 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
Novel features of FIM/NIM:
NOAA Earth System Research Laboratory - Boulder, Colorado Page 32
• 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.
Novel features of FIM/NIM:
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
10
ground 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
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
* In Namelist, Physics= ‘gfs’ or ‘grims’
Physics packagesPhysics packages
NIM 800-day aqua-planet simulationMODEL 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 ResultsPreliminary Results
Zonal Mean SSTZonal Mean SST
Hoskins et al. (1999), Tellus NIM
43
Mean Zonal Wind (3-month)Mean Zonal Wind (3-month)
perturbation pressure
meridional wind
vertical velocity
MSLP + precip
MSLP + precip
MSLP + precip
MSLP + precip
MSLP + precip
MSLP + precip
MSLP + precip
MSLP + precip
MSLP + precip
MSLP + precip
NIM aqua-planet simulationglobal precipitation map NIM precipitation
56
Microphysics sensitivities on model resolutions 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.
top related