[ncl-talk] plotting_issue

[Dennis Shea]

To convert the Cartesiian values

float jx(jx) ;

jx:long_name = "*x-coordinate in Cartesian system*" ;

jx:units = "m" ;

jx:axis = "X" ;

float iy(iy) ;

iy:long_name = "*y-coordinate in Cartesian system*" ;

iy:units = "m" ;

iy:axis = "Y" ;

float aod(time,* iy, jx*)

to georeferenced lat/lon requires *some* map information.
There is no magic tool that can surmise the necessary information.




On Thu, Jun 20, 2019 at 5:59 AM Barry Lynn via ncl-talk <ncl-talk at ucar.edu>
wrote:

> Hello:
>
> You can convert this data to a grid with lat/lon data by assigning lat/lon
> attributes to your data.  Alternatively, you could interpolate your data to
> a lat/lon grid of your choice.
>
> I can't access my programs now, so I can't send you an example.
>
> However, if your grid is curvelinear, you can use:
>
> https://www.ncl.ucar.edu/Document/Functions/Built-in/rcm2rgrid.shtml
>
> There is also
>
> https://www.ncl.ucar.edu/Document/Functions/Built-in/area_hi2lores.shtml
>
> On Thu, Jun 20, 2019 at 2:42 PM Kunal Bali via ncl-talk <ncl-talk at ucar.edu>
> wrote:
>
>> thanks for this information. but can we change cartesian grid to latlon
>> grid.
>>
>>
>>
>> ---
>> Kunal
>>
>>
>>
>>
>>
>>
>> On Wed, Jun 19, 2019 at 2:49 AM Dennis Shea <shea at ucar.edu> wrote:
>>
>>> Please look at your data.
>>>
>>> The data are on a cartesian grid.
>>> There are no latitude/longitude values.  Hence:
>>>
>>> (a) the *mp** resources are not applicable.
>>> (b) the *gsn_csm_contour_map* function is not applicable. It looks for
>>> latitude/longitude values.
>>>
>>> ===========
>>> *;;*res at trGridType        = "TriangularMesh"  ; plotting is too slow
>>> without this
>>> *;;*res at gsnAddCyclic      = False                 ; applicable for
>>> georeferenced data only
>>>
>>> *;;*res at mpMinLonF = 68.0
>>> *;;*res at mpMaxLonF = 97.5
>>> *;;*res at mpMinLatF = 5.0
>>> *;;*res at mpMaxLatF = 38.0
>>>
>>>   res at pmTickMarkDisplayMode  = "Always"
>>> *;;*res at mpGridSpacingF         =  2
>>> *;;*res at mpGridLineThicknessF   = 4.0
>>>   res at tmXBLabelFontHeightF  = 0.02
>>>   res at tmYLLabelFontHeightF  = 0.02
>>>   res at tmXBLabelFontQuality  = "High"
>>>   res at tmBorderThicknessF    = 6
>>>
>>> ;---Fix the contour levels for all plots
>>>   res at cnLevelSelectionMode = "ManualLevels"
>>>   res at cnMinLevelValF  = 0                   ; -0.1
>>>   res at cnMaxLevelValF  = 1
>>>   res at cnLevelSpacingF = 0.1
>>>   res at cnFillPalette              = "WhBlGrYeRe"
>>>
>>> *;;*plot = gsn_csm_contour_map(wks,var1,res)
>>>   plot = *gsn_csm_contour*(wks,var1,res)
>>>
>>> On Tue, Jun 18, 2019 at 12:14 PM Soma Roy <somaroy892 at gmail.com> wrote:
>>>
>>>> Send your input file.. I will inform you.
>>>>
>>>> Thanks,
>>>> Soma
>>>>
>>>> On Tue, Jun 18, 2019, 23:37 Kunal Bali <kunal.bali9 at gmail.com> wrote:
>>>>
>>>>> Dear NCL users,
>>>>>
>>>>> not able to plot,
>>>>> not receiving any error,
>>>>> script is attached
>>>>>
>>>>>
>>>>> *file summary:*
>>>>>
>>>>> dimensions:
>>>>>
>>>>> time = UNLIMITED ; // (1 currently)
>>>>>
>>>>> bnds = 2 ;
>>>>>
>>>>> jx = 261 ;
>>>>>
>>>>> iy = 217 ;
>>>>>
>>>>> variables:
>>>>>
>>>>> float time(time) ;
>>>>>
>>>>> time:standard_name = "time" ;
>>>>>
>>>>> time:long_name = "time" ;
>>>>>
>>>>> time:bounds = "time_bnds" ;
>>>>>
>>>>> time:units = "hours since 1949-12-01 00:00:00 UTC" ;
>>>>>
>>>>> time:calendar = "standard" ;
>>>>>
>>>>> time:axis = "T" ;
>>>>>
>>>>> double time_bnds(time, bnds) ;
>>>>>
>>>>> float jx(jx) ;
>>>>>
>>>>> jx:long_name = "x-coordinate in Cartesian system" ;
>>>>>
>>>>> jx:units = "m" ;
>>>>>
>>>>> jx:axis = "X" ;
>>>>>
>>>>> float iy(iy) ;
>>>>>
>>>>> iy:long_name = "y-coordinate in Cartesian system" ;
>>>>>
>>>>> iy:units = "m" ;
>>>>>
>>>>> iy:axis = "Y" ;
>>>>>
>>>>> float aod(time, iy, jx) ;
>>>>>
>>>>> aod:standard_name = "atmosphere_optical_thickness_due_to_aerosol" ;
>>>>>
>>>>> aod:long_name = "Aerosol optical thickness in the visible band" ;
>>>>>
>>>>> aod:units = "1" ;
>>>>>
>>>>> aod:cell_methods = "time: point"
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> please provide some information?
>>>>>
>>>>>
>>>>>
>>>>> ---
>>>>> Kunal
>>>>>
>>>>>
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> ncl-talk at ucar.edu
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> ncl-talk at ucar.edu
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>
> --
> Barry H. Lynn, Ph.D
> Senior Associate Scientist, Lecturer,
> The Institute of the Earth Science,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> C.E.O, Weather It Is, LTD
> Weather and Climate Focus
> http://weather-it-is.com
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190620/e0dd38a1/attachment-0001.html>