Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ...

47
Dynamics & Coupling LACE progress report 2008-2009 Filip V ´ a ˇ na [email protected] CHMI 31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 1

Transcript of Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ...

Page 1: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Dynamics & Coupling

LACE progress report

2008-2009

Filip Vana

[email protected]

CHMI

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 1

Page 2: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

New interpolators for SL

quasi-cubic spline (κ = −2) cubic Lagrange polynomial (κ = 0)

-0.60 -0.48 -0.36 -0.24 -0.12 0.12 0.24 0.36 0.48 0.60

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

10 20 30 40 50 60 70 80 90 100

grid point

10 20 30 40 50 60 70 80 90 100

-0.60 -0.48 -0.36 -0.24 -0.12 0.12 0.24 0.36 0.48 0.60

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

10 20 30 40 50 60 70 80 90 100

grid point

10 20 30 40 50 60 70 80 90 100

quadratic interpolator (κ = 1) strongly diffusive interpolator (κ = 6)

-0.60 -0.48 -0.36 -0.24 -0.12 0.12 0.24 0.36 0.48 0.60

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

10 20 30 40 50 60 70 80 90 100

grid point

10 20 30 40 50 60 70 80 90 100

-0.60 -0.48 -0.36 -0.24 -0.12 0.12 0.24 0.36 0.48 0.60

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

10 20 30 40 50 60 70 80 90 100

grid point

10 20 30 40 50 60 70 80 90 100

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 2

Page 3: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

New interpolators for SL

(κ = 0, ε = 0.00) (κ = 0, ε = 0.01) (κ = 0, ε = 0.02)

-0.60 -0.48 -0.36 -0.24 -0.12 0.12 0.24 0.36 0.48 0.60

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

10 20 30 40 50 60 70 80 90 100

grid point

10 20 30 40 50 60 70 80 90 100

-0.60 -0.48 -0.36 -0.24 -0.12 0.12 0.24 0.36 0.48 0.60

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

10 20 30 40 50 60 70 80 90 100

grid point

10 20 30 40 50 60 70 80 90 100

-0.60 -0.48 -0.36 -0.24 -0.12 0.12 0.24 0.36 0.48 0.60

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

10 20 30 40 50 60 70 80 90 100

grid point

10 20 30 40 50 60 70 80 90 100

(κ = 0, ε = 0.04) linear interpolator

-0.60 -0.48 -0.36 -0.24 -0.12 0.12 0.24 0.36 0.48 0.60

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

10 20 30 40 50 60 70 80 90 100

grid point

10 20 30 40 50 60 70 80 90 100

-0.60 -0.48 -0.36 -0.24 -0.12 0.12 0.24 0.36 0.48 0.60

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

0.0

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

heig

ht [k

m]

10 20 30 40 50 60 70 80 90 100

grid point

10 20 30 40 50 60 70 80 90 100

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 3

Page 4: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

New interpolators for SL

Experience from the real atmosphere(LACE domain, ∆x=9km)

Decreasing of κ leads to:

general improvement of wind speed and MSLpressure (conservation of mass)

other model variables improved in uppertroposphere (above 500 hPa) and in atmosphereabove 100 hPa

improvement for prognostic physical variables

detrimental effects in PBL and lower stratosphere

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 4

Page 5: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

New interpolators for SL

Experience from the real atmosphere(LACE domain, ∆x=9km)

Decreasing of κ leads to:

general improvement of wind speed and MSLpressure (conservation of mass)

other model variables improved in uppertroposphere (above 500 hPa) and in atmosphereabove 100 hPa

improvement for prognostic physical variables

detrimental effects in PBL and lower stratosphere

Decreased κ needs to be compensated by increasedhorizontal and vertical diffusion

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 4

Page 6: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

New interpolators for SL

ALADIN/CE, parallel suite (06/05/2009-08/06/2009)

Relative humidity rmse Relative humidity bias

MSLP bias wind speed bias cloudiness bias

κ = 0., εH=0.0025, εV =0. ⇒ κ = -0.6, εH,V =0.016

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 5

Page 7: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

New interpolators for SL

Conclusions

Horizontal smoothing influences KE spectra, onlyweak effect to scores

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 6

Page 8: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

New interpolators for SL

Conclusions

Horizontal smoothing influences KE spectra, onlyweak effect to scores

Vertical smoothing nearly no effect to KE spectra,strong impact to scores and mass conservation

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 6

