load "/home/server/student/homes/naraligi/contributed/gsn_code.ncl" load "/home/server/student/homes/naraligi/contributed/gsn_csm.ncl" load "/home/server/student/homes/naraligi/contributed/shea_util.ncl" load "/home/server/student/homes/naraligi/contributed/contributed.ncl" begin fil = asciiread("events_point.txt",(/37,9/),"integer") phase=fil(:,2) rfil ="~/research/composite/towspd.nc" wdfil ="~/research/composite/wspd.nc" wdfil1 ="~/research/qkscat/q_filt.nc" vcfil = addfile(wdfil,"r") tofil = addfile(rfil,"r") isofil = addfile(wdfil1,"r") phase@units = "days since 0-1-1 00:00:0.0" ymd = ut_calendar(phase,0) mon=ymd(:,1) i = ind(mon.eq.11.or.mon.eq.12.or.mon.eq.1.or. \ mon.eq.2.or.mon.eq.3.or.mon.eq.4) u1=isofil->u(phase(i),:,:) v1=isofil->v(phase(i),:,:) towsp=tofil->wspd(phase(i),:,:) vcwsp=vcfil->wspd(phase(i),:,:) edwsp=towsp-vcwsp edwsp!0="time" edwsp!1="lat" edwsp!2="lon" lat = u1&lat lon = u1&lon mnu1=dim_avg(u1(lat|:,lon|:,time|:)) mnv1=dim_avg(v1(lat|:,lon|:,time|:)) mntowsp=dim_avg(towsp(lat|:,lon|:,time|:)) mnvcwsp=dim_avg(vcwsp(lat|:,lon|:,time|:)) mnedwsp=dim_avg(edwsp(lat|:,lon|:,time|:)) mnu1!0="lat" mnv1!0="lat" mntowsp!0="lat" mnvcwsp!0="lat" mnedwsp!0="lat" mnu1!1="lon" mnv1!1="lon" mntowsp!1="lon" mnvcwsp!1="lon" mnedwsp!1="lon" print("error is after this") mnu1&lat = lat mnv1&lat = lat mntowsp&lat = lat mnvcwsp&lat = lat mnedwsp&lat = lat mnu1&lon = lon mnv1&lon = lon mntowsp&lon = lon mnvcwsp&lon = lon mnedwsp&lon = lon ;===================================================================================== wks = gsn_open_wks("ps","eddyphase3_ssum") gsn_define_colormap(wks,"gui_default") plot = new(3,graphic) res = True res@mpProjection = "Mercator" res@mpPerimOn = True res@mpFillOn = True res@mpCenterLonF = 180 res@mpLandFillColor = "black" ; set land to be gray res@gsnFrame = False ; don't advance frame res@gsnDraw = False ; don't draw res@gsnScalarContour = True res@cnFillOn = True res@cnLineLabelsOn = False res@cnLinesOn = False res@lbLabelBarOn = False res@lbLabelFontAspectF = 0.5 ;res@lbOrientation = "Vertical" ; vertical label bar ;res@pmLabelBarOrthogonalPosF = 0.18 ; move label bar closer res@mpLimitMode = "LatLon" res@mpMinLatF = -15 res@mpMaxLatF =15 res@mpMinLonF = 130 res@mpMaxLonF = 205 ;res@vpXF = 0.12 ;res@vpYF = 3.0 ;res@vpWidthF = 0.8 ;res@vpHeightF = 1.0 res@pmTickMarkDisplayMode = "Always" ;res@pmTickMarkZone = 0.5 res@gsnAddCyclic = False res@gsnSpreadColors = True res@gsnMaximize = True ;res@txFontHeightF = 0.008 ;res@vcRefMagnitudeF = 4.0 ; define vector ref mag ;res@vcRefLengthF = 0.055 ; define length of vec ref ;res@vcRefAnnoOrthogonalPosF = -1.0 ; move ref vector ;res@vcRefAnnoArrowLineColor = "black" ; change ref vector color ;res@vcRefAnnoArrowUseVecColor = False ; don't use vec color for ref ; res@vcGlyphStyle = "CurlyVector" ; turn on curley vectors ; res@vcLineArrowColor = "black" ; change vector color ; res@vcLineArrowThicknessF = 1.0 ; change vector thickness ; res@vcVectorDrawOrder = "PostDraw" ; draw vectors last ; res@vcFillArrowFillColor = "black" res@cnLevelSelectionMode = "ManualLevels" ; set manual contour levels res@cnMinLevelValF = -0.5 ; set min contour level res@cnMaxLevelValF = 0.5 ; set max contour level res@cnLevelSpacingF = 0.1 ; set contour spacing plot(0) = gsn_csm_vector_scalar_map(wks,mnu1,mnv1,mntowsp,res) plot(1) = gsn_csm_vector_scalar_map(wks,mnu1,mnv1,mnvcwsp,res) plot(2) = gsn_csm_vector_scalar_map(wks,mnu1,mnv1,mnedwsp,res) resP = True resP@gsnPanelLabelBar = True ; add common colorbar ;resP@txString = "Westerly phase (phase 7)" resP@lbLabelFontHeightF = 0.01 resP@gsnMaximize = True resP@gsnPanelBottom = 0.05 resP@gsnPanelFigureStrings= (/"a)","b)","c)"/) ; add strings to panel resP@gsnPanelFigureStringsFontHeightF = 0.010 resP@amJust = "BottomLeft" resP@gsnPanelFigureStrings= (/"a) Total","b) Vector","c) Eddies"/) ; add strings to panel resP@gsnPanelFigureStringsPerimOn = True resP@gsnPanelFigureStringsFontHeightF = 0.012 resP@gsnFrame = False ; don't advance frame gsn_panel(wks,plot,(/3,1/),resP) txres = True txres@txFontHeightF = 0.012 gsn_text_ndc(wks,"Figure : Easterly Wind Speed Annomalies for southern summer " ,0.5,0.02,txres) frame(wks) end