Difference between revisions of "Olson land map"
(New page: On this page we discuss the land map from Jennifer Olson that is currently used by the GEOS-Chem Dry deposition and Soil NOx emissions modules. == Overview == The GEOS-Chem [[Dry...) |
(→Olson land map values are zero when using GCHP compiled with GNU Fortran) |
||
(115 intermediate revisions by 4 users not shown) | |||
Line 1: | Line 1: | ||
On this page we discuss the land map from Jennifer Olson that is currently used by the GEOS-Chem [[Dry deposition]] and [[Soil NOx emissions]] modules. | On this page we discuss the land map from Jennifer Olson that is currently used by the GEOS-Chem [[Dry deposition]] and [[Soil NOx emissions]] modules. | ||
− | == | + | == Olson 2001 land map == |
− | The GEOS-Chem [[Dry deposition]] | + | The Olson 2001 land map replaces the Olson 1992 land map in GEOS-Chem for versions after v9-01-03. |
+ | |||
+ | === Overview === | ||
+ | |||
+ | The Olson 2001 land map has a native resolution of 1km x 1km, but has been binned to 0.25° x 0.25° by Eloise Marais for the use with the [[GEOS-FP|GEOS-FP met fields]] (in order to facilitate the SEAC<sup>4</sup>RS campaign). The Olson 2001 land map is now used with all Met sources in GEOS-Chem. | ||
+ | |||
+ | The Olson 2001 land map is described at [https://lta.cr.usgs.gov/GLCC this website]. | ||
+ | |||
+ | '''''[mailto:emarais@gmail.com Eloise Marais] wrote:''''' | ||
+ | |||
+ | <blockquote> | ||
+ | <p>The Olson 2001 landcover map is at a native resolution of 1km x 1km. I've identified the dominant vegetation types in each 0.25° x 0.25° degree gridbox and use this as input to GEOS-Chem.</p> | ||
+ | |||
+ | <p>The Olson 2001 landcover map also has 96 vegetation types compared with 74 for Olson 1992. For the most part vegetation types 75-96 are either not dominant vegetation types at 0.25° x 0.25° or they are crop types that I lump with other similar vegetation types (either crops or mixed forest/field vegetation). Therefore, the rebinned Olson 2001 landcover data set at 0.25° x 0.25° degrees has 74 vegetation types (adhering to the same indexing as in Olson 1992).</p> | ||
+ | |||
+ | <p>There are also new vegetation types that are defined in the Olson 2001 data set from 1-74 that were previously listed as "not used" in <tt>drydep.table</tt>. These are assigned appropriate deposition ID # and z0 values.</p> | ||
+ | |||
+ | <p>Vegetation types that are listed as "not used" in the updated <tt>drydep.table</tt> dataset are those that are not dominant at 0.25° x 0.25°, but may be present in the 1km x 1km data set.</p> | ||
+ | |||
+ | <p>The largest differences to O3 dry deposition flux for a 1-day run in June 2007 at 4x5 degree resolution include the US, and the western portion of South America.</p> | ||
+ | </blockquote> | ||
+ | |||
+ | There are actually 73 land types, not 74, in both the 1992 and 2001 Olson land maps. The table below shows the translation between the Olson 2001 and Olson 1992 land maps: | ||
+ | |||
+ | Olson 2001 Olson 1992 # in | ||
+ | LC # Description Equivalent Dry deposition | ||
+ | ============================================================================== | ||
+ | 1 Urban 1 2 | ||
+ | 2 Low Sparse Grassland 2 3 | ||
+ | 3 Coniferous Forest 3 4 | ||
+ | 4 Deciduous Conifer Forest 4 5 | ||
+ | 5 Deciduous Broadleaf Forest 5 6 | ||
+ | 6 Evergreen Broadleaf Forests 6 7 | ||
+ | 7 Tall Grasses and Shrubs 7 8 | ||
+ | 8 Bare Desert 8 9 | ||
+ | 9 Upland Tundra 9 10 | ||
+ | 10 Irrigated Grassland 10 11 | ||
+ | 11 Semi Desert 11 12 | ||
+ | 12 Glacier Ice 12 13 | ||
+ | 13 Wooded Wet Swamp 13 14 | ||
+ | 14 Inland Water 0 1 | ||
+ | 15 Sea Water 0 1 | ||
+ | 16 Shrub Evergreen 16 17 | ||
+ | 17 Shrub Deciduous 18 19 | ||
+ | 18 Mixed Forest and Field none present | ||
+ | 19 Evergreen Forest and Fields 19 20 | ||
+ | 20 Cool Rain Forest 20 21 | ||
+ | 21 Conifer Boreal Forest 21 22 | ||
+ | 22 Cool Conifer Forest 22 23 | ||
+ | 23 Cool Mixed Forest 23 24 | ||
+ | 24 Mixed Forest 24 25 | ||
+ | 25 Cool Broadleaf Forest 25 26 | ||
+ | 26 Deciduous Broadleaf Forest 26 27 | ||
+ | 27 Conifer Forest 27 28 | ||
+ | 28 Montane Tropical Forests 28 29 | ||
+ | 29 Seasonal Tropical Forest 29 30 | ||
+ | 30 Cool Crops and Towns 30 31 | ||
+ | 31 Crops and Town 31 32 | ||
+ | 32 Dry Tropical Woods 32 33 | ||
+ | 33 Tropical Rainforest 33 34 | ||
+ | 34 Tropical Degraded Forest 34 35 | ||
+ | 35 Corn and Beans Cropland 35 36 | ||
+ | 36 Rice Paddy and Field 36 37 | ||
+ | 37 Hot Irrigated Cropland 37 38 | ||
+ | 38 Cool Irrigated Cropland 38 39 | ||
+ | 39 Cold Irrigated Cropland none present | ||
+ | 40 Cool Grasses and Shrubs 40 41 | ||
+ | 41 Hot and Mild Grasses and Shrubs 41 42 | ||
+ | 42 Cold Grassland 42 43 | ||
+ | 43 Savanna (Woods) 43 44 | ||
+ | 44 Mire, Bog, Fen 44 45 | ||
+ | 45 Marsh Wetland 45 46 | ||
+ | 46 Mediterranean Scrub 46 47 | ||
+ | 47 Dry Woody Scrub 47 48 | ||
+ | 48 Dry Evergreen Woods none present | ||
+ | 49 Volcanic Rock none present | ||
+ | 50 Sand Desert none present | ||
+ | 51 Semi Desert Shrubs 51 52 | ||
+ | 52 Semi Desert Sage 52 53 | ||
+ | 53 Barren Tundra 53 54 | ||
+ | 54 Cool Southern Hemisphere Mixed Forests 54 55 | ||
+ | 55 Cool Fields and Woods 55 56 | ||
+ | 56 Forest and Field 56 57 | ||
+ | 57 Cool Forest and Field 57 58 | ||
+ | 58 Fields and Woody Savanna 58 59 | ||
+ | 59 Succulent and Thorn Scrub 59 60 | ||
+ | 60 Small Leaf Mixed Woods 60 61 | ||
+ | 61 Deciduous and Mixed Boreal Forest 61 62 | ||
+ | 62 Narrow Conifers 62 63 | ||
+ | 63 Wooded Tundra 63 64 | ||
+ | 64 Heath Scrub 64 65 | ||
+ | 65 Coastal Wetland, NW none present | ||
+ | 66 Coastal Wetland, NE none present | ||
+ | 67 Coastal Wetland, SE none present | ||
+ | 68 Coastal Wetland, SW none present | ||
+ | 69 Polar and Alpine Desert 69 70 | ||
+ | 70 Glacier Rock none present | ||
+ | 71 Salt Playas none present | ||
+ | 72 Mangrove 72 73 | ||
+ | 73 Water and Island Fringe none present | ||
+ | 74 Land, Water, and Shore (see Note 1) none present | ||
+ | 75 Land and Water, Rivers (see Note 1) none present | ||
+ | 76 Crop and Water Mixtures 36 37 | ||
+ | 77 Southern Hemisphere Conifers none present | ||
+ | 78 Southern Hemisphere Mixed Forest 32 33 | ||
+ | 79 Wet Sclerophylic Forest 26 27 | ||
+ | 80 Coastline Fringe none present | ||
+ | 81 Beaches and Dunes none present | ||
+ | 82 Sparse Dunes and Ridges none present | ||
+ | 83 Bare Coastal Dunes none present | ||
+ | 84 Residual Dunes and Beaches none present | ||
+ | 85 Compound Coastlines none present | ||
+ | 86 Rocky Cliffs and Slopes none present | ||
+ | 87 Sandy Grassland and Shrubs none present | ||
+ | 88 Bamboo none present | ||
+ | 89 Moist Eucalyptus 26 27 | ||
+ | 90 Rain Green Tropical Forest 33 34 | ||
+ | 91 Woody Savanna 43 44 | ||
+ | 92 Broadleaf Crops 29 30 | ||
+ | 93 Grass Crops 41 42 | ||
+ | 94 Crops, Grass, Shrubs 41 42 | ||
+ | 95 Evergreen Tree Crop 33 34 | ||
+ | 96 Deciduous Tree Crop 33 34 | ||
+ | |||
+ | The rebinned Olson 2001 land types at 0.25° x 0.25° resolution (73 land types total) are plotted here: | ||
+ | |||
+ | [[Image:Olson2001.png]] | ||
+ | |||
+ | --[[User:Lizzie Lundgren|Lizzie Lundgren]] ([[User talk:Lizzie Lundgren|talk]]) 17:45, 2 November 2016 (UTC) | ||
+ | |||
+ | === Implementation notes === | ||
+ | |||
+ | Adoption of the Olson 2001 land map requires the following updates to the GEOS-Chem source code & input files: | ||
+ | |||
+ | #<p>Parameter <tt>NNNTYPE</tt> in <tt>drydep_mod.F</tt> needs to be increased from 15 to 25.</p> | ||
+ | #*<p>'''UPDATE (21 May 2012):''' Parameter <tt>NNNTYPE</tt> was removed from [[GEOS-Chem v9-01-03 benchmark history#v9-01-03j|v9-01-03j]] (17 Apr 2012). It is replaced by parameter <tt>NNTYPE</tt> in <tt>CMN_SIZE_mod.F</tt>.</p> | ||
+ | #<p>Parameter <tt>NTYPE</tt> in <tt>CMN_SIZE_mod.F</tt> needs to be incresaed from 15 to 25.</p> | ||
+ | #*<p>'''UPDATE (21 May 2012):''' This was done in v9-01-03j (17 Apr 2012)</p> | ||
+ | #<p>Parameter <tt>NNNVEGTYPE</tt> can be unchanged (it currently = 74). The rebinned Olson 2001 map at 0.25° x 0.25° uses the same indexing as Olson 1992.</p> | ||
+ | #*<p>'''UPDATE (21 May 2012):''' Parameter <tt>NNNVEGTYPE</tt> was removed from v9-01-03j (17 Apr 2012). It is replaced by parameter <tt>NNVEGTYPE</tt> in <tt>CMN_SIZE_mod.F</tt>.</p> | ||
+ | #<p>A corresponding new <tt>drydep.table</tt> needs to be adopted for the Olson 2001 land map.</p> | ||
+ | #*<p>'''UPDATE (21 May 2012):''' This has been done in v9-01-03j (17 Apr 2012).</p> | ||
+ | |||
+ | Also, we will not create <tt>vegtype.global</tt> files for the Olson 2001 land map, but will instead compute the corresponding values of the <tt>IREG</tt>, <tt>IUSE</tt>, <tt>ILAND</tt>, etc. arrays on-the-fly, [[#Modifications for the grid-independent model|as described in the next section]]. | ||
+ | |||
+ | # The dominant changes in going from the Olson 1992 to the Olson 2001 include additional crops not previously defined in Olson 1992 that vary from tree crops to grass crops. Another difference is a further subdivision of torpcial type vegetation into degraded vegetation and mixed forest and fields. These differences account for the discrepancies between Olson 1992 and Olson 2001 dry deposition rates in the figures found here on the [[Leaf area indices in GEOS-Chem]] page. | ||
+ | # The most appropriate dry deposition ID (that similar to the Olson 1992 landcover types) and correspondingly dry deposition parameters are assigned to the new vegetation types that are included in Olson 2001. | ||
+ | |||
+ | --[[User:Emarais|Emarais]] 16:47, 13 April 2012 (EDT)<br> | ||
+ | --[[User:Bmy|Bob Y.]] 14:23, 21 May 2012 (EDT) | ||
+ | |||
+ | === Updates to dry deposition when using the Olson 2001 land map === | ||
+ | |||
+ | <span style="color:green">'''''This update was tested in the 1-month benchmark simulation [[GEOS-Chem v10-01 benchmark history#v10-01a|GEOS-Chem v10-01a]] and approved on 05 Feb 2014.'''''</span> | ||
+ | |||
+ | Patrick Kim implemented updates to improve dry deposition of ozone when using the Olson 2001 land map. For a full description, please see [[Dry deposition#Updates to dry deposition when using the Olson 2001 land map|this post on our ''Dry deposition'' wiki page]]. | ||
+ | |||
+ | --[[User:Melissa Payer|Melissa Sulprizio]] 17:45, 3 February 2014 (EST) | ||
== Modifications for the grid-independent model == | == Modifications for the grid-independent model == | ||
+ | |||
+ | The [[Grid-independent GEOS-Chem]] project seeks to interface the GEOS-Chem model directly into a General Circulation Model (such as NASA's GEOS-5 GCM). In order for this to happen, much of the existing GEOS-Chem source code will have to be restructured. The [[GEOS-Chem Support Team]] is actively working towards this goal, which involves both (1) building a new interface that will link GEOS-Chem to a GCM, and (2) removing legacy code that will hamper grid-independent functionality. | ||
+ | |||
+ | [[Grid-independent GEOS-Chem#Action items identified after initial test with Beijing GCM|During a recent test interfacing of GEOS-Chem with the Beijing Climate Center's CAM/GCM]], we discovered that the existing <tt>vegtype.global</tt> files (and the corresponding <tt>rdland.F</tt> routine) was not compatible with the GCM, for several reasons: | ||
+ | |||
+ | # The GCM could not read in an ASCII file. It expected input to be in netCDF format. | ||
+ | # The <tt>vegtype.global</tt> files used by GEOS-Chem were hardwired to particular grids (4° x 5°, 2° x 2.5° and the 0.5° x 0.666° nested grids). However, the GCM was using a T42 spectral grid, for which we had not created a corresponding <tt>vegtype.global</tt> file. | ||
+ | # The <tt>vegtype.global</tt> files were created by offline pre-processing. Instead, we need to employ an on-the-fly regridding algorithm to translate the Olson 0.5° x 0.5° land map to the <tt>IREG</tt>, <tt>ILAND</tt>, <tt>IUSE</tt>, etc. arrays on the GEOS-Chem grid. | ||
+ | |||
+ | To this end, we have created a new GEOS-Chem module <tt>olson_landmap_mod.F90</tt>. This module, which replaces the <tt>rdland.F</tt> routine, does the following: | ||
+ | |||
+ | # Reads in the Olson 1992 land map at 0.5° x 0.5° resolution in netCDF format, instead of ASCII | ||
+ | # Constructs the <tt>IREG</tt>, <tt>ILAND</tt>, <tt>IUSE</tt> arrays on the GEOS-Chem grid (for soil NOx emissions) | ||
+ | # Constructs the <tt>FRCLND</tt> array on the GEOS-Chem grid (for use throughout GEOS-Chem) | ||
+ | |||
+ | The <tt>olson_landmap_mod.F90</tt> was introduced to [[GEOS-Chem v9-01-03]]. It was validated with 1-month benchmark simulation [[GEOS-Chem v9-01-03 benchmark history#v9-01-03i|v9-01-03i]] (approved on 17 Apr 2012). We have also extended the functionality of <tt>olson_landmap_mod.F90</tt> so that it can also process the Olson 2001 land map as well as the Olson 1992 land map. | ||
+ | |||
+ | Also, certain ASCII input files to the dry deposition modules, such as: | ||
+ | |||
+ | # <tt>drydep.table</tt> | ||
+ | # <tt>drydep.coef</tt> | ||
+ | |||
+ | have been replaced by netCDF files in order to facilitate file I/O by external GCM's. | ||
+ | |||
+ | --[[User:Bmy|Bob Y.]] 14:28, 21 May 2012 (EDT) | ||
== References == | == References == | ||
Line 12: | Line 194: | ||
#Olson, J, ''World Ecosystems (WE1.4): Digital raster data on a 10 minute geographic 1080 x 2160 grid'', in <u>Global Ecosystems Database, version 1.0, Disc A</u>, edited by NOAA Natl. Geophys. Data Center, Boulder, Colorado, 1992. | #Olson, J, ''World Ecosystems (WE1.4): Digital raster data on a 10 minute geographic 1080 x 2160 grid'', in <u>Global Ecosystems Database, version 1.0, Disc A</u>, edited by NOAA Natl. Geophys. Data Center, Boulder, Colorado, 1992. | ||
#Wang, Y., D.J. Jacob, and J.A. Logan, ''Global simulation of tropospheric O3-NOx-hydrocarbon chemistry, 1. Model formulation'', <u>J. Geophys. Res.</u>, '''103''', D9, 10,713-10,726, 1998. | #Wang, Y., D.J. Jacob, and J.A. Logan, ''Global simulation of tropospheric O3-NOx-hydrocarbon chemistry, 1. Model formulation'', <u>J. Geophys. Res.</u>, '''103''', D9, 10,713-10,726, 1998. | ||
+ | #Olson 2011 land map reference may be found here: http://edc2.usgs.gov/glcc/globdoc2_0.php | ||
+ | |||
+ | == Previously resolved issues == | ||
+ | |||
+ | In this section we provide information about issues with the Olson land map code that have since been resolved. | ||
+ | |||
+ | === Speed up computations in subroutine COMPUTE_OLSON_LANDMAP === | ||
+ | |||
+ | <span style="color:green">'''''These updates were validated in the 1-month benchmark simulation [[GEOS-Chem v10-01 benchmark history#v10-01c|v10-01c]] and approved on 29 May 2014.'''''</span> | ||
+ | |||
+ | We have added a fix to speed up the computations in the parallel DO loop in routine <tt>COMPUTE_OLSON_LANDMAP</tt>. These computations are used to generate the arrays <tt>IREG</tt>, <tt>ILAND</tt>, and <tt>IUSE</tt>—which GEOS-Chem's legacy dry deposition module require as input—from the data in the Olson land map netCDF file. | ||
+ | |||
+ | We made the following source-code modifications: | ||
+ | |||
+ | (1) Add this section of code to the variable declaration section of subroutine <tt>COMPUTE_OLSON_LANDMAP</tt>: | ||
+ | |||
+ | ! | ||
+ | ! !DEFINED PARAMETERS: | ||
+ | ! | ||
+ | ! The following parameters are used to skip over Olson NATIVE GRID boxes | ||
+ | ! that are too far away from the GEOS-CHEM GRID BOX. This can speed up | ||
+ | ! the Olson computation by a factor of 100 or more! | ||
+ | ! | ||
+ | #if defined( GRID05x0666 ) || defined( GRID025x03125 ) | ||
+ | REAL*8, PARAMETER :: latThresh = 1d0 ! Lat threshold, nested grid | ||
+ | REAL*8, PARAMETER :: lonThresh = 1d0 ! Lon threshold, nested grid | ||
+ | #else | ||
+ | REAL*8, PARAMETER :: latThresh = 5d0 ! Lat threshold, global | ||
+ | REAL*8, PARAMETER :: lonThresh = 6d0 ! Lon threshold, global | ||
+ | #endif | ||
+ | |||
+ | (2) Add the code <tt>!$OMP SCHEDULE( DYNAMIC )</tt> to the OpenMP parallel loop declaration. This will tell each CPU to grab the next available grid box in the array. This will improve the load balancing of the loop. | ||
+ | |||
+ | !====================================================================== | ||
+ | ! Loop over all GEOS-CHEM GRID BOXES and initialize variables | ||
+ | !====================================================================== | ||
+ | !$OMP PARALLEL DO & | ||
+ | !$OMP DEFAULT( SHARED ) & | ||
+ | !$OMP PRIVATE( I, J, xedgeC_w, yedgeC_s, xedgeC_e ) & | ||
+ | !$OMP PRIVATE( yedgeC_n, dxdy4, sumArea, JJ, III ) & | ||
+ | !$OMP PRIVATE( dxdy, mapWt, II, xedge_w, yedge_s ) & | ||
+ | !$OMP PRIVATE( xedge_e, yedge_n, area, type, maxIuse ) & | ||
+ | !$OMP PRIVATE( sumIUse, uniqOlson, C, IG ) & | ||
+ | !$OMP SCHEDULE( DYNAMIC ) | ||
+ | |||
+ | (3) Modify this section of the parallel loop. The code should now look like this: | ||
+ | |||
+ | !=================================================================== | ||
+ | ! Find each NATIVE GRID BOX that fits into the GEOS-CHEM GRID BOX. | ||
+ | ! Keep track of the land types and coverage fractions. | ||
+ | !=================================================================== | ||
+ | |||
+ | ! Loop over latitudes on the NATIVE GRID | ||
+ | DO JJ = 1, J_OLSON | ||
+ | |||
+ | ! Latitude edges of this NATIVE GRID box | ||
+ | yedge_s = latedge(JJ ) ! S edge | ||
+ | yedge_n = latedge(JJ+1) ! N edge | ||
+ | |||
+ | !%%%%%% LATITUDE SHUNT TO REDUCE WALL TIME (bmy, 3/20/14) %%%%%%%%%% | ||
+ | !%%% | ||
+ | !%%% Skip further computations unless we are within LATTHRESH | ||
+ | !%%% degrees of the western edge of box (I,J). This prevents | ||
+ | !%%% excess computations and subroutine calls. | ||
+ | !%%% | ||
+ | IF ( ABS( yedge_s - yedgeC_s ) > latThresh ) CYCLE | ||
+ | !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | ||
+ | |||
+ | ! Loop over longitudes on the NATIVE GRID | ||
+ | DO III = 1, I_OLSON | ||
+ | |||
+ | ! Initialize | ||
+ | dxdy = 0e0 | ||
+ | mapWt = 0e0 | ||
+ | |||
+ | ! Find the NATIVE GRID longitude index for use below. Account for | ||
+ | ! the first GEOS-CHEM GRID box, which straddles the date line. | ||
+ | IF ( isGlobal .and. IG == 1 ) THEN | ||
+ | II = shiftLon(III) | ||
+ | ELSE | ||
+ | II = indLon(III) | ||
+ | ENDIF | ||
+ | |||
+ | ! Edges of this NATIVE GRID box | ||
+ | xedge_w = lonedge(II ) ! W edge | ||
+ | xedge_e = lonedge(II+1) ! E edge | ||
+ | |||
+ | ! Because the first GEOS-CHEM GRID BOX straddles the date line, | ||
+ | ! we have to adjust the W and E edges of the NATIVE GRID BOX to | ||
+ | ! be in monotonically increasing order. This will prevent | ||
+ | ! erronous results from being returned by GET_MAP_WT below. | ||
+ | IF ( isGlobal .and. IG == 1 .and. II >= shiftLon(1) ) THEN | ||
+ | xedge_w = xedge_w - 360e0 | ||
+ | xedge_e = xedge_e - 360e0 | ||
+ | ENDIF | ||
+ | |||
+ | !%%%%%% LONGITUDE SHUNT TO REDUCE WALL TIME (bmy, 3/20/14) %%%%%% | ||
+ | !%%% | ||
+ | !%%% Skip further computations unless we are within LONTHRESH | ||
+ | !%%% degrees of the western edge of box (I,J). This prevents | ||
+ | !%%% excess computations and subroutine calls. | ||
+ | !%%% | ||
+ | IF ( ABS( xedge_w - xedgeC_w ) > lonThresh ) CYCLE | ||
+ | !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | ||
+ | |||
+ | ! "Area" of the GEOS-CHEM GRID BOX in degrees (DLON * DLAT) | ||
+ | dxdy = ( xedge_e - xedge_w ) * ( yedge_n - yedge_s ) | ||
+ | |||
+ | ! Get the mapping weight (i.e. The fraction of the NATIVE | ||
+ | ! GRID BOX that lies w/in the GEOS-CHEM GRID BOX) | ||
+ | CALL GET_MAP_WT( xedge_w, xedge_e, xedgeC_w, xedgeC_e, & | ||
+ | yedge_s, yedge_n, yedgeC_s, yedgeC_n, & | ||
+ | mapWt ) | ||
+ | |||
+ | ... etc ... | ||
+ | |||
+ | These modifications will skip calling the subroutine <tt>GET_MAP_WT</tt> unless: | ||
+ | |||
+ | *The latitude of each Olson land map box is within a specified number of degrees (<tt>LATTHRESH</tt>) of the GEOS-Chem grid box latitude, and | ||
+ | *The longitude of each Olson land map box is within a specified number of degrees (<tt>LONTHRESH</tt>) of the GEOS-Chem grid box longitude. | ||
+ | |||
+ | By restricting when subroutine <tt>GET_MAP_WT</tt> is called in this way, we reduce a substantial amount of overhead. Profiling results with the Tuning and Analysis Utilities (TAU) have shown that this can result in a 100X speedup in execution time for <tt>COMPUTE_OLSON_LANDMAP</tt>, as shown by this figure. | ||
+ | |||
+ | [[Image:Olson_Speedup.png]] | ||
+ | |||
+ | The figure shows wall clock time (in seconds) for each subroutine on the main thread (CPU 0) of a 1-day GEOS-Chem simulation with [[GEOS-FP]] meteorology at 4° x 5° horizontal resolution. Before the fix was implemented, routine <tt>COMPUTE_OLSON_LANDMAP</tt> took about 11 seconds on the main CPU. After the fix, the time spent in the routine was reduced to 0.163 seconds, thus eliminating a major bottleneck. | ||
+ | |||
+ | --[[User:Bmy|Bob Y.]] 17:07, 30 May 2014 (EDT) | ||
+ | |||
+ | == Outstanding issues that have not yet been resolved == | ||
+ | |||
+ | === Olson land map values are read incorrectly when using GCHP compiled with GNU Fortran === | ||
+ | |||
+ | We have discovered an issue with the Olson land map that occurs when compiling [[GCHP]] with the GNU Fortran compiler. The Olson land map values are apparently being returned from disk (via ExtData) as all zeroes. A full description of the issue may be found here: https://github.com/geoschem/gchp/issues/15 | ||
− | + | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 20:06, 2 January 2019 (UTC) |
Latest revision as of 20:08, 2 January 2019
On this page we discuss the land map from Jennifer Olson that is currently used by the GEOS-Chem Dry deposition and Soil NOx emissions modules.
Contents
Olson 2001 land map
The Olson 2001 land map replaces the Olson 1992 land map in GEOS-Chem for versions after v9-01-03.
Overview
The Olson 2001 land map has a native resolution of 1km x 1km, but has been binned to 0.25° x 0.25° by Eloise Marais for the use with the GEOS-FP met fields (in order to facilitate the SEAC4RS campaign). The Olson 2001 land map is now used with all Met sources in GEOS-Chem.
The Olson 2001 land map is described at this website.
Eloise Marais wrote:
The Olson 2001 landcover map is at a native resolution of 1km x 1km. I've identified the dominant vegetation types in each 0.25° x 0.25° degree gridbox and use this as input to GEOS-Chem.
The Olson 2001 landcover map also has 96 vegetation types compared with 74 for Olson 1992. For the most part vegetation types 75-96 are either not dominant vegetation types at 0.25° x 0.25° or they are crop types that I lump with other similar vegetation types (either crops or mixed forest/field vegetation). Therefore, the rebinned Olson 2001 landcover data set at 0.25° x 0.25° degrees has 74 vegetation types (adhering to the same indexing as in Olson 1992).
There are also new vegetation types that are defined in the Olson 2001 data set from 1-74 that were previously listed as "not used" in drydep.table. These are assigned appropriate deposition ID # and z0 values.
Vegetation types that are listed as "not used" in the updated drydep.table dataset are those that are not dominant at 0.25° x 0.25°, but may be present in the 1km x 1km data set.
The largest differences to O3 dry deposition flux for a 1-day run in June 2007 at 4x5 degree resolution include the US, and the western portion of South America.
There are actually 73 land types, not 74, in both the 1992 and 2001 Olson land maps. The table below shows the translation between the Olson 2001 and Olson 1992 land maps:
Olson 2001 Olson 1992 # in LC # Description Equivalent Dry deposition ============================================================================== 1 Urban 1 2 2 Low Sparse Grassland 2 3 3 Coniferous Forest 3 4 4 Deciduous Conifer Forest 4 5 5 Deciduous Broadleaf Forest 5 6 6 Evergreen Broadleaf Forests 6 7 7 Tall Grasses and Shrubs 7 8 8 Bare Desert 8 9 9 Upland Tundra 9 10 10 Irrigated Grassland 10 11 11 Semi Desert 11 12 12 Glacier Ice 12 13 13 Wooded Wet Swamp 13 14 14 Inland Water 0 1 15 Sea Water 0 1 16 Shrub Evergreen 16 17 17 Shrub Deciduous 18 19 18 Mixed Forest and Field none present 19 Evergreen Forest and Fields 19 20 20 Cool Rain Forest 20 21 21 Conifer Boreal Forest 21 22 22 Cool Conifer Forest 22 23 23 Cool Mixed Forest 23 24 24 Mixed Forest 24 25 25 Cool Broadleaf Forest 25 26 26 Deciduous Broadleaf Forest 26 27 27 Conifer Forest 27 28 28 Montane Tropical Forests 28 29 29 Seasonal Tropical Forest 29 30 30 Cool Crops and Towns 30 31 31 Crops and Town 31 32 32 Dry Tropical Woods 32 33 33 Tropical Rainforest 33 34 34 Tropical Degraded Forest 34 35 35 Corn and Beans Cropland 35 36 36 Rice Paddy and Field 36 37 37 Hot Irrigated Cropland 37 38 38 Cool Irrigated Cropland 38 39 39 Cold Irrigated Cropland none present 40 Cool Grasses and Shrubs 40 41 41 Hot and Mild Grasses and Shrubs 41 42 42 Cold Grassland 42 43 43 Savanna (Woods) 43 44 44 Mire, Bog, Fen 44 45 45 Marsh Wetland 45 46 46 Mediterranean Scrub 46 47 47 Dry Woody Scrub 47 48 48 Dry Evergreen Woods none present 49 Volcanic Rock none present 50 Sand Desert none present 51 Semi Desert Shrubs 51 52 52 Semi Desert Sage 52 53 53 Barren Tundra 53 54 54 Cool Southern Hemisphere Mixed Forests 54 55 55 Cool Fields and Woods 55 56 56 Forest and Field 56 57 57 Cool Forest and Field 57 58 58 Fields and Woody Savanna 58 59 59 Succulent and Thorn Scrub 59 60 60 Small Leaf Mixed Woods 60 61 61 Deciduous and Mixed Boreal Forest 61 62 62 Narrow Conifers 62 63 63 Wooded Tundra 63 64 64 Heath Scrub 64 65 65 Coastal Wetland, NW none present 66 Coastal Wetland, NE none present 67 Coastal Wetland, SE none present 68 Coastal Wetland, SW none present 69 Polar and Alpine Desert 69 70 70 Glacier Rock none present 71 Salt Playas none present 72 Mangrove 72 73 73 Water and Island Fringe none present 74 Land, Water, and Shore (see Note 1) none present 75 Land and Water, Rivers (see Note 1) none present 76 Crop and Water Mixtures 36 37 77 Southern Hemisphere Conifers none present 78 Southern Hemisphere Mixed Forest 32 33 79 Wet Sclerophylic Forest 26 27 80 Coastline Fringe none present 81 Beaches and Dunes none present 82 Sparse Dunes and Ridges none present 83 Bare Coastal Dunes none present 84 Residual Dunes and Beaches none present 85 Compound Coastlines none present 86 Rocky Cliffs and Slopes none present 87 Sandy Grassland and Shrubs none present 88 Bamboo none present 89 Moist Eucalyptus 26 27 90 Rain Green Tropical Forest 33 34 91 Woody Savanna 43 44 92 Broadleaf Crops 29 30 93 Grass Crops 41 42 94 Crops, Grass, Shrubs 41 42 95 Evergreen Tree Crop 33 34 96 Deciduous Tree Crop 33 34
The rebinned Olson 2001 land types at 0.25° x 0.25° resolution (73 land types total) are plotted here:
--Lizzie Lundgren (talk) 17:45, 2 November 2016 (UTC)
Implementation notes
Adoption of the Olson 2001 land map requires the following updates to the GEOS-Chem source code & input files:
Parameter NNNTYPE in drydep_mod.F needs to be increased from 15 to 25.
UPDATE (21 May 2012): Parameter NNNTYPE was removed from v9-01-03j (17 Apr 2012). It is replaced by parameter NNTYPE in CMN_SIZE_mod.F.
Parameter NTYPE in CMN_SIZE_mod.F needs to be incresaed from 15 to 25.
UPDATE (21 May 2012): This was done in v9-01-03j (17 Apr 2012)
Parameter NNNVEGTYPE can be unchanged (it currently = 74). The rebinned Olson 2001 map at 0.25° x 0.25° uses the same indexing as Olson 1992.
UPDATE (21 May 2012): Parameter NNNVEGTYPE was removed from v9-01-03j (17 Apr 2012). It is replaced by parameter NNVEGTYPE in CMN_SIZE_mod.F.
A corresponding new drydep.table needs to be adopted for the Olson 2001 land map.
UPDATE (21 May 2012): This has been done in v9-01-03j (17 Apr 2012).
Also, we will not create vegtype.global files for the Olson 2001 land map, but will instead compute the corresponding values of the IREG, IUSE, ILAND, etc. arrays on-the-fly, as described in the next section.
- The dominant changes in going from the Olson 1992 to the Olson 2001 include additional crops not previously defined in Olson 1992 that vary from tree crops to grass crops. Another difference is a further subdivision of torpcial type vegetation into degraded vegetation and mixed forest and fields. These differences account for the discrepancies between Olson 1992 and Olson 2001 dry deposition rates in the figures found here on the Leaf area indices in GEOS-Chem page.
- The most appropriate dry deposition ID (that similar to the Olson 1992 landcover types) and correspondingly dry deposition parameters are assigned to the new vegetation types that are included in Olson 2001.
--Emarais 16:47, 13 April 2012 (EDT)
--Bob Y. 14:23, 21 May 2012 (EDT)
Updates to dry deposition when using the Olson 2001 land map
This update was tested in the 1-month benchmark simulation GEOS-Chem v10-01a and approved on 05 Feb 2014.
Patrick Kim implemented updates to improve dry deposition of ozone when using the Olson 2001 land map. For a full description, please see this post on our Dry deposition wiki page.
--Melissa Sulprizio 17:45, 3 February 2014 (EST)
Modifications for the grid-independent model
The Grid-independent GEOS-Chem project seeks to interface the GEOS-Chem model directly into a General Circulation Model (such as NASA's GEOS-5 GCM). In order for this to happen, much of the existing GEOS-Chem source code will have to be restructured. The GEOS-Chem Support Team is actively working towards this goal, which involves both (1) building a new interface that will link GEOS-Chem to a GCM, and (2) removing legacy code that will hamper grid-independent functionality.
During a recent test interfacing of GEOS-Chem with the Beijing Climate Center's CAM/GCM, we discovered that the existing vegtype.global files (and the corresponding rdland.F routine) was not compatible with the GCM, for several reasons:
- The GCM could not read in an ASCII file. It expected input to be in netCDF format.
- The vegtype.global files used by GEOS-Chem were hardwired to particular grids (4° x 5°, 2° x 2.5° and the 0.5° x 0.666° nested grids). However, the GCM was using a T42 spectral grid, for which we had not created a corresponding vegtype.global file.
- The vegtype.global files were created by offline pre-processing. Instead, we need to employ an on-the-fly regridding algorithm to translate the Olson 0.5° x 0.5° land map to the IREG, ILAND, IUSE, etc. arrays on the GEOS-Chem grid.
To this end, we have created a new GEOS-Chem module olson_landmap_mod.F90. This module, which replaces the rdland.F routine, does the following:
- Reads in the Olson 1992 land map at 0.5° x 0.5° resolution in netCDF format, instead of ASCII
- Constructs the IREG, ILAND, IUSE arrays on the GEOS-Chem grid (for soil NOx emissions)
- Constructs the FRCLND array on the GEOS-Chem grid (for use throughout GEOS-Chem)
The olson_landmap_mod.F90 was introduced to GEOS-Chem v9-01-03. It was validated with 1-month benchmark simulation v9-01-03i (approved on 17 Apr 2012). We have also extended the functionality of olson_landmap_mod.F90 so that it can also process the Olson 2001 land map as well as the Olson 1992 land map.
Also, certain ASCII input files to the dry deposition modules, such as:
- drydep.table
- drydep.coef
have been replaced by netCDF files in order to facilitate file I/O by external GCM's.
--Bob Y. 14:28, 21 May 2012 (EDT)
References
- Bey I., D. J. Jacob, R. M. Yantosca, J. A. Logan, B. Field, A. M. Fiore, Q. Li, H. Liu, L. J. Mickley, and M. Schultz, Global modeling of tropospheric chemistry with assimilated meteorology: Model description and evaluation, J. Geophys. Res., 106, 23,073-23,096, 2001
- Olson, J, World Ecosystems (WE1.4): Digital raster data on a 10 minute geographic 1080 x 2160 grid, in Global Ecosystems Database, version 1.0, Disc A, edited by NOAA Natl. Geophys. Data Center, Boulder, Colorado, 1992.
- Wang, Y., D.J. Jacob, and J.A. Logan, Global simulation of tropospheric O3-NOx-hydrocarbon chemistry, 1. Model formulation, J. Geophys. Res., 103, D9, 10,713-10,726, 1998.
- Olson 2011 land map reference may be found here: http://edc2.usgs.gov/glcc/globdoc2_0.php
Previously resolved issues
In this section we provide information about issues with the Olson land map code that have since been resolved.
Speed up computations in subroutine COMPUTE_OLSON_LANDMAP
These updates were validated in the 1-month benchmark simulation v10-01c and approved on 29 May 2014.
We have added a fix to speed up the computations in the parallel DO loop in routine COMPUTE_OLSON_LANDMAP. These computations are used to generate the arrays IREG, ILAND, and IUSE—which GEOS-Chem's legacy dry deposition module require as input—from the data in the Olson land map netCDF file.
We made the following source-code modifications:
(1) Add this section of code to the variable declaration section of subroutine COMPUTE_OLSON_LANDMAP:
! ! !DEFINED PARAMETERS: ! ! The following parameters are used to skip over Olson NATIVE GRID boxes ! that are too far away from the GEOS-CHEM GRID BOX. This can speed up ! the Olson computation by a factor of 100 or more! ! #if defined( GRID05x0666 ) || defined( GRID025x03125 ) REAL*8, PARAMETER :: latThresh = 1d0 ! Lat threshold, nested grid REAL*8, PARAMETER :: lonThresh = 1d0 ! Lon threshold, nested grid #else REAL*8, PARAMETER :: latThresh = 5d0 ! Lat threshold, global REAL*8, PARAMETER :: lonThresh = 6d0 ! Lon threshold, global #endif
(2) Add the code !$OMP SCHEDULE( DYNAMIC ) to the OpenMP parallel loop declaration. This will tell each CPU to grab the next available grid box in the array. This will improve the load balancing of the loop.
!====================================================================== ! Loop over all GEOS-CHEM GRID BOXES and initialize variables !====================================================================== !$OMP PARALLEL DO & !$OMP DEFAULT( SHARED ) & !$OMP PRIVATE( I, J, xedgeC_w, yedgeC_s, xedgeC_e ) & !$OMP PRIVATE( yedgeC_n, dxdy4, sumArea, JJ, III ) & !$OMP PRIVATE( dxdy, mapWt, II, xedge_w, yedge_s ) & !$OMP PRIVATE( xedge_e, yedge_n, area, type, maxIuse ) & !$OMP PRIVATE( sumIUse, uniqOlson, C, IG ) & !$OMP SCHEDULE( DYNAMIC )
(3) Modify this section of the parallel loop. The code should now look like this:
!=================================================================== ! Find each NATIVE GRID BOX that fits into the GEOS-CHEM GRID BOX. ! Keep track of the land types and coverage fractions. !=================================================================== ! Loop over latitudes on the NATIVE GRID DO JJ = 1, J_OLSON ! Latitude edges of this NATIVE GRID box yedge_s = latedge(JJ ) ! S edge yedge_n = latedge(JJ+1) ! N edge !%%%%%% LATITUDE SHUNT TO REDUCE WALL TIME (bmy, 3/20/14) %%%%%%%%%% !%%% !%%% Skip further computations unless we are within LATTHRESH !%%% degrees of the western edge of box (I,J). This prevents !%%% excess computations and subroutine calls. !%%% IF ( ABS( yedge_s - yedgeC_s ) > latThresh ) CYCLE !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ! Loop over longitudes on the NATIVE GRID DO III = 1, I_OLSON ! Initialize dxdy = 0e0 mapWt = 0e0 ! Find the NATIVE GRID longitude index for use below. Account for ! the first GEOS-CHEM GRID box, which straddles the date line. IF ( isGlobal .and. IG == 1 ) THEN II = shiftLon(III) ELSE II = indLon(III) ENDIF ! Edges of this NATIVE GRID box xedge_w = lonedge(II ) ! W edge xedge_e = lonedge(II+1) ! E edge ! Because the first GEOS-CHEM GRID BOX straddles the date line, ! we have to adjust the W and E edges of the NATIVE GRID BOX to ! be in monotonically increasing order. This will prevent ! erronous results from being returned by GET_MAP_WT below. IF ( isGlobal .and. IG == 1 .and. II >= shiftLon(1) ) THEN xedge_w = xedge_w - 360e0 xedge_e = xedge_e - 360e0 ENDIF !%%%%%% LONGITUDE SHUNT TO REDUCE WALL TIME (bmy, 3/20/14) %%%%%% !%%% !%%% Skip further computations unless we are within LONTHRESH !%%% degrees of the western edge of box (I,J). This prevents !%%% excess computations and subroutine calls. !%%% IF ( ABS( xedge_w - xedgeC_w ) > lonThresh ) CYCLE !%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% ! "Area" of the GEOS-CHEM GRID BOX in degrees (DLON * DLAT) dxdy = ( xedge_e - xedge_w ) * ( yedge_n - yedge_s ) ! Get the mapping weight (i.e. The fraction of the NATIVE ! GRID BOX that lies w/in the GEOS-CHEM GRID BOX) CALL GET_MAP_WT( xedge_w, xedge_e, xedgeC_w, xedgeC_e, & yedge_s, yedge_n, yedgeC_s, yedgeC_n, & mapWt ) ... etc ...
These modifications will skip calling the subroutine GET_MAP_WT unless:
- The latitude of each Olson land map box is within a specified number of degrees (LATTHRESH) of the GEOS-Chem grid box latitude, and
- The longitude of each Olson land map box is within a specified number of degrees (LONTHRESH) of the GEOS-Chem grid box longitude.
By restricting when subroutine GET_MAP_WT is called in this way, we reduce a substantial amount of overhead. Profiling results with the Tuning and Analysis Utilities (TAU) have shown that this can result in a 100X speedup in execution time for COMPUTE_OLSON_LANDMAP, as shown by this figure.
The figure shows wall clock time (in seconds) for each subroutine on the main thread (CPU 0) of a 1-day GEOS-Chem simulation with GEOS-FP meteorology at 4° x 5° horizontal resolution. Before the fix was implemented, routine COMPUTE_OLSON_LANDMAP took about 11 seconds on the main CPU. After the fix, the time spent in the routine was reduced to 0.163 seconds, thus eliminating a major bottleneck.
--Bob Y. 17:07, 30 May 2014 (EDT)
Outstanding issues that have not yet been resolved
Olson land map values are read incorrectly when using GCHP compiled with GNU Fortran
We have discovered an issue with the Olson land map that occurs when compiling GCHP with the GNU Fortran compiler. The Olson land map values are apparently being returned from disk (via ExtData) as all zeroes. A full description of the issue may be found here: https://github.com/geoschem/gchp/issues/15
--Bob Yantosca (talk) 20:06, 2 January 2019 (UTC)