[gradsusr] Help extrapolate geopt 1000 mb

[Eric Altshuler]

Hello Angelo,

You still have a slight problem in your namelist. In a namelist record, a slash (/) indicates the end of the record. Anything that comes after the / and before the start of the next record (which is indicated by an ampersand followed by a namelist record name, e.g. &datetime) is ignored by the program. In your case, in the namelist record &interp, the value of extrapolate is specified after the / indicating the end of the &interp namelist record. So your specified value of extrapolate is ignored, and its value defaults to .false. To fix the problem, simply move the line where you specify the value of extrapolate ABOVE the line containing the slash (/):

&interp
...
...
extrapolate = .true.
/
...
...

You should now get fields extrapolated below ground where necessary.

Best regards,

Eric L. Altshuler
Assistant Research Scientist
Center for Ocean-Land-Atmosphere Studies
4041 Powder Mill Road, Suite 302
Calverton, MD 20705-3106
USA

E-mail: ela at cola.iges.org
Phone: (301) 902-1257
Fax: (301) 595-9793

----- Original Message -----
From: "angelo pascualetti" <apascualetti at gmail.com>
To: gradsusr at gradsusr.org
Sent: Sunday, April 11, 2010 12:42:27 PM
Subject: [gradsusr] Help extrapolate geopt 1000 mb


Thank you very much for your answer, but I can not solve the problem. 
Change false to true and still the same problem. 
What I need is to plot the thickness I need geopotencia and the geopotential at 500 mb and 1000 mb. 
That is the problem ... I can not get the geopotential at 1000 mb. 
This is my namelist.ARWpost: 

&datetime 
start_date = '2010-04-11_00:00:00', 
end_date = '2010-04-12_00:00:00', 
interval_seconds = 21600, 
tacc = 0, 
debug_level = 0, 
/ 

&io 
io_form_input = 2, 
input_root_name = '/home/wrf/wrf3.2/WRFV3/run/wrfout_d01_2010-04-11_00:00:00' 
output_root_name = './test' 
plot = 'all_list' 
fields = 'height,geopt,theta,tc,tk,td,td2,rh,rh2,umet,vmet,pressure,u10m,v10m,wdir,wspd,wd10,ws10,slp,mcape,mcin,lcl,lfc,cape,cin,dbz,max_dbz,clfr' 
output_type = 'grads' 
mercator_defs = .true. 
/ 
split_output = .true. 
frames_per_outfile = 2 

output_type = 'grads' 
output_type = 'v5d' 

plot = 'all' 
plot = 'list' 
plot = 'all_list' 
! Below is a list of all available diagnostics 
fields = 'height,geopt,theta,tc,tk,td,td2,rh,rh2,umet,vmet,pressure,u10m,v10m,wdir,wspd,wd10,ws10,slp,mcape,mcin,lcl,lfc,cape,cin,dbz,max_dbz,clfr' 


&interp 
interp_method = 1, 
interp_levels = 1013.,1000.,975.,950.,925.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100., 
/ 
extrapolate = .true. 

interp_method = 1, ! 0 is model levels, -1 is nice height levels, 1 is user specified pressure/height 

interp_levels = 1000.,950.,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100., 
interp_levels = 0.25, 0.50, 0.75, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0, 
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr