Re: station data - map configuration

I haven't looked at the script, but if you haven't already, try setting

     res@gsnAddCyclic = False

Let me know if this fixes the problem. I need to add this to the FAQ.

Thanks,

--Mary

On Jun 11, 2010, at 8:15 AM, carina.padilha@ufpel.edu.br wrote:

> Hello All,
>
> Thanks a lot Phillips, I include the masks using mpFillDrawOrder and it works.
> Now I am trying to generate a grid to fill others areas of Brasil with data. The
> gridded data is correct, because print(grided_data) returns coherent values. But
> when I try to use draw(map) of gridded_data it is giving me a error:
> warning:ScalarFieldSetValues: coordinate array sfXArray requires 152894907
> elements: defaulting
> warning:ScalarFieldSetValues: coordinate array sfYArray requires 152894916
> elements: defaulting
>
> I think this problem is related with:
>
> res@sfXArray = lat ;(defined field of the station data)
> res@sfYArray = lon ;(defined field of the station data)
>
> Do you know how can I fix it?
>
>
> Citando Adam Phillips <asphilli@ucar.edu>:
>
>> Hello,
>> I think you cam really close with your script. I took pieces of your
>> script, modified it, used my own data file, and came up with what I
>> think you want. (Plot and script attached.)
>>
>> The key was to set mpFillDrawOrder = "PostDraw", so the map fill masks
>> every region minus the Brazil states after the contour fill plot is
>> applied.
>>
>> Hope that helps, let the group know if that isn't what you were after...
>> Adam
>>
>> On 06/09/2010 11:45 AM, carina.padilha@ufpel.edu.br wrote:
>>>
>>>
>>> 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
>>
>> --
>> __________________________________________________
>> Adam Phillips
>> asphilli@ucar.edu
>> National Center for Atmospheric Research tel: (303) 497-1726
>> Climate and Global Dynamics Division fax: (303) 497-1333
>> P.O. Box 3000
>> Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli
>>
>
>
>
>
> -------------------------------------------------
> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jun 11 08:50:47 2010