Page 9: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

New interpolators for SL

Conclusions

Horizontal smoothing influences KE spectra, onlyweak effect to scores

Vertical smoothing nearly no effect to KE spectra,strong impact to scores and mass conservation

Lagrangian cubic interpolation seems to performextremely well in the real atmosphere

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 6

Page 10: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

New interpolators for SL

Conclusions

Horizontal smoothing influences KE spectra, onlyweak effect to scores

Vertical smoothing nearly no effect to KE spectra,strong impact to scores and mass conservation

Lagrangian cubic interpolation seems to performextremely well in the real atmosphere

Still some space to re-distribute the SL diffusivityfrom the basic interpolator toward a stronger (3Dacting) diffusion scheme

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 6

Page 11: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

LACE project targeted to ∆x ≈ 4km

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 7

Page 12: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

LACE project targeted to ∆x ≈ 4km

Hydrostatic dynamics of ALADIN is simplyextensible to the NH one (see Bénard et al., 2009,QJRMS)

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 7

Page 13: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

LACE project targeted to ∆x ≈ 4km

Hydrostatic dynamics of ALADIN is simplyextensible to the NH one (see Bénard et al., 2009,QJRMS)

Benchmark the NH dynamics with respect to thehydrostatic one: eliminating the true added value ofthe NH dynamics for the real atmosphericsimulations.

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 7

Page 14: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

LACE project targeted to ∆x ≈ 4km

Hydrostatic dynamics of ALADIN is simplyextensible to the NH one (see Bénard et al., 2009,QJRMS)

Benchmark the NH dynamics with respect to thehydrostatic one: eliminating the true added value ofthe NH dynamics for the real atmosphericsimulations.

The aim is to define the resolution where it is worthto activate the NH dynamics.

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 7

Page 15: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 16: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 17: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 18: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 19: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 20: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 21: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 22: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 23: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 24: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 25: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 26: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 27: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Cold bubble test

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 8

Page 28: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Parallel suite ALADIN/CE, ∆x=9 km, L43

Identify eventual problems of NH dyn.

Learn about the systematic extra value of the NH (for that scales)

DA and LBC field are the same (hydrostatic)

Comparable model settings (geometry, physics, timestep

organization,. . . )

Dynamics as close as possible to the hydrostatic reference

Horizontal diffusion, LBC coupling (two extra variables)

SI background (T ∗, π∗)

vertical discretization

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 9

Page 29: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Impact of the SI background

0 200 400 600 800 1000timestep

1e-20

1e-19

1e-18

1e-17

1e-16

1e-15

spec

tral

norm

of

kin

etic

ener

gy [J

/kg]

NSITER=1, SIPR=100000., SITR=300.

NSITER=1, SIPR=90000., SITR=350.

NSITER=1, SIPR=90000., SITR=300.

NSITER=1, SIPR=100000., SITR=350.

NSITER=0, SIPR=100000., SITR=300. XIDT=0.05

L=10dx, H=100m, linear orographyeta coordinate, hydrostatic run

2D model, time evolution of KE

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 10

Page 30: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Vertical discretization

HYD NH

NDLNPR=0 yes -

NDLNPR=1 yes yes

VFE yes not yet

NDLNPR = 0: δl = δπl

πl

= δ ln πl = ln(

πl

πl−1

)

NDLNPR = 1: δl = δπl

πl

l−π

l−1√π

l−1

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 11

Page 31: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Vertical discretization (II.)rmse of geopotential height

default NH versus default hydr. hydr. NDLNPR=1 vs. NDLNPR=0

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 12

Page 32: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Vertical discretization (II.)rmse of geopotential height

NH versus hydr. (NDLNPR=1) VFE vs. NDLNPR=1 (both hydr.)

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 12

Page 33: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Vertical discretization (II.)rmse of geopotential height

default NH versus def. hydr. hydr. VFE vs. def. hydr.

⇒ Higher sophistication of vertical discretization has potential to improve result,

especially when provided by sufficient vertical resolution.

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 12

Page 34: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

ICI versus SI schemermse of geopotential height

NSITER = 1 vs. NSITER = 0 NSITER = 3 vs. NSITER = 1

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 13

Page 35: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

ICI versus SI scheme (II.)rmse of geopotential height

VFE (hydr.) vs. HN (NSITER = 0) VFE (hydr.) vs. NH (NSITER = 1)

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 14

Page 36: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

