最原始观测值与模拟值对比计算与绘图

最新推荐文章于 2023-12-29 16:05:59 发布

阿水水水水水

最新推荐文章于 2023-12-29 16:05:59 发布

阅读量1.1k

点赞数 2

分类专栏：脚本积累文章标签：脚本语言

本文链接：https://blog.youkuaiyun.com/qq_40256777/article/details/122155303

版权

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
file_obs="isop_obs7.csv"
file_sim="oldisopsim.csv"
f_h = asciiread("sitelist.txt",-1,"string")

lines_obs=asciiread(file_obs,-1,"string")
lines_obs@_Fillvalue="-999"
lines_sim=asciiread(file_sim,-1,"string")
delim="," 
R=new(10,float)
MAE=new(10,float)
RMSE=new(10,float)
do i=1,10
time=tofloat(str_get_field(lines_obs,1,delim))
tmp1=str_get_field(lines_obs,i+1,delim)
tmp2=str_get_field(lines_sim,i+1,delim)
allvalue=tofloat(str_get_field(lines_sim,i+1,delim))  ;把所有值都加上
tmp1@_Fillvalue&#