Running WRF with ECMWF ERA-Interim Dataset

Dataset specific topics involving WPS.

Re: Running WRF with ECMWF ERA-Interim Dataset

Postby arianna.v » Tue Apr 10, 2018 6:56 am

Hi,
I have tried your solution and it works also for me, the only problem is that I have a weird interpolation of the skin temperature on some areas: especially a coastal sea-area in north Adriatic which in July has 282K and surrounding areas (both land and sea) has 296K.

I tried changing both WPS version and simulation period, resolution and METGRID.TBL masking and interpolation.

Has anyone the same problem?
Thanks in advance!
arianna.v
 
Posts: 7
Joined: Mon Feb 20, 2017 3:09 pm

Re: Running WRF with ECMWF ERA-Interim Dataset

Postby aimeusdietger » Sun Nov 04, 2018 4:09 pm

Eowyn wrote:Hi all,

Thanks for all the contributions to this thread to date. It helped me navigate the model vs pressure level issue, , get ungrib and metgrid working, and determine that calc_ecmwf_..exe wasn't necessary if you use pressure levels. All good things.

I'm using pressure-level data, my Vtable is Vtable.ECMWF

Now, real.exe is failing with:
Missing surface temp, replaced with closest level, use_surface set to false.
Missing surface u wind, replaced with closest level, use_surface set to false.
Missing surface v wind, replaced with closest level, use_surface set to false.
PGFIO-F-219/formatted write/internal file/attempt to read/write past end of record.
In source file ../dyn_em/module_initialize_real.f90, at line number 1020

My intermediate files have TT,UU,VV at 1000hPa. So how do I get ungrib to set the surface temp, surface u wind, etc variables? I have a feeling it is in the Vtable but I don't have any experience messing with those.

Thanks in advance,
Eowyn

Update: I was able to modify METGRID.TBL to use the 1000hPa values as the surface values by following this advice from an old FAQ:
Q . If my input data has no surface information, how can I fill the surface level data with the lowest model level data?
A. If your input data does not have any surface fields (e.g. 2 m T, 10m U, V, etc.), you can simply edit metgrid/METGRID.TBL, and add a line like the following for TT, QV, UU and VV:

fill_lev=200100:TT(X)

where X will be the number of the lowest model level in your input data. For example, you can set
fill_lev=200100:TT(100000)

In this example, 1000mb temperature is filled in as surface temperature. In program REAL, you should then set use_surface = .false.


>>> Real still crashes, but for some other reason. I will edit this post if/when I solve that. EDIT 2: Fixed! I used this solution (in which a block of text is added to rrpr.f90 in ugrib/src/):
http://forum.wrfforum.com/viewtopic.php ... SFC#p16075


Hey, have you managed to solve the crash issues?? i will have to read through your post several times so as to understand how got about.. so far so good ,,thanks for the detailed information!
google.com
aimeusdietger
 
Posts: 1
Joined: Thu Nov 01, 2018 1:16 pm

Re: Running WRF with ECMWF ERA-Interim Dataset

Postby jbasconcillo » Thu Aug 29, 2019 9:23 pm

bob.bobb wrote:Thanks Jan,

your post gave me the correct hint to my misunderstanding. I hadn't realized that I had only downloaded the "surface" data from the ECWMF site. After I also downloaded the data for all "pressure levels" and put surface- and pressure-level-files in one folder, it finally worked.

One remark from my side to others having problems: I had to download the pressure-level-data on a Gaussian grid (I used N256), because ungrib wouldn't accept any others. The surface data however I could feed in on a "normal" grid... weird, but as long as it works, I will not complain :lol: .

Thanks again Jan for giving me the tip.


put surface- and pressure-level-files in one folder, it finally worked


What does it mean by this? These are two separate files, right? What does same folder exactly mean? I having the same problem.
jbasconcillo
 
Posts: 19
Joined: Thu May 16, 2019 8:56 am

Re: Running WRF with ECMWF ERA-Interim Dataset

Postby jbasconcillo » Fri Aug 30, 2019 1:54 am

bob.bobb wrote:Thanks Jan,

your post gave me the correct hint to my misunderstanding. I hadn't realized that I had only downloaded the "surface" data from the ECWMF site. After I also downloaded the data for all "pressure levels" and put surface- and pressure-level-files in one folder, it finally worked.

