[ncl-talk] Problems plotting CICE grid

[Anne]

Hi Dave,
yes, I've been plotting the 1st time step only for test purposes. Maybe I
phrased it unluckily, I call it tice because in the end I am plotting on
the t-grid, but the variable I am reading in is aice. How do you think
could I avoid reading in missing values then?
Thanks for getting back to me,
Anne

On 15 February 2018 at 17:50, David Bailey <dbailey at ucar.edu> wrote:

> 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
> <(303)%20497-1737>
> PO Box 3000
> fax  : 303-497-1700 <(303)%20497-1700>
> Boulder, CO 80307-3000
>                    http://www.cgd.ucar.edu/staff/dbailey
>

-- 
Nota automatica aggiunta dal sistema di posta.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180215/ee7dd03a/attachment.html>