210 likes | 300 Views
Future Modeling?. Hann-Ming Henry Juang RSM workshop 2011. Topics. Different IOs? To downscaling from different data? More accurate? Generalized coordinates as GFS Speed up model code Semi-Lagrangian Modulerization Fortran 99 Non dimensional code Others?.
E N D
Future Modeling? Hann-Ming Henry Juang RSM workshop 2011
Topics • Different IOs? • To downscaling from different data? • More accurate? • Generalized coordinates as GFS • Speed up model code • Semi-Lagrangian • Modulerization • Fortran 99 • Non dimensional code • Others?
A specific form for generalized hybrid where For sigma-pressure, let for all levels - decrease from 1 to 0 at low layers - increase from 0 to be pressure where B=0, then decrease to zero. for all levels For sigma-theta, let works the same, and works as in sigma-p
sigma-theta parallel M operational parallel H 2005 hurricane season
Numerical instability due to advection by the CFL condition Let i.e u=200m/s dx=1000m dt <5 sec But if we solve by then unconditional stable
A M Traditional- Upstream SL D
Starting from mid-point A M NDSL D No guessing and no iteration but one 2-D interpolation and one 2-D remapping
Instead doing following Interpolation relocation D M A remapping
For mass conservation, let’s start from continuity equation Consider 1-D and rewrite it in advection form, we have Advection form is for semi-Lagrangian, but it is not conserved if divergence is treated as force at mid-point, So divergence term should be treated with advection
Divergence term in Lagrangian sense is the change of the volume if mass is conserved, so we can write divergence form as Put it into the previous continuity equation, we have which can be seen as
We do Interpolation relocation ML MR X AR DL DR M AL
control 06h fcst specific humidity at model layer 40 nislfv
control 24h fcst specific humidity at model layer 40 nislfv
control 72h fcst specific humidity at model layer 40 nislfv
control 6hr fcst cloud water at model layer 35 nislfv
control 24hr fcst cloud water at model layer 30 nislfv
T574 Opr Eulerian TL878 NDSL T574 GC Eulerian
Day 1 Day 5 TL878 NDSL Opr T574 Eulerian
Day 1 Day 5 Opr T574 Eulerian NDSL TL878
The possible future • More generalized and/or hybrid system will be built, more accurate thermodynamics will be used. • Semi-Lagrangian will become more simple, conserving and economical code for time integration. • Gravity and acoustic waves can be resolved explicitly (show in Riemann solver talks). • Deep atmosphere or non-approximated system will be introduced for further accuracy in model dynamics.