- 积分
- 218
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-1-4
- 最后登录
- 1970-1-1
|
发表于 2021-9-19 17:46:48
|
显示全部楼层
没那么复杂。。。。。大概思路是:用NCL读取WPS处理生成的met文件,然后更改土壤湿度最后写进去,你想更改什么时段和哪个层次都能实现。脚本如:cmd = systemfunc("ls ./met_em.*.nc")
Nfile = dimsizes(cmd)
;====读取========
nc_file = addfile(cmd(i),"r")
sm000010 =nc_file ->SM000010
;======更改=======
sm000010(:,loc(1,0):loc(1,1),loc(0,0):loc(0,1)) =sm000010(:,loc(1,0):loc(1,1),loc(0,0):loc(0,1))*0.5
;========写入==========
nc_file->SM000010 = sm000010
大概就是这么个思路 ,具体脚本你自己写
|
|