[ncl-talk] Latitude units assignment

[Sourav Chatterjee]

Hi,


My original variable is as below:


    float SLA(Time, Longitude, Latitude) ;
        SLA:_FillValue = 9.96921e+36f ;
        SLA:scale_factor = 1.f ;
        SLA:coordinates = "Sea Level Anomaly Estimate" ;
        SLA:add_offset = 0.f ;
        SLA:long_name = "Sea Level Anomaly Estimate" ;
        SLA:standard_name = "Sea Level Anomaly Estimate" ;
        SLA:units = "m" ;

As the lat lon dimensions has no units I was not getting the plot properly.


I tried this:


sst = f->SLA
printVarSummary(sst)

sst!1 = "longitude"
sst!2 = "latitude"

sst&latitude at units = "degrees_north"
sst&longitude at units ="degrees_east"

printVarSummary(sst)



But still getting the error,


(0)    check_for_y_lat_coord: Warning: Data either does not contain a valid latitude coordinate array or doesn't contain one at all.
(0)    A valid latitude coordinate array should have a 'units' attribute equal to one of the following values:

Why is the units attribute is not getting assigned?


Thanks

Sourav
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171219/4b08c7b3/attachment.html>