用read_wrf_nc_change.f修改wrf的地形出现的问题

小英子

大马哈 发表于 2014-6-27 10:56
我装的3.6。3的版本

我重装了3.6.3的还是不行，我用的是pgi。
提示还是为：
[wxy@wxy em_real]$ pgf90 read_wrf_nc_change.f -L/home/wxy/modle/netcdf/lib -lnetcdf -lm \-I/home/wxy/modle/netcdf/include -Mfree -o read_wrf_nc_change
read_wrf_nc_change.f:
PGF90-S-0124-Relational expression mixes character with numeric data (read_wrf_nc_change.f: 1635)
PGF90-S-0076-Subscripts specified for non-array variable data_real (read_wrf_nc_change.f: 1636)
PGF90-S-0076-Subscripts specified for non-array variable data_real (read_wrf_nc_change.f: 1636)
PGF90-S-0076-Subscripts specified for non-array variable data_real (read_wrf_nc_change.f: 1636)
PGF90-S-0076-Subscripts specified for non-array variable data_real (read_wrf_nc_change.f: 1636)
PGF90-S-0076-Subscripts specified for non-array variable data_real (read_wrf_nc_change.f: 1636)
PGF90-S-0076-Subscripts specified for non-array variable data_real (read_wrf_nc_change.f: 1636)
PGF90-S-0072-Assignment operation illegal to external procedure data_real (read_wrf_nc_change.f: 1636)
  0 inform,   0 warnings,   8 severes, 0 fatal for user_code
PGF90-S-0302-Unmatched ENDSUBROUTINE statement (read_wrf_nc_change.f: 1685)
  0 inform,   0 warnings,   1 severes, 0 fatal for MAIN
[wxy@wxy em_real]$ ./read_wrf_nc_change.f wrfinput_d01 -EditData HGT
./read_wrf_nc_change.f:行1: Special: 未找到命令
information: ERROR: cannot open `information' (No such file or directory)
on:          ERROR: cannot open `on' (No such file or directory)
the:         ERROR: cannot open `the' (No such file or directory)
screen:      ERROR: cannot open `screen' (No such file or directory)
./read_wrf_nc_change.f:行3: Can: 未找到命令
./read_wrf_nc_change.f:行4: 未预期的符号 `(' 附近有语法错误
./read_wrf_nc_change.f:行4: `!  Can read double precision file (like WRF-Var)'

请问可以给我传一个你修改好的read_wrf_nc_change.f文件么？

大马哈

小英子 发表于 2014-6-27 23:03
我重装了3.6.3的还是不行，我用的是pgi。
提示还是为：
[wxy@wxy em_real]$ pgf90 read_wrf_nc_change. ...

不好意思 白天没有注意  我用的是官网的文件。我觉得改子程序那部分没有错  我觉得楼主可以用ncview看看数据格式是什么  我是修改的土地利用分类 按理说是整形 但是用data——int是不行的 要用real才行 是不是编译器的什么问题啊 我觉得楼主是不是换个编译器试试  
=========================
说完我又去看了看  ncview上面显示的是整形  但是我就是整形没有改出来 实形改出来了   昨天画图看着是改好了 但是没有做后续的步骤 明天运行一下 做一下对比  祝楼主找出原因  互相学习

小英子

问题解决。
是subroutine USER_CODE 修改的时候错误导致。

风浪

值得学习，谢谢楼至分享！！