Re: [errno=1101] in gsn_csm_pres_hgt

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Sat Apr 06 2013 - 15:15:25 MDT

Hi Erika,

Thanks for the files.

It turns out that the subset you are trying to plot is all missing. For debug purposes, I added this line:

printMinMax(u(10,:,{-50:50}),0)

and got this:

(0) Zonal Mean Wind: min=-32767 max=-32767

If you select a time index other than "10", like "24" (which is the first time step which didn't have all missing values), then your script works just fine.

--Mary

On Apr 4, 2013, at 12:50 AM, Erika Folova wrote:

> Hallo,
>
> Thanks Mary,
> I have already put the clean script and the data to ftp incoming:
>
> E40.mean.ncl
> E40.mean.nc
>
>
> Thanks a lot
> Erika
>
>
>
>
>
>
> On Thu, Apr 4, 2013 at 5:51 AM, Mary Haley <haley@ucar.edu> wrote:
> Erika,
>
> Unfortunately I still can't tell what the problem is. It seems like it's complaining about taking the log of your Y axis, but I think those values are okay based on your output.
>
> Can you provide a *clean* version of your script and the data?
>
> http://www.ncl.ucar.edu/report_bug.shtml#HowToFTP
>
> Thanks,
>
> --Mary
>
> On Apr 3, 2013, at 9:10 AM, Erika Folova wrote:
>
>> Hallo,
>>
>> Thanks for looking at my code.
>> Here, I put the request of printVarSummary:
>>
>> Variable: u
>> Type: float
>> Total Size: 3546048 bytes
>> 886512 values
>> Number of Dimensions: 3
>> Dimensions and sizes: [time | 528] x [level | 23] x [lat | 73]
>> Coordinates:
>> time: [508416..893376]
>> level: [1..1000]
>> lat: [90..-90]
>> Number Of Attributes: 10
>> anomaly_op_ncl : Anomalies from Annual Cycle: calcMonAnomTLLL: contributed.ncl
>> missing_value : -32767
>> long_name : U velocity
>> units : m s**-1
>> missing_value_original : -32767
>> _FillValue : -32767
>> _FillValue_original : -32767
>> wgt_runave_op_ncl : wgt_runave_n
>> average_op_ncl : dim_avg_n over dimension(s): lon
>> spatial_op_ncl : smth9; nine-pt smoother applied; p=0.25 q=-0.25 cyclic=False
>>
>>
>> (0) U velocity: min=-35.7064 max=32.2786
>>
>>
>> Variable: level (coordinate)
>> Type: integer
>> Total Size: 92 bytes
>> 23 values
>> Number of Dimensions: 1
>> Dimensions and sizes: [level | 23]
>> Coordinates:
>> Number Of Attributes: 2
>> long_name : pressure_level
>> units : millibars
>> (0) 1
>> (1) 2
>> (2) 3
>> (3) 5
>> (4) 7
>> (5) 10
>> (6) 20
>> (7) 30
>> (8) 50
>> (9) 70
>> (10) 100
>> (11) 150
>> (12) 200
>> (13) 250
>> (14) 300
>> (15) 400
>> (16) 500
>> (17) 600
>> (18) 700
>> (19) 775
>> (20) 850
>> (21) 925
>> (22) 1000
>>
>>
>> Thank you,
>> Erika Folova,
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>> On Wed, Apr 3, 2013 at 11:18 PM, Mary Haley <haley@ucar.edu> wrote:
>> Hi Erika,
>>
>> It looks like there might be something wrong with the actual data itself. I see a "printVarSummary" on "u", but this was the "u" before you deleted it and recalculated it.
>>
>> After this line:
>>
>> u= smth9_Wrap (uz, 0.25, -0.25, False)
>>
>> Add these lines and email the output back to ncl-talk:
>>
>> printVarSummary(u)
>> printMinMax(u,0)
>> print(u&level)
>>
>> If the values don't look correct, then you will need to examine the rest of your data more closely, including "filter" and "uz".
>>
>> --Mary
>>
>>
>> On Mar 31, 2013, at 4:31 AM, Erika Folova wrote:
>>
>>> Hallo,
>>>
>>> I have a very simple model output and attempted to visualize it using "gsn_csm_pres_hgt" of ncl.
>>> However, I always get this errors:
>>>
>>> warning:ContourPlotInitialize: no valid values in scalar field; ContourPlot not possible:[errno=1101]
>>> warning:LlSetValues: invalid range for Y Axis log mode; setting trYLog off
>>>
>>> Please find in the attachment my simple code, and here are the some details of the data run:
>>>
>>> u = short2flt(f->uwnd)
>>>
>>> Variable: u
>>> Type: float
>>> Total Size: 510630912 bytes
>>> 127657728 values
>>> Number of Dimensions: 4
>>> Dimensions and sizes: [time | 528] x [level | 23] x [lat | 73] x [lon | 144]
>>> Coordinates:
>>> time: [508416..893376]
>>> level: [1..1000]
>>> lat: [90..-90]
>>> lon: [ 0..357.5]
>>> Number Of Attributes: 6
>>> _FillValue_original : -32767
>>> _FillValue : -32767
>>> missing_value_original : -32767
>>> units : m s**-1
>>> long_name : U velocity
>>> missing_value : -32767
>>>
>>> ----
>>>
>>> Variable: utrend
>>> Type: float
>>> Total Size: 510630912 bytes
>>> 127657728 values
>>> Number of Dimensions: 4
>>> Dimensions and sizes: [time | 528] x [level | 23] x [lat | 73] x [lon | 144]
>>> Coordinates:
>>> time: [508416..893376]
>>> level: [1..1000]
>>> lat: [90..-90]
>>> lon: [ 0..357.5]
>>> Number Of Attributes: 7
>>> anomaly_op_ncl : Anomalies from Annual Cycle: calcMonAnomTLLL: contributed.ncl
>>> missing_value : -32767
>>> long_name : U velocity
>>> units : m s**-1
>>> missing_value_original : -32767
>>> _FillValue : -32767
>>> _FillValue_original : -32767
>>>
>>> ----
>>>
>>> Variable: filter
>>> Type: float
>>> Total Size: 510630912 bytes
>>> 127657728 values
>>> Number of Dimensions: 4
>>> Dimensions and sizes: [time | 528] x [level | 23] x [lat | 73] x [lon | 144]
>>> Coordinates:
>>> time: [508416..893376]
>>> level: [1..1000]
>>> lat: [90..-90]
>>> lon: [ 0..357.5]
>>> Number Of Attributes: 8
>>> anomaly_op_ncl : Anomalies from Annual Cycle: calcMonAnomTLLL: contributed.ncl
>>> missing_value : -32767
>>> long_name : U velocity
>>> units : m s**-1
>>> missing_value_original : -32767
>>> _FillValue : -32767
>>> _FillValue_original : -32767
>>> wgt_runave_op_ncl : wgt_runave_n
>>>
>>> ------
>>>
>>> Variable: uz
>>> Type: float
>>> Total Size: 3546048 bytes
>>> 886512 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [time | 528] x [level | 23] x [lat | 73]
>>> Coordinates:
>>> time: [508416..893376]
>>> level: [1..1000]
>>> lat: [90..-90]
>>> Number Of Attributes: 9
>>> anomaly_op_ncl : Anomalies from Annual Cycle: calcMonAnomTLLL: contributed.ncl
>>> missing_value : -32767
>>> long_name : U velocity
>>> units : m s**-1
>>> missing_value_original : -32767
>>> _FillValue : -32767
>>> _FillValue_original : -32767
>>> wgt_runave_op_ncl : wgt_runave_n
>>> average_op_ncl : dim_avg_n over dimension(s): lon
>>>
>>> ===
>>>
>>> Thanks for your help.
>>>
>>> Regards,
>>> -- Folova
>>>
>>>
>>>
>>> <test.ncl>_______________________________________________
>>> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sat Apr 6 15:15:38 2013

This archive was generated by hypermail 2.1.8 : Sun Apr 07 2013 - 21:13:30 MDT