One remark from my side to others having problems: I had to download the pressure-level-data on a Gaussian grid (I used N256), because ungrib wouldn't accept any others. The surface data however I could feed in on a "normal" grid... weird, but as long as it works, I will not complain :lol: .

Thanks again Jan for giving me the tip.



Can you please elaborate on this workaround?
jbasconcillo
 
Posts: 19
Joined: Thu May 16, 2019 8:56 am

Re: Running WRF with ECMWF ERA-Interim Dataset

Postby jbasconcillo » Fri Aug 30, 2019 7:15 am

alainaketh wrote:Hi everyone. I have successfully run a simulation using ERA-Interim data. Here are the steps:

1) Download data from http://rda.ucar.edu/datasets/ds627.0/#!description and go to "Data Access" tab at the top of the page. Please note that you need to register first in order to download the data. You are going to download pressure and surface information. Once you can access to this dataset, select the following options:
- ERA Interim atmospheric model analysis interpolated to pressure levels > Request Access
- ERA Interim atmospheric model analysis for surface > Request Access.
Be sure to download ALL the files correspondent to the dates in your simulation. Example:
I am simulating 13th January 2012 00:00-12:00.
- ERA Interim atmospheric model analysis interpolated to pressure levels > Files downloaded:
@ ei.oper.an.pl.regn128sc.2012011300
@ ei.oper.an.pl.regn128uv.2012011300
@ ei.oper.an.pl.regn128sc.2012011306
@ ei.oper.an.pl.regn128uv.2012011306
@ ei.oper.an.pl.regn128sc.2012011312
@ ei.oper.an.pl.regn128uv.2012011312

- ERA Interim atmospheric model analysis for surface > Files downloaded:
@ ei.oper.an.sfc.regn128sc.2012011300
@ ei.oper.an.sfc.regn128sc.2012011306
@ ei.oper.an.sfc.regn128sc.2012011312

This means that in your namelist.wps you should have:
&share
interval_seconds = 21600,

Vtable used is Vtable.ECMWF

2) run geogrid.
For the pressure data:
3) Type:
./link_grib.csh ../DATA/ei.oper.an.pl*
To link the GRIBFILEs to the pressure data.
4) Modify namelist.wps as follows:
&ungrib
out_format = 'WPS',
prefix = 'PRESS'
5) Run ungrib.
6) Erase the GRIBFILEs
Now repeat steps 3-5 for the surface data.
3b) Type:
./link_grib.csh ../DATA/ei.oper.an.sfc*
To link the GRIBFILEs to the surface data.
4b) Modify namelist.wps as follows:
&ungrib
out_format = 'WPS',
prefix = 'SFC'
5b) Run ungrib.

Now you have a nice set of PRESS:YYYY-MM-DD and SFC:YYY-MM-DD files.

7) Edit namelist.wps as follows:
&metgrid
fg_name = 'SFC','PRESS'
io_form_metgrid = 2
/
8) Run metgrid.

And that's it. Now the problem is that the model yields ridiculously low values of rainfall compared to those obtained using NCEP FNL initial conditions. Could someone tell me why?

Cheers,

Alaina.

EDIT: WRF keeps running happily with ECMWF data but the rainfall I obtain is still way lower than what I expect. I've noticed that when I run ungrib, the highest level at which ungrib calculates all the meteorological fields is 1000, whereas when I use NCEP data it can go up until 2013.

ungrib.log using NCEP dataset ds083.2:
PRES TT UU VV RH HGT PSFC PMSL
-------------------------------------------------------------------------------
2013.0 O O O O O O X O
2001.0 X X X X O X O X
1000.0 X X X X X

ungrib.log using ECMWF data:
PRES GEOPT HGT TT UU VV RH
-------------------------------------------------------------------------------
1000.0 X X X X X
975.0 X X X X X
950.0 X X X X X

Somebody has any ideas?

* * * * * * * *
EDIT AND SOLVED:

Check my post: viewtopic.php?p=26482#p26482


I followed these steps religiously but I cannot proceed from the wrf.exe. I am getting segmentation fault occurred error.
jbasconcillo
 
Posts: 19
Joined: Thu May 16, 2019 8:56 am

Previous

Return to Working with Various Datasets

Who is online

Users browsing this forum: No registered users and 4 guests