[ncl-talk] Problems plotting CICE grid

[David Bailey]

Can you just plot aice? What does the field look like in ncview?

Dave

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

> 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.
>



-- 

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/bc4c79ff/attachment.html>