爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
楼主: tbag

应用高精度地形数据到WRF中(ASTER,SRTM)

  [复制链接]
发表于 2018-9-14 09:54:16 | 显示全部楼层
tbag 发表于 2018-9-14 00:28
如果要尝试的话,我觉得不用改变第一个文件,可以将超过90000变成100000-100000的格式

ps: 我计算过, ...

确实,如果直接将5位数改为6位数出问题了,wrf找不到需要的ASTER地形文件,结果就是地形全部都没有。看来要试验的话估计只能按你的来了。
密码修改失败请联系微信:mofangbao
发表于 2018-9-14 09:54:28 | 显示全部楼层
tbag 发表于 2018-9-14 00:28
如果要尝试的话,我觉得不用改变第一个文件,可以将超过90000变成100000-100000的格式

ps: 我计算过, ...

确实,如果直接将5位数改为6位数出问题了,wrf找不到需要的ASTER地形文件,结果就是地形全部都没有。看来要试验的话估计只能按你的来了。
密码修改失败请联系微信:mofangbao
发表于 2018-9-14 10:29:31 | 显示全部楼层
tbag 发表于 2018-9-10 23:45
下垫面其实也没那么难,只要你知道改什么

3600 行列会有问题吗?我记得好像理论上最大是9999?但是一般好像都是1200一个文件来着
密码修改失败请联系微信:mofangbao
发表于 2018-9-14 10:29:48 | 显示全部楼层
tbag 发表于 2018-9-10 23:45
下垫面其实也没那么难,只要你知道改什么

3600 行列会有问题吗?我记得好像理论上最大是9999?但是一般好像都是1200一个文件来着
密码修改失败请联系微信:mofangbao
发表于 2018-10-9 17:02:30 | 显示全部楼层
喵^_^胖天使 发表于 2018-9-14 09:54
确实,如果直接将5位数改为6位数出问题了,wrf找不到需要的ASTER地形文件,结果就是地形全部都没有。看来 ...

请问试验成功了吗
密码修改失败请联系微信:mofangbao
发表于 2018-10-9 17:02:45 | 显示全部楼层
喵^_^胖天使 发表于 2018-9-14 09:54
确实,如果直接将5位数改为6位数出问题了,wrf找不到需要的ASTER地形文件,结果就是地形全部都没有。看来 ...

请问试验成功了吗
密码修改失败请联系微信:mofangbao
发表于 2018-10-10 11:02:52 | 显示全部楼层
本帖最后由 zjycy070206 于 2018-10-10 20:07 编辑
tbag 发表于 2017-1-13 00:09
再提一点,在hdr描述文件中,经纬度分辨是25和120,但在index文件中,经纬度一定要是最最左下角的坐标,也就 ...

楼主,请问aster的地形数据要不要转换投影方式  转wrf用的事lambert投影    还有由于没有管理员权限  没法在geog文件里添加topo_ASTER文件   在GEOGRID.TBL 里应该怎样改

密码修改失败请联系微信:mofangbao
 楼主| 发表于 2018-10-11 00:32:56 | 显示全部楼层
zjycy070206 发表于 2018-10-10 11:02
楼主,请问aster的地形数据要不要转换投影方式  转wrf用的事lambert投影    还有由于没有管理员权限  没 ...

投影的话按道理来说跟使用什么地形数据没什么关系

GEOG的问题我也遇到过,我解决的方式是自己下载了完整的wps_geog文件,然后自己Link到下载后的地方进行更改
密码修改失败请联系微信:mofangbao
发表于 2018-10-11 08:04:35 | 显示全部楼层
tbag 发表于 2018-9-14 00:17
这个问题也困扰着我说实话
我没有测试过这个,因为我感觉这个数太大会出问题,我做的是一个大岛的模拟, ...

其实可以,秘密就在index文件里
type = continuous
signed = yes
projection = regular_ll
dx = 0.00027777778
dy = 0.00027777778
known_x = 1.0
known_y = 1.0
known_lat = 24.000
known_lon = 83.000
wordsize = 2
endian = little
filename_digits=6
tile_x = 3601
tile_y = 3601
tile_z = 1
row_order = top_bottom
missing_value = 32768
units = "meters MSL"
description = "ASTER 1-sec Topography Height"
注意到了么?“filename_digits=6”
密码修改失败请联系微信:mofangbao
回复 支持 2 反对 0

使用道具 举报

发表于 2018-10-11 08:06:34 | 显示全部楼层
其实参考一下WPS_GEOG\urbfrac_nlcd2011的index就明白了,这个也是6位数的
type=continuous
projection=albers_nad83
dx=30.0
dy=30.0
known_x=1.0
known_y=104424
known_lat = 48.707394
known_lon = -130.232828
truelat1=29.5
truelat2=45.5
stdlon=-96.0
wordsize=1
scale_factor=0.01
filename_digits=6
tile_x=1990
tile_y=2748
tile_z=1
units="fraction"
description="Urban fraction derived from 30 m NLCD 2011 (22 = 50%, 23 = 90%, 24 = 95%)"
注意到了么?“filename_digits=6”
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表