ncl数据处理--变量单位换算

chongzika

    相关数据： slat = 95
      lat = 96
      lon = 144
      slon = 144
      lev = 30
      nhym = 30
      nhyi = 31
      ilev = 31
      time = 1  // unlimited
      nb2 = 2
   variables:
      double slat ( slat )
         standard_name :        latitude
         long_name :        staggered latitude
         units :        degrees_north
         axis :        Y

      double lat ( lat )
         standard_name :        latitude
         long_name :        latitude
         units :        degrees_north
         axis :        Y

      double lon ( lon )
         standard_name :        longitude
         long_name :        longitude
         units :        degrees_east
         axis :        X

float DSTSFMBL ( time, lat, lon )
         long_name :        Mobilization flux at surface
         units :        kg/m2/s
         cell_methods :        time: mean

lxlxllx89

mbl       = dstsfmbl*1000*3600*24*365
这句话运算之后 会丢失坐标信息 用copy_VarCoords 来copy一下就好饿了

longlivehj

变量属性运算过程中丢了吧！

chongzika

longlivehj 发表于 2015-5-2 09:49
变量属性运算过程中丢了吧！

怎么再添加上去呢？我用的cdo的ensmean，或者说有没有办法避免这个问题呢

longlivehj

chongzika 发表于 2015-5-2 10:13
怎么再添加上去呢？我用的cdo的ensmean，或者说有没有办法避免这个问题呢

先找到在哪儿丢的呗！

chongzika

longlivehj 发表于 2015-5-2 10:56
先找到在哪儿丢的呗！

这个要怎么找啊？我就用了cdo一个口令ensmean，10年*12个月的数据全部做平均，没有其他的操作了

longlivehj

chongzika 发表于 2015-5-2 10:58
这个要怎么找啊？我就用了cdo一个口令ensmean，10年*12个月的数据全部做平均，没有其他的操作了

不是运行area.ncl出的问题么？检查脚本里的变量啊！

chongzika

longlivehj 发表于 2015-5-2 11:02
不是运行area.ncl出的问题么？检查脚本里的变量啊！

理解错了，我以为是cdo出的问题！没感觉啥具体的错误啊，这是具体的脚本：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"
; =================================================;
begin
  f = addfile("2001-2010.nc","r")
  dstsfmbl  = f->DSTSFMBL(0,:,:)   ; read in example data [2D only here]
  mbl       = dstsfmbl*1000*3600*24*365
wks = gsn_open_wks("eps","mbl3")               ; open a ps file
gsn_define_colormap(wks,"MPL_YlOrBr")  ; choose colormap
  res  = True                                      ; plot mods desired
  res@tiMainString         = "mbl"     ; title
  res@cnFillOn             = True                  ; turn on color fill
  res@mpLandFillColor        = "white"
  res@gsnRightString        = "t/km:S:2:N:"
; =================================================;  
  res@mpMinLatF              =  20
  res@mpMaxLatF              =  60
  res@mpMinLonF              = 70
  res@mpMaxLonF              = 130
; =================================================;
  res@gsnAddCyclic  = False      ;
  res@mpDataSetName         = "Earth..4"   ; This new database contains
                                           ; divisions for other countries.
  res@mpDataBaseVersion     = "MediumRes"  ; Medium resolution database
  res@mpOutlineOn=True ; Turn on map outlines
  res@mpOutlineSpecifiers=(/"China","Xinjiang Uygur"/) ;China:states
; =================================================;
  res@pmTickMarkDisplayMode = "Always"            
plot = gsn_csm_contour_map_ce(wks, mbl, res)        ; create plot
end

chongzika

lxlxllx89 发表于 2015-5-2 17:41
mbl       = dstsfmbl*1000*3600*24*365
这句话运算之后 会丢失坐标信息 用copy_VarCoords 来copy一下就好 ...

好的，多谢赐教