ICI versus SI scheme (III.)

Surface parameters (NSITER=1 vs. NSITER=0)

bias rmse

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 15

Page 37: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Summary (for LAM, ∆x=9 km, L43)

NH vs. hydrostatic

more variability of wind in PBL (better bias,worse stde)

otherwise no significant signal at all

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 16

Page 38: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Summary (for LAM, ∆x=9 km, L43)

NH vs. hydrostatic

more variability of wind in PBL (better bias,worse stde)

otherwise no significant signal at all

ICI vs. SI

ICI generally beneficial

systematic shift of surface variables bias (drift bysurface scheme?)

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 16

Page 39: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

Summary (for LAM, ∆x=9 km, L43)

NH vs. hydrostatic

more variability of wind in PBL (better bias,worse stde)

otherwise no significant signal at all

ICI vs. SI

ICI generally beneficial

systematic shift of surface variables bias (drift bysurface scheme?)

GP-w option (LGVADW)

no significant signal at all

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 16

Page 40: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

ConclusionsNH dynamics seems to perform safely

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 17

Page 41: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

ConclusionsNH dynamics seems to perform safely

Quality of results are related to their cost:

configuration CPU increase (NEC-SX6, OpenMP)

hydr. NDLNPR=0 1.

hydr. NDLNPR=1 1.

hydr. VFE 1.15

NH SI 1.08

NH ICI (NSITER=1) 1.52

NH ICI, GP-w (NSITER=1) 1.49

NH ICI (NSITER=3) 2.32

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 17

Page 42: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

ConclusionsNH dynamics seems to perform safely

Quality of results are related to their cost:

configuration CPU increase (NEC-SX6, OpenMP)

hydr. NDLNPR=0 1.

hydr. NDLNPR=1 1.

hydr. VFE 1.15

NH SI 1.08

NH ICI (NSITER=1) 1.52

NH ICI, GP-w (NSITER=1) 1.49

NH ICI (NSITER=3) 2.32

Would be nice to have NH (ICI) with VFE. . .

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 17

Page 43: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Toward the NH dynamics

ConclusionsNH dynamics seems to perform safely

Quality of results are related to their cost:

configuration CPU increase (NEC-SX6, OpenMP)

hydr. NDLNPR=0 1.

hydr. NDLNPR=1 1.

hydr. VFE 1.15

NH SI 1.08

NH ICI (NSITER=1) 1.52

NH ICI, GP-w (NSITER=1) 1.49

NH ICI (NSITER=3) 2.32

Would be nice to have NH (ICI) with VFE. . .

To be repeated with ∆x=4.5km, L87

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 17

Page 44: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Pressure gradient term

PG scheme of Simmons & Jiabin (1990)The present form:

∇Φ + RT∇(ln p)

is replaced by:

∇Φ + RdT∇(ln p)

with T = Tv − T0 (p/p0)α

Φ = Φs + RdT0

α(p/p0)

α +∫ ps

pRdT

pdp

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 18

Page 45: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Pressure gradient term

First preliminary results (∆x=9km)2008/04/16 00UTCBase2008/04/16 18UTCValid

mma103@voodoo Fri Mar 20 14:38:23 2009 [PFALAD+18.newgp2 PFALAD+18.refgp]

PG scheme - ref scheme

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 19

Page 46: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Pressure gradient computation

Tests in 2D (HYD)

0 200 400 600 800 1000timestep

1e-21

1e-20

1e-19

1e-18

1e-17sp

ectr

al n

orm

of

kin

etic

en

erg

y

[J/k

g]

refPG

L=10dx, H=100m, linear orographysigma coordinate, hydrostatic run

Still too easy for the pressure gradient term. . .

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 20

Page 47: Dynamics&Coupling · Dynamics&Coupling LACEprogressreport 2008-2009 Filip Va´naˇ filip.vana@chmi.cz CHMI 31stEWGLAM+16thSRNWPMeetings,28thSeptember–1stOctober2009,Athens– p.

Other development

Rotated Mercator projection (including TL/AD)

More OpenMP for LAM (enabled MPI-OpenMP,intelligent scheduling,. . . )

Optimization of SL interpolations

Separation of SL quantities to be interpolated (DDH& 3D turb)

Fully elastic projection of heat

Boyd’s technique for bi-periodicity of LAM

31st EWGLAM + 16th SRNWP Meetings, 28th September – 1st October 2009, Athens – p. 21