- 积分
- 312
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-1-4
- 最后登录
- 1970-1-1
|
发表于 2022-7-31 15:24:50
|
显示全部楼层
大佬，你用过CMAQ里面的region-specific scaling of emissions using the CMAQv5.3 DESID module嘛，这个GRIDMASK files要怎么做自己的区域啊，我没看懂官方的教程。
附上官方教程：
GRIDMASK file processing for U.S. states and NOAA/NCDC climate regions:
• The Fortran code used to generate the GRIDMASK files utilized the SMOKE gridding surrogate files for surrogate categories 340 (“NLCD Land”) and 350 (“NLCD Water”). These files are named USA_340_NOFILL_12km.txt and USA_350_NOFILE_12km.txt and are available for download at ftp://newftp.epa.gov/air/emismod ... _v1.3_29nov2016.zip
• The general format of these gridding surrogate files is described at https://www.cmascenter.org/smoke ... ml/ch08s04s02.html.
In addition to the required data columns A-E of the non-header lines, the files used in the processing also had three additional columns (F -H) included for quality assurance purposes. Column F is an exclamation mark (“!”) to indicate that the following two data columns are optional and not used by SMOKE. Column G contains the total absolute value (area) of the surrogate feature (i.e. “NLCD Land” or “NLCD Water”) in the grid cell indicated by columns C and D, while column H contains the total absolute value (area) of the surrogate feature (i.e. “NLCD Land” or “NLCD Water”) in the state and county indicated by the FIPS code in column B. The ratio of column G to column H is the spatial surrogate ratio shown in column E which is used by SMOKE.
• The Fortran code utilized the information in columns B (integer FIPS code indicating state and county), C (integer indicating grid cell column), D (integer indicating grid cell row), and G (total “NLCD Land” or “NLCD Water” area in the grid cell indicated by columns C and D).
• Accounting for the offset between the 12 km Lambert Conformal grid for which the spatial surrogates were developed (grid origin x = -2736000 m, y = -2088000 m) and the 12US1 modeling grid for which the GRIDMASK file was developed (grid origin x = -2556000 m, y = -1728000 m), for a given grid cell the Fortran code computed the GRIDMASK value for a given state or climate region by dividing the sum of the “NLCD Land” and “NLCD Water” areas associated with that state or climate region through the FIPS code by the sum of the “NLCD Land” and “NLCD Water” areas associated with any FIPS code.
• For example, to compute the GRIDMASK value for “NC” for a given grid cell, the code summed all column G entries from the USA_340_NOFILL_12km.txt and USA_350_NOFILE_12km.txt files that were associated with the grid cell in question (defined by columns C and D, accounting for the offset between surrogate grid and 12US1 modeling grid) and had a FIPS code (column B) starting with “37”, and divided this value by the sum over all column G entries from the USA_340_NOFILL_12km.txt and USA_350_NOFILE_12km.txt files that were associated with the grid cell in question regardless of FIPS code.
• For the NOAA climate region GRIDMASK, states were assigned to climate regions based on the Karl and Koss (1984) definition of climate regions available at https://www.ncdc.noaa.gov/monito ... climate-regions.php (note that Washington, D.C. is not included in any of the climate regions on the website but was included with the “Northeast” region for the generation of this GRIDMASK file):
o CLIMATE_REGION_1: Northwest (OR, WA, ID)
o CLIMATE_REGION_2: West (CA, NV)
o CLIMATE_REGION_3: West North Central (MT, WY, ND, SD, NE)
o CLIMATE_REGION_4: Southwest (UT, AZ, NM, CO)
o CLIMATE_REGION_5: South (KS, OK, TX, LA, AR, MS)
o CLIMATE_REGION_6: Central (MO, IL, IN, KY, TN, OH, WV)
o CLIMATE_REGION_7: East North Central (MN, IA, WI, MI)
o CLIMATE_REGION_8: Northeast (MD, DE, NJ, PA, NY, CT, RI, MA, VT, NH, ME) + Washington, D.C.
o CLIMATE_REGION_9: Southeast (VA, NC, SC, GA, AL, GA)
Thomas R. Karl and Walter James Koss, 1984: "Regional and National Monthly, Seasonal, and Annual Temperature Weighted by Area, 1895-1983." Historical Climatology Series 4-3, National Climatic Data Center, Asheville, NC, 38 pp |
|
[复制链接]
窥视卡
雷达卡
显身卡
楼主
