根据CFSR逐小时再分析数据计算某经纬位置日值序列

hongkang

＝＝＝＝＝
　问题
＝＝＝＝＝
利用ＣＤＯ、ｗｇｒｉｂ２等工具，根据CFSR逐小时再分析数据（ｇｒｉｂ２格式）计算某经纬位置日值序列

＝＝＝＝＝
单个文件/时间处理命令示例
＝＝＝＝＝
＃全球资料裁剪到中国区域
wgrib2 apcp.cdas1.201609.grb2 -set_grib_type c2 -small_grib 109:110 35:37 apcp.cdas1.201609.region.grb2
＃按时间转为ｎｃ格式
wgrib2  tmp2m.cdas1.201609.region.grb2 -netcdf tmp2m_201609.ncwgrib2  apcp.cdas1.201609.region.grb2  -d 1 -netcdf 201609.nc　　　　　＃第一个时次wgrib2  wnd10m.cdas1.201609.region.grb2 -d 1 -netcdf u_201609.nc　　　＃第一个时次wgrib2  wnd10m.cdas1.201609.region.grb2 -d 2 -netcdf v_201609.nc　　　＃第二个时次
＃不同时间文件合并为一个文件并转世界时
cdo mergetime tmp2m*nc  tmp2m.nccdo shifttime,+8hours tmp2m.nc tmp2m_BTC.nc
＃计算积温
cdo daymean tmp2m.nccdo addc,-273.15 tmp2m.nc tmp2m_Celsius.nc
＃日平均风速（u、v 在同一个文件）
cdo expr,'wnd_spd=sqrt(UGRD_10maboveground*UGRD_10maboveground+VGRD_10maboveground*VGRD_10maboveground)' wnd10m.cdas1.201609.region.nc wnd_spd10m.cdas1.201609.region.nc
＃日平均风速（ｕ、ｖ不在同一个文件）
cdo chname,uwnd,ws -sqrt -add -sqr -selname,uwnd ifile1 -sqr -selname,vwnd ifile2 ofile


＝＝＝＝＝＝
多文件批量处理
＝＝＝＝＝＝
计算积温（global_accumulated_temperature_calculate.py）
计算日值（DataFromCFSR.py）

李白

好东西！感谢楼主

爱听雨的青蛙

楼主厉害啊，学习学习！

浪子

谢谢楼主分享

速冻小黄鱼

非常感谢！都很有用！

ybhldl

感谢分享，赚钱赚钱赚钱，为了下载心仪的东东！

rabin_xu

逐小时的CFSR数据从哪里下载啊？

hongkang

rabin_xu 发表于 2017-12-21 21:34
逐小时的CFSR数据从哪里下载啊？

#! /bin/csh -f
# Usage: Download CFSR reanalysis data to specific dir
# Args: $1-password $2-username $3-year $4-mon $5-dir

set passwd = $1
set username = $2
set year = $3
set i = $4
set dir = $5

set num_chars = `echo "$passwd" |awk '{print length($0)}'`
if ($num_chars == 0) then
  echo "You need to set your password before you can continue"
  echo "  see the documentation in the script"
  exit
endif
@ num = 1
set newpass = ""
while ($num <= $num_chars)
  set c = `echo "$passwd" |cut -b{$num}-{$num}`
  if ("$c" == "&") then
    set c = "%26";
  else
    if ("$c" == "?") then
      set c = "%3F"
    else
      if ("$c" == "=") then
        set c = "%3D"
      endif
    endif
  endif
  set newpass = "$newpass$c"
  @ num ++
end
set passwd = "$newpass"
#
set cert_opt = ""
wget $cert_opt -O auth_status.rda.ucar.edu --save-cookies auth.rda.ucar.edu.$$ --post-data="email=$username&passwd=$passwd&action=login" https://rda.ucar.edu/cgi-bin/login

## current year
## set year = `date +%Y`
set opts = "-P $dir/$year -N"
set mons = (01 02 03 04 05 06 07 08 09 10 11 12)

wget $cert_opt $opts --load-cookies auth.rda.ucar.edu.$$ rda.ucar.edu/data/ds094.1/$year/tmax.cdas1.$year$mons[$i].grb2
wget $cert_opt $opts --load-cookies auth.rda.ucar.edu.$$ rda.ucar.edu/data/ds094.1/$year/tmp2m.cdas1.$year$mons[$i].grb2
wget $cert_opt $opts --load-cookies auth.rda.ucar.edu.$$ rda.ucar.edu/data/ds094.1/$year/tmin.cdas1.$year$mons[$i].grb2
wget $cert_opt $opts --load-cookies auth.rda.ucar.edu.$$ rda.ucar.edu/data/ds094.1/$year/rh2m.cdas1.$year$mons[$i].grb2
wget $cert_opt $opts --load-cookies auth.rda.ucar.edu.$$ rda.ucar.edu/data/ds094.1/$year/apcp.cdas1.$year$mons[$i].grb2
wget $cert_opt $opts --load-cookies auth.rda.ucar.edu.$$ rda.ucar.edu/data/ds094.1/$year/wnd10m.cdas1.$year$mons[$i].grb2
wget $cert_opt $opts --load-cookies auth.rda.ucar.edu.$$ rda.ucar.edu/data/ds094.1/$year/dswsfc.cdas1.$year$mons[$i].grb2

hongkang

rabin_xu 发表于 2017-12-21 21:34
逐小时的CFSR数据从哪里下载啊？

wget -O auth_status.rda.ucar.edu --save-cookies auth.rda.ucar.edu.$$ --post-data="email=***&passwd=***=login" https://rda.ucar.edu/cgi-bin/login

year='2019'
month='03'
wget $cert_opt $opts --load-cookies auth.rda.ucar.edu.$$ rda.ucar.edu/data/ds094.1/$year/tmax.cdas1.$year$month.grb2

932686787

请问这是用cdo吗？裁剪就一行命令？