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