Re: [ferret_users] Inf, NAN and **** output from rect_to_curv function

[Stephen Guimond <sguimond@xxxxxxx>]

Hi Ansley,

My file was in binary format read in with the /format=stream option in ferret so this must explain the discrepancy in data type (r8 vs. r4) and thus the problem with the rect_to_curv function.  Does ferret only convert to single precision when the data file is NetCDF?

Steve

----- Original Message -----
From: Ansley Manke <Ansley.B.Manke@xxxxxxxx>
Date: Tuesday, September 2, 2008 10:19 am
Subject: Re: [ferret_users] Inf, NAN and **** output from rect_to_curv function
To: guimond@xxxxxxxxxxxxx
Cc: oar.pmel.ferret_users@xxxxxxxx

> Hi Steve,
> Thanks for reporting back. If the data is coming in from netCDF 
> files, 
> Ferret should convert to single precision on input (except for 
> coordinate axis values, that is 1-D axes, not curvilinear 
> coordinate 
> data which is read as variable fields and is in single precision ). 
> So 
> I'm not sure what may have been happening to make this happen, but 
> it's 
> good to know you're off and running.
> 
> Ansley
> 
> Stephen Guimond wrote:
> > I believe I found an answer to my own question.  My input data to 
> the 
> > rect_to_curv function was in double precision (read in with type 
> r8) 
> > and the output coordinates were in single precision therefore the 
> > output must have been confused with the type and returned inf, 
> NAN and 
> > **** values.  All I did was output the data to a file and then 
> read it 
> > back in with single precision and the call to rect_to_curv worked.
> >
> > ----- Original Message -----
> > From: Stephen Guimond <sguimond@xxxxxxx>
> > Date: Monday, September 1, 2008 4:31 am
> > Subject: [ferret_users] Inf, NAN and **** output from 
> rect_to_curv 
> > function
> > To: oar.pmel.ferret_users@xxxxxxxx
> >
> > > Hello,
> > >
> > > I am trying to grid a rectilinear grid onto a curvilinear grid but
> > > just get a bunch of inf, NAN and *** for my output.  The grid 
> of my
> > > input variable (called ntest):
> > >
> > > name       axis              # pts   start                end
> > > NX        LONGITUDE          200mi   126.12W(233.89)     
> > > 113.89W(246.11) NY        LATITUDE           200 i   14.253N    
>   
> > >      25.746N
> > > ZAX       Z (KM)              71 i   0.035583             15.116
> > >
> > > The grid of my destination curvilinear coordinates (same 2-D
> > > dimensions for lat and lon):
> > >
> > > name       axis              # pts   start                end
> > > X1        X                  180 r   1                    180
> > > Y1        Y                  200 r   1                    200
> > >
> > > The values of the 2-D lat and lon can be seen from the attached
> > > plot, which shows the general rectilinear grid and an overlay of
> > > the curvilinear grid.  Here is my call to the function:
> > >
> > > let testing = rect_to_curv(ntest[d=1], clon, clat, 2)
> > >
> > > Again, this gives me inf, NAN and *** for output.  Not sure if the
> > > destination coords need to be defined on a lat/lon axis or if they
> > > can be just a listing of lat/lon in degrees as I have it here. 
> > > This wasn't clear from the documentation.
> > >
> > > Thanks,
> > > Steve
> > > =======================================================
> > > Stephen R. Guimond
> > > Graduate Research Assistant
> > > Florida State University
> > > Center for Ocean-Atmospheric Prediction Studies (COAPS)
> > > =======================================================
> > >
> > >
> > >
> > >
> >
> > =======================================================
> > Stephen R. Guimond
> > Graduate Research Assistant
> > Florida State University
> > Center for Ocean-Atmospheric Prediction Studies (COAPS)
> > =======================================================
> 
> 
> 

=======================================================
Stephen R. Guimond
Graduate Research Assistant
Florida State University
Center for Ocean-Atmospheric Prediction Studies (COAPS)
=======================================================