HWRF Nesting. Xuejin Zhang AOML/HRD & UM/CIMAS

Transcription

1 HWRF Nesting Xuejin Zhang AOML/HRD & UM/CIMAS Acknoedgements Coaborators: EMC Hurricane Project (Lead: M. Avicha); DTC (Lead: K. Neman); HRD Modeing Group (Lead: Gopaakrishnan) Computer resource: Jet supercomputer HFIP support

2 HFIP Vision and Goas ( ) VISION Organize the hurricane community to dramaticay improve numerica forecast guidance to NHC in 5-10 years. GOALS Reduce numerica forecast errors in track and intensity by 20% in 5 years, 50% in 10 years Extend forecast guidance to 7 days ith ski comparabe to 5 days at project inception Increase probabiity of predicting rapid intensification (RI) at day 1 to 90% and 60% at day 5 Improve products on storm surge prediction 1

3 Horizonta-Tempora Scaes of Time 1 mo Micro γ scae Atmospheric Processes Micro β scae Micro α scae Meso γ scae Meso β scae Meso α scae Macro β scae Barocinic aves Macro α scae Standing aves Utra ong aves 1 day 1 hour 1 min 1 sec Horizonta Length Pumes Roughness Turbuence Dust devis Thermas Wakes Tornados Deep convection Short gravity aves Fronts (from cross front to aong front) Nocturna o eve jet Squa ines Inertia aves Thunderstorms Coud custers I.G.W. Lake & Mtn Cear air disturbances turbuence Urban effects Tropica cycone (vortex to convective scae) Topography disturbance Tida aves 20 m 200 m 2 km 20 km 200 km 2000 km km Adopted from Oranski (1975) 2

4 Scientific Objectives Preserve across-scae processes on TC genesis, intensifying, decaying, and andfa processes ithin an integrated modeing system to represent the fu-scae spectrum of atmospheric aves and study on muti-scae interactions e.g. TC-terrain interaction, and andfa processes and QPF etc. Enhance resoved resoution that can represent TC inner core physics and can predict TC dynamics ( 3 km or ess) Non-hydrostatic mode becomes required Physics schemes shoud be suitabe to the high-resoution mode 3

5 Mode Deveopment Objectives Taior a too that is operationay feasibe and transferabe ithin the operationa resource capabiity Address the vision and goas of HFIP Meet the NHC operationa requirements Deveop mode capabiities required by the advancements on observing and data assimiation technoogy 4

6 HWRF Grid System Horizonta: Rotated Latitude-Longitude E-Grid Vertica: σ-p hybrid Nesting ØCoincident grid points at start and end points ØFixed nest ratio: 1:3 ØMoving one grid each time 5

7 Rotated Latitude-Longitude Domain Reguar Latitude-Longitude map background Rotated Latitude-Longitude map background Credit: Zavisa Janjic 6

8 Mass points Veocity points E-GRID Structure i = IMS IMS+1 IME-1 i = IME j = JME JME - 1 JME - 2 JMS+2 JMS+1 j = JMS 7

9 Sigma-Pressure Hybrid Coordinate Credit: Zavisa Janjic 8

10 Idea of Moving Nest Moving Nest Setup Domain starting point at mass point Grids aigning at mass points Fine grids adjusting mass Leading edge initiaization Feedback to coarse grids Domain movement Innermost domain fooing the intermediate domain in the current operationa HWRF Mass points (H) Wind points (V) 9

11 Nesting Scheme and Feedback 1 TC one TC to 1 Parae TC one Parae TC to To-ay interaction To-ay interaction The operationa HWRF 10

12 Nest Motion Soution MSLP or vorticity aone is not enough Ne method is neary 100% successfu Ony probems are hen nest movement interva is too sma for storm speed. Rare: decaying extra-tropica systems 11

13 Vortex Tracking Nests Track nine smoothed fieds: Vorticity 10 m, 850 mb, 700 mb Wind minimum 10 m, 850 mb, 700 mb Height 850 mb, 700 mb Membrane MSLP (deveoped by Hui-Ya Chuang at EMC) Discard fieds that are far from the average Fina average is a ne ocation 12

14 Membrane MSLP Project atmosphere on pressure eves Extrapoate virtua temperature on pressure surfaces Smooth atmosphere Integrate to get P(z=0) 13

15 Mass Adjustment Coarse-res Domain Fine-res Domain Hybrid p-σ Coordinate og(p) space Interpoation Hybrid p-σ Coordinate Directy interpoated from parent hybrid to nested hybrid (sigma eves don t match!) Ø The height, temperature and moisture fieds from the parent domain eves are interpoated horizontay onto nested grids on the same hybrid eves in the parent domain. Ø The refined meteoroogica fieds from the parent domain are verticay interpoated onto hybrid eves in the nested domain by using the high-resoution topography pseudo hydrostatic mass baancing. 14

16 HWRF NESTED BOUNDARY CONDITIONS j=jme j=jme-1 j=jme-2 i=1 i=2 i=ime-1 j=1 i=ime Ø Externa Data: Prescribed from the parent after mass adjustment of hydrostatic variabes/direct interpoation of other variabes Ø Bending Zone: Simpe 4 point averaging Ø Dynamic Interface: Mode integration starts from 3 rd ro/coumn Ø Leading edge: fooing same procedure of initia condition j=3 j=2 15

17 Feedback 9 point average Weighted equay beteen the parent and nest average Same vertica mass adjustment procedure 16

