日志
海平面气压场
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
wks_name="2019to112_windbab_surface"
wk = "pdf"
wks = gsn_open_wks(wk,wks_name)
gsn_define_colormap(wks,"matlab_jet")
filename_1 = "./fnl_data/fnl_20190101_00_00.grib2"
filename_2 = "./fnl_data/fnl_20190101_06_00.grib2"
filename_3 = "./fnl_data/fnl_20190101_12_00.grib2"
filename_4 = "./fnl_data/fnl_20190101_18_00.grib2"
filename_5 = "./fnl_data/fnl_20190102_00_00.grib2"
filename_6 = "./fnl_data/fnl_20190102_06_00.grib2"
filename_7 = "./fnl_data/fnl_20190102_12_00.grib2"
filename_8 = "./fnl_data/fnl_20190102_18_00.grib2"
filename_9 = "./fnl_data/fnl_20190103_00_00.grib2"
filename_10 = "./fnl_data/fnl_20190103_06_00.grib2"
filename_11 = "./fnl_data/fnl_20190103_12_00.grib2"
filename_12 = "./fnl_data/fnl_20190103_18_00.grib2"
filename_13 = "./fnl_data/fnl_20190104_00_00.grib2"
filename_14 = "./fnl_data/fnl_20190104_06_00.grib2"
filename_15 = "./fnl_data/fnl_20190104_12_00.grib2"
filename_16 = "./fnl_data/fnl_20190104_18_00.grib2"
filename_17 = "./fnl_data/fnl_20190105_00_00.grib2"
filename_18 = "./fnl_data/fnl_20190105_06_00.grib2"
filename_19 = "./fnl_data/fnl_20190105_12_00.grib2"
filename_20 = "./fnl_data/fnl_20190105_18_00.grib2"
filename_21 = "./fnl_data/fnl_20190106_00_00.grib2"
filename_22 = "./fnl_data/fnl_20190106_06_00.grib2"
filename_23 = "./fnl_data/fnl_20190106_12_00.grib2"
filename_24 = "./fnl_data/fnl_20190106_18_00.grib2"
filename_25 = "./fnl_data/fnl_20190107_00_00.grib2"
filename_26 = "./fnl_data/fnl_20190107_06_00.grib2"
filename_27 = "./fnl_data/fnl_20190107_12_00.grib2"
filename_28 = "./fnl_data/fnl_20190107_18_00.grib2"
filename_29 = "./fnl_data/fnl_20190108_00_00.grib2"
filename_30 = "./fnl_data/fnl_20190108_06_00.grib2"
filename_31 = "./fnl_data/fnl_20190108_12_00.grib2"
filename_32 = "./fnl_data/fnl_20190108_18_00.grib2"
filename_33 = "./fnl_data/fnl_20190109_00_00.grib2"
filename_34 = "./fnl_data/fnl_20190109_06_00.grib2"
filename_35 = "./fnl_data/fnl_20190109_12_00.grib2"
filename_36 = "./fnl_data/fnl_20190109_18_00.grib2"
filename_37 = "./fnl_data/fnl_20190110_00_00.grib2"
filename_38 = "./fnl_data/fnl_20190110_06_00.grib2"
filename_39 = "./fnl_data/fnl_20190110_12_00.grib2"
filename_40 = "./fnl_data/fnl_20190110_18_00.grib2"
filename_41 = "./fnl_data/fnl_20190111_00_00.grib2"
filename_42 = "./fnl_data/fnl_20190111_06_00.grib2"
filename_43 = "./fnl_data/fnl_20190111_12_00.grib2"
filename_44 = "./fnl_data/fnl_20190111_18_00.grib2"
filename_45 = "./fnl_data/fnl_20190112_00_00.grib2"
filename_46 = "./fnl_data/fnl_20190112_06_00.grib2"
filename_47 = "./fnl_data/fnl_20190112_12_00.grib2"
filename_48 = "./fnl_data/fnl_20190112_18_00.grib2"
f_1 = addfile(filename_1, "r")
f_2 = addfile(filename_2, "r")
f_3 = addfile(filename_3, "r")
f_4 = addfile(filename_4, "r")
f_5 = addfile(filename_5, "r")
f_6 = addfile(filename_6, "r")
f_7 = addfile(filename_7, "r")
f_8 = addfile(filename_8, "r")
f_9 = addfile(filename_9, "r")
f_10 = addfile(filename_10, "r")
f_11 = addfile(filename_11, "r")
f_12 = addfile(filename_12, "r")
f_13 = addfile(filename_13, "r")
f_14 = addfile(filename_14, "r")
f_15 = addfile(filename_15, "r")
f_16 = addfile(filename_16, "r")
f_17 = addfile(filename_17, "r")
f_18 = addfile(filename_18, "r")
f_19 = addfile(filename_19, "r")
f_20 = addfile(filename_20, "r")
f_21 = addfile(filename_21, "r")
f_22 = addfile(filename_22, "r")
f_23 = addfile(filename_23, "r")
f_24 = addfile(filename_24, "r")
f_25 = addfile(filename_25, "r")
f_26 = addfile(filename_26, "r")
f_27 = addfile(filename_27, "r")
f_28 = addfile(filename_28, "r")
f_29 = addfile(filename_29, "r")
f_30 = addfile(filename_30, "r")
f_31 = addfile(filename_31, "r")
f_32 = addfile(filename_32, "r")
f_33 = addfile(filename_33, "r")
f_34 = addfile(filename_34, "r")
f_35 = addfile(filename_35, "r")
f_36 = addfile(filename_36, "r")
f_37 = addfile(filename_37, "r")
f_38 = addfile(filename_38, "r")
f_39 = addfile(filename_39, "r")
f_40 = addfile(filename_40, "r")
f_41 = addfile(filename_41, "r")
f_42 = addfile(filename_42, "r")
f_43 = addfile(filename_43, "r")
f_44 = addfile(filename_44, "r")
f_45 = addfile(filename_45, "r")
f_46 = addfile(filename_46, "r")
f_47 = addfile(filename_47, "r")
f_48 = addfile(filename_48, "r")
lev_1 = f_1->MSLET_P0_L101_GLL0
lev_2 = f_2->MSLET_P0_L101_GLL0
lev_3 = f_3->MSLET_P0_L101_GLL0
lev_4 = f_4->MSLET_P0_L101_GLL0
lev_5 = f_5->MSLET_P0_L101_GLL0
lev_6 = f_6->MSLET_P0_L101_GLL0
lev_7 = f_7->MSLET_P0_L101_GLL0
lev_8 = f_8->MSLET_P0_L101_GLL0
lev_9 = f_9->MSLET_P0_L101_GLL0
lev_10 = f_10->MSLET_P0_L101_GLL0
lev_11 = f_11->MSLET_P0_L101_GLL0
lev_12 = f_12->MSLET_P0_L101_GLL0
lev_13 = f_13->MSLET_P0_L101_GLL0
lev_14 = f_14->MSLET_P0_L101_GLL0
lev_15 = f_15->MSLET_P0_L101_GLL0
lev_16 = f_16->MSLET_P0_L101_GLL0
lev_17 = f_17->MSLET_P0_L101_GLL0
lev_18 = f_18->MSLET_P0_L101_GLL0
lev_19 = f_19->MSLET_P0_L101_GLL0
lev_20 = f_20->MSLET_P0_L101_GLL0
lev_21 = f_21->MSLET_P0_L101_GLL0
lev_22 = f_22->MSLET_P0_L101_GLL0
lev_23 = f_23->MSLET_P0_L101_GLL0
lev_24 = f_24->MSLET_P0_L101_GLL0
lev_25 = f_25->MSLET_P0_L101_GLL0
lev_26 = f_26->MSLET_P0_L101_GLL0
lev_27 = f_27->MSLET_P0_L101_GLL0
lev_28 = f_28->MSLET_P0_L101_GLL0
lev_29 = f_29->MSLET_P0_L101_GLL0
lev_30 = f_30->MSLET_P0_L101_GLL0
lev_31 = f_31->MSLET_P0_L101_GLL0
lev_32 = f_32->MSLET_P0_L101_GLL0
lev_33 = f_33->MSLET_P0_L101_GLL0
lev_34 = f_34->MSLET_P0_L101_GLL0
lev_35 = f_35->MSLET_P0_L101_GLL0
lev_36 = f_36->MSLET_P0_L101_GLL0
lev_37 = f_37->MSLET_P0_L101_GLL0
lev_38 = f_38->MSLET_P0_L101_GLL0
lev_39 = f_39->MSLET_P0_L101_GLL0
lev_40 = f_40->MSLET_P0_L101_GLL0
lev_41 = f_41->MSLET_P0_L101_GLL0
lev_42 = f_42->MSLET_P0_L101_GLL0
lev_43 = f_43->MSLET_P0_L101_GLL0
lev_44 = f_44->MSLET_P0_L101_GLL0
lev_45 = f_45->MSLET_P0_L101_GLL0
lev_46 = f_46->MSLET_P0_L101_GLL0
lev_47 = f_47->MSLET_P0_L101_GLL0
lev_48 = f_48->MSLET_P0_L101_GLL0
lev_avg = (lev_1 + lev_2 + lev_3 + lev_4 + lev_5 + lev_6 + \
lev_7 + lev_8 + lev_9 + lev_10 + lev_11 + lev_12 + \
lev_13 + lev_14 + lev_15 + lev_16 + lev_17 + lev_18 + \
lev_19 + lev_20 + lev_21 + lev_22 + lev_23 + lev_24 + \
lev_25 + lev_26 + lev_27 + lev_28 + lev_29 + lev_30 + \
lev_31 + lev_32 + lev_33 + lev_34 + lev_35 + lev_36 + \
lev_37 + lev_38 + lev_39 + lev_40 + lev_41 + lev_42 + \
lev_43 + lev_44 + lev_45 + lev_46 + lev_47 + lev_48) / 48
DATADir = "./fnl_data/"
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190101*")
f = addfiles(FILES,"r")
U=f[:]->UGRD_P0_L100_GLL0
U500_1 = U(30,:,:)
U500_2 = U(61,:,:)
U500_3 = U(92,:,:)
U500_4 = U(123,:,:)
U500_sum1 = (U500_1+U500_2+U500_3+U500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190102*")
f = addfiles(FILES,"r")
U=f[:]->UGRD_P0_L100_GLL0
U500_1 = U(30,:,:)
U500_2 = U(61,:,:)
U500_3 = U(92,:,:)
U500_4 = U(123,:,:)
U500_sum2 = (U500_1+U500_2+U500_3+U500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190103*")
f = addfiles(FILES,"r")
U=f[:]->UGRD_P0_L100_GLL0
U500_1 = U(30,:,:)
U500_2 = U(61,:,:)
U500_3 = U(92,:,:)
U500_4 = U(123,:,:)
U500_sum3 = (U500_1+U500_2+U500_3+U500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190104*")
f = addfiles(FILES,"r")
U=f[:]->UGRD_P0_L100_GLL0
U500_1 = U(30,:,:)
U500_2 = U(61,:,:)
U500_3 = U(92,:,:)
U500_4 = U(123,:,:)
U500_sum4 = (U500_1+U500_2+U500_3+U500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190105*")
f = addfiles(FILES,"r")
U=f[:]->UGRD_P0_L100_GLL0
U500_1 = U(30,:,:)
U500_2 = U(61,:,:)
U500_3 = U(92,:,:)
U500_4 = U(123,:,:)
U500_sum5 = (U500_1+U500_2+U500_3+U500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190106*")
f = addfiles(FILES,"r")
U=f[:]->UGRD_P0_L100_GLL0
U500_1 = U(30,:,:)
U500_2 = U(61,:,:)
U500_3 = U(92,:,:)
U500_4 = U(123,:,:)
U500_sum6 = (U500_1+U500_2+U500_3+U500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190107*")
f = addfiles(FILES,"r")
U=f[:]->UGRD_P0_L100_GLL0
U500_1 = U(30,:,:)
U500_2 = U(61,:,:)
U500_3 = U(92,:,:)
U500_4 = U(123,:,:)
U500_sum7 = (U500_1+U500_2+U500_3+U500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190108*")
f = addfiles(FILES,"r")
U=f[:]->UGRD_P0_L100_GLL0
U500_1 = U(30,:,:)
U500_2 = U(61,:,:)
U500_3 = U(92,:,:)
U500_4 = U(123,:,:)
U500_sum8 = (U500_1+U500_2+U500_3+U500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190109*")
f = addfiles(FILES,"r")
U=f[:]->UGRD_P0_L100_GLL0
U500_1 = U(30,:,:)
U500_2 = U(61,:,:)
U500_3 = U(92,:,:)
U500_4 = U(123,:,:)
U500_sum9 = (U500_1+U500_2+U500_3+U500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190110*")
f = addfiles(FILES,"r")
U=f[:]->UGRD_P0_L100_GLL0
U500_1 = U(30,:,:)
U500_2 = U(61,:,:)
U500_3 = U(92,:,:)
U500_4 = U(123,:,:)
U500_sum10 = (U500_1+U500_2+U500_3+U500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190111*")
f = addfiles(FILES,"r")
U=f[:]->UGRD_P0_L100_GLL0
U500_1 = U(30,:,:)
U500_2 = U(61,:,:)
U500_3 = U(92,:,:)
U500_4 = U(123,:,:)
U500_sum11 = (U500_1+U500_2+U500_3+U500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190112*")
f = addfiles(FILES,"r")
U=f[:]->UGRD_P0_L100_GLL0
U500_1 = U(30,:,:)
U500_2 = U(61,:,:)
U500_3 = U(92,:,:)
U500_4 = U(123,:,:)
U500_sum12 = (U500_1+U500_2+U500_3+U500_4)
U500_sum_all1 = U500_sum1 + U500_sum2 + U500_sum3 + U500_sum4 + U500_sum5 + U500_sum6 + U500_sum7 + U500_sum8 + U500_sum9 + U500_sum10 + U500_sum11 + U500_sum12
U500_ave_1 = U500_sum_all1/48 ;;48个fnl文件
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190101*")
f = addfiles(FILES,"r")
V=f[:]->VGRD_P0_L100_GLL0
V500_1 = V(30,:,:)
V500_2 = V(61,:,:)
V500_3 = V(92,:,:)
V500_4 = V(123,:,:)
V500_sum1 = (V500_1+V500_2+V500_3+V500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190102*")
f = addfiles(FILES,"r")
V=f[:]->VGRD_P0_L100_GLL0
V500_1 = V(30,:,:)
V500_2 = V(61,:,:)
V500_3 = V(92,:,:)
V500_4 = V(123,:,:)
V500_sum2 = (V500_1+V500_2+V500_3+V500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190103*")
f = addfiles(FILES,"r")
V=f[:]->VGRD_P0_L100_GLL0
V500_1 = V(30,:,:)
V500_2 = V(61,:,:)
V500_3 = V(92,:,:)
V500_4 = V(123,:,:)
V500_sum3 = (V500_1+V500_2+V500_3+V500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190104*")
f = addfiles(FILES,"r")
V=f[:]->VGRD_P0_L100_GLL0
V500_1 = V(30,:,:)
V500_2 = V(61,:,:)
V500_3 = V(92,:,:)
V500_4 = V(123,:,:)
V500_sum4 = (V500_1+V500_2+V500_3+V500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190105*")
f = addfiles(FILES,"r")
V=f[:]->VGRD_P0_L100_GLL0
V500_1 = V(30,:,:)
V500_2 = V(61,:,:)
V500_3 = V(92,:,:)
V500_4 = V(123,:,:)
V500_sum5 = (V500_1+V500_2+V500_3+V500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190106*")
f = addfiles(FILES,"r")
V=f[:]->VGRD_P0_L100_GLL0
V500_1 = V(30,:,:)
V500_2 = V(61,:,:)
V500_3 = V(92,:,:)
V500_4 = V(123,:,:)
V500_sum6 = (V500_1+V500_2+V500_3+V500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190107*")
f = addfiles(FILES,"r")
V=f[:]->VGRD_P0_L100_GLL0
V500_1 = V(30,:,:)
V500_2 = V(61,:,:)
V500_3 = V(92,:,:)
V500_4 = V(123,:,:)
V500_sum7 = (V500_1+V500_2+V500_3+V500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190108*")
f = addfiles(FILES,"r")
V=f[:]->VGRD_P0_L100_GLL0
V500_1 = V(30,:,:)
V500_2 = V(61,:,:)
V500_3 = V(92,:,:)
V500_4 = V(123,:,:)
V500_sum8 = (V500_1+V500_2+V500_3+V500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190109*")
f = addfiles(FILES,"r")
V=f[:]->VGRD_P0_L100_GLL0
V500_1 = V(30,:,:)
V500_2 = V(61,:,:)
V500_3 = V(92,:,:)
V500_4 = V(123,:,:)
V500_sum9 = (V500_1+V500_2+V500_3+V500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190110*")
f = addfiles(FILES,"r")
V=f[:]->VGRD_P0_L100_GLL0
V500_1 = V(30,:,:)
V500_2 = V(61,:,:)
V500_3 = V(92,:,:)
V500_4 = V(123,:,:)
V500_sum10 = (V500_1+V500_2+V500_3+V500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190111*")
f = addfiles(FILES,"r")
V=f[:]->VGRD_P0_L100_GLL0
V500_1 = V(30,:,:)
V500_2 = V(61,:,:)
V500_3 = V(92,:,:)
V500_4 = V(123,:,:)
V500_sum11 = (V500_1+V500_2+V500_3+V500_4)
FILES = systemfunc (" ls -1 " + DATADir + "fnl_20190112*")
f = addfiles(FILES,"r")
V=f[:]->VGRD_P0_L100_GLL0
V500_1 = V(30,:,:)
V500_2 = V(61,:,:)
V500_3 = V(92,:,:)
V500_4 = V(123,:,:)
V500_sum12 = (V500_1+V500_2+V500_3+V500_4)
V500_sum_all1 = V500_sum1 + V500_sum2 + V500_sum3 + V500_sum4 + V500_sum5 + V500_sum6 + V500_sum7 + V500_sum8 + V500_sum9 + V500_sum10 + V500_sum11 + V500_sum12
V500_ave_1 = V500_sum_all1/48 ;;48个fnl文件
;;;;;复制属性的文件
m = addfile("./fnl_data/fnl_20190102_18_00.grib2","r")
sea_m = m->MSLET_P0_L101_GLL0
u_m = m->UGRD_P0_L100_GLL0(13,:,:)
v_m = m->VGRD_P0_L100_GLL0(13,:,:)
lev_avg_hpa = lev_avg/100
;;;复制属性
copy_VarMeta(sea_m, lev_avg_hpa)
copy_VarMeta(u_m, U500_ave_1)
copy_VarMeta(v_m, V500_ave_1)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
sealev = lev_avg_hpa
UU = U500_ave_1
VV = V500_ave_1
;print(sealev)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
plot = new(1,graphic)
res = True
res@gsnFrame = False
res@gsnDraw = False
res@mpMinLatF = 0.
res@mpMaxLatF = 70.
res@mpMinLonF = 60.
res@mpMaxLonF = 150.
res@mpDataSetName = "./Earth..4"
res@mpDataBaseVersion = "MediumRes" ; or "Ncarg4_1"
res@mpAreaMaskingOn = True
res@mpMaskAreaSpecifiers = (/"China"/)
res@mpOutlineSpecifiers = (/"China","China:Provinces"/)
;res@mpOutlineSpecifiers = (/"China"/)
res@mpFillOn = False ;;;;不填充地图
res@tmXTLabelsOn = False
res@tmYROn = False
res@tmXTOn = False
res@gsnLeftString = ""
res@gsnRightString = ""
res@pmTickMarkDisplayMode="Always"
res@tmXBLabelFontHeightF = 0.015
res@tmXBLabelFontThicknessF = 0.03
res@tmYLLabelFontHeightF = 0.015
res@tmYLLabelFontThicknessF = 0.03
;res@gsnTitleString = "500hPa"
res@cnFillOn = True
res@cnLinesOn = False
res@cnLineLabelsOn = False
res@cnInfoLabelOn = False
res@cnLevelSelectionMode = "ManualLevels"
res@lbOrientation = "vertical"
; res@cnLevelSelectionMode = "ExplicitLevels"
; res@cnMinLevelValF = 4850
; res@cnMaxLevelValF = 5910
; res@cnLevelSpacingF = 80
res@gsnSpreadColors = True
res@lbBoxLinesOn = False
res@lbLabelAutoStride= True
res@lbLabelBarOn = False ;;;关闭小图色标
res@lbLabelFontThicknessF = 20. ;;;font
res@lbLabelFontAspectF = 0.5
res@lbTitleString = "gpm"
res@lbTitlePosition = "Right"
res@pmLabelBarOrthogonalPosF = 0.15
res@pmLabelBarParallelPosF = 0.5
res@pmLabelBarSide = "Left"
res@lbLabelFontHeightF = 0.005
res@cnFillOn = False
res@cnLinesOn = True
res@cnLineLabelsOn = True
res@cnLineLabelFont = "times-roman"
res@cnLineLabelFontHeightF = 0.015 ;;标签字体大小
res@cnLineLabelInterval = 1
res@cnLineThicknessF = 2.5 ;;line的粗细
res@cnMissingValPerimOn = True
res@cnLineLabelBackgroundColor = 0
res@cnInfoLabelOn = False
;res@cnLevelSelectionMode = "ExplicitLevels"
;res@cnLevels = (/516,520,524,528,532,536,,,,,/)
;res@cnLevels = (/496,500,504,508,512,516,520,524,528,532,536,540,544,548,552,556,560,564,568,572,576,580,584,588/)
;;;修改x y轴字体大小
res@tmXBLabelFontHeightF = 0.013
res@tmYLLabelFontHeightF = 0.013
res@tmXBLabelFont = "times-roman"
res@tmYLLabelFont = "times-roman"
;;加上°N
res@pmTickMarkDisplayMode = "Always" ; Turn on map tickmarks
;res@lbTitleString = "500hPa"
;res@tiXAxisString="gpm"
; res@gsnRBString = "gpm" ;invalid
; res@gsnPanelFigureStrings = "gpm"
; res@tmXTLabelsOn= True
; res@tiXAxisString="gpm"
; res@tiXAxisPosition = "Right"
; res@tiXAxis = True
;;;画矢量
vec_res = True
vec_res@gsnDraw=False
vec_res@gsnFrame=False
vec_res@vcMinDistanceF = 0.025
vec_res@vcRefLengthF = 0.025 ;;参考矢量长度
vec_res@vcMinFracLengthF = 0.01 ;;长度
vec_res@vcMinMagnitudeF = 0.05 ;;less than 10,delete
vec_res@vcMaxMagnitudeF = 50 ;;more than 10,delete
vec_res@vcMinDistanceF = 0.03 ;;距离
vec_res@vcLineArrowThicknessF = 1.5 ;;粗细
;vec_res@vcGlyphStyle = "CurlyVector" ;;矢量箭头形状
;vec_res@vcGlyphStyle = "LineArrow"
vec_res@vcGlyphStyle = "WindBarb"
;;包括LineArrow，FillArrow，CurlyVector，WindBarb
;;LineArrow vcLineArrowHeadMinSizeF=0.01 ;控制箭头的大小
;;LineArrow
vec_res@vcRefAnnoOn = True
vec_res@vcRefMagnitudeF = 4 ;;参考矢量的值
vec_res@vcRefAnnoString1On = True ;;;false关掉指示风/箭头上面标注
vec_res@vcRefAnnoString1 ="4 m/s"
;vec_res@vcLabelFontHeightF = 5
;res2@tmXBLabelFont = 26 ;; Set font to times-roman
;res2@tmYLLabelFont = 26
;res2@lbLabelFont = 26
;res2@lbLabelFontHeightF = 0.013 ;;font size
vec_res@vcRefAnnoFontHeightF = 0.01
vec_res@vcRefAnnoFont = 26
vec_res@vcRefAnnoString2On = True ;箭头下面的标注
vec_res@vcRefAnnoString2 = ""
;vec_res@vcRefAnnoOrthogonalPosF = -1.305 ;标签位置，负为向上移，正为向下
vec_res@vcRefAnnoOrthogonalPosF = -1.09
vec_res@vcRefAnnoParallelPosF = 1 ;标签位置，负为左移，正为右移
vec_res@vcVectorDrawOrder = "Postdraw"
vec_res@gsnLeftString = ""
vec_res@gsnRightString = ""
plot(0) = gsn_csm_contour_map(wks,sealev,res)
wind0 = gsn_csm_vector(wks,UU,VV,vec_res)
overlay(plot(0), wind0)
pres = True
;;;;;;; pres@txString = "500hPa"
;pres@txFontHeightF = 0.02
pres@gsnPanelLabelBar = True
pres@pmLabelBarOrthogonalPosF = -0.015
pres@pmLabelBarWidthF = 0.5
pres@lbLabelFontHeightF = 0.012
pres@gsnPanelFigureStrings = (/"(d)"/)
pres@gsnPanelFigureStringsFontHeightF = 0.026
pres@gsnPanelFigureStringsFont = "times-roman"
;pres@gsnPanelFigureStringsHeightF = 0.015
pres@amJust = "TopLeft"
gsn_panel(wks,plot,(/1,1/),pres)
pres@lbOrientation = "vertical"
pres@lbLabelFontHeightF =.012
end