Defining sic&lat and sic&lon

From: 진현근 <narshyas2_at_nyahnyahspammersnyahnyah>
Date: Fri May 17 2013 - 14:32:36 MDT

thank you for your help
but I have another problem.
getting your advice, I added some sentences into my code
but just got error messages.

fatal:(lat) is not a named dimension in variable (sic).
fatal:["Execute.c":7743]:Execute: Error occurred at or near line 49 in file

>>>>>>>>>>>>>>>>>>>>>>>>>>>>
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")*
*lsmask = lm->LSMASK
lsm = landsea_mask(lsmask,sic&lat,sic&lon) ; create land sea mask for
sic array
sic = mask(sic,conform(sic,lsm,(/1,2/)).ge.1,False) ; conform 2D lsm
array to 3D array, mask out land values*

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

         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 14:32:49 2013

This archive was generated by hypermail 2.1.8 : Thu May 30 2013 - 11:38:10 MDT