- 积分
- 24
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-10-17
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
请高人指点
我想输出类似于这样的nc文件
但是总是不对,我的code是这样的
; predefine coordinate variables (预先定义坐标变量)
nlat = dimsizes(demlat)
nlon = dimsizes(demlon)
dimlatSizes = (/nlat/) ;######改
dimlonSizes = (/nlat/)
dimphSizes = (/117,148/)
;;;;改成自己的网格设置*******(e_sn-1,e_we-1)*********************
dimlatUnlim = (/False/)
dimlonUnlim = (/False/)
dimphUnlim = (/False,False/)
; explicit dimension naming or getvardims can be used (明确维名称,或者可以用getvardims命令)
;filedimdef( ncf, dim_names, dim_sizes, dimUnlim )
dim_names = (/ "col","row"/)
filedimdef(difout,"lat",dimlatSizes,dimlatUnlim)
filedimdef(difout,"lon",dimlonSizes,dimlonUnlim)
filedimdef(difout,dim_names,dimphSizes,dimphUnlim)
;var_names = (/"difphase"/)
;var_types = (/"float"/)
var_nameslat1D = (/ "lat" /)
var_nameslon1D = (/ "lon"/)
var_names2D = (/ "difphase" /)
varvar_types1D = (/ "float" /)
varvar_types2D = (/ "float", "float" /)
filevardef( difout, var_nameslat1D, varvar_types1D, (/ "lat"/) ) ;input lat 1D
filevardef( difout, var_nameslon1D, varvar_types1D, (/ "lon"/) )
filevardef( difout, var_names2D, varvar_types2D, (/ "col", "row" /) )
;filevardef(fout,var_names,var_types,(/"lat","lon","phase"/))
lat@long_time = "latitude"
lat@units = "degrees_north"
lon@long_time = "longitude"
lon@units = "degrees_east"
difphase@long_time = "diff phase"
difphase@units = "arc degree"
;;do k=1,dim(0)-1
;demlat1=demlat
;demlon1=demlon
;h1 = phase2
; ; end do
filevarattdef(difout,"lat",demlat)
filevarattdef(difout,"lon",demlon)
filevarattdef(difout,"difphase",h1)
; ;predefine each variable’s attributes (预先定义每个变量的属性)
;setfileoption(fout,”SuppressDefineMode”,False) ; optional
; output values only (输出值)
fout->lat = (/demlat/)
fout->lon = (/demlon/)
fout->difphase = (/h1/)
先行谢过
|
-
目标文件
|