请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 5351|回复: 1

[经验总结] 利用shell+ncl脚本循环处理数据并输出

[复制链接]

新浪微博达人勋

发表于 2022-5-27 10:53:51 | 显示全部楼层 |阅读模式

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

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

x
NCL专辑 | 提高NCL代码的运行效率的各种方法

方法参考的是这一篇的内容,并成功实现

我需要对我所有的数据进行纬向平均,得到的数据应该为三维time,lat,lev,
如果直接用cdo中的zonmean,lon这一维度还存在,所以在网上找到了shell和ncl脚本配合使用输出文件

shell脚本:(在ubuntu下,使用 /. 调用shell脚本,我的脚本是zon.sh)
  1. export datadirectory='/azong/data/'         #文件路径
  2. declare -i numOFfile=20                       #文件总数,整形
  3. #开始循环
  4. for ((i=0;i<$numOFfile;i=i+1))
  5. do
  6.    export fileNO=${i}      #循环一次更新声明一次文件编号,由NCL中的getenv函数读取
  7.    ncl zon.ncl    #你要循环运行的脚本
  8. done
复制代码
ncl脚本:(文件名称为zon.ncl)
  1. begin
  2.   filedirectory = getenv("datadirectory")  ;调用文件目录
  3.   fs = systemfunc("ls "+filedirectory+"*")   ;列出文件
  4.   i = stringtoint(getenv("fileNO"))    ;读文件号,转成整型

  5.     ;读文件内容
  6.   f := addfile(fs(i),"r")   
  7.   u = f->uas      ;赋值
  8.   uas = dim_avg_n_Wrap(u,2)   ;计算平均

  9. fout=addfile("zon"+fs(i), "c") ;表示给输出数据加前缀zon
  10. fout->uas=uas                  

  11. end
复制代码



密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2023-3-2 18:14:03 | 显示全部楼层
多谢分享 学习一下{:eb502:}{:eb502:}
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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