metgrid.exe problems

Having problems running the WRF Preprocessing Sysmte (WPS)?

Re: metgrid.exe problems

Postby AchternDiek » Mon Jul 02, 2012 4:49 am

Hi kurkutes,

Please, check the following two issues:

During ungrib.exe: do you get a confirmation that PSFC and PMSL (and all other fields of interest) are read and extracted successfully?
Did you double check that your entire model domain is covered by your input data?
AchternDiek
 
Posts: 147
Joined: Wed Jun 06, 2012 6:10 am
Location: Norway

Re: metgrid.exe problems

Postby jan » Mon Jul 02, 2012 4:40 pm

kurkutes,
is your input data on a global grid? From wrfhelp:
We only support regular lat/lon grid data, and if it is Gaussian, it needs to be global. We do not support reduced Gaussian grid. If it is not pressure level data, you'll also need to run the utility program calc_ecmwf_p.exe.

And it looks as if PRES is missing in your Vtable.
Are you sure, that all needed variables are read by ungrib in a correct way? On the other hand, it could also be, that you need to specify the missing values of your input data in METGRID.TBL, then metgrid does know how to deal with these points otherwise there might be the warning you have.
jan
 
Posts: 174
Joined: Wed Jun 15, 2011 3:48 am
Location: Kiel, Germany, GEOMAR

Re: metgrid.exe problems

Postby kurkutes » Wed Jul 04, 2012 12:50 pm

Hi AchternDiek and Jan,
Yeah i have checked the ungrib.log, and it didnt report any error.
Well, i have checked that all the values covers the required domain.
And the data is GEM- Regional data, its not global data.

I have been trying to sort out this problem since long, but i didnt get success,

PLease help me guys,

Thanking Yours,
sopan
kurkutes
 
Posts: 7
Joined: Mon May 07, 2012 11:27 am

Re: metgrid.exe problems

Postby AchternDiek » Thu Jul 05, 2012 3:31 am

Hi,

No error in ungrib.exe is a good start.
But did you really see a confirmation that the pressure fields were extracted correctly?

If you don't know how to check that, you could post a section of your ungrib.log.

Basically, for every timestep all fields are listed. Then a 'X' indicates that a variable was read successfully, while a '0' indicates that it was not.

Good Luck,
AchternDiek
AchternDiek
 
Posts: 147
Joined: Wed Jun 06, 2012 6:10 am
Location: Norway

Re: metgrid.exe problems

Postby jan » Thu Jul 05, 2012 5:09 am

Hi sopan,
can you provide me your GEM data on a ftp server somewhere? I could try to read them with your namelist settings. This seems easier for me.
jan
 
Posts: 174
Joined: Wed Jun 15, 2011 3:48 am
Location: Kiel, Germany, GEOMAR

Re: metgrid.exe problems

Postby kurkutes » Thu Jul 05, 2012 3:09 pm

Hi AchternDiek and Jan,
The ftp for downloading the data is
http://dd.weatheroffice.gc.ca/model_gem_regional/high_resolution/grib2/06/
There are separate grib2 fils for each variable t each level. And the data is with 15km resolution. You can see all the information about variables and map here
http://www.weatheroffice.gc.ca/grib/grib2_reg_HR_e.html

Hi AchternDiek
This is the part of my ungrib.exe .
Please see the Vtable and namelist.wps

Inventory for date = 2012-06-07 06:00:00

PRES TT UU VV RH HGT SPECHUMD PRES PRMSL SOILHGT SNOW SKINTEMP
-------------------------------------------------------------------------------
2013.0 O O O O O O O X O O O
2001.0 X X X X O O X O O O O
1015.0 X X X X X X
1000.0 X X X X X X
985.0 X X X X X X
970.0 X X X X X X
950.0 X X X X X X
925.0 X X X X X X
900.0 X X X X X X
875.0 X X X X X X
850.0 X X X X X X
800.0 X X X X X X
750.0 X X X X X X
700.0 X X X X X X
650.0 X X X X X X
600.0 X X X X X X
550.0 X X X X X X
500.0 X X X X X X
450.0 X X X X X X
400.0 X X X X X X
350.0 X X X X X X
300.0 X X X X X X
275.0 X X X X X X
250.0 X X X X X X
225.0 X X X X X X
200.0 X X X X X X
175.0 X X X X X X
150.0 X X X X X X
100.0 X X X X X X
50.0 X X X X X X
-------------------------------------------------------------------------------
Subroutine DATINT: Interpolating 3-d files to fill in any missing data...
Looking for data at time 2012-06-07_00
Found file: FILE:2012-06-07_00
Looking for data at time 2012-06-07_03
Found file: FILE:2012-06-07_03
Looking for data at time 2012-06-07_06
Found file: FILE:2012-06-07_06
End Subroutine DATINT.

**********
Deleting temporary files created by ungrib...
**********

Deleting file: ./PFILE:2012-06-07_00
Deleting file: ./PFILE:2012-06-07_03
Deleting file: ./PFILE:2012-06-07_06

**********
Done deleting temporary files.
**********

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Successful completion of ungrib. !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Here is the Vtable ,

