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

小英子

使用天目介绍的方法：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=4273

用vi命令修改此程序中的subroutine USER_CODE这一
改为：
!-------------------------------------------------------------------------------------------
subroutine USER_CODE (data_real,data_dp_real,data_int,dim1,dim2,dim3,var)
if(var=='HGT') then
    data_real(160:200,125:155,1:1)=data_real(160:200,125:155,1:1)*0.5
    end if
    end subroutine USER_CODE                     


!------------------------------READ FIRST--------------------------------------------------  


运行出错：
[wxy@wxy em_real]$ pgf90 read_wrf_nc_change.f -L/home/wxy/wxy/netcdf/lib -lnetcdf -lm \-I/home/wxy/wxy/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: 1595)
PGF90-S-0076-Subscripts specified for non-array variable data_real (read_wrf_nc_change.f: 1596)
PGF90-S-0076-Subscripts specified for non-array variable data_real (read_wrf_nc_change.f: 1596)
PGF90-S-0076-Subscripts specified for non-array variable data_real (read_wrf_nc_change.f: 1596)
PGF90-S-0076-Subscripts specified for non-array variable data_real (read_wrf_nc_change.f: 1596)
PGF90-S-0076-Subscripts specified for non-array variable data_real (read_wrf_nc_change.f: 1596)
PGF90-S-0076-Subscripts specified for non-array variable data_real (read_wrf_nc_change.f: 1596)
PGF90-S-0072-Assignment operation illegal to external procedure data_real (read_wrf_nc_change.f: 1596)
  0 inform,   0 warnings,   8 severes, 0 fatal for user_code
PGF90-S-0302-Unmatched ENDSUBROUTINE statement (read_wrf_nc_change.f: 1638)
  0 inform,   0 warnings,   1 severes, 0 fatal for MAIN
[wxy@wxy em_real]$ chmod +x read_wrf_nc_change.f
[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)'
[wxy@wxy em_real]$

flying鹰

./read_wrf_nc_change.f wrfinput_d01 -EditData HGT错在这里了
改为./read_wrf_nc_change wrfinput_d01 -EditData HGT试试

小英子

flying鹰 发表于 2014-6-26 21:44
./read_wrf_nc_change.f wrfinput_d01 -EditData HGT错在这里了
改为./read_wrf_nc_change wrfinput_d01 - ...

[wxy@wxy em_real]$ ./read_wrf_nc_change wrfinput_d01 -EditData HGT
bash: ./read_wrf_nc_change: 没有那个文件或目录
这样好像也不行！
read_wrf_nc是从官网上下载的。
请问是不是
subroutine USER_CODE (data_real,data_dp_real,data_int,dim1,dim2,dim3,var)
if(var=='HGT') then
    data_real(160:200,125:155,1:1)=data_real(160:200,125:155,1:1)*0.5
    end if
    end subroutine USER_CODE
修改错了？

flying鹰

哦哦，是你的pgf90 read_wrf_nc_change.f -L/home/wxy/wxy/netcdf/lib -lnetcdf -lm \-I/home/wxy/wxy/netcdf/include -Mfree -o read_wrf_nc_change没成功
你先试试把data_real(160:200,125:155,1:1)=data_real(160:200,125:155,1:1)*0.5改为data_real=data_real*0.5然后运行
pgf90 read_wrf_nc_change.f -L/home/wxy/wxy/netcdf/lib -lnetcdf -lm  -I/home/wxy/wxy/netcdf/include -Mfree -o read_wrf_nc_change试试

小英子

flying鹰 发表于 2014-6-26 21:55
哦哦，是你的pgf90 read_wrf_nc_change.f -L/home/wxy/wxy/netcdf/lib -lnetcdf -lm \-I/home/wxy/wxy/netc ...

还是一样的提示错误，换为官网上没改动的subroutine USER_CODE，也是出现相同的错误！

flying鹰

小英子 发表于 2014-6-26 22:03
还是一样的提示错误，换为官网上没改动的subroutine USER_CODE，也是出现相同的错误！

你检查下netcdf安装路径～

flying鹰

试试pgf90 read_wrf_nc_change.f -L/home/wxy/model/netcdf/lib -lnetcdf -lm  -I/home/wxy/model/netcdf/include -Mfree -o read_wrf_nc_change

大马哈

楼主 您好 我刚刚也在运行这个程序  我按照fortran中的文件 写了如下 命令
================================================
ifort  read_wrf_nc.f90 -L/home/users/jshu/zyc/wrftest/netcdf/lib -lnetcdf  -lm  -I/home/users/jshu/zyc/wrftest/netcdf/include -Mfree  -o wrfread3 &>1.log
==============================================
其中有些 option不是很理解 但是想依葫芦画瓢 路径对就应该ok

但是报错中出现
===========================================================
/home/users/jshu/zyc/wrftest/netcdf-4.1.1/libncdap3/oc/http.c:44: undefined reference to `curl_easy_setopt'
21 /home/users/jshu/zyc/wrftest/netcdf-4.1.1/libncdap3/oc/http.c:79: undefined reference to `curl_easy_strerror'
22 /home/users/jshu/zyc/wrftest/netcdf-4.1.1/libncdap3/oc/http.c:49: undefined reference to `curl_easy_setopt'
23 /home/users/jshu/zyc/wrftest/netcdf-4.1.1/libncdap3/oc/http.c:54: undefined reference to `curl_easy_setopt'
24 /home/users/jshu/zyc/wrftest/netcdf-4.1.1/libncdap3/oc/http.c:59: undefined reference to `curl_easy_setopt'
25 /home/users/jshu/zyc/wrftest/netcdf-4.1.1/libncdap3/oc/http.c:63: undefined reference to `curl_easy_perform'
=============================================================
也就是说 编译f文件的时候 他在【/home/users/jshu/zyc/wrftest/netcdf-4.1.1/】文件夹下寻找 这是*。gz解压后的文件  
不知道楼主能否帮助解答    看到其他有类似的报错是netcdf未安装好 但是我按照全路径运行时可行


=======================
早上重新装了一个低级版本的可以了

小英子

大马哈 发表于 2014-6-26 23:05
楼主 您好 我刚刚也在运行这个程序  我按照fortran中的文件 写了如下 命令
============================= ...

你好，请问你是装的哪个版本的netcdf？

大马哈

小英子 发表于 2014-6-27 10:34
你好，请问你是装的哪个版本的netcdf？

我装的3.6。3的版本