[ncl-talk] gsn_add_polygon function

[Guilherme Martins]

Hi Alan,

Thanks a lot!

It works!

Guilherme.

*Guilherme Martins*
https://sites.google.com/site/jgmsantos


2015-03-01 15:50 GMT-03:00 Alan Brammer <abrammer at albany.edu>:

> oh if your data only has 12 time dimensions then you can possibly switch
> out y1&time at lines 77/78 with an array of 1-12.
>
> The below would  rely on the nmes always being 12, or I'm sure you can
> make it more robust, based on the needs.
>
> dummytime = ispan(1,12,1)
> do k=0,nmes-1
>    yp(k)              = y1(k) + dx
>    xp(k)              = dummytime(k)
>    xp(sigma*nmes-1-k) = dummytime(k)
>    yp(sigma*nmes-1-k) = y1(k) - dx
> end do
>
>
> On Sun, Mar 1, 2015 at 1:43 PM, Guilherme Martins <jgmsantos at gmail.com>
> wrote:
>
>> Hi Alan,
>>
>> I'll take a look in your sugestion. Thanks.
>>
>> The informations about my file using the ncdump
>>
>> netcdf fluxos.ERAI {
>> dimensions:
>> lon = 1 ;
>> lat = 1 ;
>> time = UNLIMITED ; // (12 currently)
>> variables:
>> double lon(lon) ;
>> lon:standard_name = "longitude" ;
>> lon:long_name = "longitude" ;
>> lon:units = "degrees_east" ;
>> lon:axis = "X" ;
>> double lat(lat) ;
>> lat:standard_name = "latitude" ;
>> lat:long_name = "latitude" ;
>> lat:units = "degrees_north" ;
>> lat:axis = "Y" ;
>> double time(time) ;
>> time:standard_name = "time" ;
>> time:units = "days since 1979-01-01 12:00:00" ;
>> time:calendar = "standard" ;
>> float fnorte(time, lat, lon) ;
>> fnorte:code = 3 ;
>> float fsul(time, lat, lon) ;
>> fsul:code = 5 ;
>> float fleste(time, lat, lon) ;
>> fleste:code = 2 ;
>> float foeste(time, lat, lon) ;
>> foeste:code = 4 ;
>> float fbalanco(time, lat, lon) ;
>> fbalanco:code = 1 ;
>>
>> // global attributes:
>> :CDI = "Climate Data Interface version 1.6.5 (
>> http://code.zmaw.de/projects/cdi)" ;
>> :Conventions = "CF-1.4" ;
>> :history = "Sun Mar 01 01:16:49 2015: cdo -O -s merge norte.nc sul.nc
>> leste.nc oeste.nc balanco.nc fluxos.ERAI.nc\n",
>> "Sun Mar 01 01:16:44 2015: cdo -O -s -r -ymonmean -monmean
>> -setcalendar,standard -settaxis,1979-01-01,12:00:00,1day -setcode,1
>> -chname,var1,fbalanco balanco.tmp.nc balanco.nc\n",
>> "Sun Mar 01 01:16:44 2015: cdo -O -s -f nc input,r1x1 balanco.tmp.nc" ;
>> :CDO = "Climate Data Operators version 1.6.5 (
>> http://code.zmaw.de/projects/cdo)" ;
>> }
>>
>> Guilherme.
>>
>> *Guilherme Martins*
>> https://sites.google.com/site/jgmsantos
>>
>>
>> 2015-03-01 15:25 GMT-03:00 Alan Brammer <abrammer at albany.edu>:
>>
>> Hard to know without the file information.
>>> It looks like your plotting the data on ax axis of  1-12.  The x points
>>> of your polygon therefore need to be with respect to this 1-12 axis.
>>> Whereas you're defining them based on the time coordinate of your data
>>> which I presume is not 1-12 currently.
>>>
>>> Either scale the x coordinates of the polygon to match the 1-12 axis or
>>> change the plot to be with respect to y1&time.
>>>
>>> Throw a print(xp)  between line 81 and the end and if the values are not
>>> between 1-12 then the polygon is being plotted but somewhere else not
>>> within the plot bounds.
>>>
>>> Hopefully that makes sense.
>>>
>>>
>>> ~Alan.
>>>
>>>
>>> On Sun, Mar 1, 2015 at 12:08 PM, Guilherme Martins <jgmsantos at gmail.com>
>>> wrote:
>>>
>>>> Hi NCL users!
>>>>
>>>> I'm trying to insert shaded area in my dataset observed using the
>>>> gsn_add_polygon function, but I'm not getting. I don't get errors and the
>>>> shaded area don't appear.
>>>>
>>>> My files are in attached.
>>>>
>>>> I'm using ncl 6.2.1.
>>>>
>>>> Any help is very appreciated.
>>>>
>>>> *Guilherme Martins*
>>>> https://sites.google.com/site/jgmsantos
>>>>
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150301/ddf03f0e/attachment.html