[ncl-talk] wrf_user_vert_interp and pressure fields

[Trostel, John M.]

When I extract the pressure (p) from a wrf_out file and print out the values along the x-axis for all heights, I get reasonable values at the bottom of the column.
Code: -------

press = wrf_user_getvar(a, "p", it) ; pressure
opt=True
opt at fout="pt_array.txt"
write_matrix(press(:,0,:), "2000f12.3",opt)

Code: -------

Output: ---------
100172.508  100170.219  100169.672  100169.789  100171.039  100171.391  100170.969  100170.711  100170.172  100170.289  100170.516  100170.922  100171.945  100173.680  100176.086  100179.180  100181.867  100185.047  100190.961  100197.172  100202.883  100209.180  100216.094  100223.086  100230.359  100238.453  100245.367  100252.328  100258.961  100264.789  100268.188  100270.227  100272.344  100272.016  100271.617  100271.117  100271.031  100269.219  100266.398  100263.109  100258.648  100254.227  100250.789  100246.883  100243.375  100239.750  100236.336  100233.117  100229.617  100227.500  100226.031  100223.445  100220.094  100217.211  100213.883  100211.539  100208.453  100205.266  100201.211  100197.305  100192.141  100188.383  100184.867  100181.508  100179.055  100177.297  100177.180  100177.172  100178.836  100181.391  100185.430
Output: ---------

When I take that pressure field and apply wrf_user_vert_interp, I get output values that are divided by 100 for the first values of the lowest level.

Code: -------

opts at field_type = "p"
press_interp  = wrf_user_vert_interp( a,press,vert_coord,height_levels,opts)
opt=True
opt at fout="qewr_p_interp_array.txt"
write_matrix(press_interp(:,0,:), "2000f12.3",opt)

Code: -------

Output: ---------

1004.600    1004.593    1004.586    1004.579    1004.571    1004.563    1004.555    1004.716    1004.539    1004.531    1004.523    1004.516    1004.508    1004.501    1004.494    1004.487    1004.482    1004.476    1004.470    1004.465    1004.460    1004.455    1004.451    1004.447    1004.442    1004.438    1004.435    1004.431    1004.430    1004.419    1004.412    1004.407    1004.448    1004.460    1004.477    1004.498    1004.524    1004.538    1004.552    1004.460    1004.446    1004.434    1004.421    1004.406    1004.393    1004.377    1004.380    1004.378    1004.518    1004.547    1004.636    1004.691    1004.531    1004.779  100475.281  100446.828  100466.047  100466.289  100464.000  100467.219  100465.859  100446.922  100457.930  100459.453  100461.820  100464.930  100463.008  100460.797  100468.297  100450.555  100451.336  100449.453  100457.414  100457.813  100448.898  100453.891  100459.102  100460.977  100462.531  100452.063  100451.789

Output: ---------


This doesn't happen for other data...


--
John Trostel
Director - Severe Storms Research Center
Georgia Tech Research Institute
Atlanta, GA 30332-0857
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180926/80494bfb/attachment.html>