爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 941|回复: 0

[经验总结] NCL读取POP海洋模式数据并插值到经纬网格

[复制链接]

新浪微博达人勋

发表于 2023-10-23 10:18:21 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

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)
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表