[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