Re: Defining sic&lat and sic&lon

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Fri May 17 2013 - 15:44:29 MDT

Did you look at the variable "sic" ?

    printVarSummary(sic)

On 5/17/13 2:32 PM, 진현근 wrote:
> 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
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri May 17 15:44:38 2013

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