应用高精度地形数据到WRF中（ASTER,SRTM）

月夜猫妖

参考：http://www2.mmm.ucar.edu/wrf/users/docs/user_guide_v4/v4.0/users_guide_chap3.html#_Description_of_index
“34. FILENAME_DIGITS : An integer specifying the number of digits used in the names of data tiles. Possible values are 5 or 6. Default value is 5.”
明明白白得写了这个参数是用来调节文件名长度的，默认情况是5，也就是不写的时候是5位的

月夜猫妖

手动重命名？不存在的，送上一段小脚本

1. i=1;tmp1=0;tmp2=0; for x in ASTGTM2_N24E*; do tmp1=$[($i-1)*3601+1];tmp2=$[$i*3601];new=$(printf "%06d-" ${tmp1});new=$new$(printf "%06d.000001-003601" ${tmp2});gdal_translate -of ENVI -co INTERLEAVE=BSQ ${x} ${new};let i=i+1; done;

i=1;tmp1=0;tmp2=0; for x in ASTGTM2_N24E*; do tmp1=$[($i-1)*3601+1];tmp2=$[$i*3601];new=$(printf "%06d-" ${tmp1});new=$new$(printf "%06d.000001-003601" ${tmp2});gdal_translate -of ENVI -co INTERLEAVE=BSQ ${x} ${new};let i=i+1; done;

tbag

月夜猫妖 发表于 2018-10-11 08:10
参考：http://www2.mmm.ucar.edu/wrf/users/docs/user_guide_v4/v4.0/users_guide_chap3.html#_Description ...

可以的可以的~有了这个就拯救了我的科研一部分~多谢多谢~这个选项在index中确实没注意过

tbag

月夜猫妖 发表于 2018-10-11 08:24
手动重命名？不存在的，送上一段小脚本

这个命名脚本倒是没问题，但是前提的建立在已经添加了“空白”的tile文件把这个所有的范围都给覆盖到~

不过还是很方便~

月夜猫妖

tbag 发表于 2018-10-11 10:35
这个命名脚本倒是没问题，但是前提的建立在已经添加了“空白”的tile文件把这个所有的范围都给覆盖到~

...

技术不算太好只能一次一次一条一条地改，没法做到一次性全部搞定

喵^_^胖天使

感谢，原来还有这个参数选项。

zjycy070206

Timing for main: time 2017-06-06_00:01:05 on domain   2:    0.59290 elapsed seconds
d02 2017-06-06_00:01:05            2 points exceeded cfl=2 in domain d02 at time 2017-06-06_00:01:05 hours
d02 2017-06-06_00:01:05  MAX AT i,j,k:            6          10           3 vert_cfl,w,d(eta)=   2.758896       13.95187      1.2999952E-02
d02 2017-06-06_00:01:05            4 points exceeded cfl=2 in domain d02 at time 2017-06-06_00:01:05 hours
d02 2017-06-06_00:01:05  MAX AT i,j,k:            7          10           3 vert_cfl,w,d(eta)=   3.007503      -124.0469      1.2999952E-02
d02 2017-06-06_00:01:05            8 points exceeded cfl=2 in domain d02 at time 2017-06-06_00:01:05 hours
d02 2017-06-06_00:01:05  MAX AT i,j,k:            6          10           3 vert_cfl,w,d(eta)=   4.099665      -5.231970      1.2999952E-02
forrtl: severe (174): SIGSEGV, segmentation fault occurred

我使用地形数据后，转wrf总是报这种错误

tbag

zjycy070206 发表于 2018-10-18 20:54
Timing for main: time 2017-06-06_00:01:05 on domain   2:    0.59290 elapsed seconds
d02 2017-06-06_ ...

信息太少了昂

zjycy070206

tbag 发表于 2018-10-19 00:22
信息太少了昂

我可能发现问题在哪了   geogrid.exe的时候  topo_ASTER/index 中 unrecognized option in entry 16

是 row_order=top_bottom 这个代码识别不了 不知道该怎么办

zjycy070206

zjycy070206 发表于 2018-10-19 14:54
我可能发现问题在哪了   geogrid.exe的时候  topo_ASTER/index 中 unrecognized option in entry 16

...

结果问题是filename_digits=6这个代码识别不了     WPS要3.8以上才行