Re: Contour folding and attribute errors

Hi Justin,

Are the errors coming from the printVarSummary? I can't really tell
from looking at the output you included below. This error message
occurs if you have an existing attribute and try to reassign it to a
value that's a type that can't be coerced. For example, the code:

res = True
res_at_mpMinLatF = 10
res_at_mpMinLatF = 10.

Will cause the same fatal attribute error. Can you personally email me
your script and indicate where you are seeing the fatal errors? It
would be even better if you could also provide a data file.

Thanks,

--Mary

On May 20, 2008, at 3:50 PM, Justin Glisan wrote:

> Hello group,
>
> As a new person to NCL, I've run into a problem plotting a polar
> stereographic plot. Attached, you will find an output image with a
> distortion
> from the Prime Meridian to the NPole. It almost looks as if the data
> is folded
> in these place.
>
> Below is header information from the two input files (MM5 and
> ERA40) and then printVarSummary from the variables in play. Can anyone
> offer any
> suggestions? For one of the variable I'm getting multiple attribute
> errors and
> I'm not sure how to address them.
>
> Thank you,
>
> Justin
>
> __________________________________________________________
> netcdf 8007M (PAN-ARCTIC MM5)
> dimensions:
> icross = 94 ;
> jcross = 127 ;
> time = UNLIMITED ; // (1 currently)
> variables:
> double latitcrs(icross, jcross) ;
> latitcrs:long_name = "Latitude" ;
> latitcrs:standard_name = "Lat" ;
> latitcrs:units = "degrees_north" ;
> double longicrs(icross, jcross) ;
> longicrs:long_name = "Longitude" ;
> longicrs:standard_name = "Lon" ;
> longicrs:units = "degrees_east" ;
> double psealvlc(time, icross, jcross) ;
> psealvlc:long_name = "Sea Level Pressure" ;
> psealvlc:standard_name = "MSLP" ;
> psealvlc:units = "Pa" ;
> double time(time) ;
> time:long_name = "time" ;
> time:standard_name = "time" ;
> time:calendar = "gregorian" ;
> time:units = "hours since 1979-09-01 00:00:0.0"
>
> __________________________________________________________
> netcdf 8007E (ERA40)
> dimensions:
> longitude = 144 ;
> latitude = 73 ;
> time = UNLIMITED ; // (1 currently)
> variables:
> float longitude(longitude) ;
> longitude:units = "degrees_east" ;
> longitude:long_name = "longitude" ;
> float latitude(latitude) ;
> latitude:units = "degrees_north" ;
> latitude:long_name = "latitude" ;
> int time(time) ;
> time:units = "hours since 1900-01-01 00:00:0.0" ;
> time:long_name = "time" ;
> double msl(time, latitude, longitude) ;
> msl:_FillValue = -32767s ;
> msl:missing_value = -32767. ;
> msl:units = "Pa" ;
> msl:long_name = "Mean sea level pressure" ;
> __________________________________________________________
> PrintVarSummay:
>
>
> Variable: MSLP_MM5
> Type: double
> Total Size: 95504 bytes
> 11938 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 1] x [icross | 94] x [jcross | 127]
> Coordinates:
> time: [7666.5..7666.5]
> Number Of Attributes: 3
> units : Pa
> standard_name : MSLP
> long_name : Sea Level Pressure
> fatal:Attribute assignment type mismatch
> fatal:Attribute assignment type mismatch
>
> __________________________________________________________
> Variable: MSLP_ERA
> Type: double
> Total Size: 84096 bytes
> 10512 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 1] x [latitude | 73] x [longitude |
> 144]
> Coordinates:
> time: [705993..705993]
> latitude: [90..-90]
> longitude: [ 0..357.5]
> Number Of Attributes: 4
> long_name : Mean sea level pressure
> units : Pa
> missing_value : -32767
> _FillValue : -32767
> fatal:Attribute assignment type mismatch
>
> __________________________________________________________
> Variable: MSLP_ERA_Flip (ERA data flipped from dec. lat. to asc. lat.)
> Type: double
> Total Size: 84096 bytes
> 10512 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 1] x [latitude | 73] x [longitude |
> 144]
> Coordinates:
> time: [705993..705993]
> latitude: [-90..90]
> longitude: [ 0..357.5]
> Number Of Attributes: 4
> _FillValue : -32767
> missing_value : -32767
> units : Pa
> long_name : Mean sea level pressure
> fatal:Attribute assignment type mismatch
> fatal:Attribute assignment type mismatch
> fatal:Attribute assignment type mismatch
> fatal:Attribute assignment type mismatch
> fatal:Attribute assignment type mismatch
>
> __________________________________________________________
> Variable: newgrid (REGRID to PanARC MM5 Grid)
> Type: double
> Total Size: 95504 bytes
> 11938 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 1] x [lat | 94] x [lon | 127]
> Coordinates:
> time: [705993..705993]
> lat: [ -90..89.99999463558197]
> lon: [ 0..357.1653642654419]
> Number Of Attributes: 4
> long_name : Mean sea level pressure
> units : Pa
> missing_value : -32767
> _FillValue : -32767
>
> __________________________________________________________
> Variable: diff (Difference between MM5 and ERA40)
> Type: double
> Total Size: 95504 bytes
> 11938 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 1] x [lat | 94] x [lon | 127]
> Coordinates:
> time: [705993..705993]
> lat: [ -90..89.99999463558197]
> lon: [ 0..357.1653642654419]
> Number Of Attributes: 4
> _FillValue : -32767
> missing_value : -32767
> units : Pa
>
>
> ___________________________
> Justin Glisan
> Doctoral Research Assistant
> Geological and Atmospheric Sciences
> Iowa State University
>
> <bias_test.png>_______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue May 20 2008 - 20:30:55 MDT