Contour folding and attribute errors

From: Justin Glisan <glisanj_at_nyahnyahspammersnyahnyah>
Date: Tue, 20 May 2008 16:50:28 -0500 (CDT)

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

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

bias_test.png
Received on Tue May 20 2008 - 15:50:28 MDT

This archive was generated by hypermail 2.2.0 : Wed May 21 2008 - 17:47:28 MDT