1 / 28

Conservative Ca scade R emapping between S pherical grids ( CaRS )

NCAR, SIParCS Student Internship Presentations, Friday 10 th August, 2007. Conservative Ca scade R emapping between S pherical grids ( CaRS ). Arunasalam Rahunanthan Department of Mathematics, University of Wyoming. Introduction. Need for CaRS For a better higher order method.

laurel
Download Presentation

Conservative Ca scade R emapping between S pherical grids ( CaRS )

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. NCAR, SIParCS Student Internship Presentations, Friday 10th August, 2007. Conservative Cascade Remapping between Spherical grids (CaRS) Arunasalam Rahunanthan Department of Mathematics, University of Wyoming.

  2. Introduction • Need for CaRS • For a better higher order method. • For the coupling between different model components data. • Goal • Higher order accurate interpolation of field variables from one spherical grid to another without violating conservation and monotonicity.

  3. Implementation • SCRIP – a conservative remapping procedure on the sphere (Jones, Mon.Wea.Rev, 1999) • Advantage - Great geometric flexibility and capable of handling different types of spherical grids. • Disadvantage – Low order method. • CaRS – a cascade remapping between cubed-sphere grids and the RLL grids • CaRS – area based (Lauritzen & Nair, Mon.Wea.Rev, July 2007) • CaRS performance compared with SCRIP • CaRS – length based • CaRS length based approach compared with CaRS area based approach

  4. Cascade interpolation • Algorithm based on cascade interpolation method developed for semi-Lagrangian advection schemes • A two dimensional interpolation problem split into two one dimensional problems • First sweep - Interpolation from ‘o’ to ‘x’ along the source grid lines • Second sweep – The resulting field interpolated from ‘x’ to • Allows for high-order sub grid cell construction and advanced monotone filters

  5. Cubed-sphere latitudes and longitudes • The entire cubed-sphere grid reconstructed with a family of horizontal and vertical grid lines. • The cubed-sphere latitudes as vertically stacked closed curves (squared patterns) • The cubed-sphere latitudes belts intersected by a set of cubed-sphere longitudes

  6. Gnomonic projection and cubed-sphere grid Gnomonic projection North pole South pole • Great circle arcs are straight lines on gnomonic projection

  7. First sweep Source grid

  8. First sweep Source grid

  9. First sweep Source grid

  10. First sweep Source grid

  11. First sweep Target grid (intermediate grid) Source grid

  12. First sweep Target grid (intermediate grid) Source grid

  13. First sweep Target grid (intermediate grid) Source grid

  14. First sweep Target grid (intermediate grid) Source grid North pole North pole South pole South pole

  15. Second sweep Source grid (intermediate grid)

  16. Second sweep Source grid (intermediate grid)

  17. Second sweep Source grid (intermediate grid)

  18. Second sweep Source grid (intermediate grid) Target grid (cubed-sphere grid)

  19. Second sweep Source grid (intermediate grid) Target grid (cubed-sphere grid)

  20. Second sweep Source grid (intermediate grid) Target grid (cubed-sphere grid)

  21. Fields used for testing CaRS

  22. Results – Mapping from a coarse RLL grid to a fine cubed-sphere grid (a) Field :: Y3216 Reconstruction :: PCM Nc = 130,Nλ= 128, Nµ=64 Length based CaRS Area based CaRS

  23. Results – Mapping from a fine cubed-sphere grid to a coarse RLL grid (b) Field :: Vortex, Reconstruction :: PPM Nc = 130,Nλ= 128, Nµ=64 Length based CaRS Area based CaRS

  24. Results – Mapping between a coarse RLL grid and a fine cubed-sphere grid Nc = 130,Nλ= 128, Nµ=64 (a) Field :: Y3216 ,Mapping :: RLL grid  Cubed-sphere grid (b) Field :: Vortex, Mapping :: Cubed-sphere grid  RLL grid

  25. Conclusions • The length based CaRS is true representation of the one-dimensional remapping and it could be easily employed in more complicated problems. • We need to sacrifice little bit of accuracy while employing the length based CaRS compared to the area based CaRS.

  26. Future works RLL grid RLL grid - longitude RLL grid - latitude Cubed- sphere grid :: top panel • By projecting the spherical grids on a Tangent plane, the lengths will be measured for the cascade remapping. • Area preserving projection. • Both grid lines in unified co-ordinate system.

  27. Appendix A - Remapping between 1D grids • Given source grid and cell average values on source grid remap to target grid. • Reconstruction of sub-grid scale distribution with mass-conservation and monotonicity as constraints. • Piecewise Constant Method (PCoM) • Piecewise Linear Method (PLM) • Piecewise Parabolic Method (PPM) • Piecewise Cubic Method (PCM)

  28. Appendix B - Length measured on different projections Nc = 33,Nλ= 128, Nµ=64 Field = Y3216 , Reconstruction :: PCM First sweep – based on length measured on (λ,µ=sinθ) plane Second sweep – based on length measured on (λ,µ=sinθ) plane First sweep – based on length measured on (λ,θ) plane Second sweep – based on length measured on (λ,θ) plane First sweep – based on arc length Second sweep – based on arc length

More Related