Help with Skew-T charts

[HathewayWill]

Hello everyone,

Working on a new Skew-T chart that will read in WRF output and create model skew-t.

So I'm having some strange issues.

First off I am following and modifying these tutorials.
https://unidata.github.io/MetPy/latest/examples/plots/Skew-T_Layout.html

Now I have run the 18z wrf and looked at the 0z Model skew-t vs the real skew-t taken at KDRT

Actual Observation:

python code that is in the correct format according to tutorial but wrong chart is produced.

Python code that is wrong according to the tutoiral but the chart is correct.

Now I have attached the two different python codes here:

Skew-T_Issues.zip

The issue seems to be with the lat lon functions in my python code:

This is the section I think it is happening:

• Correct format wrong chart

# Get variables needed for Skew-T
    p1 = wrf.getvar(ncfile, "pres", timeidx=0, units='hPa')
    T1 = wrf.getvar(ncfile, "temp", timeidx=0, units='degC')
    Td1 = wrf.getvar(ncfile, "td", timeidx=0, units='degC')
    u1 = wrf.getvar(ncfile, "ua", timeidx=0, units='kt')
    v1 = wrf.getvar(ncfile, "va", timeidx=0, units='kt')

    p = p1[:, x_y[0], x_y[1]] * units.hPa
    T = T1[:, x_y[0], x_y[1]] * units.degC
    Td = Td1[:, x_y[0], x_y[1]] * units.degC
    u = u1[:, x_y[0], x_y[1]] * units("knots")
    v = v1[:, x_y[0], x_y[1]] * units("knots")

• Wrong Format right chart

 # Get variables needed for Skew-T
    p1 = wrf.getvar(ncfile, "pres", timeidx=0, units='hPa')
    T1 = wrf.getvar(ncfile, "temp", timeidx=0, units='degC')
    Td1 = wrf.getvar(ncfile, "td", timeidx=0, units='degC')
    u1 = wrf.getvar(ncfile, "ua", timeidx=0, units='kt')
    v1 = wrf.getvar(ncfile, "va", timeidx=0, units='kt')

    p = p1[:, x_y[1], x_y[0]] * units.hPa
    T = T1[:, x_y[1], x_y[0]] * units.degC
    Td = Td1[:, x_y[1], x_y[0]] * units.degC
    u = u1[:, x_y[1], x_y[0]] * units("knots")
    v = v1[:, x_y[1], x_y[0]] * units("knots")

Any ideas?

[dopplershift]

I don't understand what makes you say the bottom code is "wrong format" if it produces the correct chart?

As we indicated in the email support, everything seems to be working correctly based on how the wrf-python API is working.

[HathewayWill]

sorry for the confusion.

[HathewayWill]

The issue was that I thought the netcdf data was in lat, long format but it was in long, lat format.