Re: regridding of ECMWF to WRF grid

From: BasitAli Khan <BasitAli.Khan_at_nyahnyahspammersnyahnyah>
Date: Fri Feb 01 2013 - 05:39:29 MST

Thanks Dennis and Dave. Yes, lonFlip did the job, and it all came up nice
and clean.

Thanks again.

----
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
Skype name: basit.a.khan
On 1/29/13 9:10 PM, "Dennis Shea" <shea@ucar.edu> wrote:
>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_20060
>>427_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
>>
>
________________________________
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
Received on Fri Feb 1 05:39:54 2013

This archive was generated by hypermail 2.1.8 : Fri Feb 01 2013 - 15:31:04 MST