[ncl-talk] wrf_user_vert_interp and pressure fields

[Bill Ladwig]

Hi John,

You left some information out of your snippets showing what you set as your
'vert_coord' and 'height_levels' arguments to wrf_user_vert_interp. It
looks like it tried to convert part of your data to hPa, but without seeing
how you called the function, I'm not able to see what's going on. If you
set your 'vert_coord' to pressure, it wants your 'height_levels' to be in
units of hPa, but if you're trying to interpolate pressure as your variable
type and interpolate to pressure in hPa, it might be going off the rails
(and I'm not sure that's what you really wanted to do).

Bill

On Wed, Sep 26, 2018 at 11:07 AM Trostel, John M. <
john.trostel at gtri.gatech.edu> wrote:

> 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
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180926/c0bdb876/attachment.html>