[img]GRIB1| Level| From | To | metgrid | metgrid| metgrid |GRIB2|GRIB2|GRIB2|GRIB2|
Param| Type |Level1|Level2| Name | Units | Description |Discp|Catgy|Param|Level|
-----+------+------+------+----------+---------+-----------------------------------------+-----------------------+
11 | 100 | * | | TT | K | Temperature | 0 | 0 | 0 | 100 |
33 | 100 | * | | UU | m s-1 | U | 0 | 2 | 2 | 100 |
34 | 100 | * | | VV | m s-1 | V | 0 | 2 | 3 | 100 |
52 | 100 | * | | RH | % | Relative Humidity | 0 | 1 | 1 | 100 |
7 | 100 | * | | HGT | m | Height | 0 | 3 | 5 | 100 |
11 | 105 | 2 | | TT | K | Temperature at 2 m | 0 | 0 | 0 | 103 |
| 105 | 2 | | RH | % | Relative Humidity at 2 m (n/a in p file)| 0 | 1 | 0 | 103 |
51 | 105 | 2 | | SPECHUMD | kg kg-1 | | 0 | 1 | 0 | 100 |
33 | 105 | 10 | | UU | m s-1 | U at 10 m | 0 | 2 | 2 | 103 |
34 | 105 | 10 | | VV | m s-1 | V at 10 m | 0 | 2 | 3 | 103 |
1 | 1 | 0 | | PRES | Pa | Surface Pressure | 0 | 3 | 0 | 1 |
2 | 102 | 0 | | PRMSL | Pa | Sea-level Pressure | 0 | 3 | 1 | 101 |
7 | 1 | 0 | | SOILHGT | m | Terrain field of source analysis | 2 | 0 | 7 | 1 |
65 | 1 | 0 | | SNOW | kg m-2 | Water equivalent snow depth | 0 | 1 | 13 | 1 |
11 | 1 | 0 | | SKINTEMP | K | Temperature at ground sfc (from s file) | 0 | 0 | 0 | 1 |
-----+------+------+------+----------+---------+-----------------------------------------+-----------------------+[/img]

In the previous last replies i had also posted my namelist.wps.
I am really appreciates you guys are helping me out to solve my problems.

Thanking Yours,
With Warm regards,
kurkutes
kurkutes
 
Posts: 7
Joined: Mon May 07, 2012 11:27 am

Re: metgrid.exe problems

Postby AchternDiek » Tue Jul 10, 2012 4:34 am

Hi kurkutes,

It seems that your pressure fields are extracted correctly during ungrib.exe.
If your grib data really covers your entire domain and everything runs smoothly during geogrid.exe,
I am running out of ideas for your issue.

Unfortunately I do not have the time to try your files with my WRF.
But perhaps jan does?

Good luck,
AchternDiek
AchternDiek
 
Posts: 147
Joined: Wed Jun 06, 2012 6:10 am
Location: Norway

Re: metgrid.exe problems

Postby jan » Tue Jul 10, 2012 5:23 am

You do not have a land-sea-mask field in your Vtable, which could be a potential problem for ungrib while interpolating data. You should add this field. The other fields do look ok.
Unfortunately I will not have the time to download all the required data from the http site you gave me, this will take too much of my sparse working time.
And I am still not sure if the data you want to use (reg GEM) really works, since wrfhelp says it only works with regular lat/lon or global gaussian data sets.

Required fields:
WRF requires a number of fields as input.

3D Data (e.g. data on pressure levels)
Temperature
U and V components of Wind
Geopotential Height
Relative Humidity
(the code can calculate RH if Specific Humidity is available;
this is controlled in the Vtable)

2D Data
Surface Pressure
Mean Sea Level Pressure
Skin Temperature
2-meter Temperature
2-meter Relative or Specific Humidity
10-meter U and V components of wind

If any masked field is ingested, then a LANDSEA field is recommended
If the WRF model is going to be run with the Noah LSM model, then at least 2 levels of Soil Temperature and Soil Moisture are required.
SST is needed if you plan on updating SST during the model run (often done for very long simulations).
Soil Height is recommended but not required. If Soil Height is supplied, then SOIL Temperatures and TSK can be adjusted to the WRF model terrain height.
Water equivalent snow depth (SNOW) is a nice field to have, but not required.
SEAICE is good to have for a high latitude winter case, but it is not required.

One more point:
You do use "&mod_levs" in your namelist.wps, which is for the program mod_levs.exe. I have never tried this tool. But try ungrib and metgrid without that option! The changes are not large anyway.
jan
 
Posts: 174
Joined: Wed Jun 15, 2011 3:48 am
Location: Kiel, Germany, GEOMAR

Re: metgrid.exe problems

Postby kurkutes » Tue Jul 10, 2012 10:49 am

Hi AchternDiek and Jan,
Thank you very much for your kind responses,
I am not sure if i have made any mistakes while setting up the domains. I have also tried WRFDomain WIzard to set up the domains, but still running out of luck. The GEM-CMC data which i am using as an input to WRF model has stereographic map projection.
And frankly speaking i do not have much experience about how to set up the domain using this projection.
As Jan said that wrfhelp says it only works with regular lat/lon or global gaussian data sets, i am also not sure whether i can use Reg-GEM CMC data-sets as an input to RF model, if we can then how to use it ??
I dont know how to handle this problem.
Well, i will keep trying to solve this problem, i will be very thankful to you guys if in the meanwhile you come up with any solution.
Thanking Yours,

Kurkutes
kurkutes
 
Posts: 7
Joined: Mon May 07, 2012 11:27 am

Re: metgrid.exe problems

Postby jan » Fri Jul 13, 2012 4:01 pm

Information from wrfhelp:
WPS is able to process Polar-stereographic data. Usually WPS will recognize the grids defined by NCEP. If your GRIB file has different definition of grid, then WPS may not recognize it.

How could you run geogrid.exe with
Code: Select all
map_proj = 'polar stereographic',

"polar" would be correct.
Have you tried to inlcude the land-sea mask?
jan
 
Posts: 174
Joined: Wed Jun 15, 2011 3:48 am
Location: Kiel, Germany, GEOMAR

PreviousNext

Return to Runtime Problems

Who is online

Users browsing this forum: No registered users and 1 guest