Loading…
The mass-conservative cell-integrated semi-Lagrangian advection scheme on the sphere
A mass-conservative cell-integrated semi-Lagrangian (CISL) scheme is presented and tested for 2D transport on the sphere. The total mass is conserved exactly and the mass of each individual grid cell is conserved in general. The scheme is based on a general scheme developed by Machenhauer and Olk th...
Saved in:
Published in: | Monthly weather review 2002-03, Vol.130 (3), p.649-667 |
---|---|
Main Authors: | , |
Format: | Article |
Language: | English |
Subjects: | |
Citations: | Items that this one cites |
Online Access: | Get full text |
Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
cited_by | |
---|---|
cites | cdi_FETCH-LOGICAL-c409t-7cdd6fe934581317300598c9602c909a3482e4bc41ad57867fc9a1deb256dbc83 |
container_end_page | 667 |
container_issue | 3 |
container_start_page | 649 |
container_title | Monthly weather review |
container_volume | 130 |
creator | NAIR, Ramachandran D MACHENHAUER, Bennert |
description | A mass-conservative cell-integrated semi-Lagrangian (CISL) scheme is presented and tested for 2D transport on the sphere. The total mass is conserved exactly and the mass of each individual grid cell is conserved in general. The scheme is based on a general scheme developed by Machenhauer and Olk that has increased cost effectiveness without loss of accuracy, compared to the CISL scheme of Rancic. A regular latitude-longitude grid is used on the sphere and upstream trajectories from the corner points of the regular grid cells (the Eulerian cells) define the corner points of the departure cells. The sides in these so-called Lagrangian cells are generally defined as straight lines in a ( lambda , mu ) plane, where lambda is the longitude and mu is the sine of the latitude. The mass distribution within each Eulerian grid cell is defined by quasi-biparabolic functions, which are used to integrate analytically the mass in each Lagrangian computational cell. The auxiliary computational cells are polygons with each side parallel to the coordinate axis. Also, the computational cells have the same area as the Lagrangian cells they approximate. They were introduced in order to simplify the analytical integrals of mass. Near the poles, the east and west sides of certain Lagrangian cells cannot be approximated by straight lines in the ( lambda , mu ) plane, and are instead represented by straight lines in polar tangent plane coordinates. Each of the latitudinal belts of Lagrangian cells in the polar caps are split up into several latitudinal belts of subcells, which can be approximated by computational cells as in the case of cells closer to the equator. One latitudinal belt in each hemisphere, which encloses the Eulerian pole (singular belt), is treated in a special way. First the total mass in the singular belt is determined and then it is redistributed among the cells in the belt using weights determined by a traditional SL scheme at the midpoints of the cells. By this procedure the total mass is still conserved while the conservation is only approximately maintained for the individual cells in the singular belt. These special treatments in the polar caps fit well into the general structure of the code and can be implemented with minor modifications in the code used for the rest of the sphere. Compared to two other conservative advection schemes implemented on the sphere the CISL scheme used here was found to be competitive in terms of accuracy for the same resolut |
doi_str_mv | 10.1175/1520-0493(2002)130<0649:tmccis>2.0.co;2 |
format | article |
fullrecord | <record><control><sourceid>proquest_cross</sourceid><recordid>TN_cdi_proquest_miscellaneous_18294134</recordid><sourceformat>XML</sourceformat><sourcesystem>PC</sourcesystem><sourcerecordid>18294134</sourcerecordid><originalsourceid>FETCH-LOGICAL-c409t-7cdd6fe934581317300598c9602c909a3482e4bc41ad57867fc9a1deb256dbc83</originalsourceid><addsrcrecordid>eNpFkFtLwzAUgIMoOKf_oS-KPnQ7ufQSFUGKzsFkD9bnkKWnW6SX2XQD_70pG_qUE_j4zuEjZEphQmkSTWnEIAQh-S0DYHeUwyPEQt73tTHWPbEJTEz7wE7I6I88JSOPJqHnxDm5cO4LAOJYsBHJ8w0GtXYuNG3jsNvr3u4xMFhVoW16XHe6xyJwWNtwof2vWVvdBLrYo-lt2wTObLDGwE-9N7ntBju8JGelrhxeHd8x-Xx9ybO3cLGczbPnRWgEyD5MTFHEJUouopRymnCASKZGxsCMBKm5SBmKlRFUF1GSxklppKYFrlgUFyuT8jG5OXi3Xfu9Q9er2rrhdN1gu3OKpkwKyoUHZwfQdK1zHZZq29ladz-KghqqqqGVGlqpoaryVdVQVeXvWTb_UEyBypaKedP1caV2RlelD-Kr_-u4kGkUS_4LPGt8Ug</addsrcrecordid><sourcetype>Aggregation Database</sourcetype><iscdi>true</iscdi><recordtype>article</recordtype><pqid>18294134</pqid></control><display><type>article</type><title>The mass-conservative cell-integrated semi-Lagrangian advection scheme on the sphere</title><source>EZB Electronic Journals Library</source><creator>NAIR, Ramachandran D ; MACHENHAUER, Bennert</creator><creatorcontrib>NAIR, Ramachandran D ; MACHENHAUER, Bennert</creatorcontrib><description>A mass-conservative cell-integrated semi-Lagrangian (CISL) scheme is presented and tested for 2D transport on the sphere. The total mass is conserved exactly and the mass of each individual grid cell is conserved in general. The scheme is based on a general scheme developed by Machenhauer and Olk that has increased cost effectiveness without loss of accuracy, compared to the CISL scheme of Rancic. A regular latitude-longitude grid is used on the sphere and upstream trajectories from the corner points of the regular grid cells (the Eulerian cells) define the corner points of the departure cells. The sides in these so-called Lagrangian cells are generally defined as straight lines in a ( lambda , mu ) plane, where lambda is the longitude and mu is the sine of the latitude. The mass distribution within each Eulerian grid cell is defined by quasi-biparabolic functions, which are used to integrate analytically the mass in each Lagrangian computational cell. The auxiliary computational cells are polygons with each side parallel to the coordinate axis. Also, the computational cells have the same area as the Lagrangian cells they approximate. They were introduced in order to simplify the analytical integrals of mass. Near the poles, the east and west sides of certain Lagrangian cells cannot be approximated by straight lines in the ( lambda , mu ) plane, and are instead represented by straight lines in polar tangent plane coordinates. Each of the latitudinal belts of Lagrangian cells in the polar caps are split up into several latitudinal belts of subcells, which can be approximated by computational cells as in the case of cells closer to the equator. One latitudinal belt in each hemisphere, which encloses the Eulerian pole (singular belt), is treated in a special way. First the total mass in the singular belt is determined and then it is redistributed among the cells in the belt using weights determined by a traditional SL scheme at the midpoints of the cells. By this procedure the total mass is still conserved while the conservation is only approximately maintained for the individual cells in the singular belt. These special treatments in the polar caps fit well into the general structure of the code and can be implemented with minor modifications in the code used for the rest of the sphere. Compared to two other conservative advection schemes implemented on the sphere the CISL scheme used here was found to be competitive in terms of accuracy for the same resolution. In addition the CISL scheme has the advantage over these schemes that it is applicable for Courant numbers larger than one. In plane geometry the scheme of Rancic had an overhead factor of 2.5 in CPU time compared to a traditional bicubic semi-Lagrangian scheme. This factor is reduced to 1.1 for the Machenhauer and Olk scheme on the plane while on the sphere the factor is found to be 1.28 for the present scheme. This overhead seems to be a reasonable price to pay for increased accuracy and exact mass conservation.</description><identifier>ISSN: 0027-0644</identifier><identifier>EISSN: 1520-0493</identifier><identifier>DOI: 10.1175/1520-0493(2002)130<0649:tmccis>2.0.co;2</identifier><identifier>CODEN: MWREAB</identifier><language>eng</language><publisher>Boston, MA: American Meteorological Society</publisher><subject>Earth, ocean, space ; Exact sciences and technology ; External geophysics ; Geophysics. Techniques, methods, instrumentation and models</subject><ispartof>Monthly weather review, 2002-03, Vol.130 (3), p.649-667</ispartof><rights>2002 INIST-CNRS</rights><lds50>peer_reviewed</lds50><oa>free_for_read</oa><woscitedreferencessubscribed>false</woscitedreferencessubscribed><cites>FETCH-LOGICAL-c409t-7cdd6fe934581317300598c9602c909a3482e4bc41ad57867fc9a1deb256dbc83</cites></display><links><openurl>$$Topenurl_article</openurl><openurlfulltext>$$Topenurlfull_article</openurlfulltext><thumbnail>$$Tsyndetics_thumb_exl</thumbnail><link.rule.ids>314,780,784,27924,27925</link.rule.ids><backlink>$$Uhttp://pascal-francis.inist.fr/vibad/index.php?action=getRecordDetail&idt=13498569$$DView record in Pascal Francis$$Hfree_for_read</backlink></links><search><creatorcontrib>NAIR, Ramachandran D</creatorcontrib><creatorcontrib>MACHENHAUER, Bennert</creatorcontrib><title>The mass-conservative cell-integrated semi-Lagrangian advection scheme on the sphere</title><title>Monthly weather review</title><description>A mass-conservative cell-integrated semi-Lagrangian (CISL) scheme is presented and tested for 2D transport on the sphere. The total mass is conserved exactly and the mass of each individual grid cell is conserved in general. The scheme is based on a general scheme developed by Machenhauer and Olk that has increased cost effectiveness without loss of accuracy, compared to the CISL scheme of Rancic. A regular latitude-longitude grid is used on the sphere and upstream trajectories from the corner points of the regular grid cells (the Eulerian cells) define the corner points of the departure cells. The sides in these so-called Lagrangian cells are generally defined as straight lines in a ( lambda , mu ) plane, where lambda is the longitude and mu is the sine of the latitude. The mass distribution within each Eulerian grid cell is defined by quasi-biparabolic functions, which are used to integrate analytically the mass in each Lagrangian computational cell. The auxiliary computational cells are polygons with each side parallel to the coordinate axis. Also, the computational cells have the same area as the Lagrangian cells they approximate. They were introduced in order to simplify the analytical integrals of mass. Near the poles, the east and west sides of certain Lagrangian cells cannot be approximated by straight lines in the ( lambda , mu ) plane, and are instead represented by straight lines in polar tangent plane coordinates. Each of the latitudinal belts of Lagrangian cells in the polar caps are split up into several latitudinal belts of subcells, which can be approximated by computational cells as in the case of cells closer to the equator. One latitudinal belt in each hemisphere, which encloses the Eulerian pole (singular belt), is treated in a special way. First the total mass in the singular belt is determined and then it is redistributed among the cells in the belt using weights determined by a traditional SL scheme at the midpoints of the cells. By this procedure the total mass is still conserved while the conservation is only approximately maintained for the individual cells in the singular belt. These special treatments in the polar caps fit well into the general structure of the code and can be implemented with minor modifications in the code used for the rest of the sphere. Compared to two other conservative advection schemes implemented on the sphere the CISL scheme used here was found to be competitive in terms of accuracy for the same resolution. In addition the CISL scheme has the advantage over these schemes that it is applicable for Courant numbers larger than one. In plane geometry the scheme of Rancic had an overhead factor of 2.5 in CPU time compared to a traditional bicubic semi-Lagrangian scheme. This factor is reduced to 1.1 for the Machenhauer and Olk scheme on the plane while on the sphere the factor is found to be 1.28 for the present scheme. This overhead seems to be a reasonable price to pay for increased accuracy and exact mass conservation.</description><subject>Earth, ocean, space</subject><subject>Exact sciences and technology</subject><subject>External geophysics</subject><subject>Geophysics. Techniques, methods, instrumentation and models</subject><issn>0027-0644</issn><issn>1520-0493</issn><fulltext>true</fulltext><rsrctype>article</rsrctype><creationdate>2002</creationdate><recordtype>article</recordtype><recordid>eNpFkFtLwzAUgIMoOKf_oS-KPnQ7ufQSFUGKzsFkD9bnkKWnW6SX2XQD_70pG_qUE_j4zuEjZEphQmkSTWnEIAQh-S0DYHeUwyPEQt73tTHWPbEJTEz7wE7I6I88JSOPJqHnxDm5cO4LAOJYsBHJ8w0GtXYuNG3jsNvr3u4xMFhVoW16XHe6xyJwWNtwof2vWVvdBLrYo-lt2wTObLDGwE-9N7ntBju8JGelrhxeHd8x-Xx9ybO3cLGczbPnRWgEyD5MTFHEJUouopRymnCASKZGxsCMBKm5SBmKlRFUF1GSxklppKYFrlgUFyuT8jG5OXi3Xfu9Q9er2rrhdN1gu3OKpkwKyoUHZwfQdK1zHZZq29ladz-KghqqqqGVGlqpoaryVdVQVeXvWTb_UEyBypaKedP1caV2RlelD-Kr_-u4kGkUS_4LPGt8Ug</recordid><startdate>20020301</startdate><enddate>20020301</enddate><creator>NAIR, Ramachandran D</creator><creator>MACHENHAUER, Bennert</creator><general>American Meteorological Society</general><scope>IQODW</scope><scope>AAYXX</scope><scope>CITATION</scope><scope>7TG</scope><scope>KL.</scope></search><sort><creationdate>20020301</creationdate><title>The mass-conservative cell-integrated semi-Lagrangian advection scheme on the sphere</title><author>NAIR, Ramachandran D ; MACHENHAUER, Bennert</author></sort><facets><frbrtype>5</frbrtype><frbrgroupid>cdi_FETCH-LOGICAL-c409t-7cdd6fe934581317300598c9602c909a3482e4bc41ad57867fc9a1deb256dbc83</frbrgroupid><rsrctype>articles</rsrctype><prefilter>articles</prefilter><language>eng</language><creationdate>2002</creationdate><topic>Earth, ocean, space</topic><topic>Exact sciences and technology</topic><topic>External geophysics</topic><topic>Geophysics. Techniques, methods, instrumentation and models</topic><toplevel>peer_reviewed</toplevel><toplevel>online_resources</toplevel><creatorcontrib>NAIR, Ramachandran D</creatorcontrib><creatorcontrib>MACHENHAUER, Bennert</creatorcontrib><collection>Pascal-Francis</collection><collection>CrossRef</collection><collection>Meteorological & Geoastrophysical Abstracts</collection><collection>Meteorological & Geoastrophysical Abstracts - Academic</collection><jtitle>Monthly weather review</jtitle></facets><delivery><delcategory>Remote Search Resource</delcategory><fulltext>fulltext</fulltext></delivery><addata><au>NAIR, Ramachandran D</au><au>MACHENHAUER, Bennert</au><format>journal</format><genre>article</genre><ristype>JOUR</ristype><atitle>The mass-conservative cell-integrated semi-Lagrangian advection scheme on the sphere</atitle><jtitle>Monthly weather review</jtitle><date>2002-03-01</date><risdate>2002</risdate><volume>130</volume><issue>3</issue><spage>649</spage><epage>667</epage><pages>649-667</pages><issn>0027-0644</issn><eissn>1520-0493</eissn><coden>MWREAB</coden><abstract>A mass-conservative cell-integrated semi-Lagrangian (CISL) scheme is presented and tested for 2D transport on the sphere. The total mass is conserved exactly and the mass of each individual grid cell is conserved in general. The scheme is based on a general scheme developed by Machenhauer and Olk that has increased cost effectiveness without loss of accuracy, compared to the CISL scheme of Rancic. A regular latitude-longitude grid is used on the sphere and upstream trajectories from the corner points of the regular grid cells (the Eulerian cells) define the corner points of the departure cells. The sides in these so-called Lagrangian cells are generally defined as straight lines in a ( lambda , mu ) plane, where lambda is the longitude and mu is the sine of the latitude. The mass distribution within each Eulerian grid cell is defined by quasi-biparabolic functions, which are used to integrate analytically the mass in each Lagrangian computational cell. The auxiliary computational cells are polygons with each side parallel to the coordinate axis. Also, the computational cells have the same area as the Lagrangian cells they approximate. They were introduced in order to simplify the analytical integrals of mass. Near the poles, the east and west sides of certain Lagrangian cells cannot be approximated by straight lines in the ( lambda , mu ) plane, and are instead represented by straight lines in polar tangent plane coordinates. Each of the latitudinal belts of Lagrangian cells in the polar caps are split up into several latitudinal belts of subcells, which can be approximated by computational cells as in the case of cells closer to the equator. One latitudinal belt in each hemisphere, which encloses the Eulerian pole (singular belt), is treated in a special way. First the total mass in the singular belt is determined and then it is redistributed among the cells in the belt using weights determined by a traditional SL scheme at the midpoints of the cells. By this procedure the total mass is still conserved while the conservation is only approximately maintained for the individual cells in the singular belt. These special treatments in the polar caps fit well into the general structure of the code and can be implemented with minor modifications in the code used for the rest of the sphere. Compared to two other conservative advection schemes implemented on the sphere the CISL scheme used here was found to be competitive in terms of accuracy for the same resolution. In addition the CISL scheme has the advantage over these schemes that it is applicable for Courant numbers larger than one. In plane geometry the scheme of Rancic had an overhead factor of 2.5 in CPU time compared to a traditional bicubic semi-Lagrangian scheme. This factor is reduced to 1.1 for the Machenhauer and Olk scheme on the plane while on the sphere the factor is found to be 1.28 for the present scheme. This overhead seems to be a reasonable price to pay for increased accuracy and exact mass conservation.</abstract><cop>Boston, MA</cop><pub>American Meteorological Society</pub><doi>10.1175/1520-0493(2002)130<0649:tmccis>2.0.co;2</doi><tpages>19</tpages><oa>free_for_read</oa></addata></record> |
fulltext | fulltext |
identifier | ISSN: 0027-0644 |
ispartof | Monthly weather review, 2002-03, Vol.130 (3), p.649-667 |
issn | 0027-0644 1520-0493 |
language | eng |
recordid | cdi_proquest_miscellaneous_18294134 |
source | EZB Electronic Journals Library |
subjects | Earth, ocean, space Exact sciences and technology External geophysics Geophysics. Techniques, methods, instrumentation and models |
title | The mass-conservative cell-integrated semi-Lagrangian advection scheme on the sphere |
url | http://sfxeu10.hosted.exlibrisgroup.com/loughborough?ctx_ver=Z39.88-2004&ctx_enc=info:ofi/enc:UTF-8&ctx_tim=2025-01-07T17%3A59%3A50IST&url_ver=Z39.88-2004&url_ctx_fmt=infofi/fmt:kev:mtx:ctx&rfr_id=info:sid/primo.exlibrisgroup.com:primo3-Article-proquest_cross&rft_val_fmt=info:ofi/fmt:kev:mtx:journal&rft.genre=article&rft.atitle=The%20mass-conservative%20cell-integrated%20semi-Lagrangian%20advection%20scheme%20on%20the%20sphere&rft.jtitle=Monthly%20weather%20review&rft.au=NAIR,%20Ramachandran%20D&rft.date=2002-03-01&rft.volume=130&rft.issue=3&rft.spage=649&rft.epage=667&rft.pages=649-667&rft.issn=0027-0644&rft.eissn=1520-0493&rft.coden=MWREAB&rft_id=info:doi/10.1175/1520-0493(2002)130%3C0649:tmccis%3E2.0.co;2&rft_dat=%3Cproquest_cross%3E18294134%3C/proquest_cross%3E%3Cgrp_id%3Ecdi_FETCH-LOGICAL-c409t-7cdd6fe934581317300598c9602c909a3482e4bc41ad57867fc9a1deb256dbc83%3C/grp_id%3E%3Coa%3E%3C/oa%3E%3Curl%3E%3C/url%3E&rft_id=info:oai/&rft_pqid=18294134&rft_id=info:pmid/&rfr_iscdi=true |