- 积分
- 319
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-4-13
- 最后登录
- 1970-1-1
|
发表于 2018-12-11 20:17:53
|
显示全部楼层
请教天目大神:我用30米精度的GIS土地利用数据替换USGS数据,按照国内外论坛里的教程,在GIS里重投影、重采样、转ASCII等过程之后,又用write_geogrid.c、write_geogrid.o以及别人写的.f90文件生成了00001-05323.00001-05077这样的文件,然后按照wrf user form里面的教程,写index和修改GEOGRID.TBL, 在运行geogrid.exe时花的时间很长,约30分钟,生成的log文件很大,查看log文件,里面出现大量的“WARNING: In source tile /home/bloom/DATA/geog/wuhan2010/00001-05323.00001-05077, point (1, 1) has an invalid category of 241”这样的信息。我没管它,继续到real.exe时出现报错“error in the grid%tsk
i,j= 1 1
grid%landmask= 0.0000000E+00
grid%tsk, grid%sst, grid%tmn= 0.0000000E+00 0.0000000E+00 0.0000000E+00
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 2920
grid%tsk unreasonable”。
想问一下天目大神,我这是哪个环节出了问题呢?是GIS转换的问题,还是index或者GEOGRID.TBL没有写对?拜托拜托,卡死在这里了。
我的ascii头文件如下:
ncols 5323
nrows 5077
xllcorner 113.55347408316
yllcorner 29.859405496465
cellsize 0.00031751915813163
NODATA_value -9999
我的index文件如下:
type=categorical
category_min=1
category_max=28
projection=regular_ll
dx=0.00031751915813163
dy=0.00031751915813163
known_x=1.0
known_y=1.0
known_lat=29.859405496465
known_lon=113.55347408316
wordsize=1
tile_x=5323
tile_y=5077
tile_z=1
units="category"
description="24-category USGS landuse-from gis"
iswater=16
islake=28
isurban=1
isice=24
isoilwater=14
mminlu="USGS"
我的GEOGRID.TBL.ARW文件如下(只修改了LANDUSEF部分):
===============================
name=LANDUSEF
priority=2
dest_type=categorical
z_dim_name=land_cat
landmask_water = default:16
interp_option = 1s:nearest_neighbor
rel_path= 1s:wuhan2010/
===============================
name=LANDUSEF
priority=1
dest_type=categorical
z_dim_name=land_cat
dominant = LU_INDEX
landmask_water = modis_lakes:17,21 # Calculate a landmask from this field
landmask_water = nlcd2006_9s:17 # Calculate a landmask from this field
landmask_water = nlcd2006_30s:17 # Calculate a landmask from this field
landmask_water = nlcd2011_9s:17 # Calculate a landmask from this field
landmask_water = nlcd2006:17 # Calculate a landmask from this field
landmask_water = ssib_10m:16 # Calculate a landmask from this field
landmask_water = ssib_5m:16 # Calculate a landmask from this field
landmask_water = modis_15s:17 # Calculate a landmask from this field
landmask_water = modis_30s:17 # Calculate a landmask from this field
landmask_water = modis_30s_lake:17,21 # Calculate a landmask from this field
landmask_water = usgs_30s:16 # Calculate a landmask from this field
landmask_water = usgs_lakes:16,28 # Calculate a landmask from this field
landmask_water = usgs_2m:16 # Calculate a landmask from this field
landmask_water = usgs_5m:16 # Calculate a landmask from this field
landmask_water = usgs_10m:16 # Calculate a landmask from this field
landmask_water = default:17,21 # Calculate a landmask from this field
interp_option = modis_lakes:nearest_neighbor
interp_option = nlcd2006_9s:average_gcell(0.0)
interp_option = nlcd2006_30s:average_gcell(0.0)
interp_option = nlcd2011_9s:average_gcell(0.0)
interp_option = nlcd2006:nearest_neighbor
interp_option = ssib_10m:four_pt
interp_option = ssib_5m:four_pt
interp_option = modis_15s:nearest_neighbor
interp_option = modis_30s:nearest_neighbor
interp_option = modis_30s_lake:nearest_neighbor
interp_option = usgs_30s:nearest_neighbor
interp_option = usgs_lakes:nearest_neighbor
interp_option = usgs_2m:four_pt
interp_option = usgs_5m:four_pt
interp_option = usgs_10m:four_pt
interp_option = default:nearest_neighbor
rel_path = modis_lakes:modis_landuse_21class_30s/
rel_path = nlcd2006_9s:nlcd2006_ll_9s/
rel_path = nlcd2006_30s:nlcd2006_ll_30s/
rel_path = nlcd2011_9s:nlcd2011_ll_9s/
rel_path = nlcd2006:nlcd2006_ll_30s/
rel_path = ssib_10m:ssib_landuse_10m/
rel_path = ssib_5m:ssib_landuse_5m/
rel_path = modis_15s:modis_landuse_20class_15s/
rel_path = modis_30s:modis_landuse_20class_30s/
rel_path = modis_30s_lake:modis_landuse_20class_30s_with_lakes/
rel_path = usgs_30s:landuse_30s/
rel_path = usgs_lakes:landuse_30s_with_lakes/
rel_path = usgs_2m:landuse_2m/
rel_path = usgs_5m:landuse_5m/
rel_path = usgs_10m:landuse_10m/
rel_path = default:modis_landuse_20class_30s_with_lakes/
=============================== |
|