关于WRF下垫面制作的问题

arch_bloom

请教论坛里的各位大神：
我用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/
===============================

arch_bloom

在论坛里看过天目、兰溪等大神的帖子，但还是没找出我的问题。实在痛苦。。。

APTX4869

很早以前改过，也出现同样警告，不过run成功了

arch_bloom

APTX4869 发表于 2018-12-12 09:19
很早以前改过，也出现同样警告，不过run成功了

我的除了有warnning, 在运行real.exe时还出现了FATAL CALLED，所以没办法运行WRF.EXE了！谢谢你的留言。

APTX4869

arch_bloom 发表于 2018-12-12 09:40
我的除了有warnning, 在运行real.exe时还出现了FATAL CALLED，所以没办法运行WRF.EXE了！谢谢你的留言。

你可以先检查这个错误是不是替换的问题？

arch_bloom

APTX4869 发表于 2018-12-12 15:09
你可以先检查这个错误是不是替换的问题？

麻烦你能否说的稍微详细一点，是GEOGRID.TBL.ARW里面的条目有问题吗？另外，你在GIS里重投影的时候，选的是哪一个坐标系里的WGS_84？我发现地理坐标系和投影坐标系里都有WGS_84，我选的是地理坐标系里的。

弱水三千

感觉是数据替换出的问题，（1,1）LU_INDEX咋是241，检查一下。另外想请教一下，GIS土地利用替换USGS时，需要重分类吗？

arch_bloom

弱水三千 发表于 2018-12-14 11:28
感觉是数据替换出的问题，（1,1）LU_INDEX咋是241，检查一下。另外想请教一下，GIS土地利用替换USGS时，需 ...

谢谢你的回复！我们研究室买的GIS数据的分类和USGS的不一样，所以需要按照USGS的分类方法进行重分类。请问一下我这个问题该怎么修改呢？再次感谢！

arch_bloom

arch_bloom 发表于 2018-12-14 13:10
谢谢你的回复！我们研究室买的GIS数据的分类和USGS的不一样，所以需要按照USGS的分类方法进行重分类。请 ...

重分类的方法在常鸣博士的一篇文章里有提到。https://www.guokr.com/blog/441223/

弱水三千

arch_bloom 发表于 2018-12-14 13:10
谢谢你的回复！我们研究室买的GIS数据的分类和USGS的不一样，所以需要按照USGS的分类方法进行重分类。请 ...

检查一下你的LU_INDEX值呢！