18 HWRF nameist.input: nest specific! &domains max_dom = 3, grid_id = 1, 2, 3, parent_id = 0, 1, 2, &time_contro start_year = 2008, 2008, 2008, start_month = 09, 09, 09, start_day = 09, 09, 09, start_hour = 00, 00, 00, start_minute = 00, 00, 00, start_second = 00, 00, 00, end_year = 2008, 2008, 2008, end_month = 09, 09, 09, end_day = 09, 09, 09, end_hour = 06, 06, 06, end_minute = 00, 00, 00, end_second = 00, 00, 00,./frame/LOGICAL FUNCTION nests_to_open CALL n_get_max_dom ( 1, max_dom ) DO nestid = 2, max_dom IF (.NOT. active_domain( nestid ) ) THEN CALL n_get_parent_id ( nestid, parent_id ) IF ( parent_id.eq. parent%id ) THEN CALL n_get_start_year (nestid,s_yr) CALL n_get_end_year (nestid,e_yr) CALL n_get_start_month (nestid,s_mm) CALL n_get_end_month (nestid,e_mm) CALL n_get_start_day (nestid,s_dd) CALL n_get_end_day (nestid,e_dd) CALL n_get_start_hour (nestid,s_h) CALL n_get_end_hour (nestid,e_h) CALL n_get_start_minute (nestid,s_m) CALL n_get_end_minute (nestid,e_m) CALL n_get_start_second (nestid,s_s) CALL n_get_end_second (nestid,e_s). ENDDO 17

19 &domains time_step = 54, time_step_fract_num = 0, time_step_fract_den = 1, max_dom = 3, s_e = 1, 1, 1, e_e = 216, 60, 60, s_sn = 1, 1, 1, e_sn = 432, 100, 100, s_vert = 1, 1, 1 e_vert = 43, 43, 43, parent_grid_ratio = 1, 3, 3, parent_time_step_ratio = 1, 3, 3, dx =.18, 0.06, 0.02, dy =.18, 0.06, 0.02, Domain configuration for the nest./frame/module modue_integrate RECURSIVE SUBROUTINE integrate DO WHILE ( nests_to_open( grid, nestid, kid ) )! nestid is index into mode_config_rec (modue_configure) of the grid! to be opened; kid is index into an open sot in grid's ist of chidren a_nest_as_opened =.true. CALL med_pre_nest_initia ( grid, nestid, config_fags ) CALL aoc_and_configure_domain ( domain_id = nestid, & #ifdef DM_PARALLEL active_this_task = domain_active_this_task( nestid ), & #endif grid = ne_nest, & parent = grid, & kid = kid ) CALL Setup_Timekeeping (ne_nest) CALL med_nest_initia ( grid, ne_nest, config_fags ) IF ( grid%active_this_task ) THEN IF ( grid%dfi_stage == DFI_STARTFWD ) THEN CALL rf_dfi_startfd_init(ne_nest) ENDIF IF (couper_on) CALL cp_defdomain( ne_nest ) ENDIF! active_this_task END DO 18

20 HWRF: ISTART, JSTART AND GRID MOTION &domains grid_id = 1, 2, 3, i_parent_start = 0, ISTART1, ISTART2, j_parent_start = 0, JSTART1, JSTART2, num_moves = -99, scorner_dynamic.f can automaticay determine the southest corner index: istart?, jstart?, ithin the parent domain (coarseresoution) for the chid domain (fine-resoution). Option: num_moves = -99 for automatic grid motion, specificay for Hurricanes; For ocating the initia grid based on storm center, e are providing a utiity../hrfutiities/vortex_init/hrf_set_ijstart/scorner_dynamic.f. This code is compied as part of the hrf-utiities component of HWRF. 19

21 Registry/Registry.NMM Interpoation Routines share/interp_fcn.f Three cases: upscae(up/u), donscae(don/d), boundary forcing(bdy/f) Four methods: nearest neighbor (Near), binary copy (Copy), mass adjustment (Mass), veocity (Ve) state rea u ijkb dyn_nmm 1 v i01rh02u=(upve)d=(donve)f=(bdyve) state rea v ijkb dyn_nmm 1 v i01rh02u=(upve)d=(donve)f=(bdyve) State rea f_ice ikj dyn_nmm 1 rhd=(donmassikj:@eextrap,0.0)u=(upmassikj:@eextrap,0.0) state rea qv ijkfbt moist 1 m rhu=(upmass:@ecopy,0.0), = extrapoation method (beo ground): copy oest mode = extrapoate using constant 5.5 beo = ineary extrapoate to constant at 1030 mb 20

22 Summary The HWRF nesting as deveoped for hurricane forecast appication. The saient features are: Mode configuration is easiy controed through nameist.input Teescopic nests foo the moving system Mutipe nest moving options are avaiabe Mass baance beteen coarse grids to fine grids (donscae and feedback) are directy adjusted on σ eves Mutipe up- and don-scae and forcing agorithms are avaiabe Reference Zhang, X., S.G. Gopaakrishnan, S. Trahan, T.S. Quirino, Q. Liu, Z. Zhang, G. Aaka, and and V. Taapragada. Representing mutipe scaes in the Hurricane Weather Research and Forecasting modeing system: Design of mutipe sets of movabe mutieve nesting and the basin-scae HWRF forecast verification. Weather and Forecasting, 31(6): , doi: /waf-d