- 积分
- 47
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-11-5
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2014-11-30 21:27:39
|
显示全部楼层
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
in = addfile("atmos.nc","r")
ts = in->TS(0,:,:)
oro = in->ORO(0,:,:)
;************************************************
; use mask function to mask out land and then ocean data
; ocean=0,land=1,sea_ice=2
;************************************************
land_only = ts ; trick to keep cv's and atts
ocean_only = ts
land_only = mask(ts,oro,1) ; Return the values of ts where oro=1
ocean_only = mask(ts,oro,0) ; Return the values of ts where oro=0
|
|