- 积分
- 102
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-2-15
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2021-2-25 23:35:12
|
显示全部楼层
begin
filein=addfile("D:\CN05.1\2019\CN05.1_Pre_2019_month_025x025.nc","r")
t=filein->time
t@units="hours since 1961-01-01 00:00:00"
date=cd_calendar(t,0)
year=date(:,0)
month=date(:,1)
js=filein->pre
js@_FillValue = js@missing_value
jsa=new((/163,283/),float)
jsa=(js(3,:,:)+js(4,:,:)+js(5,:,:))
js1=jsa(69:109,161:201)
lon=filein->lon(161:201)
lon@units = "degrees_north"
lat=filein->lat(69:109)
lat@units = "degrees_east"
fill_areas=(/"China:states"/)
wks = gsn_open_wks("png","premam")
gsn_define_colormap(wks,"rainbow")
res=Trueres@mpOutlineSpecifiers = (/"China:states","Taiwan"/) ;China:states
res@mpDataBaseVersion="MediumRes"
res@mpDataSetName="Earth..4" ;包含中国省界
res@cnFillOn = True ; turn on color fill
res@cnLinesOn = False ; turn off the contour lines
res@cnLineLabelsOn = False ; turn the line labels off
;js1=jsa(69:109,161:201)
res@tiMainString="HUABEIpre_MAM"
plot=gsn_csm_contour_map(wks,js1,res)
end
这是脚本 |
|