[ncl-talk] coordinate problem
Md. Jalal Uddin
dmjalal90 at gmail.com
Tue Sep 11 20:32:27 MDT 2018
Thank you for your help! It works now.
On Tue, Sep 11, 2018 at 11:54 PM Adam Phillips <asphilli at ucar.edu> wrote:
> 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>
>
--
*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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180912/d734a552/attachment-0001.html>
More information about the ncl-talk
mailing list