WRF修改土壤湿度，土壤湿度从外部输入怎么实现？

shrily521

根据天目神眉大大的http://bbs.06climate.com/forum.php?mod=viewthread&tid=4273来修改wrfinput_d0，修改里面的土壤湿度，
用vi命令打开read_wrf_nc_change.f90修改此程序中的subroutine USER_CODE这一块
例：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
想用wrfout文件中的土壤湿度（已经用grads转成dat数据）来修改这部分。不知道应该怎么写？
我把这块改成了
open(9999,file='/gpfshome/zhangsw/lishyi/Program/351/WRFV3/modify_initial_field/smois.2001.06.d01.6.dat',status='unknown',form='unformatted',access='direct')
read(9999,rec=1) data_real1
if(var=='SMOIS')then
data_real=data_real1
endif
运行./read_wrf_nc wrfinput_d01 -EditData SMOIS
出现
Changing variable: SMOIS for time 2001-06-01_00:00:00
Segmentation fault (core dumped)
报错~

求助~求指正~谢谢