登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
[color=rgba(0, 0, 0, 0.9)]Step1: 下载相关插值权重文件; [color=rgba(0, 0, 0, 0.9)]单个文件下载:
svn export https://svn-ccsm-inputdata.cgd.u ... _bilin_da_100716.nc注意文件包含了:插值函数PopLatLon popLatLonV 需要的参数信息[color=rgba(0, 0, 0, 0.9)]批量文件下载(通常不需要):
svn export https://svn-ccsm-inputdata.cgd.ucar.edu/trunk/ncl
[color=rgba(0, 0, 0, 0.9)]Step2: 将权重文件加载到环境变量;
[color=rgba(0, 0, 0, 0.9)]首先,将文件保存在个人文件目录下如:weights
[color=rgba(0, 0, 0, 0.9)]然后,加载到环境变量:
[color=rgba(0, 0, 0, 0.9)]vi ~/.bashrc
export NCL_POP_REMAP=/mnt/c/Users/Heng/Desktop/CESM2/NCl/pop2/weights/[color=rgba(0, 0, 0, 0.9)]source ~/.bashrc [color=rgba(0, 0, 0, 0.9)]Step3: NCL脚本读取POP模式数据并转换到经纬网格,注意加载popRemap.ncl文件。
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/popRemap.ncl"
list_of_files=systemfunc("ls /mnt/h/CESM_data_output/G1850ECO-f19_g16/*.nc")file_list=addfiles(list_of_files, "r") ListSetType(file_list, "cat") IRON_FLUX=file_list[:]->IRON_FLUX
IRON_FLUX:=PopLatLon(IRON_FLUX, "gx1v6", "1x1d", "bilin", "da", "100716")
IRON_FLUX_avg:=IRON_FLUX(0:4,:,:)do n=0,4,1i=0IRON_FLUX_avg(i,:,:)=dim_avg_n_Wrap(IRON_FLUX(i:i+4,:,:), 0)i=i+12end doIRON_FLUX_avg:=dim_avg_n_Wrap(IRON_FLUX_avg, 0)
outfile_name="/mnt/c/Users/Heng/Desktop/CESM2/NCl/pop2/dust_carbon/G1850ECO-f19_g19-photo_TOT_zint_100m.nc"system("rm -rf " + outfile_name)outfile=addfile(outfile_name, "c")outfile->IRON_FLUX=IRON_FLUX_avg
printVarSummary(IRON_FLUX_avg)
|