[ncl-talk] WRF interpolating (amendment)

[Ehsan Taghizadeh]

 Yes, exactly.
I want to interpolate WRF data on station locations, which have been included in "stations.dat" file.The problem is that some of the stations in this file don't locate in WRF domain, but most of them get values. If the station doesn't locate in the domain, how it could get interpolated values. For example, 40899 and 40890 are out of domain 2, and they got bad value "9969209968386869046778552952102584320.0", correctly. On the other hand, 40745, also is out of domain 2, but it got value! Logically, It must have "9969209968386869046778552952102584320.0", too.
SincerelyEhsan    On Wednesday, September 11, 2019, 07:16:56 PM GMT+4:30, Barry Lynn <barry.h.lynn at gmail.com> wrote:  
 
 Hello:
DId you mean that you want to interpolate WRF data to station locations?
On Wed, Sep 11, 2019 at 3:08 PM Ehsan Taghizadeh via ncl-talk <ncl-talk at ucar.edu> wrote:

(Pardon me for incomplete previous email)Dears,I hope this email finds you well.I have a script for interpolating WRF data on some stations, using "rcm2points". The input file for stations includes stations which some of them don't locate on wrf out domain 2. I thought these stations must get missing value or so. But all stations, except a few ones, got values!May I ask you how it could possible?I've attached the script (wrfncl.ncl), the input file for stations (stations.dat), an image showing domain 2 (wrfd2.png), also an image showing all stations locations (stations.png). The wrfout file is about 2G. Should I put this file in ftp?Some information of my system is as below:$ uname -a
Linux localhost.localdomain 3.10.0-514.2.2.el7.x86_64 #1 SMP Tue Dec 6 23:06:41 UTC 2016 x86_64 x86_64 x86_64 GNU/Linux
$ ncl -V
6.6.2
Also, about some variables:Variable: wlatType: floatTotal Size: 383616 bytes            95904 valuesNumber of Dimensions: 2Dimensions and sizes:   [south_north | 333] x [west_east | 288]Coordinates:Number Of Attributes: 6  FieldType :   104  MemoryOrder : XY  description : LATITUDE, SOUTH IS NEGATIVE  units :       degree_north  stagger :  coordinates : XLONG XLAT(0)     LATITUDE, SOUTH IS NEGATIVE (degree_north) : min=31.2974   max=40.4085(0)     LONGITUDE, WEST IS NEGATIVE (degree_east) : min=43.4514   max=53.5448
Variable: TC2Type: floatTotal Size: 8055936 bytes            2013984 valuesNumber of Dimensions: 3Dimensions and sizes:   [21] x [333] x [288]Coordinates:Number Of Attributes: 1  units :       C
Variable: TC2_stType: floatTotal Size: 32172 bytes            8043 valuesNumber of Dimensions: 2Dimensions and sizes:   [21] x [383]Coordinates:Number Of Attributes: 1  _FillValue :  9.96921e+36
Variable: SMType: floatTotal Size: 8055936 bytes            2013984 valuesNumber of Dimensions: 3Dimensions and sizes:   [Time | 21] x [south_north | 333] x [west_east | 288]Coordinates:Number Of Attributes: 6  FieldType :   104  MemoryOrder : XYZ  description : SOIL MOISTURE  units :       m3 m-3  stagger :     Z  coordinates : XLONG XLAT XTIME
Variable: SM_stType: floatTotal Size: 32172 bytes            8043 valuesNumber of Dimensions: 2Dimensions and sizes:   [21] x [383]Coordinates:Number Of Attributes: 1  _FillValue :  9.96921e+36
I'll be thankful if I hear from you.

SincerelyEhsan
_______________________________________________
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.DSenior 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

  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190911/f4844f00/attachment.html>