爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 7563|回复: 4

POP数据处理!

[复制链接]
发表于 2014-11-5 09:38:25 | 显示全部楼层 |阅读模式

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

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

x
程序如下:
;************************************************
; pop2lat_1.ncl
;************************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"  
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/popRemap.ncl"  
;************************************************
begin
  in  = addfile("case01.pop.h.nmonth.1908-01.nc","r")
  ice = in->SST
  printVarSummary(ice)
  print(in)
;************************************************
; convert pop to a 1x1 degree grid
;************************************************
  remap = PopLatLon(ice,"gx1v6","1x1d","bilin","da","696451")
  printVarSummary(remap)
;************************************************
; create plots
;************************************************
  wks = gsn_open_wks("x11","pop2lat")             ; open ps file
  gsn_define_colormap(wks,"WhBlGrYeRe")         ; choose colormap
;
; This will not be necessary in V6.1.0 and later. Named colors can
; be used without having to first add them to the color map.
;
  i = NhlNewColor(wks,0.7,0.7,0.7)              ; add gray to colormap  
  res                      = True
  res@gsnPolar             = "NH"               ; choose hemisphere  
  res@cnFillOn             = True               ; turn on color
  res@cnLinesOn            = False              ; no contour lines  
  res@mpLandFillColor      = "grey"             ; color of land
  res@mpMinLatF            = 50                 ; adjust map limit
  res@cnLevelSelectionMode = "ManualLevels"     ; set manual contour levels
  res@cnMinLevelValF       =  0.                ; set min contour level
  res@cnMaxLevelValF       =  5.                ; set max contour level
  res@cnLevelSpacingF      =  0.5               ; set contour spacing  
  res@gsnSpreadColors     = True                ; use full colormap
  res@gsnSpreadColorStart = 2
  res@gsnSpreadColorEnd   = 97  
  res@tiMainString        = "gx1v3 to 1x1d"  
; because this is ice data, which has a gap in the tropics, we need to
; explicitly pass the range of the data to plot.  Since we are coloring
; the vectors, this range should also match the MinLatF above, since the
; range for the colors is chosen over the full data passed, and not the
; map limits.
  plot = gsn_csm_contour_map_polar(wks,remap(0,{50:90},:),res)
end

报错如下:
(0)         
ls: cannot access map_*: No such file or directory
(0)        create_rmpPopFileName: file not recognized/found: map_gx1v6_to_1x1d_bilin_da_696451.nc
(0)        create_rmpPopFileName: dir searched: ./

求大神指教啊!
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-11-5 09:45:37 | 显示全部楼层
求大神指教啊
密码修改失败请联系微信:mofangbao
发表于 2014-11-5 09:52:04 | 显示全部楼层
这我也想知道
密码修改失败请联系微信:mofangbao
 楼主| 发表于 2014-11-10 10:25:42 | 显示全部楼层
这个问题已经解决!传一个脚本给大家一起分析!

ESMF_regrid_24.ncl

5.56 KB, 下载次数: 56, 下载积分: 金钱 -5

密码修改失败请联系微信:mofangbao
发表于 2016-12-13 09:22:19 | 显示全部楼层
我也提示是这个错误 是哪里有问题呀 应该怎么改
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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