- 积分
- 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/
=============================== |
|
[复制链接]
窥视卡
雷达卡
显身卡