station data - map configuration

Hi

I am using station data from some cities of south of Brazil. They are stored
into a .csv file, where the first field is longitude, the second field is
latitude and the next field is the variable. I made a script that fill the area
from the point of the data, that means, it is filled in part of south of Brazil
but also it is filled over part of the other countries. I want to see a figure
that fill south of Brazil and don't plot the variable over the other countries
of South America, nor over the Atlantic Ocean. In fact, I want to fill just the
south state of Brazil. I want something like this example:
http://www.ncl.ucar.edu/Applications/Images/mask_8_lg.png. Is there a way to
make this using ncl? Perhaps converting
this file to grid and using the grid file? The figure generated by the script is
attached.
Above, the script that I made:

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"

begin

  dados = "dados_2010_04.csv"

  lines = asciiread(dados,-1,"string")

campo1 = str_get_field(lines,1,";")

lon = stringtofloat(campo1)

campo2 = str_get_field(lines,2,";")

lat = stringtofloat(campo2)

campo3 = str_get_field(lines,3,";")

dec1 = stringtofloat(campo3)

wks = gsn_open_wks("ps","figura1")

gsn_define_colormap(wks,"BlRe")

res = True

res@gsnFrame = False

res@gsnMaximize = True

res@gsnPaperOrientation = "portrait"

res@cnFillOn = True

res@cnFillColors = (/70,66,60,43,40,36,31,28,17,14/)

res@cnLevelSelectionMode = "ExplicitLevels"

res@cnLevels = (/20,30,40,50,60,70,80,90,100,150/)

;;;;;;;;;;;;Title

; tabela de acentos

atilde = "a:H-13V2F35:D:FV-2H3:"

cedil = "c:H-13F35:K:FH2:"

ecirc = "e:H-13V2F35:C:FV-2H3:"

degree = ":S:o:N:"

res@tiMainString = "Precipita"+cedil+atilde+"o 1"+degree+" dec"+ecirc+"ndio de
abril de 2010"

res@sfXArray = lon

res@sfYArray = lat

res@mpProjection = "LambertConformal"

res@mpLambertParallel1F = -30

res@mpLambertParallel2F = -20

res@mpLambertMeridianF = -45

res@mpLimitMode = "LatLon"

res@mpMinLatF = -34

res@mpMaxLatF = -27

res@mpMinLonF = -58

res@mpMaxLonF = -49

res@mpDataSetName = "Earth..4"

res@mpOutlineOn = True

res@mpDataBaseVersion = "MediumRes"

res@mpOutlineSpecifiers = (/"Brazil:states"/)

res@mpFillAreaSpecifiers = (/"Land", "Brazil:states","Water"/)

res@mpSpecifiedFillColors = (/0,-1,49/)

res@mpOutlineBoundarySets = "AllBoundaries"

map = gsn_csm_contour_map(wks,dec1,res)

draw(map)

end

-------------------------------------------------
This mail sent through IMP: http://horde.org/imp/

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk