[ncl-talk] Problem with CFS initial condition data plot in NCL

[Mary Haley]

I think you are close.  The lat/lon array is actually 410 x 720, and not
720 x 410, so you want something like this:

nlon = getfilevardimsizes (fin,"xaxis_1")
nlat = getfilevardimsizes (fin,"yaxis_1")

; "Conform" the 1D lat/lon arrays to be 2D arrays of nlat x nlon in size
lat2d = conform_dims((/nlat,nlon/),fin->yaxis_1,0)
lon2d = conform_dims((/nlat,nlon/),fin->xaxis_1,1)

gen_res = True
gen_res at sfXArray = lon2d
gen_res at sfYArray = lat2d

On Tue, Aug 4, 2015 at 8:18 AM, Vimal Koul <koulvimal18 at gmail.com> wrote:

> Hi Nicolas,
>
> Thanks for your help.
>
> But now, the following error messages are seen:
>
> warning:ScalarFieldSetValues: 2d coordinate array sfXArray has an
> incorrect dimension size: defaulting sfXArray
> warning:ScalarFieldSetValues: 2d coordinate array sfYArray has an
> incorrect dimension size: defaulting sfYArray
> warning:ContourPlotDraw: out of range coordinates encountered; standard
> AreaFill rendering method may be unreliable;
>  consider setting the resource trGridType to "TriangularMesh" if
> coordinates contain missing values
> fatal:ContourPlotDraw: Workspace reallocation would exceed maximum size
> 100000000
> fatal:ContourPlotDraw: draw error
> fatal:ContourPlotDraw: draw error
> fatal:PlotManagerDraw: error in plot draw
> fatal:_NhlPlotManagerDraw: Draw error
>
> Please note that I don't have the grid resolution information except nlat,
> nlon and depth levels.
>
> Modified script is attached.
>
>
> Regards,
> Vimal
>
>
> Sincerely
> Vimal Koul
>
> Graduate Student
> Department of Atmospheric and Space Science
> Savitribai Phule Pune University
> India
>
>
>
> On 4 August 2015 at 16:41, Nicolas GASNIER <ngprod41 at gmail.com> wrote:
>
>> Hi,
>>
>> With the xaxis_1 and yaxis_1 variables you should try to build a lat and
>> lon 2d arrays matching your data grid dimensions. Then use these 2d
>> arrays as sfXArray and sfYArray plotting ressource.
>>
>> Something like this (untested) :
>>
>> xaxis_1 = fin->xaxis_1
>> yaxis_1 = fin->yaxis_1
>>
>> lon2d = new((/720, 410/), float)
>> do i=0,719
>>      lon2d(i, :) = xaxis_1
>> end do
>>
>> lat2d = new((/720, 410/), float)
>> do i=0,409
>>      lat2d(:, i) = yaxis_1
>> end do
>>
>>
>> gen_res = True
>> gen_res at sfXArray = lon2d
>> gen_res at sfYArray = lat2d
>>
>> Then pass that ressource to your plotting function.
>>
>> Hope that helps.
>>
>> Nicolas Gasnier
>>
>>
>> Le 04/08/2015 08:23, ncl-talk-request at ucar.edu a écrit :
>> > Date: Tue, 4 Aug 2015 11:53:07 +0530
>> > From: Vimal Koul<koulvimal18 at gmail.com>
>> > Subject: [ncl-talk] Problem with CFS initial condition data plot in
>> >       NCL
>> > To:ncl-talk at ucar.edu
>> > Message-ID:
>> >       <
>> CAJVe0MCbueZr-1EZHk6eOD1VULML0XEOMaGchvsM5rczLjjFSQ at mail.gmail.com>
>> > Content-Type: text/plain; charset="utf-8"
>> >
>> > Dear Users,
>> >
>> > I am trying to plot an ocean data file which is a part of CFS High
>> > Resolution Initial Conditions (http://nomads.ncdc.noaa.gov/data.php).
>> >
>> > File: ocnanl.gdas.1985050100.ocean_temp_salt.res.nc (uploaded to ncl
>> ftp
>> > link) or download from:https://copy.com/7bVgOvDfQKDjcg07
>> >
>> > Unfortunately, printVarSummary shows that the variable "temp" I'm
>> trying to
>> > plot does not have coordinate information, however, nlat, nlon are
>> > provided. So I used fspan (and tried other functions also) to attach
>> > coordinates, but the resulting plot is messed up. The plot and ncl
>> script
>> > are attached.
>> >
>> > Can you please help me in defining the lat long coordinate information
>> for
>> > this type of data?
>> >
>> > I am completely new to NCL, infact just started last week.
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> 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/20150805/fe5dae78/attachment.html