Running WRF with Copernicus Data Service ERA5 pressure level

Dataset specific topics involving WPS.

Running WRF with Copernicus Data Service ERA5 pressure level

Postby zhuohaifeng » Thu Oct 11, 2018 8:47 pm

Hi,
I usually use " Vtable -> ungrib/Variable_Tables/Vtable.ERA-interim.pl " to do ungrib with ECMWF's ERA5 data before, and it works very well. But after the ERA5 migrate to Copernocus Data Service since this September 2018, the Vtable.ERA-interim.pl doesn't work at all. I did a test with the CDS ERA5 dataset with Vtable.ERA-interim.pl, only when I download grib data globally it goes well, it does not work while I set the 'area' keyword. The problem is that I don't want to download the whole global area data, I just want to download an area with high resolution.
Does anyone could give me a new Vtable.ERA5.pl to me for running WRF with ERA5 pressure levels?



Here is the problem I have.
While I do ungrib.exe with the CDS ERA5 pressure level data now. The data value are all 0.00 while writing in WPS format using ungrib.exe/WRF if setting " 'area': '90/0/-90/180' " for downloading python code.

1, I use get_CDS_ERA5_2wrf.py to download ERA5 data;
{{
[#43#haifeng@zhuohf ~/2018/CDS_ERA5/test]$vim get_CDS_ERA5_2wrf.py --> attachmented
[#44#haifeng@zhuohf ~/2018/CDS_ERA5/test]$python3 get_CDS_ERA5_2wrf.py
2018-10-08 14:17:07,576 INFO Sending request to https://cds.climate.copernicus.eu/api/v ... gle-levels
2018-10-08 14:17:08,482 INFO Request is queued
2018-10-08 14:17:09,707 INFO Request is running
2018-10-08 14:17:17,524 INFO Request is completed
2018-10-08 14:17:17,524 INFO Downloading http://136.156.132.110/cache-compute-00 ... 47e65.grib to ERA5-sf-2017-01-01.grib (63.6M)
2018-10-08 14:21:13,123 INFO Download rate 276.3K/s
2018-10-08 14:21:13,123 INFO Sending request to https://cds.climate.copernicus.eu/api/v ... ure-levels
2018-10-08 14:21:14,037 INFO Request is queued
2018-10-08 14:21:15,263 INFO Request is running
2018-10-08 14:21:47,820 INFO Request is completed
2018-10-08 14:21:47,820 INFO Downloading http://136.156.132.110/cache-compute-00 ... e674b.grib to ERA5-pl-2017-01-01.grib (489.6M)
2018-10-08 14:28:45,002 INFO Download rate 1.2M/s
2018-10-08 14:28:45,003 INFO Sending request to https://cds.climate.copernicus.eu/api/v ... gle-levels
2018-10-08 14:28:46,483 INFO Request is queued
2018-10-08 14:28:47,774 INFO Request is running
2018-10-08 14:28:55,753 INFO Request is completed
2018-10-08 14:28:56,139 INFO Downloading http://136.156.132.105/cache-compute-00 ... 85893.grib to ERA5-sf-2017-01-02.grib (63.6M)
}}

2, I use g1print.exe to got a file named g1print.log;
{{
[#11#haifeng@zhuohf ~/WRFV3_9_model/WPS/util]$./g1print.exe ~/2018/CDS_ERA5/test/ERA5-sf-2017-01-01.grib > g1print.log
[#12#haifeng@zhuohf ~/WRFV3_9_model/WPS/util]$vim g1print.log --> attachmented
}}

3, I use ncl_filedump got a file named ncl_filedump.log.
{{
[#10#haifeng@zhuohf ~/2018/CDS_ERA5/test]$ncl_filedump ERA5-sf-2017-01-01.grib > ncl_filedump.log
[#11#haifeng@zhuohf ~/2018/CDS_ERA5/test]$vim ncl_filedump.log
[#12#haifeng@zhuohf ~/2018/CDS_ERA5/test]$ncl_filedump -v 10U_GDS0_SFC ERA5-sf-2017-01-01.grib > ncl_filedump.log
[#13#haifeng@zhuohf ~/2018/CDS_ERA5/test]$vim ncl_filedump.log -->attachmented
}}
4, then do ungrib.exe, got a file named ungrib.log.
{{
[#50#haifeng@zhuohf ~/WRFV3_9_model/WPS]$./link_grib.csh ~/2018/CDS_ERA5/test/ERA5-sf-2017-01-01.grib
[#51#haifeng@zhuohf ~/WRFV3_9_model/WPS]$vim namelist.wps
[#52#haifeng@zhuohf ~/WRFV3_9_model/WPS]$ls -ls
Total 42268
4 drwxr-xr-x. 2 haifeng haifeng 4096 Sep 30 14:06 arch
24 drwxrwxr-x. 2 haifeng haifeng 24576 Sep 30 14:06 bohai_ERA5_0p125_3km
4 -rwxr-xr-x. 1 haifeng haifeng 1765 Sep 30 14:06 clean
8 -rwxr-xr-x. 1 haifeng haifeng 4796 Sep 30 14:06 compile
12 -rwxr-xr-x. 1 haifeng haifeng 11514 Sep 30 14:06 configure
4 -rw-rw-r--. 1 haifeng haifeng 3527 Sep 30 16:43 configure.wps
41908 -rw-rw-r--. 1 haifeng haifeng 42909776 Sep 30 14:06 geo_em.d01.nc
4 drwxr-xr-x. 4 haifeng haifeng 4096 Sep 30 17:27 geogrid
0 lrwxrwxrwx. 1 haifeng haifeng 23 Sep 30 17:27 geogrid.exe -> geogrid/src/geogrid.exe
12 -rw-rw-r--. 1 haifeng haifeng 9303 Sep 30 14:06 geogrid.log
0 lrwxrwxrwx. 1 haifeng haifeng 56 Oct 8 14:31 GRIBFILE.AAA -> /home/haifeng/2018/CDS_ERA5/test/ERA5-sf-2017-01-01.grib
4 -rwxr-xr-x. 1 haifeng haifeng 1331 Sep 30 14:00 link_grib.csh
4 drwxr-xr-x. 3 haifeng haifeng 4096 Sep 30 17:27 metgrid
0 lrwxrwxrwx. 1 haifeng haifeng 23 Sep 30 17:27 metgrid.exe -> metgrid/src/metgrid.exe
244 -rw-rw-r--. 1 haifeng haifeng 247677 Sep 30 14:06 metgrid.log
4 -rw-r--r--. 1 haifeng haifeng 1485 Oct 8 14:31 namelist.wps
4 -rw-r--r--. 1 haifeng haifeng 2887 Sep 30 14:06 namelist.wps.all_options
4 -rwxr-xr-x. 1 haifeng haifeng 2079 Sep 30 14:00 namelist.wps.fire
4 -rw-r--r--. 1 haifeng haifeng 1810 Sep 30 14:00 namelist.wps.global
4 -rw-r--r--. 1 haifeng haifeng 656 Sep 30 14:00 namelist.wps.nmm
8 -rw-r--r--. 1 haifeng haifeng 6337 Sep 30 14:06 README
4 drwxr-xr-x. 4 haifeng haifeng 4096 Sep 30 17:27 ungrib
0 lrwxrwxrwx. 1 haifeng haifeng 21 Sep 30 17:27 ungrib.exe -> ungrib/src/ungrib.exe
4 drwxr-xr-x. 3 haifeng haifeng 4096 Oct 8 14:23 util
0 lrwxrwxrwx. 1 haifeng haifeng 44 Oct 8 14:12 Vtable -> ungrib/Variable_Tables/Vtable.ERA-interim.pl
[#53#haifeng@zhuohf ~/WRFV3_9_model/WPS]$vim namelist.wps
[#57#haifeng@zhuohf ~/WRFV3_9_model/WPS]$./ungrib.exe > zzz
c_open.c: No such file or directory
Name of source model =>ECMWF
Name of source model =>ECMWF
Name of source model =>ECMWF
Name of source model =>ECMWF
[#58#haifeng@zhuohf ~/WRFV3_9_model/WPS]$vim namelist.wps
[#59#haifeng@zhuohf ~/WRFV3_9_model/WPS]$vim ungrib.log -->attachmented
}}
ungrib.log:

2018-10-08 14:33:13.426 --- Inventory for date = 2017-01-01 00:00:00
2018-10-08 14:33:13.426 --- PRES GEOPT HGT TT UU VV RH DEWPT LANDSEA SOILGEO SOILHGT PSFC PMSL SKINTEMP SEAICE SST SNOW_DEN SNOW_EC SNOW SNOWH ST000007 ST007028 ST028100 ST100289 SM000007 SM007028 SM028100 SM100289
2018-10-08 14:33:13.426 --- -------------------------------------------------
2018-10-08 14:33:13.426 --- 200100 O O X X X O X X O O X X X O X X X O O X X X X X X X X
2018-10-08 14:33:13.426 --- -------------------------------------------------
2018-10-08 14:33:13.426 --- DEBUG: writing in WPS format iunit = 13, map%igrid = 0
2018-10-08 14:33:13.428 --- DEBUG: hdate = 2017-01-01_00:00:00, xfcst = 0.000000
2018-10-08 14:33:13.428 --- DEBUG: map%source = ECMWF, field = TT, units = K
2018-10-08 14:33:13.428 --- DEBUG: Desc = Temperature, level = 200100.000000
2018-10-08 14:33:13.428 --- DEBUG: map%nx = 721, map%ny = 481
2018-10-08 14:33:13.428 --- DEBUG: max value = 0.000000 , min value = 0.000000

Again, I try to not use 'area' keyword by " # 'area': '90/0/-30/180', "
1, I use get_CDS_ERA5_2wrf.py ( which is not use the 'area' keyword ) to download ERA5 data;
{{
[#45#haifeng@zhuohf ~/2018/CDS_ERA5/test]$mkdir 20170101
[#46#haifeng@zhuohf ~/2018/CDS_ERA5/test]$mv ERA5-* 20170101/
[#47#haifeng@zhuohf ~/2018/CDS_ERA5/test]$ls
20170101 checkProcess.sh get_CDS_ERA5_2wrf.py ncl_filedump.log
[#48#haifeng@zhuohf ~/2018/CDS_ERA5/test]$cp get_CDS_ERA5_2wrf.py get_CDS_ERA5_2wrf_global.py
[#49#haifeng@zhuohf ~/2018/CDS_ERA5/test]$vi get_CDS_ERA5_2wrf_global.py -->attachmented
[#50#haifeng@zhuohf ~/2018/CDS_ERA5/test]$python3 get_CDS_ERA5_2wrf_global.py
2018-10-08 14:41:38,120 INFO Sending request to https://cds.climate.copernicus.eu/api/v ... gle-levels
2018-10-08 14:41:39,248 INFO Request is queued
2018-10-08 14:41:42,256 INFO Request is running
2018-10-08 14:41:48,401 INFO Request is completed
2018-10-08 14:41:48,402 INFO Downloading http://136.156.132.110/cache-compute-00 ... 407fd.grib to ERA5-sf-2017-01-01.grib (191.6M)
2018-10-08 14:45:52,885 INFO Download rate 802.5K/s
2018-10-08 14:45:52,885 INFO Sending request to https://cds.climate.copernicus.eu/api/v ... ure-levels
2018-10-08 14:45:53,960 INFO Request is queued
2018-10-08 14:45:55,209 INFO Request is running
2018-10-08 14:46:27,940 INFO Request is completed
2018-10-08 14:46:27,941 INFO Downloading http://136.156.132.105/cache-compute-00 ... 81f5e.grib to ERA5-pl-2017-01-01.grib (1.4G)
}}

2, then do ungrib.exe, got a file named ungrib.log.global
{{
[#65#haifeng@zhuohf ~/WRFV3_9_model/WPS]$./link_grib.csh ~/2018/CDS_ERA5/test/ERA5-sf-2017-01-01.grib
[#66#haifeng@zhuohf ~/WRFV3_9_model/WPS]$ls -ls
总用量 145388
4 drwxr-xr-x. 2 haifeng haifeng 4096 9月 30 14:06 arch
24 drwxrwxr-x. 2 haifeng haifeng 24576 9月 30 14:06 bohai_ERA5_0p125_3km
4 -rwxr-xr-x. 1 haifeng haifeng 1765 9月 30 14:06 clean
8 -rwxr-xr-x. 1 haifeng haifeng 4796 9月 30 14:06 compile
12 -rwxr-xr-x. 1 haifeng haifeng 11514 9月 30 14:06 configure
4 -rw-rw-r--. 1 haifeng haifeng 3527 9月 30 16:43 configure.wps
25744 -rw-rw-r--. 1 haifeng haifeng 26361284 10月 8 14:33 FILE-sf:2017-01-01_00
25744 -rw-rw-r--. 1 haifeng haifeng 26361284 10月 8 14:33 FILE-sf:2017-01-01_06
25744 -rw-rw-r--. 1 haifeng haifeng 26361284 10月 8 14:33 FILE-sf:2017-01-01_12
25744 -rw-rw-r--. 1 haifeng haifeng 26361284 10月 8 14:33 FILE-sf:2017-01-01_18
41908 -rw-rw-r--. 1 haifeng haifeng 42909776 9月 30 14:06 geo_em.d01.nc
4 drwxr-xr-x. 4 haifeng haifeng 4096 9月 30 17:27 geogrid
0 lrwxrwxrwx. 1 haifeng haifeng 23 9月 30 17:27 geogrid.exe -> geogrid/src/geogrid.exe
12 -rw-rw-r--. 1 haifeng haifeng 9303 9月 30 14:06 geogrid.log
0 lrwxrwxrwx. 1 haifeng haifeng 56 10月 8 14:47 GRIBFILE.AAA -> /home/haifeng/2018/CDS_ERA5/test/ERA5-sf-2017-01-01.grib
4 -rwxr-xr-x. 1 haifeng haifeng 1331 9月 30 14:00 link_grib.csh
4 drwxr-xr-x. 3 haifeng haifeng 4096 9月 30 17:27 metgrid
0 lrwxrwxrwx. 1 haifeng haifeng 23 9月 30 17:27 metgrid.exe -> metgrid/src/metgrid.exe
244 -rw-rw-r--. 1 haifeng haifeng 247677 9月 30 14:06 metgrid.log
4 -rw-r--r--. 1 haifeng haifeng 1485 10月 8 14:33 namelist.wps
4 -rw-r--r--. 1 haifeng haifeng 2887 9月 30 14:06 namelist.wps.all_options
4 -rwxr-xr-x. 1 haifeng haifeng 2079 9月 30 14:00 namelist.wps.fire
4 -rw-r--r--. 1 haifeng haifeng 1810 9月 30 14:00 namelist.wps.global
4 -rw-r--r--. 1 haifeng haifeng 656 9月 30 14:00 namelist.wps.nmm
8 -rw-r--r--. 1 haifeng haifeng 6337 9月 30 14:06 README
4 drwxr-xr-x. 4 haifeng haifeng 4096 9月 30 17:27 ungrib
0 lrwxrwxrwx. 1 haifeng haifeng 21 9月 30 17:27 ungrib.exe -> ungrib/src/ungrib.exe
132 -rw-rw-r--. 1 haifeng haifeng 134608 10月 8 14:33 ungrib.log.old
4 drwxr-xr-x. 3 haifeng haifeng 4096 10月 8 14:23 util
0 lrwxrwxrwx. 1 haifeng haifeng 44 10月 8 14:12 Vtable -> ungrib/Variable_Tables/Vtable.ERA-interim.pl
12 -rw-rw-r--. 1 haifeng haifeng 9588 10月 8 14:33 zzz
[#67#haifeng@zhuohf ~/WRFV3_9_model/WPS]$./ungrib.exe > zzz
c_open.c: No such file or directory
Name of source model =>ECMWF
Name of source model =>ECMWF
Name of source model =>ECMWF
Name of source model =>ECMWF
[#68#haifeng@zhuohf ~/WRFV3_9_model/WPS]$vim ungrib.log
[#69#haifeng@zhuohf ~/WRFV3_9_model/WPS]$mv ungrib.log ungrib.log.global
}}
ungrib.log.global:
2018-10-08 14:48:14.104 --- DEBUG: Calling output: 2017-01-01_00
2018-10-08 14:48:14.104 --- Inventory for date = 2017-01-01 00:00:00
2018-10-08 14:48:14.104 --- PRES GEOPT HGT TT UU VV RH DEWPT LANDSEA SOILGEO SOILHGT PSFC PMSL SKINTEMP SEAICE SST SNOW_DEN SNOW_EC SNOW SNOWH ST000007 ST007028 ST028100 ST100289 SM000007 SM007028 SM028100 SM100289
2018-10-08 14:48:14.104 --- -------------------------------------------------
2018-10-08 14:48:14.104 --- 200100 O O X X X O X X O O X X X O X X X O O X X X X X X X X
2018-10-08 14:48:14.104 --- -------------------------------------------------
2018-10-08 14:48:14.104 --- DEBUG: writing in WPS format iunit = 13, map%igrid = 0
2018-10-08 14:48:14.111 --- DEBUG: hdate = 2017-01-01_00:00:00, xfcst = 0.000000
2018-10-08 14:48:14.111 --- DEBUG: map%source = ECMWF, field = TT, units = K
2018-10-08 14:48:14.111 --- DEBUG: Desc = Temperature, level = 200100.000000
2018-10-08 14:48:14.111 --- DEBUG: map%nx = 1440, map%ny = 721
2018-10-08 14:48:14.112 --- DEBUG: max value = 315.139069 , min value = 217.137115

Above is in single level, and the pressure levels data is as the same as that.
zhuohaifeng
 
Posts: 1
Joined: Wed Oct 10, 2018 9:36 pm

Re: Running WRF with Copernicus Data Service ERA5 pressure l

Postby dcvz » Fri Oct 12, 2018 8:11 pm

WPS does not support sub-setted (or regional) Gaussian grids. It only works for global Gaussian grids. It's not a matter of the Vtable, but in the code itself. I don't know if it's possible to modify the WPS source code to make it work, but I suspect that if it were simple it would have already been implemented.
dcvz
 
Posts: 176
Joined: Tue Apr 15, 2008 12:02 am


Return to Working with Various Datasets

Who is online

Users browsing this forum: No registered users and 4 guests