Re: isosurface plot problem

Dennis

Thanks for the clarification. That is what I'd thought, but I've looked
a this for quite awhile and I'm confused by the results. On only a
couple of pot outputs do I get anything approximating coverage of the
entire study area. I'd have thought that, with
cnRes@cnLevelSelectionMode = "AutomaticLevels"
I should see a general coverage. Is it possible that the granularity of
my model data is insufficient to provide adequate interpolation?

Thanks, Gerry

On 11/19/2011 06:09 PM, Dennis Shea wrote:
> The *_n functions were developed so that 'dimension reordering'
> would not be necessary. For
> modelTemp(time,zlev,latitude,longitude)
> 0 1 2 3
> The correct useage is
>
> isot = int2p_n_Wrap(modelTemp, depth, tlev, 0,1)
>
> D
> On 11/19/11 1:13 PM, donna Cote wrote:
>> Gerry, I'm know I'm guessing on what you are trying to do.
>>
>> I'm reading this (see quote or web page:
>> http://www.ncl.ucar.edu/Document/Functions/Contributed/int2p_n_Wrap.shtml),
>> and I don't think "depth"'s right-most dimension is the level.
>>
>> "modelTemp" has the coordinates ordered as:
>> [time | 99] x [zlev | 22] x [latitude | 261] x [longitude | 401]
>>
>> and this function wants (I think)
>> [time | 99] x [latitude | 261] x [longitude | 401] x [zlev | 22]
>>
>> The same note applies to the depth and tlev vars.
>>
>> isot = int2p_n_Wrap(modelTemp, depth, tlev, 0,1)
>> ^pin^ ^xin^ ^pout^
>>
>> It looks like you have to reorder the dimensions, which requires that
>> all the dimensions be named. Ref:
>> http://www.ncl.ucar.edu/Document/Language/reorder.shtml
>>
>> -Donna
>>
>> On 11/17/11 4:54 PM, Gerald Creager wrote:
>> > Dennis, Mary,
>> >
>> > I've uploaded createIsoTest.ncl and an example PNG (Model1-20C.png) to
>> > your incoming directory. My model .nc file was apparently too big but
>> > can be snagged from
>> > http://mesonet.tamu.edu/iasnfs_forecast_2011111112_2012021712.nc
>> >
>> > I was expecting to see an isosurface over the whole Gulf of Mexico
>> > region. This would be for the 20C isotherm, if I'm not completely
>> > delusional.
>> >
>> > Any help or suggestions you could offer would be appreciated.
>> >
>> > Thanks, Gerry
>>
>>
>> On 11/19/11 1:00 PM, ncl-talk-request@ucar.edu wrote:
>>> Date: Fri, 18 Nov 2011 17:10:16 -0600
>>> From: Gerald Creager<gerry.creager@tamu.edu>
>>> Subject: Re: isosurface plot problem
>>> To: Rick Brownrigg<brownrig@ucar.edu>
>>> Cc:"ncl-talk@ucar.edu"<ncl-talk@ucar.edu>
>>>
>>> Howdy, Rick
>>>
>>> On 11/18/2011 04:42 PM, Rick Brownrigg wrote:
>>>>> Hi Gerry,
>>>>>
>>>>> I don't have answers for you, but that is a most curious result. The colors shown in the plot you attached aren't in the palette you specified (gsn_define_colormap(wks,"BlAqGrYeOrRe") ). In fact, they are from the "default" palette, and appear to start with index 4 (as per cnRes@gsnSpreadColorStart = 4), and are for the most part every other color from the palette, with a few exceptions. Its as though the palette is not getting set to what you specified. Are you getting any errors/warnings about that?
>>> No, I'm getting no errors on the pallete. I'm getting a warning on
>>> int2p_n, but that's all.
>>>
>>>>> Also, please excuse my ignorance of the science, but the range of depth intervals is from 40 to 150m; is it possible parts of the 20C isotherm are at depths outside that range?
>>> They shouldn't be, unless I've done something wrong. Depth should be:
>>> 0,
>>> 20,40,60,80,100,150,200,250,300,400,500,700,800,900,1000,1500,2000,2500,3000,3500
>>> or, 22 levels that are NOT evenly spaced. I anticipate the 20C isotherm
>>> would lie in the forst 4-6 levels but for really deep loop currents, it
>>> can be a bit deeper...
>>>
>>> Please help me out (I've stared at this 'til I'm cross-eyed): Where am I
>>> limiting the levels?
>>>
>>> gerry
>>>
>>>>> FWIW...
>>>>> Rick
>>>>>
>>>>> On Nov 18, 2011, at 11:06 AM, Gerald Creager wrote:
>>>>>
>>>>>>> I've uploaded createIsoTest.ncl and an example PNG (Model1-20C.png) to your incoming directory. My model .nc file was apparently too big but can be snagged from
>>>>>>> http://mesonet.tamu.edu/iasnfs_forecast_2011111112_2012021712.nc
>>>>>>>
>>>>>>> I was expecting to see an isosurface over the whole Gulf of Mexico region. This would be for the 20C isotherm, if I'm not completely delusional.
>>>>>>>
>>>>>>> Any help or suggestions you could offer would be appreciated.
>>>>>>>
>>>>>>> Thanks, Gerry
>> _______________________________________________
>> 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

-- 
Gerry Creager -- gerry.creager@tamu.edu
Texas Mesonet -- AATLT, Texas A&M University
Cell: 979.229.5301 Office: 979.458.4020 FAX: 979.862.3983
Office: 7607 Eastmark Drive, Suite 112, College Station, TX 77840
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk