[ncl-talk] Does someone have a suggestion how might solve this?

Barry Lynn barry.h.lynn at gmail.com
Tue Mar 7 08:29:47 MST 2017


Thank you for giving it a go.
---------- Forwarded message ----------
From: Barry Lynn <barry.h.lynn at gmail.com>
Date: Mon, Mar 6, 2017 at 6:55 PM
Subject: An NCL Mystery (regarding writing to wrfinput_d03)
To: ncl-talk <ncl-talk at ucar.edu>


Hello:

I wrote a program (feed_data.ncl -- with a bit of help from Dennis, Mary,
etc) to read an ascii file and write the data to the variable "l_obs" in
the wrfinput_d03.nc file.

When I read the asci file and then write the data to wrfinput_d03.nc, and
then plot the data (using l_obs_map.ncl) it looks absolutely fine
(attached).

However, when I access the data within the WRF itself it is not located at
the correct indices.  For instance, points that should be at about latitude
30.95, longitude 36.6 are at latitude 28, longitude 35).

I realize that there is a switching of indices from the NCL to the WRF, but
I assumed that this was done correctly/automatically​
 lobs_output_1.3km.txt
<https://drive.google.com/file/d/0B3E3GjUeKpBdMTJGOHd0R00ya0E/view?usp=drive_web>
​ as the size of the array in wrfinput_d03 is correct.

Just to reiterate, I can read the data with ncl from wrfinput_d03.nc and
map it (and it corresponds to what I see using a different mapping
program).  However, when I access the array in the WRF, as
l_obs(i,k_lev,j), where k_lev is the time increment, it is not in the
correct location.

Thank you for your suggestions.

Barry




-- 
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/20170307/ac875cec/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: feed_data.ncl
Type: application/octet-stream
Size: 3265 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170307/ac875cec/attachment-0002.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: l_obs_map.pdf
Type: application/pdf
Size: 3515354 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170307/ac875cec/attachment-0001.pdf 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: l_obs_map.ncl
Type: application/octet-stream
Size: 2662 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170307/ac875cec/attachment-0003.obj 


More information about the ncl-talk mailing list