[ncl-talk] Warning: The 'lon2d' attribute must either be the same dimension
Barry Lynn
barry.h.lynn at gmail.com
Thu Jan 5 18:38:43 MST 2017
Hi:
I think you meant:
; snowh_rate_dom = snowh_rate(y_start:y_end,x_start:x_end)
snowh_rate_dom = (/snowh_rate(y_start:y_end,x_start:x_end*)/)*
printVarSummary(snowh_rate_dom)
printVarSummary(snowh_rate)
lat_dom = xlat2d(y_start:y_end,x_start:x_end)
lon_dom = xlon2d(y_start:y_end,x_start:x_end)
snowh_rate_*dom*@lat2d= lat_dom
snowh_rate_*dom*@lon2d = lon_dom
Thank you! -- I had adapted this code from a rectangular grid, and that's
why it didn't work (I needed to strip out the old curvelinear coordinates,
as you noted).
I will try to attach both for future reference, after making some
additional changes to the placement/size of the maps.
On Fri, Jan 6, 2017 at 2:21 AM, Mary Haley <haley at ucar.edu> wrote:
> Thanks for cleaning up the script.
>
> I was able to get it to run and I see what the problem is.
>
> You are subscripting snowh_rate with:
>
> snowh_rate_dom = snowh_rate(y_start:y_end,x_start:x_end)
>
> and then trying to attach coordinate arrays to it with:
>
> lat_dom = xlat2d(y_start:y_end,x_start:x_end)
> lon_dom = xlon2d(y_start:y_end,x_start:x_end)
> snowh_rate_dom!0 = "lat_dom"
> snowh_rate_dom!1 = "lon_dom"
> snowh_rate_dom at lat_dom = lat_dom
> snowh_rate_dom at lon_dom = lon_dom
>
> This is not going to work, because snowh_rate_dom is still a curvilinear
> grid (i.e. a grid defined by 2D lat / lon coordinates), so you still need
> to do the "@lat2d" and "@lon2d" attachments.
>
> Furthermore, when when you do this:
>
> snowh_rate_dom = snowh_rate(y_start:y_end,x_start:x_end)
>
> the original lat2d and lon2d attributes from snowh_rate are copied over,
> and you end up with a lat2d and lon2d that are the wrong size. That's why
> you are getting these two errors:
>
> 0) is_valid_latlon2d_attr: Warning: The 'lat2d' attribute must either be
> the same dimension sizes as the data, or one element larger in both
> directions. Your data will most likely not be overlaid on the map
> correctly.
> (0) is_valid_latlon2d_attr: Warning: The 'lon2d' attribute must either be
> the same dimension sizes as the data, or one element larger in both
> directions. Your data will most likely not be overlaid on the map
> correctly.
>
> You almost have the right idea.
>
> First, when you subscript snowh_rate, make sure the lat2d and lon2d don't
> get copied over by using the special (/ and /) syntax:
>
> snowh_rate_dom = (/snowh_rate(y_start:y_end,x_start:x_end/)) ;
> enclosing a variable with (/ and /) will strip all metadata
>
> Now you need to attach new lat2d and lon2d to snowh_rate_dom that are
> subscripted the same way. You already have this with lat_dom and lon_dom,
> so just use those:
>
> snowh_rate at lat2d = lat_dom
> snowh_rate at lon2d = lon_dom
>
> Make sure you remove these lines:
>
> snowh_rate_dom!0 = "lat_dom"
> snowh_rate_dom!1 = "lon_dom"
> snowh_rate_dom at lat_dom = lat_dom
> snowh_rate_dom at lon_dom = lon_dom
>
> --Mary
>
>
>
>
>
>
> On Thu, Jan 5, 2017 at 10:26 AM, Barry Lynn <barry.h.lynn at gmail.com>
> wrote:
>
>> Hi Mary:
>>
>> This reads the files that I previously attached.
>>
>> I left the wrf calls, but put them in a logical if statement.
>>
>> I did not remove the code that sets the label bars, makes nice maps.
>>
>> If I need to do this, please let me know.
>>
>> The code now gives these errors (below), but I think you can ignore them.
>> (I am not sure how to manually add this information, for this "clean"
>> version.
>>
>> It still gives the one about the dimensions.
>>
>> (Thanks.)
>>
>> 0) is_valid_latlon2d_attr: Warning: The 'lat2d' attribute must either be
>> the same dimension sizes as the data, or one element larger in both
>> directions. Your data will most likely not be overlaid on the map
>> correctly.
>> (0) check_for_y_lat_coord: Warning: Data either does not contain a valid
>> latitude coordinate array or doesn't contain one at all.
>> (0) A valid latitude coordinate array should have a 'units' 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) is_valid_latlon2d_attr: Warning: The 'lon2d' attribute must either
>> be the same dimension sizes as the data, or one element larger in both
>> directions. Your data will most likely not be overlaid on the map
>> correctly.
>> (0) check_for_lon_coord: Warning: Data either does not contain a valid
>> longitude coordinate array or doesn't contain one at all.
>> (0) A valid longitude coordinate array should have a 'units' attribute
>> equal to one of the following values:
>> (0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
>> 'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' '
>>
>> On Thu, Jan 5, 2017 at 7:11 PM, Mary Haley <haley at ucar.edu> wrote:
>>
>>> Barry,
>>>
>>> Offline for a bit. I still can't run your script because it's assuming
>>> that I have a WRF output file, but the "file_for_ncltalk.nc" file you
>>> gave me is not a WRF output file. I
>>> t's very time-consuming for me to try to figure out how to change your
>>> script to work with this file, so please, send me a *clean* script that
>>> works with the files you've sent me so far.
>>>
>>> Thanks,
>>>
>>> --Mary
>>>
>>>
>>>
>>> On Thu, Jan 5, 2017 at 9:58 AM, Barry Lynn <barry.h.lynn at gmail.com>
>>> wrote:
>>>
>>>> Hello:
>>>>
>>>> I have attached the wrf_times.file.
>>>>
>>>> However, I added res at gsnAddCyclic = False
>>>>
>>>> and sres at gsnAddCyclic = False, but the error remains.
>>>>
>>>> The problem pertains not to the mapping of the terrain ("res"), but of
>>>> the snowh_rate_dom variable.
>>>>
>>>> Since, I am plotting a subdomain of terrain by setting:
>>>>
>>>> res at mpFillOn = False
>>>> res at mpMaxLatF = lat_end ; specify the plot domain
>>>> res at mpMinLatF = lat_beg
>>>> res at mpMinLonF = lon_beg ;
>>>> res at mpMaxLonF = lon_end
>>>>
>>>> I need to excerpt snowh_rate_dom from snowh_rate.
>>>>
>>>> I do this as follows, and this is where the error message pertains to.
>>>>
>>>> (Note, if you try to run the program, you have to read in the uploaded
>>>> nc file, not wrf_output, and set ntimes = 0.)
>>>>
>>>> i_loc = 1
>>>> j_loc = 1
>>>> GET_IJ::get_ij(xlat2d,xlon2d,lat_beg,lon_beg,i_loc,j_loc,dim
>>>> s2d(0),dims2d(1))
>>>> x_start = i_loc
>>>> y_start = j_loc
>>>>
>>>> print("x_start = " + x_start)
>>>> print("y_start = " + y_start)
>>>> GET_IJ::get_ij(xlat2d,xlon2d,lat_end,lon_end,i_loc,j_loc,dim
>>>> s2d(0),dims2d(1))
>>>> x_end = i_loc
>>>> y_end = j_loc
>>>> print("x_end = " + x_end)
>>>> print("y_end = " + y_end)
>>>> snowh_rate_dom = snowh_rate(y_start:y_end,x_start:x_end)
>>>> printVarSummary(snowh_rate_dom)
>>>> printVarSummary(snowh_rate)
>>>>
>>>> lat_dom = xlat2d(y_start:y_end,x_start:x_end)
>>>> lon_dom = xlon2d(y_start:y_end,x_start:x_end)
>>>> snowh_rate_dom!0 = "lat_dom"
>>>> snowh_rate_dom!1 = "lon_dom"
>>>> snowh_rate_dom at lat_dom = lat_dom
>>>> snowh_rate_dom at lon_dom = lon_dom
>>>>
>>>>
>>>>
>>>> On Thu, Jan 5, 2017 at 5:45 PM, Mary Haley <haley at ucar.edu> wrote:
>>>>
>>>>> Barry,
>>>>>
>>>>> I was unable to run the script because I was missing an ASCII file.
>>>>>
>>>>> I think the problem is that NCL is trying to add a longitude cyclic
>>>>> point, and you don't want to do this for WRF data. What's happening is that
>>>>> NCL is taking your 270 x 270 data array and turning it into a 270 x 271
>>>>> array, but the lat2d and lon2d are untouched (because gsnAddCyclic doesn't
>>>>> touch these), and hence you end up with two arrays of different sizes.
>>>>>
>>>>> Try setting:
>>>>>
>>>>> res at gsnAddCyclic = False
>>>>>
>>>>> --Mary
>>>>>
>>>>>
>>>>> On Thu, Jan 5, 2017 at 1:42 AM, Barry Lynn <barry.h.lynn at gmail.com>
>>>>> wrote:
>>>>>
>>>>>> Hello:
>>>>>>
>>>>>> I thought that I have done what was indicated on this page, but I
>>>>>> continue to get an error in dimensioning the second variable. I am trying
>>>>>> to plot a subset of the first variable (no problem), and then extract the
>>>>>> same subset from the second variable. I've done it before, but
>>>>>> it doesn't work here.
>>>>>>
>>>>>> I been at this about 4 or 5 hours, tried different things, including
>>>>>> adapting the suggestions here:
>>>>>>
>>>>>> https://www.ncl.ucar.edu/Support/talk_archives/2012/2697.html
>>>>>>
>>>>>> I uploaded an extraction of the original WRF file, in case this is
>>>>>> need.
>>>>>>
>>>>>> Thank you.
>>>>>>
>>>>>> ftp> put wrfsnow.ncl
>>>>>> local: wrfsnow.ncl remote: wrfsnow.ncl
>>>>>> 229 Entering Extended Passive Mode (|||49155|).
>>>>>> 150 Ok to send data.
>>>>>> 100% |********************************************************************|
>>>>>> 7521 3.06 MiB/s 00:00 ETA
>>>>>> 226 Transfer complete.
>>>>>> 7521 bytes sent in 00:00 (11.93 KiB/s)
>>>>>> ftp> put wrfsnow.out.file
>>>>>> local: wrfsnow.out.file remote: wrfsnow.out.file
>>>>>> 229 Entering Extended Passive Mode (|||49161|).
>>>>>> 150 Ok to send data.
>>>>>> 100% |********************************************************************|
>>>>>> 2566 8.21 MiB/s 00:00 ETA
>>>>>> 226 Transfer complete.
>>>>>> 2566 bytes sent in 00:00 (6.15 KiB/s)
>>>>>> ftp> put file_for_ncltalk.nc
>>>>>> local: file_for_ncltalk.nc remote: file_for_ncltalk.nc
>>>>>> 229 Entering Extended Passive Mode (|||49179|).
>>>>>> 150 Ok to send data.
>>>>>> 100% |********************************************************************|
>>>>>> 2279 KiB 278.71 KiB/s 00:00 ETA
>>>>>> 226 Transfer complete.
>>>>>> 2333728 bytes sent in 00:08 (260.45 KiB/s)
>>>>>>
>>>>>>
>>>>>> --
>>>>>> Barry H. Lynn, Ph.D
>>>>>> Senior 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
>>>>>> US 914 432 3108 <(914)%20432-3108>
>>>>>>
>>>>>> _______________________________________________
>>>>>> 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 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
>>>> US 914 432 3108 <(914)%20432-3108>
>>>>
>>>
>>>
>>
>>
>> --
>> Barry H. Lynn, Ph.D
>> Senior 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
>> US 914 432 3108 <(914)%20432-3108>
>>
>
>
--
Barry H. Lynn, Ph.D
Senior 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
US 914 432 3108
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170106/3af0870a/attachment.html
More information about the ncl-talk
mailing list