[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