Re: fatal:mask: dimension sizes of parameter 0 and parameter 1 do not match

Hi Ciara,

On line 17:

    data_mask=mask(data,shpfile,opt)

I think you want to use function shapefile_mask_data(). This is an
external function that can be obtained by looking at example #14 at
the link below; that example is also one of performing masking in the
fashion you are trying to do:

     http://www.ncl.ucar.edu/Applications/shapefiles.shtml

I hope that helps…

Rick

On Wed, 5 Mar 2014 12:40:17 +0000
  Ciara O'Hara <cohara1985@hotmail.com> wrote:
> Hi,
>
> I am running the script below, and getting the following error
>message:
>
> ------------------------------------------------------------------------------------------------------------------------
> fatal:mask: dimension sizes of parameter 0 and parameter 1 do not
>match
> fatal:["Execute.c":8128]:Execute: Error occurred at or near line 17
>in file /home/cohara/scripts/rainfall/for_map/contour_map_script.ncl
> ------------------------------------------------------------------------------------------------------------------------
>
> I am trying to use a shapefile to mask out the sea, so that my plot
>has nice clean edges around the coast. I understand that it is a
>latitude and longitude dimension size problem, but I'm not sure how
>to fix it. I've run print_Var_Summary on one of my data variables
>(they are all the same except for the date attached), here is the
>output:
>
> ------------------------------------------------
> ncl 7> printVarSummary(data)
>
> Variable: data
> Type: float
> Total Size: 188136 bytes
> 47034 values
> Number of Dimensions: 2
> Dimensions and sizes: [lat | 234] x [lon | 201]
> Coordinates:
> lat: [ 50..56.99]
> lon: [ -11.. -5]
> Number Of Attributes: 4
> height : 0
> time : 19810930
> table : 1
> _FillValue : -9e+33
> ------------------------------------------------
>
> Any suggestions would be appreciated,
>
> Best Regards,
> Ciara
>
> ---------------------------------
> script:
> ---------------------------------
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
>
> begin
> n=0
> i=1981
> print("Hi")
> y=systemfunc("ls
>/home/cohara/RainfallData/for_map/timsum_*.nc")
> x=addfiles(y,"r")
> do while(i.le.2013)
> data=x[n]->var61(0,0,:,:)
> shpfile="/home/cohara/IRL_adm/IRL_adm0.shp"
> opt=True
> opt@return_mask=True
> data_mask=mask(data,shpfile,opt)
> data_land_mask=where(data_mask.eq.1,data,data@_FillValue)
> copy_VarMeta(data,data_land_mask)
> plres=True
> plres@gsEdgesOn=True
> plres@gsLineColor="black"
> plres@gsEdgeColor="black"
> resources=True
> resources@gsnMaximize=True
> resources@gsnDraw=False
> resources@gsnFrame=False
> resources@tiMainString="Total Annual Rainfall " +
>sprinti("%0.4i",i)
> resources@cnMonoLineColor=False ;Turn of the drawing
>of contour lines in 1 colour
> resources@tiMainFontColor="black"
> resources@tiMainFont="times-roman"
> resources@cnFillOn=True ;Turn on contour line
>fill
> resources@cnMonoFillPattern=True ;Turn on single fill
>pattern (default)
> resources@cnMonoFillColor=False ;(default)
> resources@cnMonoLineColor=False ;(default)
> resources@cnLinesOn=False ;(default)
> resources@tmYLOn=False ;Removes tick marks
>from left of y axis
> resources@tmYROn=False
> resources@tmXTOn=False
> resources@tmXBOn=False
> resources@cnLineLabelsOn=False ;Removes contour line
>labels
> resources@cnInfoLabelOn=False ;Removes the
>informational label
> resources@pmLabelBarDisplayMode="Always";Turns on the label
>bar
> resources@lbLabelBarOn=True
> resources@lbBoxLinesOn=False ;Removes lines from
>between the label bar colours
> resources@lbTitleOn=True ;Turns on label bar
>title
> resources@lbTitleString="Rainfall (mm)"
> resources@lbPerimOn=False ;Removes perimeter
>from label bar
> resources@lbTitleFontColor="black"
> resources@lbTitleFont="times-roman"
> resources@lbOrientation="Horizontal"
> resources@lbLabelPosition="Bottom"
> resources@lbTitlePosition="Top"
> ;resources@cnExplicitLegendLabelsOn=True
> ;resources@lgLabelStrings=(/"4","5","6","7","8","9","10","11","12","13","14"/)
> ;resources@cnLineLabelStrings=resources@lbLabelStrings
> resources@lbTitleDirection="Across"
> resources@lbTitleFontHeightF=0.02
> resources@lbLabelFont="times-roman"
> resources@lbLabelFontColor="black"
> ;resources@cnLevelSelectionMode="ManualLevels" ;Set the
>label bar levels manually
> ;resources@cnLevelSpacingF=0.2
> resources@cnLevelSelectionMode="ExplicitLevels"
> resources@cnLevels=(/0,200,400,600,800,1000,1200,1400,1600,1800,2000,2200,2400,2600,2800,3000,3200,3400,3600,3800,4000/)
> resources@mpProjection="Mercator" ;Map Projection
> resources@gsnAddCyclic=False
> resources@mpLimitMode="LatLon" ;The viewable portion
>of the map determined by latitude and longitude
> resources@mpMinLatF=51
> resources@mpMaxLatF=55.7
> resources@mpMinLonF=-11
> resources@mpMaxLonF=-5.2
> resources@mpGridAndLimbOn=True
> resources@mpDataSetName="Earth..4" ;The data set used
>(most up to date)
> resources@mpDataBaseVersion="MediumRes"
> resources@mpDataResolution="Finest" ;Map render
>resolution
> resources@mpOutlineOn=False
> resources@mpFillOn=False
> ;resources@vpWidthF=1
> ;resources@vpHeightF=0.7
> wtype="png"
> ;wtype@wkWidth=2500
> ;wtype@wkHeight=2500
> ;sfiles=(/"/home/cohara/IRL_adm/IRL_adm0","/home/cohara/GBR_adm/GBR_adm0"/)+".shp"
> xwks=gsn_open_wks(wtype,"Total_" + sprinti("%0.4i",i))
> gsn_define_colormap(xwks,"precip4_11lev")
> plot=gsn_csm_contour_map(xwks,data_land_mask,resources)
> poly0=gsn_add_shapefile_polylines(xwks,plot,"/home/cohara/IRL_adm/IRL_adm1.shp",plres)
> poly1=gsn_add_shapefile_polylines(xwks,plot,"/home/cohara/GBR_adm/GBR_adm2.shp",plres)
> poly2=gsn_add_shapefile_polylines(xwks,plot,"/home/cohara/IRL_adm/IRL_adm0.shp",plres)
> draw(plot)
> frame(xwks)
> n=n+1
> i=i+1
> end do
> end
>
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Mar 5 08:37:49 2014