Re: isosurface plot problem

Dennis

Changes made, but the plots continue to not show much. Here's the
printVarSummary output:

ncl 0> load "/root/NCL/createIsoTest.ncl"
warning:int2p_n: 3267792 input array(s) contained all missing data. No
interpolation performed on these arrays
(0) ===

Variable: modelTemp
Type: float
Total Size: 911806632 bytes
             227951658 values
Number of Dimensions: 4
Dimensions and sizes: [time | 99] x [zlev | 22] x [latitude | 261] x
[longitude | 401]
Coordinates:
             time: [ 0..98]
             zlev: [ 0..3500]
             latitude: [18..31]
             longitude: [-98..-78]
Number Of Attributes: 5
   name : sea_water_potential_temperature
   long_name : potential temperature
   units : deg C
   missing_value : -9999
   _FillValue : -9999
(0) potential temperature: min=2.993 max=29.832
(0) ===

Variable: depth
Type: float
Total Size: 911806632 bytes
             227951658 values
Number of Dimensions: 4
Dimensions and sizes: [time | 99] x [zlev | 22] x [latitude | 261] x
[longitude | 401]
Coordinates:
             time: [ 0..98]
             zlev: [ 0..3500]
             latitude: [18..31]
             longitude: [-98..-78]
Number Of Attributes: 1
   units : m
(0) min=0 max=3500
(0) ===

Variable: tlev
Type: float
Total Size: 4 bytes
             1 values
Number of Dimensions: 1
Dimensions and sizes: [tlev | 1]
Coordinates:
Number Of Attributes: 1
   units : deg C
(0) min=20 max=20
(0) ===

Resultant plots are in http://mesonet.tamu.edu/NCL

Thanks, Gerry

On 11/21/2011 02:58 PM, Dennis Shea wrote:
> Gerry,
>
> It is best to send to ncl-talk.
>
> Your code looks good. After
>
> isot = int2p_n_Wrap(mmodelTempodelTemp, depth, tlev,0,1
>
> add the following to show that things [eg, units] are as u expect.
>
> print("===")
> printVarSummary(modelTemp)
> printMinMax(modelTemp,0)
> print("===")
> printVarSummary(depth)
> printMinMax(depth,0)
> print("===")
> printVarSummary(tlev)
> printMinMax(tlev,0)
> print("===")
>
> Comment out
>
> cnRes@cnLevelSelectionMode = "ManualLevels" ; Manually define levels
> cnRes@cnMinLevelValF = 40 ; Minimum depth (meters)
> cnRes@cnMaxLevelValF = 150 ; Maximum depth (meters)
> ;cnRes@cnLevelSpacingF = 5.
>
>
> and let ncl choose the limits
>
> On 11/21/2011 10:16 AM, Gerald Creager wrote:
>> Look at the pngs in http://mesonet.tamu.edu/NCL for an example of what
>> I'm getting...
>>
>> gc
>>
>> 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