load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
DataDir1 = "/data1/loucx/COPY/CNAQ/data/cctm/"
flist1 = systemfunc("ls "+DataDir1+"CONC.*_d03")
Sfile1 = addfiles(flist1, "r")
ListSetType (Sfile1, "join")
isop=Sfile1[:]->ISOP
o3=Sfile1[:]->O3
wrffile=addfile("/data1/loucx/COPY/WRF/DATA/WRF/wrfout_d03_2020-08-01_00:00:00","r")
lon2d=wrffile->XLONG(0,3:176,3:170)
lat2d=wrffile->XLAT(0,3:176,3:170)
f_h = asciiread("sitelist.csv",-1,"string")
delim=","
lon=tofloat(str_get_field(f_h,5,delim))
lat=tofloat(str_get_field(f_h,6,delim))
idnum=str_get_field(f_h,1,delim)
simISOP=rcm2points(lat2d, lon2d, isop, lat, lon, 0)
simISOP!0="day"
simISOP!1="hour"
simISOP!2="lay"
simISOP!3="site"
simO3=rcm2points(lat2d, lon2d, o3, lat, lon, 0)
simO3!0="day"
simO3!1="hour"
simO3!2="lay"
simO3!3="site"
record_simISOP=simISOP(site|:,day|:,hour|:,lay|:)
record_simO3=simO3(site|:,day|:,hour|:,lay|:)
;printVarSummary(record_simO3)
do i=1,53,1
write_table("提取O3-ISOP_"+idnum(i-1)+".csv","w",[/record_simO3(i-1,:,0:23,0)*48000/22.4,record_simISOP(i-1,:,0:23,0)*1000/],"%8.4f,%8.4f")
end do
更新提取数据(单位统一)
最新推荐文章于 2024-06-13 14:46:26 发布