[Wrf-users] ECMWF and WRF3

[Daniel van Dijke]

Hi Ravan,

What's the error you have with the ECMWF model level data? If you show me I might be able to help you.
Since I don't have any problem using it. 

Cheers,

Daniël

PS To Jean-Louis, You don't need to modify calc_ecmwf_p.F, just use mod_lev.exe to remove the unwanted levels computed by calc_ecmwf_p. Make sure to set your "mod_levs" in your namelist:

&mod_levs
 press_pa = 18.0, 39.0, ....... 90.0, 91.0
/

-----Oorspronkelijk bericht-----
Van: wrf-users-bounces at ucar.edu [mailto:wrf-users-bounces at ucar.edu] Namens wrfhelp
Verzonden: dinsdag 24 juni 2008 3:36
Aan: Jean-Louis Monge
CC: wrf-users at ucar.edu
Onderwerp: Re: [Wrf-users] ECMWF and WRF3

If Ravan can give these suggestions a try and let us know, it would be  
helpful. Let us know.
Thanks for helping out,

wrfhelp

On Jun 23, 2008, at 11:59 AM, Jean-Louis Monge wrote:

> On Monday 23 June 2008 17:48:13 Ravan Ahmadov wrote:
>> Dear all,
>>
>> Is someone able to run successfully WPS/WRF3 with ECMWF model level
>> data? The data which worked in the previous versions of WRF/WPS don't
>> work with the latest ones.
>>
>
> Hi Ravan,
>
> It's a long story, but you can do it. Be patient.
> Here is a summary of my own fights against wrf-3 and ECMWF model  
> levels.
> I hope I did not forget too many things.
>
> Fisrst, there is a risk of segfault in geogrid, depending on the  
> compiler you
> use. It was already present in version 2, so you can skip if you are  
> happy
> with it. Here is my diff in WPS/geogrid/src/source_data_module.F :
> ************************
> @@ -632 +633,2 @@
> -            if (ie - is <= 0) ie = 128
> +            if ((ie - is <= 0) .and. (is <= 128)) then
> +               ie = 128
> @@ -633,0 +636 @@
> +            end if
> @@ -663 +669,2 @@
> -            if (ie - is <= 0) ie = 128
> +            if ((ie - is <= 0) .and. (is <= 128)) then
> +               ie = 128
> @@ -664,0 +672 @@
> +            end if
> **************************
>
> Then, the way they compute PSFC in real.F in WRFV3 is different. To  
> make it
> compatible with ECMWF model levels, you have to force the use of  
> subroutine
> sfcprs2 in module_interp_real.F. This is done by adding this line to
> the "domains" section of namelist.input :
>    sfcp_to_sfcp = .true.
>
> Next, if you download lnsp instead of sp from Mars, it works but you  
> have to
> modify ungrib.F :
> *****************************
> @@ -278 +278,2 @@
> -           if (((field == "SST").or.(field == "SKINTEMP")) .and. &
> +           if (((field == "SST").or.(field == "SKINTEMP").or. &
> +                (field == "LOGSFP")) .and.                    & !  
> modif jlm
> *****************************
> Finally, I have modified calc_ecmwf_p.F to make it possible to use  
> only a
> subset of the 96  ECMWF model levels. It saves disk space and  
> download time.
> If some people are interested, I will mail them the modified code.
>
> Last : do not use "ucmcall = 1,1,1" in the "physics" section of  
> namelist.input
> if you use two-way nesting with "input_from_file  
> = .true., .true., .true."
> For me, it fills all the fields with NaNs at the first simulation  
> hour.
>
> I hope this helps.
>
> -- 
> Jean-Louis Monge
> Ingénieur-conseil
> 09 50 51 14 43
> jeanlouis.monge at free.fr

_______________________________________________
Wrf-users mailing list
Wrf-users at ucar.edu
http://mailman.ucar.edu/mailman/listinfo/wrf-users


-- 

This e-mail is from Meteo Consult B.V., a MeteoGroup company. For more information, see http://www.weer.nl/gebruiksvoorwaarden.

This e-mail may contain confidential information. Only the addressee is permitted to read, copy, distribute or otherwise use this e-mail or any attachments. If you have received it in error, please contact the sender immediately. Any opinion expressed in this e-mail is personal to the sender and may not reflect the opinion of MeteoGroup.

Any e-mail reply to this address may be subject to interception or monitoring for operational reasons or for lawful business practices.