delta_SLP= SLP-SLP_f copy_VarMeta(u_850,delta_u) copy_VarMeta(v_850,delta_v) copy_VarMeta(SLP,delta_SLP) delta_u@units="m/s" delta_v@units="m/s" delta_SLP@units="hPa" ;画图 wks_type = "png" wks_type@wkHeight = 2400 wks_type@wkWidth = 2400 wks = gsn_open_wks(wks_type...
CLDTOT_f= Fire_T*100 copy_VarMeta(Fire_H,CLDHGH_f) copy_VarMeta(Fire_M,CLDMED_f) copy_VarMeta(Fire_L,CLDLOW_f) copy_VarMeta(Fire_T,CLDTOT_f) CLDHGH_f@units="percent" CLDMED_f@units="percent" CLDLOW_f@units="percent" CLDTOT_f@units="percent" ;delta delta_H = CLDHGH-CLDHG...
pm_grid&lon = lon_grid 还有曲面网格属性的赋值,一般直接用copy_VarMeta函数直接从模式结果的变量里复制属性就可以了。但有时候也会进行一些改动,需要手动赋值。与标准网格数据不同的是,经纬度是二维数组,用@进行赋值 pm_cmaq!0 = "time" pm_cmaq!1 = "lon2d" pm_cmaq@lon2d = lon_cmaq pm_cmaq!2 =...
yyfac=yyyymm_to_yyyyfrac(yyyymm,0.0);sst1=dtrend_msg_n(yyfac,sst,True,False,0)sst1=dtrend_quadratic_msg_n(sst,True,False,0)copy_VarMeta(sst,sst1)clStrt=198101clLast=201012iClmStrt=ind(yyyymm.eq.clStrt)iClmLast=ind(yyyymm.eq.clLast)xClm=clmMonTLL(sst1(iClmStrt:iClmLast,:,:))...
"r") air2m = short2flt(f_air2m->air(rec_s:rec_e,0,{-90:90},:)) ; T at 2m ;;;DJF 平均& 异常;;; ;JFM季节平均,实际是12/1/2月三个月平均,因为从1979年12月开始截取 copy_VarMeta(sst(0,:,:),sst_DJF(0,:,:)) sst_DJF!0 = "year" sst_DJF&year=year sst_ano = dim...
·计算中可能造成元数据的丢失,可以使用copy_VarMeta、copy_VarCoords和copy_VarAtts函数复制元数据,或在计算之前以y=x形式预定义目标变量。 ·数组运算中<和>具有裁剪功能,例如,sst=sst>sice将数组sst中小于sice的元素替换成sice,这样结合@_FillValue属性就可以将小于等于阈值sice的值定义为无效数据。 20-20 10-...
; 去除海温趋势及年循环 yyfac = yyyymm_to_yyyyfrac(yyyymm, 0.0) ; sst1 = dtrend_msg_n(yyfac, sst, True, False,0) sst1 = dtrend_quadratic_msg_n(sst, True, False,0) copy_VarMeta(sst, sst1) clStrt = 198101 clLast = 201012 iClmStrt = ind(yyyymm.eq.clStrt) iClmLast = ind...
copy_VarMeta(v3d_template, ws_agl_locs) ws_agl_locs&height = hgt_out ws_agl_locs&reftime = reftime ws_agl_locs&leadtime = leadtime ws_agl_locs@description = "wind speed" ws_agl_locs@units = "m/s"if vars@speed filevardef(fout, "SPEED",typeof(ws_agl_locs),dims_3d) ...
u = where(u.eq."nan", u@_FillValue, false_value) ; ---change the time type to integer not int64 time0 = f->valid_time time = toint(time0) copy_VarMeta(time0, time) ; printVarSummary(time) timeARR = cd_calendar(time, 3) Computing...
使用dim_avg会失去T 的元数据: Tzon = dim_avg (T) 使用dim_avg_Wrap会保留元数据: Tzon = dim_avg_Wrap (T) 可用的wrappers一览表: /Document/Functions/Contributed/ copy_VarAtts, copy_VarCoords, copy_VarMeta, copyatt等函数由用户提供,用于显式复制坐 标变量或属性。 5.6.2 类型转换 能够将一...