[ncl-talk] Problems plotting CICE grid

[David Bailey]

Hi Anne,

I can try to help with this. What if you just plot "aice" as it is? Your
variables tice and var_dtrend are three dimensional. Are you just plotting
a single time slice? I think you are introducing some missing values with
your var_dtrend computation.

Dave

On Thu, Feb 15, 2018 at 9:30 AM, Anne <anne.seidenglanz at unive.it> wrote:

> Hello,
> I have the problem that I can perfectly plot *tice* from CICE output
> (based on *aice,* I have used one of the special ncl scripts for this
> purpose), but not anymore after I detrended it (dtrend; along with it using
> both *copy_VarMeta* and *copy_VarCoord* of the original in order to keep
> the metadata). The plot (attached) looks distorted.
>
>
> I think the error is in the *lat2d* and *lon2d*, when doing
> printVarSummary:
>
> Variable: *tice*
> Type: float
> Total Size: 846028800 bytes
>             211507200 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 144] x [y | 1020] x [x | 1440]
> Coordinates:
>             time: [727779..736325]
> Number Of Attributes: 8
>   long_name : ice area  (aggregate)
>   units : %
>   coordinates : TLON TLAT
>   _FillValue : 1e+30
>   missing_value : 1e+30
>   comment : none
>   cell_methods : time: mean
>   time_rep : averaged
>
> and printVarSummary od *var_dtrend:*
>
> Variable: *var_dtrend*
> Type: float
> Total Size: 846028800 bytes
>             211507200 values
> Number of Dimensions: 3
> Dimensions and sizes: [144] x [1020] x [1440]
> Coordinates:
> Number Of Attributes: 10
>   missing_value : 1e+30
> *  lat2d : <ARRAY of 1468800 elements>*
> *  lon2d : <ARRAY of 1468800 elements>*
>   long_name : ice area  (aggregate)
>   units : %
>   coordinates : TLON TLAT
>   comment : none
>   cell_methods : time: mean
>   time_rep : averaged
>   _FillValue : 1e+30
>
> lat2d and lon2d are extra listed here again, and I don't know why. I have
> also tried to add
> var_dtrend at lon2d = tlon2d
> var_dtrend at lat2d = tlat2d
>
> but this didn't solve the issue. I also attached the script.
>
> Is the simplest thing just to interpolate it anyway?
>
> Thanks, any help appreciated!
> Anne
>
> Nota automatica aggiunta dal sistema di posta.
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>


-- 

David A Bailey
 email: dbailey_at_ucar.edu
National Center for Atmospheric Research   phone: 303-497-1737
PO Box 3000
fax  : 303-497-1700
Boulder, CO 80307-3000
                   http://www.cgd.ucar.edu/staff/dbailey
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180215/769cc3a0/attachment.html>