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

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 05:40:30 2014