Removing for inland data

Hi,
Using Sea Ice Concentration data,
I am trying to find total SIC area(values in the non-land area),
but I couldn't removing inland data
my code as follows

Thank you for your help in advance

>>>>>>>>>>>>>>>>>>>>>>>>>>>>
begin
        year = 2006

        f = addfile("SON_"+year+".nc","r") ;#

        lon = f->longitude
        lat = f->latitude
        time = f->time
        sic = f->sic

        ; float lon(lon) ;

                lon@units = "degrees_east" ;

        ; float lat(lat) ;

                lat@units = "degrees_north" ;

        ; float sic(time, lat, lon) ;

                sic@long_name = "Sea-ice concentration" ;

                sic@units = "%" ;

                sic@grid_name = "linear" ;

                sic@grid_type = "uniform" ;

                sic@level_description = "Earth surface" ;

                sic@time_statistic = "average" ;

        ; sic@_FillValue = 1.e+20f ;

        ; lm = addfile("$NCARG_ROOT/lib/ncarg/data/cdf/landsea.nc", "r") ;
landsea mask
        ; lsmask = lm->LSMASK
        ; shpfile = "/home/hkjin/seaice/GSHHS_shp/c/GSHHS_c_L1.shp"
        ; land_mask = create_mask_from_shapefile(sic,shpfile) ;#

        ; lsm = landsea_mask(lsmask,lat,lon)

;==============

         wks = gsn_open_wks("x11","sic_"+year+"_SON") ;#

         res = True

         gsn_define_colormap(wks,"seaice_2")

        ;********************** map config.

          res@mpOutlineOn = True ; turn off outline

          res@mpFillOn = True ; turn off map fill

          res@mpFillDrawOrder = "Predraw"

          res@tiMainString = "Seasonal Sea Ice Concentration :
"+year+", SON" ;#

          res@gsnLeftString = "+year+" ; line#

          res@cnFillOn = True ; turn on color fill

          res@cnFillDrawOrder = "Predraw"

        res@gsnPolar = "NH" ; "NH" or "SH"

        res@mpMinLatF = 45 ; specify min lat

        plot = gsn_csm_contour_map_polar(wks,sic(0,:,:),res)

end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri May 17 09:15:44 2013