Re: regridding of ECMWF to WRF grid

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Tue Jan 29 2013 - 11:10:09 MST

The "_" are missing values denoted via the _FillValue attribute.

THE ERA are (nominally) 0->360; WRF are -180 to 180

Use 'lonFlip' located in contributed.ncl. This works on global grids only.

    t2m = era_file->2T_GDS0_SFC ; (:,:,:) ; get T2m (var 167)
    printVarSummary)
    t2m = lonFlip(t2m) ; make -180 to +180
    printVarSummary(t2m)

On 01/29/2013 12:57 AM, BasitAli Khan wrote:
> Thanks David and Denis,
>
> I tried the code, it ran fine and interpolated the ECMWF data to WRF
> grid but dumped "_" in most of the interpolated grid domain. Well, '_'
> is not a default missing value, i checked ecmwf variable before
> interpolation, it does not contain any '_' or missing data but after
> interpolation xgrd had the data in the following form:
>
> T2M =
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
> _, _,
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, 305.5426, 305.4818,
> 305.4208, 305.3597, 305.2985, 305.2372, 305.1757, 305.114, 305.0522,
> 304.9902, 304.9281, 304.8658, 304.8033, 304.7407, 304.6779, 304.615,
> 304.5568, 304.6154, 304.6741, 304.733, 304.7921, 304.8513, 304.9108,
> 304.9704, 305.0303, 305.0903, 305.1505, 305.2109, 305.2716, 305.3325,
> 305.3935, 305.4547, 305.5163, 305.5607, 305.5962, 305.6317, 305.6672,
> 305.7028, 305.7383, 305.7739, 305.8095, 305.8451, 305.8807, 305.9164,
> 305.952, 305.9877, 306.0234, 306.0592, 306.0949, 306.111, 306.0729,
> 306.0346, 305.9962, 305.9576, 305.9188, 305.8799, 305.8408, 305.8015,
> 305.7621, 305.7225, 305.6827, 305.6427, 305.6025, 305.5622, 305.5216,
> 305.4809, 305.4523, 305.4257, 305.3994, 305.3732, 305.3473, 305.3217,
> 305.2962, 305.271, 305.246, 305.2213, 305.1969, 305.1726, 305.1487,
> 305.125, 305.1015, 305.0783, 305.0619, 305.0533, 305.0447, 305.0362,
> 305.0279, 305.0196, 305.0115, 305.0035, 304.9956, 304.9878, 304.9802,
> 304.9726, 304.9651, 304.9579, 304.9507, 304.9436, 304.9363, 304.9231,
> 304.9098, 304.8964, 304.8828, 304.8692, 304.8555, 304.8416, 304.8277,
> 304.8136, 304.7994, 304.7852, 304.7708, 304.7563, 304.7416, 304.7269,
> 304.712, 304.7005, 304.6913, 304.6819, 304.6724, 304.6627, 304.6528,
> 304.6427, 304.6324, 304.6219, 304.6111, 304.6002, 304.5891, 304.5777,
> 304.5662, 304.5544, 304.5424, 304.5343, 304.5432, 304.5517, 304.56,
> 304.568, 304.5757, 304.583, 304.5901, 304.5968, 304.6033, 304.6094,
> 304.6153, 304.6208, 304.6259, 304.6308, 304.6354, 304.6396, 304.6326,
> 304.6216, 304.6103, 304.5988, 304.5869, 304.5749, 304.5625, 304.5499,
> 304.537, 304.5238, 304.5104, 304.4966, 304.4826, 304.4683, 304.4537,
> 304.4389, 304.4266, 304.422, 304.4177, 304.4138, 304.4102, 304.4071,
> 304.4043, 304.4019, 304.4, 304.3984, 304.3973, 304.3965, 304.3961,
> 304.3962, 304.3967, 304.3976, 304.3988, 304.4434, 304.4978, 304.5516,
> 304.6049, 304.6577, 304.7099, 304.7616, 304.8128, 304.8634, 304.9135,
> 304.963, 305.012, 305.0605, 305.1083, 305.1556, 305.2024, 305.2323,
> 305.2268, 305.2212, 305.2155, 305.2096, 305.2037, 305.1976, 305.1914,
> 305.1851, 305.1787, 305.1721, 305.1656, 305.1588, 305.1519, 305.1449,
> 305.1378, 305.1305, 305.1208, 305.1106, 305.1005, 305.0904, 305.0804,
> 305.0704,
>
>
> I am using the following code:
>
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
> begin
>
> wrf_grid = addfile("/project/k14/basit/wrf-sims/wrf_164/20060520.nc","r")
>
> lat2d = wrf_grid->XLAT(0,:,:)
> lon2d = wrf_grid->XLONG(0,:,:)
>
> era_file =
> addfile("/project/k14/basit/wrf-sims/wrf_164/ei_oper_an_sfc_075x075_20060427_20060620.nc","r")
>
> t2m = era_file->2T_GDS0_SFC(:,:,:) ; get T2m (var 167)
> g0_lat_1 = era_file->g0_lat_1(:)
> g0_lon_2 = era_file->g0_lon_2(:)
>
>
> ; invert latitudes for rgrid2rcm
> t2m = t2m(:,::-1,:)
> g0_lat_1 = g0_lat_1(::-1)
>
> printVarSummary(t2m)
> printVarSummary(g0_lat_1)
> printVarSummary(g0_lon_2)
>
> f = addfile("t2m_erai-wrfgridX.nc","c") ; Create blank file to write to
>
> g0_lat_1!0= "south_north"
> g0_lon_2!0= "west_east"
>
> g0_lat_1@units="degree_north"
> g0_lon_2@units="degree_east"
>
> xgrd = rgrid2rcm_Wrap(t2m&g0_lat_1, t2m&g0_lon_2, t2m, lat2d, lon2d, 0)
>
> xgrd!0 = "Time" ; restore meta data
> xgrd!1 = "south_north"
> xgrd!2 = "west_east"
>
> printVarSummary(xgrd)
>
> filedimdef(f,"time",-1,True)
>
> f->T2M = xgrd(:,:,:)
>
> end
>
>
> One thing that i noticed is that the longitude array g0_lon_2 is from
> 0 to 360 which is different from WRF. I tried to change 0 - 360
> longitude array to make it consistent with wrf by applying the
> following, but this also did not have any effect on the outcome.
> dimsLon = dimsizes(g0_lon_2)
> g0_lon_2X = new((/dimsLon/),float)
> do i = 0, dimsLon-1, 1
> if(g0_lon_2(i) .le. 180) then
>
> g0_lon_2X(i) = g0_lon_2(i)
> else
> g0_lon_2X(i) = g0_lon_2(i) - 360
> end if
> end do
>
> delete(g0_lon_2)
> g0_lon_2 = g0_lon_2X
>
>
> So, I guess, the problem lies somewhere else. I would greatly appreciate
> if you could suggest some tip/solution. I've not yet tried ESMF regrid
> method suggested by Denis Shea though.
>
> Cheers,
>
>
> ----
>
> Basit A. Khan, Ph.D.
>
> Postdoctoral Research Fellow
>
> Division of Physical Sciences & Engineering
>
> Office# 3204, Level 3, Building 1,
>
> King Abdullah University of Science & Technology
>
> 4700 King Abdullah Blvd, Box 2753, Thuwal 23955 6900,
>
> Kingdom of Saudi Arabia.
>
> Office: +966(0)2 808 0276, Mobile: +966(0)5 0860 3617
>
> E-mail: basitali.khan@kaust.edu.sa <mailto:basitali.khan@kaust.edu.sa>
>
> Skype name: basit.a.khan
>
>
> From: David Jones <jonesd647@yahoo.ca <mailto:jonesd647@yahoo.ca>>
> Reply-To: David Jones <jonesd647@yahoo.ca <mailto:jonesd647@yahoo.ca>>
> Date: Sun, 27 Jan 2013 11:39:39 -0800
> To: "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>" <ncl-talk@ucar.edu
> <mailto:ncl-talk@ucar.edu>>, KAUST <basitali.khan@kaust.edu.sa
> <mailto:basitali.khan@kaust.edu.sa>>
> Subject: Re: regridding of ECMWF to WRF grid
>
> You could also try using the NCL function "rgrid2rcm" like in the script
> below.
> I've used this for ERA-I data before.
>
> Dave.
>
> ======================================================
>
> ; interpolation script
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>
> begin
>
> wrf_grid = addfile("wrfout_d01_2001-12-24_00:00:00.nc","r")
>
> lat2d = wrf_grid->XLAT(0,:,:)
> lon2d = wrf_grid->XLONG(0,:,:)
>
> files = asciiread("erai_files.asc",-1,"string") ; -1 for all
> print(files)
> erai_grid = addfiles(files,"r")
> t2m = erai_grid[:]->2T_GDS4_SFC(:,:,:)
> printVarSummary(t2m)
>
> t2m = t2m(:,::-1,:) ; INVERT LATITUDES FOR rgrid2rcm
> printVarSummary(t2m)
>
> g4_lat_1 = erai_grid[:]->g4_lat_1(:)
> printVarSummary(g4_lat_1)
> g4_lat_1 = g4_lat_1(::-1)
> printVarSummary(g4_lat_1)
> g4_lon_2 = erai_grid[:]->g4_lon_2(:)
>
> f = addfile("t2m_erai-wrfgrid.nc","c") ; Create blank file to write to
>
> ;xgrd = rgrid2rcm(x&lat, x&lon, x, lat2d, lon2d, 0)
> xgrd = rgrid2rcm_Wrap(t2m&g4_lat_1, t2m&g4_lon_2, t2m, lat2d, lon2d, 0)
>
> xgrd!0 = "Time" ; restore meta data
> xgrd!1 = "south_north"
> xgrd!2 = "west_east"
>
> printVarSummary(xgrd)
>
> filedimdef(f,"time",-1,True)
>
> f->T2M = xgrd(:,:,:)
>
> end
>
> =====================================================
>
>
> ------------------------------------------------------------------------
> *From:* "ncl-talk-request@ucar.edu <mailto:ncl-talk-request@ucar.edu>"
> <ncl-talk-request@ucar.edu <mailto:ncl-talk-request@ucar.edu>>
> *To:* ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>
> *Sent:* Sunday, January 27, 2013 8:00:04 PM
> *Subject:* ncl-talk Digest, Vol 110, Issue 33
>
> Send ncl-talk mailing list submissions to
> ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> or, via email, send a message with subject or body 'help' to
> ncl-talk-request@ucar.edu <mailto:ncl-talk-request@ucar.edu>
>
> You can reach the person managing the list at
> ncl-talk-owner@ucar.edu <mailto:ncl-talk-owner@ucar.edu>
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of ncl-talk digest..."
>
>
> Today's Topics:
>
> 1. Re: How to define the rec_type in the function of fbinrecread
> to read a special binary file (Guangshan Chen)
> 2. Re: How to define the rec_type in the function of fbinrecread
> to read a special binary file (shea@ucar.edu <mailto:shea@ucar.edu>)
> 3. Re: How to define the rec_type in the function of fbinrecread
> to read a special binary file (Guangshan Chen)
> 4. regridding of ECMWF to WRF grid (BasitAli Khan)
> 5. Re: regridding of ECMWF to WRF grid (BasitAli Khan)
> 6. Re: regridding of ECMWF to WRF grid (shea@ucar.edu
> <mailto:shea@ucar.edu>)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Sat, 26 Jan 2013 19:19:22 -0600
> From: Guangshan Chen <gchen9@gmail.com <mailto:gchen9@gmail.com>>
> Subject: Re: How to define the rec_type in the function of
> fbinrecread to read a special binary file
> To: Dennis Shea <shea@ucar.edu <mailto:shea@ucar.edu>>
> Cc: "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> NCL"
> <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>>
> Message-ID: <EB1B55CD-DCAD-4D64-964E-B7A743E45D54@gmail.com
> <mailto:EB1B55CD-DCAD-4D64-964E-B7A743E45D54@gmail.com>>
> Content-Type: text/plain; charset=iso-8859-1
>
> Dear Dennis and Dave,
>
> Thanks for the suggestion.
>
> It does not work.
>
> Hopeful NCL can read read multiple types at a time one day!
>
> Guangshan
>
> On Jan 25, 2013, at 2:31 PM, Dennis Shea <shea@ucar.edu
> <mailto:shea@ucar.edu>> wrote:
>
>> NCL's fbin*read functions can only read one type at a time/
>>
>> fili = "example01"
>> ii = fbinrecread(fili, 0, 3, "integer")
>> rr = fbinrecread(fili, 0, 3, "float" )
>>
>> a = ii(0)
>> b = rr(1)
>> c = rr(2)
>>
>> Good luck
>>
>> On 01/25/2013 09:17 AM, Guangshan Chen wrote:
>>> Dear all,
>>>
>>> I have a binary file. It is written as the following in the fortran code:
>>>
>>> integer a
>>> real b,c
>>>
>>> a=10
>>> b=20.0
>>> c=30.0
>>>
>>> open(11, file="example01", form="unformatted")
>>> write(11) a, b, c
>>>
>>>
>>> There is only record in this binary file. If I would like to use fbinrecread to read
>>> finrecread("example01", 0, 3, ????). How to define the rec_type?
>>>
>>> a is integer, b and c are real.
>>>
>>> Thanks for any help.
>>>
>>> Guangshan
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>>http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>
>
>
>
> ------------------------------
>
> Message: 2
> Date: Sat, 26 Jan 2013 19:38:13 -0700
> From: shea@ucar.edu <mailto:shea@ucar.edu>
> Subject: Re: How to define the rec_type in the function of
> fbinrecread to read a special binary file
> To: "Guangshan Chen" <gchen9@gmail.com <mailto:gchen9@gmail.com>>
> Cc: "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> NCL"
> <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>>
> Message-ID:
> <7a1ed3f797131e150d70b790c8c4bd1f.squirrel@webmail.cgd.ucar.edu
> <mailto:7a1ed3f797131e150d70b790c8c4bd1f.squirrel@webmail.cgd.ucar.edu>>
> Content-Type: text/plain;charset=iso-8859-1
>
> The outlined approch has worked on multiple machines and operating systems.
>
> Perhaps the fortran sequential file was created on a 'big-endian' machine
> and you are reading on a little-endian machine or vice versa.
>
> There is no automatic way to test. It is the user's responsibility to use
> thr sppropriate 'setfileoption' .
>
>> Dear Dennis and Dave,
>>
>> Thanks for the suggestion.
>>
>> It does not work.
>>
>> Hopeful NCL can read read multiple types at a time one day!
>>
>> Guangshan
>>
>> On Jan 25, 2013, at 2:31 PM, Dennis Shea <shea@ucar.edu <mailto:shea@ucar.edu>> wrote:
>>
>>> NCL's fbin*read functions can only read one type at a time/
>>>
>>> fili = "example01"
>>> ii = fbinrecread(fili, 0, 3, "integer")
>>> rr = fbinrecread(fili, 0, 3, "float" )
>>>
>>> a = ii(0)
>>> b = rr(1)
>>> c = rr(2)
>>>
>>> Good luck
>>>
>>> On 01/25/2013 09:17 AM, Guangshan Chen wrote:
>>>> Dear all,
>>>>
>>>> I have a binary file. It is written as the following in the fortran
>>>> code:
>>>>
>>>> integer a
>>>> real b,c
>>>>
>>>> a=10
>>>> b=20.0
>>>> c=30.0
>>>>
>>>> open(11, file="example01", form="unformatted")
>>>> write(11) a, b, c
>>>>
>>>>
>>>> There is only record in this binary file. If I would like to use
>>>> fbinrecread to read
>>>> finrecread("example01", 0, 3, ????). How to define the rec_type?
>>>>
>>>> a is integer, b and c are real.
>>>>
>>>> Thanks for any help.
>>>>
>>>> Guangshan
>>>>
>>>> _______________________________________________
>>>> ncl-talk mailing list
>>>> List instructions, subscriber options, unsubscribe:
>>>>http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>
>>>
>>
>>
>
>
>
>
> ------------------------------
>
> Message: 3
> Date: Sat, 26 Jan 2013 20:41:21 -0600
> From: Guangshan Chen <gchen9@gmail.com <mailto:gchen9@gmail.com>>
> Subject: Re: How to define the rec_type in the function of
> fbinrecread to read a special binary file
> To: shea@ucar.edu <mailto:shea@ucar.edu>
> Cc: "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu> NCL"
> <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>>
> Message-ID: <42A91B9A-015E-4323-8044-32133DEC2CE1@gmail.com
> <mailto:42A91B9A-015E-4323-8044-32133DEC2CE1@gmail.com>>
> Content-Type: text/plain; charset=iso-8859-1
>
> Dear Dennis,
>
> Thanks.
>
> It was created on a "big-endian" machine. I have set the setfileoption
> to read big-endian.
>
> Guangshan
>
> On Jan 26, 2013, at 8:38 PM, shea@ucar.edu <mailto:shea@ucar.edu> wrote:
>
>> The outlined approch has worked on multiple machines and operating systems.
>>
>> Perhaps the fortran sequential file was created on a 'big-endian' machine
>> and you are reading on a little-endian machine or vice versa.
>>
>> There is no automatic way to test. It is the user's responsibility to use
>> thr sppropriate 'setfileoption' .
>>
>>> Dear Dennis and Dave,
>>>
>>> Thanks for the suggestion.
>>>
>>> It does not work.
>>>
>>> Hopeful NCL can read read multiple types at a time one day!
>>>
>>> Guangshan
>>>
>>> On Jan 25, 2013, at 2:31 PM, Dennis Shea <shea@ucar.edu <mailto:shea@ucar.edu>> wrote:
>>>
>>>> NCL's fbin*read functions can only read one type at a time/
>>>>
>>>> fili = "example01"
>>>> ii = fbinrecread(fili, 0, 3, "integer")
>>>> rr = fbinrecread(fili, 0, 3, "float" )
>>>>
>>>> a = ii(0)
>>>> b = rr(1)
>>>> c = rr(2)
>>>>
>>>> Good luck
>>>>
>>>> On 01/25/2013 09:17 AM, Guangshan Chen wrote:
>>>>> Dear all,
>>>>>
>>>>> I have a binary file. It is written as the following in the fortran
>>>>> code:
>>>>>
>>>>> integer a
>>>>> real b,c
>>>>>
>>>>> a=10
>>>>> b=20.0
>>>>> c=30.0
>>>>>
>>>>> open(11, file="example01", form="unformatted")
>>>>> write(11) a, b, c
>>>>>
>>>>>
>>>>> There is only record in this binary file. If I would like to use
>>>>> fbinrecread to read
>>>>> finrecread("example01", 0, 3, ????). How to define the rec_type?
>>>>>
>>>>> a is integer, b and c are real.
>>>>>
>>>>> Thanks for any help.
>>>>>
>>>>> Guangshan
>>>>>
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> List instructions, subscriber options, unsubscribe:
>>>>>http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>>
>>>
>>>
>>
>>
>
>
>
> ------------------------------
>
> Message: 4
> Date: Sun, 27 Jan 2013 09:25:08 +0000
> From: BasitAli Khan <BasitAli.Khan@kaust.edu.sa
> <mailto:BasitAli.Khan@kaust.edu.sa>>
> Subject: regridding of ECMWF to WRF grid
> To: "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>" <ncl-talk@ucar.edu
> <mailto:ncl-talk@ucar.edu>>
> Message-ID:
> <5B0AEC9F228B0E4093880B97E72BBDCE9D1710@WTHEMXND03.KAUST.EDU.SA
> <mailto:5B0AEC9F228B0E4093880B97E72BBDCE9D1710@WTHEMXND03.KAUST.EDU.SA>>
> Content-Type: text/plain; charset="us-ascii"
>
> Hello,
>
> I want to regrid Era Interim Reanalysis data (0.75 x 0.75 deg) to my
> high resolution WRF gird (3 km). Just wondering which
> regridding/interpolation function is best to use for this purpose. I
> know there are heaps of functions and regridding scripts available but
> could not find anything that suits my requirement. WRF has a curvilinear
> grid while Era-Interim Re-analysis (if i am not wrong) uses Gaussian
> grid. I would greatly appreciate if anyone have the script and would
> like to share it.
>
> Thanks in advance
> Best regards,
> basit
>
> ________________________________
>
> This message and its contents including attachments are intended solely
> for the original recipient. If you are not the intended recipient or
> have received this message in error, please notify me immediately and
> delete this message from your computer system. Any unauthorized use or
> distribution is prohibited. Please consider the environment before
> printing this email.
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL:
> http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20130127/5f4bdea3/attachment.html
>
>
> ------------------------------
>
> Message: 5
> Date: Sun, 27 Jan 2013 11:23:25 +0000
> From: BasitAli Khan <BasitAli.Khan@kaust.edu.sa
> <mailto:BasitAli.Khan@kaust.edu.sa>>
> Subject: Re: regridding of ECMWF to WRF grid
> To: "bjfisel@gmail.com <mailto:bjfisel@gmail.com>" <bjfisel@gmail.com
> <mailto:bjfisel@gmail.com>>, "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>"
> <ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>>
> Message-ID:
> <5B0AEC9F228B0E4093880B97E72BBDCE9D1724@WTHEMXND03.KAUST.EDU.SA
> <mailto:5B0AEC9F228B0E4093880B97E72BBDCE9D1724@WTHEMXND03.KAUST.EDU.SA>>
> Content-Type: text/plain; charset="windows-1252"
>
> Hi,
> That's right, its generally an easy fix, however, unfortunately WPS does
> not interpolate all the Era-Interim variables to WRF grid e.g T2, U10
> and V10. To output these vars i probably would need to make changes in
> V-table and/or at some other places and may also need to recompile WPS.
> I am not sure if i can take values of the first level of TT, UU and VV
> from met_em_d0* with any degree of confidence but i would appreciate if
> you could share your experience on this.
>
> Cheers,
> ----
> Basit A. Khan, Ph.D.
> Postdoctoral Research Fellow
> Division of Physical Sciences & Engineering
> Office# 3204, Level 3, Building 1,
> King Abdullah University of Science & Technology
> 4700 King Abdullah Blvd, Box 2753, Thuwal 23955 ?6900,
> Kingdom of Saudi Arabia.
>
> Office: +966(0)2 808 0276, Mobile: +966(0)5 0860 3617
> E-mail: basitali.khan@kaust.edu.sa
> <mailto:basitali.khan@kaust.edu.sa><mailto:basitali.khan@kaust.edu.sa
> <mailto:basitali.khan@kaust.edu.sa>>
> Skype name: basit.a.khan
>
> From: "bjfisel@gmail.com
> <mailto:bjfisel@gmail.com><mailto:bjfisel@gmail.com
> <mailto:bjfisel@gmail.com>>" <bjfisel@gmail.com
> <mailto:bjfisel@gmail.com><mailto:bjfisel@gmail.com
> <mailto:bjfisel@gmail.com>>>
> Date: Sun, 27 Jan 2013 05:09:10 -0600
> To: KAUST <basitali.khan@kaust.edu.sa
> <mailto:basitali.khan@kaust.edu.sa><mailto:basitali.khan@kaust.edu.sa
> <mailto:basitali.khan@kaust.edu.sa>>>
> Subject: Re: regridding of ECMWF to WRF grid
>
> Hi,
>
> You should ingest the ECMWF data into the WRF Preprocessing System
> (WPS), which will put the data on the WRF 3-km grid. There is no need to
> regrid.
>
>
> On Jan 27, 2013, at 3:25, BasitAli Khan <BasitAli.Khan@kaust.edu.sa
> <mailto:BasitAli.Khan@kaust.edu.sa><mailto:BasitAli.Khan@kaust.edu.sa
> <mailto:BasitAli.Khan@kaust.edu.sa>>> wrote:
>
> Hello,
>
> I want to regrid Era Interim Reanalysis data (0.75 x 0.75 deg) to my
> high resolution WRF gird (3 km). Just wondering which
> regridding/interpolation function is best to use for this purpose. I
> know there are heaps of functions and regridding scripts available but
> could not find anything that suits my requirement. WRF has a curvilinear
> grid while Era-Interim Re-analysis (if i am not wrong) uses Gaussian
> grid. I would greatly appreciate if anyone have the script and would
> like to share it.
>
> Thanks in advance
> Best regards,
> basit
>
> ________________________________
>
> This message and its contents including attachments are intended solely
> for the original recipient. If you are not the intended recipient or
> have received this message in error, please notify me immediately and
> delete this message from your computer system. Any unauthorized use or
> distribution is prohibited. Please consider the environment before
> printing this email.
> _______________________________________________
> ncl-talk mailing list
> 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/20130127/fed20ba9/attachment.html
>
>
> ------------------------------
>
> Message: 6
> Date: Sun, 27 Jan 2013 08:28:00 -0700
> From: shea@ucar.edu <mailto:shea@ucar.edu>
> Subject: Re: regridding of ECMWF to WRF grid
> To: "BasitAli Khan" <BasitAli.Khan@kaust.edu.sa
> <mailto:BasitAli.Khan@kaust.edu.sa>>
> Cc: "ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>" <ncl-talk@ucar.edu
> <mailto:ncl-talk@ucar.edu>>
> Message-ID:
> <f83aabad438e5ab23a9bb29a3d4a54ea.squirrel@webmail.cgd.ucar.edu
> <mailto:f83aabad438e5ab23a9bb29a3d4a54ea.squirrel@webmail.cgd.ucar.edu>>
> Content-Type: text/plain;charset=iso-8859-1
>
> rectilinear [ECMWF] to curvilinear [WRF]: You could use ESMF regrid
> example 9 and (say) example 20 as a start. Generating the weight file may
> take some time but, once the weight file is generated, the actual
> regridding is quite fast.
>
>> Hi,
>> That's right, its generally an easy fix, however, unfortunately WPS does
>> not interpolate all the Era-Interim variables to WRF grid e.g T2, U10 and
>> V10. To output these vars i probably would need to make changes in V-table
>> and/or at some other places and may also need to recompile WPS. I am not
>> sure if i can take values of the first level of TT, UU and VV from
>> met_em_d0* with any degree of confidence but i would appreciate if you
>> could share your experience on this.
>>
>> Cheers,
>> ----
>> Basit A. Khan, Ph.D.
>> Postdoctoral Research Fellow
>> Division of Physical Sciences & Engineering
>> Office# 3204, Level 3, Building 1,
>> King Abdullah University of Science & Technology
>> 4700 King Abdullah Blvd, Box 2753, Thuwal 23955 ?6900,
>> Kingdom of Saudi Arabia.
>>
>> Office: +966(0)2 808 0276, Mobile: +966(0)5 0860 3617
>> E-mail:basitali.khan@kaust.edu.sa
> <mailto:basitali.khan@kaust.edu.sa><mailto:basitali.khan@kaust.edu.sa
> <mailto:basitali.khan@kaust.edu.sa>>
>> Skype name: basit.a.khan
>>
>> From: "bjfisel@gmail.com <mailto:bjfisel@gmail.com><mailto:bjfisel@gmail.com
> <mailto:bjfisel@gmail.com>>"
>> <bjfisel@gmail.com <mailto:bjfisel@gmail.com><mailto:bjfisel@gmail.com
> <mailto:bjfisel@gmail.com>>>
>> Date: Sun, 27 Jan 2013 05:09:10 -0600
>> To: KAUST <basitali.khan@kaust.edu.sa
> <mailto:basitali.khan@kaust.edu.sa><mailto:basitali.khan@kaust.edu.sa
> <mailto:basitali.khan@kaust.edu.sa>>>
>> Subject: Re: regridding of ECMWF to WRF grid
>>
>> Hi,
>>
>> You should ingest the ECMWF data into the WRF Preprocessing System (WPS),
>> which will put the data on the WRF 3-km grid. There is no need to regrid.
>>
>>
>> On Jan 27, 2013, at 3:25, BasitAli Khan
>> <BasitAli.Khan@kaust.edu.sa
> <mailto:BasitAli.Khan@kaust.edu.sa><mailto:BasitAli.Khan@kaust.edu.sa
> <mailto:BasitAli.Khan@kaust.edu.sa>>> wrote:
>>
>> Hello,
>>
>> I want to regrid Era Interim Reanalysis data (0.75 x 0.75 deg) to my
>> high resolution WRF gird (3 km). Just wondering which
>> regridding/interpolation function is best to use for this purpose. I know
>> there are heaps of functions and regridding scripts available but could
>> not find anything that suits my requirement. WRF has a curvilinear grid
>> while Era-Interim Re-analysis (if i am not wrong) uses Gaussian grid. I
>> would greatly appreciate if anyone have the script and would like to share
>> it.
>>
>> Thanks in advance
>> Best regards,
>> basit
>>
>> ________________________________
>>
>> This message and its contents including attachments are intended solely
>> for the original recipient. If you are not the intended recipient or have
>> received this message in error, please notify me immediately and delete
>> this message from your computer system. Any unauthorized use or
>> distribution is prohibited. Please consider the environment before
>> printing this email.
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>>http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>>http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>
>
>
> ------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk@ucar.edu <mailto:ncl-talk@ucar.edu>
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
> End of ncl-talk Digest, Vol 110, Issue 33
> *****************************************
>
>
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jan 29 11:10:22 2013

This archive was generated by hypermail 2.1.8 : Tue Jan 29 2013 - 22:44:26 MST