NCL读取POP海洋模式数据并插值到经纬网格

@Xiaodie

[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)