[ncl-talk] coordinate problem

Adam Phillips asphilli at ucar.edu
Tue Sep 11 09:53:34 MDT 2018


Hi Jalal,
That you for sending a clean easy-to-read script. It looks like your
mean_prc array is dimensioned  lon x lat. gsn_csm_contour_map expects the
data to be arranged as lat x lon. Change this:
plot = gsn_csm_contour_map(wks, mean_prc,res)
to this:
plot = gsn_csm_contour_map(wks, mean_prc(lat|:,lon|:),res)

Also you might want to check that when you used coordinate subscripting
that you are subscripting the correct dimensions:
 prc = f->pre(:,{-5:30},{40:100})  ; (line 5)

Hope that helps. If you have further questions please let ncl-talk know.
Adam



On Mon, Sep 10, 2018 at 8:03 PM Md. Jalal Uddin <dmjalal90 at gmail.com> wrote:

> Hi sir,
> Please find the results that you mentioned. The clean script is attached
> herewith.
>
> Variable: prc
> Type: double
> Total Size: 316780800 bytes
>             39597600 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 7071] x [longitude | 140] x [latitude | 40]
> Coordinates:
>             time: [ 730..6939]
>             longitude: [-4.875..29.875]
>             latitude: [40.125..49.875]
> Number Of Attributes: 4
>   long_name : precipitation
>   units : mm
>   _FillValue : -nan
>   missing_value : -nan
>
> Variable: latitude (coordinate)
> Type: double
> Total Size: 320 bytes
>             40 values
> Number of Dimensions: 1
> Dimensions and sizes: [latitude | 40]
> Coordinates:
> Number Of Attributes: 4
>   standard_name : latitude
>   long_name : latitude
>   units : degrees_north
>   axis : Y
>
> Variable: longitude (coordinate)
> Type: double
> Total Size: 1120 bytes
>             140 values
> Number of Dimensions: 1
> Dimensions and sizes: [longitude | 140]
> Coordinates:
> Number Of Attributes: 4
>   standard_name : longitude
>   long_name : longitude
>   units : degrees_east
>   axis : X
>
>
> On Tue, Sep 11, 2018 at 1:23 AM Adam Phillips <asphilli at ucar.edu> wrote:
>
>> Hi Jalal,
>> Can you run this part of your script and send the results to ncl-talk?
>> Also, please include a clean, complete script so we can see how you are
>> creating the plot.
>> f = addfile("non_cy.nc", "r")
>>  prc = f->pre(:,{-5:30},{40:100})
>>  prc at long_name = "precipitation"
>>  prc at units = "mm"
>>  prc&longitude at units = "degrees_east"
>>  prc&latitude at units = "degrees_north"
>>  printVarSummary(prc)
>>  printVarSummary(prc&latitude)
>>  printVarSummary(prc&longitude)
>>
>> As always please respond to ncl-talk.
>> Adam
>>
>>
>>
>>
>>
>>
>>
>> On Sun, Sep 9, 2018 at 10:11 PM Md. Jalal Uddin <dmjalal90 at gmail.com>
>> wrote:
>>
>>> Hi all,
>>> I am getting the following errors:
>>> (0) check_for_y_lat_coord: Warning: Data either does not contain
>>> (0) a valid latitude coordinate array or doesn't contain one at all.
>>> (0) A valid latitude coordinate array should have a 'units'
>>> (0) attribute equal to one of the following values:
>>> (0)     'degrees_north' 'degrees-north' 'degree_north' 'degrees north'
>>> 'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
>>> (0) check_for_lon_coord: Warning: Data either does not contain
>>> (0) a valid longitude coordinate array or doesn't contain one at all.
>>> (0) A valid longitude coordinate array should have a 'units'
>>> (0) attribute equal to one of the following values:
>>> (0)     'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
>>> 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
>>>
>>> I tried to fix these problems by:
>>> f = addfile("non_cy.nc", "r")
>>>  prc = f->pre(:,{-5:30},{40:100})
>>>  prc at long_name = "precipitation"
>>>  prc at units = "mm"
>>>  prc&longitude at units = "degrees_east"
>>>  prc&latitude at units = "degrees_north"
>>>  ;printVarSummary(prc)
>>>
>>>  mean_prc = dim_avg_n_Wrap(prc,0)
>>>  mean_prc at long_name = "precipitation"
>>>  mean_prc at units = "mm"
>>>  mean_prc&longitude at units = "degrees_east"
>>>  mean_prc&latitude at units = "degrees_north"
>>>
>>> However, I am getting those errors again. How can I solve these
>>> problems? Please find the script herewith. I used TRMM daily rainfall data.
>>>
>>> Best Regards
>>> Jalal
>>> --
>>> *Md. Jalal Uddin*
>>> MSc in Applied Meteorology (English Language)
>>> Nanjing University of Information, Science and Technology, China
>>> B.Sc. in Disaster Management (Hons.)
>>> Patuakhali Science and Technology University, Bangladesh.
>>> Cell: +8613260859092, +8801792052662
>>> Web: www.dmjalal90.weebly.com
>>> Facebook: jalal.hossen.39
>>> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
>>> Twitter: dmjalal90
>>> Skype: dmjalal90
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>
>>
>> --
>> Adam Phillips
>> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
>> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>>
>> <http://www.cgd.ucar.edu/staff/asphilli>
>>
>
>
> --
> *Md. Jalal Uddin*
> MSc in Applied Meteorology (English Language)
> Nanjing University of Information, Science and Technology, China
> B.Sc. in Disaster Management (Hons.)
> Patuakhali Science and Technology University, Bangladesh.
> Cell: +8613260859092, +8801792052662
> Web: www.dmjalal90.weebly.com
> Facebook: jalal.hossen.39
> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
> Twitter: dmjalal90
> Skype: dmjalal90
>


-- 
Adam Phillips
Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/   303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180911/bbfc9276/attachment.html>


More information about the ncl-talk mailing list