[Met_help] [rt.rap.ucar.edu #100132] History for unexpected results from point_stat with NEAREST GRIDPOINT method
John Halley Gotway via RT
met_help at ucar.edu
Mon Jul 12 11:23:00 MDT 2021
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear Met-Help,
I am testing out the point_stat tool (V10.0.0) for 1 station pair, and trying to see if it is grabbing the correct value of windspeed from the nearest gridpoint in my input forecast netcdf grid. In my ascii obs text I purposefully put the station lat/lon as one of the gridpoints in my forecast field for lat lon indexes (starting with 0 in the ncl fashion) of 218,140. The lat/lon that this corresponds to is (37.61275, -122.3618). The 10-m windspeed at this point is 3.157048 m/s, however the point_stat*_mpr.txt results of the matched pair has the windspeed as 3.57623 m/s, which is not correct. In fact, none of the surrounding points have that value either. Why is point stat giving me a different value? I don't think that it would be doing any kind of vertical interpolation, but not sure. My expectation is that it should not be. I attached my PointStatConfig file, and the other info I provided below that may help.
Thanks,
A.J. Eiserloh - Data Scientist
PG&E - Meteorology Systems and Analytics
Applied Technology Services
3400 Crow Canyon Rd., San Ramon, CA 94583-1393
925-307-4492
* Here is my pre ascii2nc ascii input columns for the point:
ADPSFC KSFO 20210301_000000 37.61275 -122.3618 0.0 32 0 10.0 NA 11.318018760716
* The *mpr.txt results for the one matched pair:
OBS_LAT OBS_LON OBS_LVL OBS_ELV FCST OBS OBS_QC CLIMO_MEAN CLIMO_STDEV CLIMO_CDF
37.61275 -122.3618 NA 10 3.57623 11.31802 NA NA NA NA
* ncdump of my netcdf file:
netcdf wrfsfc_d03_20210301_0000_f00000-cfcomp {
dimensions:
Time = UNLIMITED ; // (1 currently)
south_north = 480 ;
west_east = 396 ;
DateStrLen = 19 ;
soil_layers_stag = 4 ;
variables:
float WSPD_10M(Time, south_north, west_east) ;
WSPD_10M:_FillValue = 1.e+20f ;
WSPD_10M:FieldType = 104 ;
WSPD_10M:MemoryOrder = "XY " ;
WSPD_10M:coordinates = "lon lat" ;
WSPD_10M:description = "U at 10 M" ;
WSPD_10M:level = "10 m" ;
WSPD_10M:long_name = "Wind Speed at 10 m" ;
WSPD_10M:stagger = "" ;
WSPD_10M:standard_name = "wind_speed" ;
WSPD_10M:units = "m/s" ;
char Times(Time, DateStrLen) ;
float XLAT(Time, south_north, west_east) ;
XLAT:FieldType = 104 ;
XLAT:MemoryOrder = "XY " ;
XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
XLAT:units = "degrees_north" ;
XLAT:stagger = "" ;
XLAT:coordinates = "XLONG XLAT" ;
XLAT:long_name = "Latitude" ;
XLAT:standard_name = "latitude" ;
float XLONG(Time, south_north, west_east) ;
XLONG:FieldType = 104 ;
XLONG:MemoryOrder = "XY " ;
XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;
XLONG:units = "degrees_east" ;
XLONG:stagger = "" ;
XLONG:coordinates = "XLONG XLAT" ;
XLONG:long_name = "Longitude" ;
XLONG:standard_name = "longitude" ;
float Q2(Time, south_north, west_east) ;
Q2:FieldType = 104 ;
Q2:MemoryOrder = "XY " ;
Q2:description = "QV at 2 M" ;
Q2:units = "kg kg-1" ;
Q2:stagger = "" ;
Q2:coordinates = "lon lat" ;
Q2:long_name = "Specific Humidity at 2 m" ;
Q2:standard_name = "specific_humidity" ;
Q2:level = "2 m" ;
Q2:_FillValue = 1.e+20f ;
float T2(Time, south_north, west_east) ;
T2:FieldType = 104 ;
T2:MemoryOrder = "XY " ;
T2:description = "TEMP at 2 M" ;
T2:units = "K" ;
T2:stagger = "" ;
T2:coordinates = "lon lat" ;
T2:long_name = "Air Temperature at 2 m" ;
T2:standard_name = "air_temperature" ;
T2:level = "2 m" ;
T2:_FillValue = 1.e+20f ;
float PSFC(Time, south_north, west_east) ;
PSFC:FieldType = 104 ;
PSFC:MemoryOrder = "XY " ;
PSFC:description = "SFC PRESSURE" ;
PSFC:units = "Pa" ;
PSFC:stagger = "" ;
PSFC:coordinates = "lon lat" ;
PSFC:long_name = "Air Pressure at the Surface" ;
PSFC:standard_name = "surface_pressure" ;
PSFC:level = "surface" ;
PSFC:_FillValue = 1.e+20f ;
float U10(Time, south_north, west_east) ;
U10:FieldType = 104 ;
U10:MemoryOrder = "XY " ;
U10:description = "U at 10 M" ;
U10:units = "m/s" ;
U10:stagger = "" ;
U10:coordinates = "lon lat" ;
U10:long_name = "u-Component of wind at 10 m (Earth)" ;
U10:standard_name = "eastward_wind" ;
U10:level = "10 m" ;
U10:_FillValue = 1.e+20f ;
float V10(Time, south_north, west_east) ;
V10:FieldType = 104 ;
V10:MemoryOrder = "XY " ;
V10:description = "V at 10 M" ;
V10:units = "m/s" ;
V10:stagger = "" ;
V10:coordinates = "lon lat" ;
V10:long_name = "v-Component of wind at 10 m (Earth)" ;
V10:standard_name = "northward_wind" ;
V10:level = "10 m" ;
V10:_FillValue = 1.e+20f ;
double time(Time) ;
time:standard_name = "time" ;
time:long_name = "Time" ;
time:calendar = "standard" ;
time:units = "hours since 1970-01-01 00:00:00" ;
float lat(south_north, west_east) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float lon(south_north, west_east) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
// global attributes:
:TITLE = " OUTPUT FROM WRF V4.1.1 MODEL" ;
:START_DATE = "2021-03-01_00:00:00" ;
:WEST-EAST_GRID_DIMENSION = 397 ;
:SOUTH-NORTH_GRID_DIMENSION = 481 ;
:BOTTOM-TOP_GRID_DIMENSION = 51 ;
:DX = 2000.f ;
:DY = 2000.f ;
:AERCU_OPT = 0 ;
:AERCU_FCT = 1.f ;
:IDEAL_CASE = 0 ;
:DIFF_6TH_SLOPEOPT = 0 ;
:AUTO_LEVELS_OPT = 2 ;
:DIFF_6TH_THRESH = 0.1f ;
:DZBOT = 50.f ;
:DZSTRETCH_S = 1.3f ;
:DZSTRETCH_U = 1.1f ;
:GRIDTYPE = "C" ;
:DIFF_OPT = 2 ;
:KM_OPT = 4 ;
:DAMP_OPT = 3 ;
:DAMPCOEF = 0.2f ;
:KHDIF = 0.f ;
:KVDIF = 0.f ;
:MP_PHYSICS = 8 ;
:RA_LW_PHYSICS = 4 ;
:RA_SW_PHYSICS = 4 ;
:SF_SFCLAY_PHYSICS = 5 ;
:SF_SURFACE_PHYSICS = 4 ;
:BL_PBL_PHYSICS = 6 ;
:CU_PHYSICS = 0 ;
:SF_LAKE_PHYSICS = 0 ;
:SURFACE_INPUT_SOURCE = 3 ;
:SST_UPDATE = 0 ;
:GRID_FDDA = 0 ;
:GFDDA_INTERVAL_M = 0 ;
:GFDDA_END_H = 0 ;
:GRID_SFDDA = 0 ;
:SGFDDA_INTERVAL_M = 0 ;
:SGFDDA_END_H = 0 ;
:HYPSOMETRIC_OPT = 2 ;
:USE_THETA_M = 1 ;
:GWD_OPT = 0 ;
:SF_URBAN_PHYSICS = 0 ;
:SF_SURFACE_MOSAIC = 0 ;
:SF_OCEAN_PHYSICS = 0 ;
:WEST-EAST_PATCH_START_UNSTAG = 1 ;
:WEST-EAST_PATCH_END_UNSTAG = 396 ;
:WEST-EAST_PATCH_START_STAG = 1 ;
:WEST-EAST_PATCH_END_STAG = 397 ;
:SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
:SOUTH-NORTH_PATCH_END_UNSTAG = 480 ;
:SOUTH-NORTH_PATCH_START_STAG = 1 ;
:SOUTH-NORTH_PATCH_END_STAG = 481 ;
:BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
:BOTTOM-TOP_PATCH_END_UNSTAG = 50 ;
:BOTTOM-TOP_PATCH_START_STAG = 1 ;
:BOTTOM-TOP_PATCH_END_STAG = 51 ;
:GRID_ID = 3 ;
:PARENT_ID = 2 ;
:I_PARENT_START = 93 ;
:J_PARENT_START = 75 ;
:PARENT_GRID_RATIO = 3 ;
:DT = 8.f ;
:CEN_LAT = 38.00961f ;
:CEN_LON = -121.0455f ;
:TRUELAT1 = 31.f ;
:TRUELAT2 = 45.38f ;
:MOAD_CEN_LAT = 38.19999f ;
:STAND_LON = -121.08f ;
:POLE_LAT = 90.f ;
:POLE_LON = 0.f ;
:GMT = 0.f ;
:JULYR = 2021 ;
:JULDAY = 60 ;
:MAP_PROJ = 1 ;
:MAP_PROJ_CHAR = "Lambert Conformal" ;
:MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
:NUM_LAND_CAT = 21 ;
:ISWATER = 17 ;
:ISLAKE = 21 ;
:ISICE = 15 ;
:ISURBAN = 13 ;
:ISOILWATER = 14 ;
:HYBRID_OPT = 2 ;
:ETAC = 0.2f ;
:file_type = "NETCDF_NCCF" ;
:title = "wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
:Conventions = "CF 1.6, Standard Name Table v19" ;
:source = "wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
:notes = "" ;
:created_by = "AJEB" ;
:institution = "PG&E" ;
:creation_date = "Mon Jun 7 17:59:23 UTC 2021" ;
:history = "Mon Jun 7 17:59:26 2021: ncap2 -O -s T2F=float(T2F) wrfsfc_d03_20210301_0000_f00000-cfcomp.nc wrfsfc_d03_20210301_0000_f00000-cfcomp.nc\n",
"Mon Jun 7 17:59:26 2021: ncap2 -O -s T2F=((T2-273.15)*(9.0/5.0))+32.0 wrfsfc_d03_20210301_0000_f00000-cfcomp.nc wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
:NCO = "netCDF Operators version 4.9.8 (Homepage = http://nco.sf.net, Code = http://github.com/nco/nco)" ;
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: unexpected results from point_stat with NEAREST GRIDPOINT method
From: John Halley Gotway
Time: Mon Jun 07 15:47:21 2021
AJ,
I see that you're not getting the expected results from the nearest
neighbor interpolation method in point-stat. There's nothing I can see
in
the Point-Stat config file you sent that is an obvious problem.
If I had the actual model/obs data available to me, I would start by
running the plot_data_plane tool to plot the 10-m wind data from your
NetCDF file.
*plot_data_plane wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
<http://wrfsfc_d03_20210301_0000_f00000-cfcomp.nc> WSPD_10M.ps
'name="WSPD_10M"; level="(0,*,*)";'*
And I'd check to make sure that the data isn't upside-down or
backwards. If
that looks good, I might also run the plot_point_obs tool. You can
actually
plot WSPD_10M field from above as a base image and then overlay your
point.
Check to make sure that point is really in the right spot on the
grid...
confirm that sign of the longitude is correct... etc.
If that all checks out, next I'd run point_stat through a debugger to
see
line-by-line what's going on.
But I'd need sample data to debug in this way. If you'd like to
provide
that, you could post it to our anonymous FTP site following these
instructions:
https://dtcenter.org/community-code/model-evaluation-tools-met/met-
help-desk#ftp
Thanks,
John Halley Gotway
On Mon, Jun 7, 2021 at 12:45 PM Eiserloh Jr., A.J. via RT
<met_help at ucar.edu>
wrote:
>
> Mon Jun 07 12:45:45 2021: Request 100132 was acted upon.
> Transaction: Ticket created by AJEB at pge.com
> Queue: met_help
> Subject: unexpected results from point_stat with NEAREST
GRIDPOINT
> method
> Owner: Nobody
> Requestors: AJEB at pge.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100132 >
>
>
> Dear Met-Help,
>
> I am testing out the point_stat tool (V10.0.0) for 1 station pair,
and
> trying to see if it is grabbing the correct value of windspeed from
the
> nearest gridpoint in my input forecast netcdf grid. In my ascii obs
text I
> purposefully put the station lat/lon as one of the gridpoints in my
> forecast field for lat lon indexes (starting with 0 in the ncl
fashion) of
> 218,140. The lat/lon that this corresponds to is (37.61275,
-122.3618). The
> 10-m windspeed at this point is 3.157048 m/s, however the
> point_stat*_mpr.txt results of the matched pair has the windspeed
as
> 3.57623 m/s, which is not correct. In fact, none of the surrounding
points
> have that value either. Why is point stat giving me a different
value? I
> don't think that it would be doing any kind of vertical
interpolation, but
> not sure. My expectation is that it should not be. I attached my
> PointStatConfig file, and the other info I provided below that may
help.
>
> Thanks,
>
> A.J. Eiserloh - Data Scientist
> PG&E - Meteorology Systems and Analytics
> Applied Technology Services
> 3400 Crow Canyon Rd., San Ramon, CA 94583-1393
> 925-307-4492
>
>
>
>
> * Here is my pre ascii2nc ascii input columns for the point:
>
> ADPSFC KSFO 20210301_000000 37.61275 -122.3618 0.0 32 0 10.0 NA
> 11.318018760716
>
>
> * The *mpr.txt results for the one matched pair:
>
> OBS_LAT OBS_LON OBS_LVL OBS_ELV FCST
> OBS OBS_QC CLIMO_MEAN CLIMO_STDEV CLIMO_CDF
> 37.61275 -122.3618 NA 10
> 3.57623 11.31802 NA NA
> NA NA
>
>
> * ncdump of my netcdf file:
>
> netcdf wrfsfc_d03_20210301_0000_f00000-cfcomp {
> dimensions:
> Time = UNLIMITED ; // (1 currently)
> south_north = 480 ;
> west_east = 396 ;
> DateStrLen = 19 ;
> soil_layers_stag = 4 ;
> variables:
> float WSPD_10M(Time, south_north, west_east) ;
> WSPD_10M:_FillValue = 1.e+20f ;
> WSPD_10M:FieldType = 104 ;
> WSPD_10M:MemoryOrder = "XY " ;
> WSPD_10M:coordinates = "lon lat" ;
> WSPD_10M:description = "U at 10 M" ;
> WSPD_10M:level = "10 m" ;
> WSPD_10M:long_name = "Wind Speed at 10 m" ;
> WSPD_10M:stagger = "" ;
> WSPD_10M:standard_name = "wind_speed" ;
> WSPD_10M:units = "m/s" ;
> char Times(Time, DateStrLen) ;
> float XLAT(Time, south_north, west_east) ;
> XLAT:FieldType = 104 ;
> XLAT:MemoryOrder = "XY " ;
> XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
> XLAT:units = "degrees_north" ;
> XLAT:stagger = "" ;
> XLAT:coordinates = "XLONG XLAT" ;
> XLAT:long_name = "Latitude" ;
> XLAT:standard_name = "latitude" ;
> float XLONG(Time, south_north, west_east) ;
> XLONG:FieldType = 104 ;
> XLONG:MemoryOrder = "XY " ;
> XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;
> XLONG:units = "degrees_east" ;
> XLONG:stagger = "" ;
> XLONG:coordinates = "XLONG XLAT" ;
> XLONG:long_name = "Longitude" ;
> XLONG:standard_name = "longitude" ;
> float Q2(Time, south_north, west_east) ;
> Q2:FieldType = 104 ;
> Q2:MemoryOrder = "XY " ;
> Q2:description = "QV at 2 M" ;
> Q2:units = "kg kg-1" ;
> Q2:stagger = "" ;
> Q2:coordinates = "lon lat" ;
> Q2:long_name = "Specific Humidity at 2 m" ;
> Q2:standard_name = "specific_humidity" ;
> Q2:level = "2 m" ;
> Q2:_FillValue = 1.e+20f ;
> float T2(Time, south_north, west_east) ;
> T2:FieldType = 104 ;
> T2:MemoryOrder = "XY " ;
> T2:description = "TEMP at 2 M" ;
> T2:units = "K" ;
> T2:stagger = "" ;
> T2:coordinates = "lon lat" ;
> T2:long_name = "Air Temperature at 2 m" ;
> T2:standard_name = "air_temperature" ;
> T2:level = "2 m" ;
> T2:_FillValue = 1.e+20f ;
> float PSFC(Time, south_north, west_east) ;
> PSFC:FieldType = 104 ;
> PSFC:MemoryOrder = "XY " ;
> PSFC:description = "SFC PRESSURE" ;
> PSFC:units = "Pa" ;
> PSFC:stagger = "" ;
> PSFC:coordinates = "lon lat" ;
> PSFC:long_name = "Air Pressure at the Surface" ;
> PSFC:standard_name = "surface_pressure" ;
> PSFC:level = "surface" ;
> PSFC:_FillValue = 1.e+20f ;
> float U10(Time, south_north, west_east) ;
> U10:FieldType = 104 ;
> U10:MemoryOrder = "XY " ;
> U10:description = "U at 10 M" ;
> U10:units = "m/s" ;
> U10:stagger = "" ;
> U10:coordinates = "lon lat" ;
> U10:long_name = "u-Component of wind at 10 m
(Earth)" ;
> U10:standard_name = "eastward_wind" ;
> U10:level = "10 m" ;
> U10:_FillValue = 1.e+20f ;
> float V10(Time, south_north, west_east) ;
> V10:FieldType = 104 ;
> V10:MemoryOrder = "XY " ;
> V10:description = "V at 10 M" ;
> V10:units = "m/s" ;
> V10:stagger = "" ;
> V10:coordinates = "lon lat" ;
> V10:long_name = "v-Component of wind at 10 m
(Earth)" ;
> V10:standard_name = "northward_wind" ;
> V10:level = "10 m" ;
> V10:_FillValue = 1.e+20f ;
> double time(Time) ;
> time:standard_name = "time" ;
> time:long_name = "Time" ;
> time:calendar = "standard" ;
> time:units = "hours since 1970-01-01 00:00:00" ;
> float lat(south_north, west_east) ;
> lat:long_name = "Latitude" ;
> lat:standard_name = "latitude" ;
> lat:units = "degrees_north" ;
> float lon(south_north, west_east) ;
> lon:long_name = "Longitude" ;
> lon:standard_name = "longitude" ;
> lon:units = "degrees_east" ;
>
> // global attributes:
> :TITLE = " OUTPUT FROM WRF V4.1.1 MODEL" ;
> :START_DATE = "2021-03-01_00:00:00" ;
> :WEST-EAST_GRID_DIMENSION = 397 ;
> :SOUTH-NORTH_GRID_DIMENSION = 481 ;
> :BOTTOM-TOP_GRID_DIMENSION = 51 ;
> :DX = 2000.f ;
> :DY = 2000.f ;
> :AERCU_OPT = 0 ;
> :AERCU_FCT = 1.f ;
> :IDEAL_CASE = 0 ;
> :DIFF_6TH_SLOPEOPT = 0 ;
> :AUTO_LEVELS_OPT = 2 ;
> :DIFF_6TH_THRESH = 0.1f ;
> :DZBOT = 50.f ;
> :DZSTRETCH_S = 1.3f ;
> :DZSTRETCH_U = 1.1f ;
> :GRIDTYPE = "C" ;
> :DIFF_OPT = 2 ;
> :KM_OPT = 4 ;
> :DAMP_OPT = 3 ;
> :DAMPCOEF = 0.2f ;
> :KHDIF = 0.f ;
> :KVDIF = 0.f ;
> :MP_PHYSICS = 8 ;
> :RA_LW_PHYSICS = 4 ;
> :RA_SW_PHYSICS = 4 ;
> :SF_SFCLAY_PHYSICS = 5 ;
> :SF_SURFACE_PHYSICS = 4 ;
> :BL_PBL_PHYSICS = 6 ;
> :CU_PHYSICS = 0 ;
> :SF_LAKE_PHYSICS = 0 ;
> :SURFACE_INPUT_SOURCE = 3 ;
> :SST_UPDATE = 0 ;
> :GRID_FDDA = 0 ;
> :GFDDA_INTERVAL_M = 0 ;
> :GFDDA_END_H = 0 ;
> :GRID_SFDDA = 0 ;
> :SGFDDA_INTERVAL_M = 0 ;
> :SGFDDA_END_H = 0 ;
> :HYPSOMETRIC_OPT = 2 ;
> :USE_THETA_M = 1 ;
> :GWD_OPT = 0 ;
> :SF_URBAN_PHYSICS = 0 ;
> :SF_SURFACE_MOSAIC = 0 ;
> :SF_OCEAN_PHYSICS = 0 ;
> :WEST-EAST_PATCH_START_UNSTAG = 1 ;
> :WEST-EAST_PATCH_END_UNSTAG = 396 ;
> :WEST-EAST_PATCH_START_STAG = 1 ;
> :WEST-EAST_PATCH_END_STAG = 397 ;
> :SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
> :SOUTH-NORTH_PATCH_END_UNSTAG = 480 ;
> :SOUTH-NORTH_PATCH_START_STAG = 1 ;
> :SOUTH-NORTH_PATCH_END_STAG = 481 ;
> :BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
> :BOTTOM-TOP_PATCH_END_UNSTAG = 50 ;
> :BOTTOM-TOP_PATCH_START_STAG = 1 ;
> :BOTTOM-TOP_PATCH_END_STAG = 51 ;
> :GRID_ID = 3 ;
> :PARENT_ID = 2 ;
> :I_PARENT_START = 93 ;
> :J_PARENT_START = 75 ;
> :PARENT_GRID_RATIO = 3 ;
> :DT = 8.f ;
> :CEN_LAT = 38.00961f ;
> :CEN_LON = -121.0455f ;
> :TRUELAT1 = 31.f ;
> :TRUELAT2 = 45.38f ;
> :MOAD_CEN_LAT = 38.19999f ;
> :STAND_LON = -121.08f ;
> :POLE_LAT = 90.f ;
> :POLE_LON = 0.f ;
> :GMT = 0.f ;
> :JULYR = 2021 ;
> :JULDAY = 60 ;
> :MAP_PROJ = 1 ;
> :MAP_PROJ_CHAR = "Lambert Conformal" ;
> :MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
> :NUM_LAND_CAT = 21 ;
> :ISWATER = 17 ;
> :ISLAKE = 21 ;
> :ISICE = 15 ;
> :ISURBAN = 13 ;
> :ISOILWATER = 14 ;
> :HYBRID_OPT = 2 ;
> :ETAC = 0.2f ;
> :file_type = "NETCDF_NCCF" ;
> :title = "wrfsfc_d03_20210301_0000_f00000-cfcomp.nc"
;
> :Conventions = "CF 1.6, Standard Name Table v19" ;
> :source = "wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc" ;
> :notes = "" ;
> :created_by = "AJEB" ;
> :institution = "PG&E" ;
> :creation_date = "Mon Jun 7 17:59:23 UTC 2021" ;
> :history = "Mon Jun 7 17:59:26 2021: ncap2 -O -s
> T2F=float(T2F) wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
> wrfsfc_d03_20210301_0000_f00000-cfcomp.nc\n",
> "Mon Jun 7 17:59:26 2021: ncap2 -O -s
> T2F=((T2-273.15)*(9.0/5.0))+32.0 wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc
> wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
> :NCO = "netCDF Operators version 4.9.8 (Homepage =
> http://nco.sf.net, Code = http://github.com/nco/nco)" ;
>
>
>
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #100132] unexpected results from point_stat with NEAREST GRIDPOINT method
From: Eiserloh Jr., A.J.
Time: Tue Jun 08 17:03:59 2021
Hi John,
I am unable to ftp to send data to the ftp address unless I get a
firewall exception request from my company's IT, and that will take
some time. I tried plotting with the plot_data_plane command but get
another error that is not helpful.
Thanks,
-AJ
[root at 698df1a69b34 nccompli_test]# plot_data_plane
wrfsfc_d03_20210301_0000_f00000-cfcomp.nc out.ps WSPD_10M
ERROR :
ERROR : yyerror() -> syntax error in file "/tmp/met_config_273_0"
ERROR :
ERROR : line = 2
ERROR :
ERROR : column = 8
ERROR :
ERROR : text = "(nul)"
ERROR :
ERROR :
ERROR : (nul)
ERROR : (nul)
ERROR :
-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Monday, June 7, 2021 2:47 PM
To: Eiserloh Jr., A.J. <AJEB at pge.com>
Subject: Re: [rt.rap.ucar.edu #100132] unexpected results from
point_stat with NEAREST GRIDPOINT method
*****CAUTION: This email was sent from an EXTERNAL source. Think
before clicking links or opening attachments.*****
AJ,
I see that you're not getting the expected results from the nearest
neighbor interpolation method in point-stat. There's nothing I can see
in the Point-Stat config file you sent that is an obvious problem.
If I had the actual model/obs data available to me, I would start by
running the plot_data_plane tool to plot the 10-m wind data from your
NetCDF file.
*plot_data_plane wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
<https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwrfsfc_d03_20210301_0000_f00000-
cfcomp.nc%2F&data=04%7C01%7CAJEB%40pge.com%7Cae1ed37b041e47bfa3a708d929fdd547%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C637586992454384982%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=U3uNNDnM%2FqUKIC6UpHBuSjEeM%2BbtEQaP5v%2Bvuly%2FxiU%3D&reserved=0>
WSPD_10M.ps 'name="WSPD_10M"; level="(0,*,*)";'*
And I'd check to make sure that the data isn't upside-down or
backwards. If that looks good, I might also run the plot_point_obs
tool. You can actually plot WSPD_10M field from above as a base image
and then overlay your point.
Check to make sure that point is really in the right spot on the
grid...
confirm that sign of the longitude is correct... etc.
If that all checks out, next I'd run point_stat through a debugger to
see line-by-line what's going on.
But I'd need sample data to debug in this way. If you'd like to
provide that, you could post it to our anonymous FTP site following
these
instructions:
https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdtcenter.org%2Fcommunity-
code%2Fmodel-evaluation-tools-met%2Fmet-help-
desk%23ftp&data=04%7C01%7CAJEB%40pge.com%7Cae1ed37b041e47bfa3a708d929fdd547%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C637586992454384982%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=CqU3jPPMQuJrjK2vQoPEHnTUo9Erb3hQvddTVrB3LC8%3D&reserved=0
Thanks,
John Halley Gotway
On Mon, Jun 7, 2021 at 12:45 PM Eiserloh Jr., A.J. via RT
<met_help at ucar.edu>
wrote:
>
> Mon Jun 07 12:45:45 2021: Request 100132 was acted upon.
> Transaction: Ticket created by AJEB at pge.com
> Queue: met_help
> Subject: unexpected results from point_stat with NEAREST
> GRIDPOINT method
> Owner: Nobody
> Requestors: AJEB at pge.com
> Status: new
> Ticket <URL:
>
https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Frt.r
>
ap.ucar.edu%2Frt%2FTicket%2FDisplay.html%3Fid%3D100132&data=04%7C0
>
1%7CAJEB%40pge.com%7Cae1ed37b041e47bfa3a708d929fdd547%7C44ae661aece641
>
aabc967c2c85a08941%7C0%7C0%7C637586992454384982%7CUnknown%7CTWFpbGZsb3
>
d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7
>
C1000&sdata=LuRAF64SY4mRnzYWzWw5YN8%2FKHVDJTkwzWj%2BolIhDmA%3D&
> ;reserved=0 >
>
>
> Dear Met-Help,
>
> I am testing out the point_stat tool (V10.0.0) for 1 station pair,
and
> trying to see if it is grabbing the correct value of windspeed from
> the nearest gridpoint in my input forecast netcdf grid. In my ascii
> obs text I purposefully put the station lat/lon as one of the
> gridpoints in my forecast field for lat lon indexes (starting with 0
> in the ncl fashion) of 218,140. The lat/lon that this corresponds to
> is (37.61275, -122.3618). The 10-m windspeed at this point is
3.157048
> m/s, however the point_stat*_mpr.txt results of the matched pair
has
> the windspeed as
> 3.57623 m/s, which is not correct. In fact, none of the surrounding
> points have that value either. Why is point stat giving me a
different
> value? I don't think that it would be doing any kind of vertical
> interpolation, but not sure. My expectation is that it should not
be.
> I attached my PointStatConfig file, and the other info I provided
below that may help.
>
> Thanks,
>
> A.J. Eiserloh - Data Scientist
> PG&E - Meteorology Systems and Analytics Applied Technology Services
> 3400 Crow Canyon Rd., San Ramon, CA 94583-1393
> 925-307-4492
>
>
>
>
> * Here is my pre ascii2nc ascii input columns for the point:
>
> ADPSFC KSFO 20210301_000000 37.61275 -122.3618 0.0 32 0 10.0 NA
> 11.318018760716
>
>
> * The *mpr.txt results for the one matched pair:
>
> OBS_LAT OBS_LON OBS_LVL OBS_ELV FCST
> OBS OBS_QC CLIMO_MEAN CLIMO_STDEV CLIMO_CDF
> 37.61275 -122.3618 NA 10
> 3.57623 11.31802 NA NA
> NA NA
>
>
> * ncdump of my netcdf file:
>
> netcdf wrfsfc_d03_20210301_0000_f00000-cfcomp {
> dimensions:
> Time = UNLIMITED ; // (1 currently)
> south_north = 480 ;
> west_east = 396 ;
> DateStrLen = 19 ;
> soil_layers_stag = 4 ;
> variables:
> float WSPD_10M(Time, south_north, west_east) ;
> WSPD_10M:_FillValue = 1.e+20f ;
> WSPD_10M:FieldType = 104 ;
> WSPD_10M:MemoryOrder = "XY " ;
> WSPD_10M:coordinates = "lon lat" ;
> WSPD_10M:description = "U at 10 M" ;
> WSPD_10M:level = "10 m" ;
> WSPD_10M:long_name = "Wind Speed at 10 m" ;
> WSPD_10M:stagger = "" ;
> WSPD_10M:standard_name = "wind_speed" ;
> WSPD_10M:units = "m/s" ;
> char Times(Time, DateStrLen) ;
> float XLAT(Time, south_north, west_east) ;
> XLAT:FieldType = 104 ;
> XLAT:MemoryOrder = "XY " ;
> XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
> XLAT:units = "degrees_north" ;
> XLAT:stagger = "" ;
> XLAT:coordinates = "XLONG XLAT" ;
> XLAT:long_name = "Latitude" ;
> XLAT:standard_name = "latitude" ;
> float XLONG(Time, south_north, west_east) ;
> XLONG:FieldType = 104 ;
> XLONG:MemoryOrder = "XY " ;
> XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;
> XLONG:units = "degrees_east" ;
> XLONG:stagger = "" ;
> XLONG:coordinates = "XLONG XLAT" ;
> XLONG:long_name = "Longitude" ;
> XLONG:standard_name = "longitude" ;
> float Q2(Time, south_north, west_east) ;
> Q2:FieldType = 104 ;
> Q2:MemoryOrder = "XY " ;
> Q2:description = "QV at 2 M" ;
> Q2:units = "kg kg-1" ;
> Q2:stagger = "" ;
> Q2:coordinates = "lon lat" ;
> Q2:long_name = "Specific Humidity at 2 m" ;
> Q2:standard_name = "specific_humidity" ;
> Q2:level = "2 m" ;
> Q2:_FillValue = 1.e+20f ;
> float T2(Time, south_north, west_east) ;
> T2:FieldType = 104 ;
> T2:MemoryOrder = "XY " ;
> T2:description = "TEMP at 2 M" ;
> T2:units = "K" ;
> T2:stagger = "" ;
> T2:coordinates = "lon lat" ;
> T2:long_name = "Air Temperature at 2 m" ;
> T2:standard_name = "air_temperature" ;
> T2:level = "2 m" ;
> T2:_FillValue = 1.e+20f ;
> float PSFC(Time, south_north, west_east) ;
> PSFC:FieldType = 104 ;
> PSFC:MemoryOrder = "XY " ;
> PSFC:description = "SFC PRESSURE" ;
> PSFC:units = "Pa" ;
> PSFC:stagger = "" ;
> PSFC:coordinates = "lon lat" ;
> PSFC:long_name = "Air Pressure at the Surface" ;
> PSFC:standard_name = "surface_pressure" ;
> PSFC:level = "surface" ;
> PSFC:_FillValue = 1.e+20f ;
> float U10(Time, south_north, west_east) ;
> U10:FieldType = 104 ;
> U10:MemoryOrder = "XY " ;
> U10:description = "U at 10 M" ;
> U10:units = "m/s" ;
> U10:stagger = "" ;
> U10:coordinates = "lon lat" ;
> U10:long_name = "u-Component of wind at 10 m
(Earth)" ;
> U10:standard_name = "eastward_wind" ;
> U10:level = "10 m" ;
> U10:_FillValue = 1.e+20f ;
> float V10(Time, south_north, west_east) ;
> V10:FieldType = 104 ;
> V10:MemoryOrder = "XY " ;
> V10:description = "V at 10 M" ;
> V10:units = "m/s" ;
> V10:stagger = "" ;
> V10:coordinates = "lon lat" ;
> V10:long_name = "v-Component of wind at 10 m
(Earth)" ;
> V10:standard_name = "northward_wind" ;
> V10:level = "10 m" ;
> V10:_FillValue = 1.e+20f ;
> double time(Time) ;
> time:standard_name = "time" ;
> time:long_name = "Time" ;
> time:calendar = "standard" ;
> time:units = "hours since 1970-01-01 00:00:00" ;
> float lat(south_north, west_east) ;
> lat:long_name = "Latitude" ;
> lat:standard_name = "latitude" ;
> lat:units = "degrees_north" ;
> float lon(south_north, west_east) ;
> lon:long_name = "Longitude" ;
> lon:standard_name = "longitude" ;
> lon:units = "degrees_east" ;
>
> // global attributes:
> :TITLE = " OUTPUT FROM WRF V4.1.1 MODEL" ;
> :START_DATE = "2021-03-01_00:00:00" ;
> :WEST-EAST_GRID_DIMENSION = 397 ;
> :SOUTH-NORTH_GRID_DIMENSION = 481 ;
> :BOTTOM-TOP_GRID_DIMENSION = 51 ;
> :DX = 2000.f ;
> :DY = 2000.f ;
> :AERCU_OPT = 0 ;
> :AERCU_FCT = 1.f ;
> :IDEAL_CASE = 0 ;
> :DIFF_6TH_SLOPEOPT = 0 ;
> :AUTO_LEVELS_OPT = 2 ;
> :DIFF_6TH_THRESH = 0.1f ;
> :DZBOT = 50.f ;
> :DZSTRETCH_S = 1.3f ;
> :DZSTRETCH_U = 1.1f ;
> :GRIDTYPE = "C" ;
> :DIFF_OPT = 2 ;
> :KM_OPT = 4 ;
> :DAMP_OPT = 3 ;
> :DAMPCOEF = 0.2f ;
> :KHDIF = 0.f ;
> :KVDIF = 0.f ;
> :MP_PHYSICS = 8 ;
> :RA_LW_PHYSICS = 4 ;
> :RA_SW_PHYSICS = 4 ;
> :SF_SFCLAY_PHYSICS = 5 ;
> :SF_SURFACE_PHYSICS = 4 ;
> :BL_PBL_PHYSICS = 6 ;
> :CU_PHYSICS = 0 ;
> :SF_LAKE_PHYSICS = 0 ;
> :SURFACE_INPUT_SOURCE = 3 ;
> :SST_UPDATE = 0 ;
> :GRID_FDDA = 0 ;
> :GFDDA_INTERVAL_M = 0 ;
> :GFDDA_END_H = 0 ;
> :GRID_SFDDA = 0 ;
> :SGFDDA_INTERVAL_M = 0 ;
> :SGFDDA_END_H = 0 ;
> :HYPSOMETRIC_OPT = 2 ;
> :USE_THETA_M = 1 ;
> :GWD_OPT = 0 ;
> :SF_URBAN_PHYSICS = 0 ;
> :SF_SURFACE_MOSAIC = 0 ;
> :SF_OCEAN_PHYSICS = 0 ;
> :WEST-EAST_PATCH_START_UNSTAG = 1 ;
> :WEST-EAST_PATCH_END_UNSTAG = 396 ;
> :WEST-EAST_PATCH_START_STAG = 1 ;
> :WEST-EAST_PATCH_END_STAG = 397 ;
> :SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
> :SOUTH-NORTH_PATCH_END_UNSTAG = 480 ;
> :SOUTH-NORTH_PATCH_START_STAG = 1 ;
> :SOUTH-NORTH_PATCH_END_STAG = 481 ;
> :BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
> :BOTTOM-TOP_PATCH_END_UNSTAG = 50 ;
> :BOTTOM-TOP_PATCH_START_STAG = 1 ;
> :BOTTOM-TOP_PATCH_END_STAG = 51 ;
> :GRID_ID = 3 ;
> :PARENT_ID = 2 ;
> :I_PARENT_START = 93 ;
> :J_PARENT_START = 75 ;
> :PARENT_GRID_RATIO = 3 ;
> :DT = 8.f ;
> :CEN_LAT = 38.00961f ;
> :CEN_LON = -121.0455f ;
> :TRUELAT1 = 31.f ;
> :TRUELAT2 = 45.38f ;
> :MOAD_CEN_LAT = 38.19999f ;
> :STAND_LON = -121.08f ;
> :POLE_LAT = 90.f ;
> :POLE_LON = 0.f ;
> :GMT = 0.f ;
> :JULYR = 2021 ;
> :JULDAY = 60 ;
> :MAP_PROJ = 1 ;
> :MAP_PROJ_CHAR = "Lambert Conformal" ;
> :MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
> :NUM_LAND_CAT = 21 ;
> :ISWATER = 17 ;
> :ISLAKE = 21 ;
> :ISICE = 15 ;
> :ISURBAN = 13 ;
> :ISOILWATER = 14 ;
> :HYBRID_OPT = 2 ;
> :ETAC = 0.2f ;
> :file_type = "NETCDF_NCCF" ;
> :title = "wrfsfc_d03_20210301_0000_f00000-cfcomp.nc"
;
> :Conventions = "CF 1.6, Standard Name Table v19" ;
> :source = "wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc" ;
> :notes = "" ;
> :created_by = "AJEB" ;
> :institution = "PG&E" ;
> :creation_date = "Mon Jun 7 17:59:23 UTC 2021" ;
> :history = "Mon Jun 7 17:59:26 2021: ncap2 -O -s
> T2F=float(T2F) wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
> wrfsfc_d03_20210301_0000_f00000-cfcomp.nc\n",
> "Mon Jun 7 17:59:26 2021: ncap2 -O -s
> T2F=((T2-273.15)*(9.0/5.0))+32.0 wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc
> wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
> :NCO = "netCDF Operators version 4.9.8 (Homepage =
>
https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fnco.sf.net%2F&data=04%7C01%7CAJEB%40pge.com%7Cae1ed37b041e47bfa3a708d929fdd547%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C637586992454384982%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=7Ltw5nmecCqII5PAhzrftQ51SLq6r0wmSDIud9Yb%2Fas%3D&reserved=0,
Code =
https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fgithub.com%2Fnco%2Fnco&data=04%7C01%7CAJEB%40pge.com%7Cae1ed37b041e47bfa3a708d929fdd547%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C637586992454384982%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=ej%2BZ8pSY5CBB3UjOMNyzXgxxhACXPCuizxKShx%2BVDp4%3D&reserved=0)"
;
>
>
>
------------------------------------------------
Subject: unexpected results from point_stat with NEAREST GRIDPOINT method
From: John Halley Gotway
Time: Tue Jun 08 17:29:45 2021
AJ,
Please try this command:
plot_data_plane wrfsfc_d03_20210301_0000_f00000-cfcomp.nc out.ps
'name="WSPD_10M";
level="(0,*,*)";'
That last option is like a mini-inline configuration file that tells
the
tool what data to read from that input file.
John
On Tue, Jun 8, 2021 at 5:04 PM Eiserloh Jr., A.J. via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100132 >
>
> Hi John,
>
> I am unable to ftp to send data to the ftp address unless I get a
firewall
> exception request from my company's IT, and that will take some
time. I
> tried plotting with the plot_data_plane command but get another
error that
> is not helpful.
>
> Thanks,
> -AJ
>
> [root at 698df1a69b34 nccompli_test]# plot_data_plane
> wrfsfc_d03_20210301_0000_f00000-cfcomp.nc out.ps WSPD_10M
> ERROR :
> ERROR : yyerror() -> syntax error in file "/tmp/met_config_273_0"
> ERROR :
> ERROR : line = 2
> ERROR :
> ERROR : column = 8
> ERROR :
> ERROR : text = "(nul)"
> ERROR :
> ERROR :
> ERROR : (nul)
> ERROR : (nul)
> ERROR :
>
>
>
> -----Original Message-----
> From: John Halley Gotway via RT <met_help at ucar.edu>
> Sent: Monday, June 7, 2021 2:47 PM
> To: Eiserloh Jr., A.J. <AJEB at pge.com>
> Subject: Re: [rt.rap.ucar.edu #100132] unexpected results from
point_stat
> with NEAREST GRIDPOINT method
>
> *****CAUTION: This email was sent from an EXTERNAL source. Think
before
> clicking links or opening attachments.*****
>
> AJ,
>
> I see that you're not getting the expected results from the nearest
> neighbor interpolation method in point-stat. There's nothing I can
see in
> the Point-Stat config file you sent that is an obvious problem.
>
> If I had the actual model/obs data available to me, I would start by
> running the plot_data_plane tool to plot the 10-m wind data from
your
> NetCDF file.
>
> *plot_data_plane wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
> <
>
https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwrfsfc_d03_20210301_0000_f00000-
cfcomp.nc%2F&data=04%7C01%7CAJEB%40pge.com%7Cae1ed37b041e47bfa3a708d929fdd547%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C637586992454384982%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=U3uNNDnM%2FqUKIC6UpHBuSjEeM%2BbtEQaP5v%2Bvuly%2FxiU%3D&reserved=0>
> WSPD_10M.ps 'name="WSPD_10M"; level="(0,*,*)";'*
>
> And I'd check to make sure that the data isn't upside-down or
backwards.
> If that looks good, I might also run the plot_point_obs tool. You
can
> actually plot WSPD_10M field from above as a base image and then
overlay
> your point.
> Check to make sure that point is really in the right spot on the
grid...
> confirm that sign of the longitude is correct... etc.
>
> If that all checks out, next I'd run point_stat through a debugger
to see
> line-by-line what's going on.
>
> But I'd need sample data to debug in this way. If you'd like to
provide
> that, you could post it to our anonymous FTP site following these
> instructions:
>
>
https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdtcenter.org%2Fcommunity-
code%2Fmodel-evaluation-tools-met%2Fmet-help-
desk%23ftp&data=04%7C01%7CAJEB%40pge.com%7Cae1ed37b041e47bfa3a708d929fdd547%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C637586992454384982%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=CqU3jPPMQuJrjK2vQoPEHnTUo9Erb3hQvddTVrB3LC8%3D&reserved=0
>
> Thanks,
> John Halley Gotway
>
> On Mon, Jun 7, 2021 at 12:45 PM Eiserloh Jr., A.J. via RT <
> met_help at ucar.edu>
> wrote:
>
> >
> > Mon Jun 07 12:45:45 2021: Request 100132 was acted upon.
> > Transaction: Ticket created by AJEB at pge.com
> > Queue: met_help
> > Subject: unexpected results from point_stat with NEAREST
> > GRIDPOINT method
> > Owner: Nobody
> > Requestors: AJEB at pge.com
> > Status: new
> > Ticket <URL:
> >
https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Frt.r
> >
ap.ucar.edu%2Frt%2FTicket%2FDisplay.html%3Fid%3D100132&data=04%7C0
> >
1%7CAJEB%40pge.com%7Cae1ed37b041e47bfa3a708d929fdd547%7C44ae661aece641
> >
aabc967c2c85a08941%7C0%7C0%7C637586992454384982%7CUnknown%7CTWFpbGZsb3
> >
d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7
> >
C1000&sdata=LuRAF64SY4mRnzYWzWw5YN8%2FKHVDJTkwzWj%2BolIhDmA%3D&
> > ;reserved=0 >
> >
> >
> > Dear Met-Help,
> >
> > I am testing out the point_stat tool (V10.0.0) for 1 station pair,
and
> > trying to see if it is grabbing the correct value of windspeed
from
> > the nearest gridpoint in my input forecast netcdf grid. In my
ascii
> > obs text I purposefully put the station lat/lon as one of the
> > gridpoints in my forecast field for lat lon indexes (starting with
0
> > in the ncl fashion) of 218,140. The lat/lon that this corresponds
to
> > is (37.61275, -122.3618). The 10-m windspeed at this point is
3.157048
> > m/s, however the point_stat*_mpr.txt results of the matched pair
has
> > the windspeed as
> > 3.57623 m/s, which is not correct. In fact, none of the
surrounding
> > points have that value either. Why is point stat giving me a
different
> > value? I don't think that it would be doing any kind of vertical
> > interpolation, but not sure. My expectation is that it should not
be.
> > I attached my PointStatConfig file, and the other info I provided
below
> that may help.
> >
> > Thanks,
> >
> > A.J. Eiserloh - Data Scientist
> > PG&E - Meteorology Systems and Analytics Applied Technology
Services
> > 3400 Crow Canyon Rd., San Ramon, CA 94583-1393
> > 925-307-4492
> >
> >
> >
> >
> > * Here is my pre ascii2nc ascii input columns for the point:
> >
> > ADPSFC KSFO 20210301_000000 37.61275 -122.3618 0.0 32 0 10.0 NA
> > 11.318018760716
> >
> >
> > * The *mpr.txt results for the one matched pair:
> >
> > OBS_LAT OBS_LON OBS_LVL OBS_ELV FCST
> > OBS OBS_QC CLIMO_MEAN CLIMO_STDEV CLIMO_CDF
> > 37.61275 -122.3618 NA 10
> > 3.57623 11.31802 NA NA
> > NA NA
> >
> >
> > * ncdump of my netcdf file:
> >
> > netcdf wrfsfc_d03_20210301_0000_f00000-cfcomp {
> > dimensions:
> > Time = UNLIMITED ; // (1 currently)
> > south_north = 480 ;
> > west_east = 396 ;
> > DateStrLen = 19 ;
> > soil_layers_stag = 4 ;
> > variables:
> > float WSPD_10M(Time, south_north, west_east) ;
> > WSPD_10M:_FillValue = 1.e+20f ;
> > WSPD_10M:FieldType = 104 ;
> > WSPD_10M:MemoryOrder = "XY " ;
> > WSPD_10M:coordinates = "lon lat" ;
> > WSPD_10M:description = "U at 10 M" ;
> > WSPD_10M:level = "10 m" ;
> > WSPD_10M:long_name = "Wind Speed at 10 m" ;
> > WSPD_10M:stagger = "" ;
> > WSPD_10M:standard_name = "wind_speed" ;
> > WSPD_10M:units = "m/s" ;
> > char Times(Time, DateStrLen) ;
> > float XLAT(Time, south_north, west_east) ;
> > XLAT:FieldType = 104 ;
> > XLAT:MemoryOrder = "XY " ;
> > XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
> > XLAT:units = "degrees_north" ;
> > XLAT:stagger = "" ;
> > XLAT:coordinates = "XLONG XLAT" ;
> > XLAT:long_name = "Latitude" ;
> > XLAT:standard_name = "latitude" ;
> > float XLONG(Time, south_north, west_east) ;
> > XLONG:FieldType = 104 ;
> > XLONG:MemoryOrder = "XY " ;
> > XLONG:description = "LONGITUDE, WEST IS NEGATIVE"
;
> > XLONG:units = "degrees_east" ;
> > XLONG:stagger = "" ;
> > XLONG:coordinates = "XLONG XLAT" ;
> > XLONG:long_name = "Longitude" ;
> > XLONG:standard_name = "longitude" ;
> > float Q2(Time, south_north, west_east) ;
> > Q2:FieldType = 104 ;
> > Q2:MemoryOrder = "XY " ;
> > Q2:description = "QV at 2 M" ;
> > Q2:units = "kg kg-1" ;
> > Q2:stagger = "" ;
> > Q2:coordinates = "lon lat" ;
> > Q2:long_name = "Specific Humidity at 2 m" ;
> > Q2:standard_name = "specific_humidity" ;
> > Q2:level = "2 m" ;
> > Q2:_FillValue = 1.e+20f ;
> > float T2(Time, south_north, west_east) ;
> > T2:FieldType = 104 ;
> > T2:MemoryOrder = "XY " ;
> > T2:description = "TEMP at 2 M" ;
> > T2:units = "K" ;
> > T2:stagger = "" ;
> > T2:coordinates = "lon lat" ;
> > T2:long_name = "Air Temperature at 2 m" ;
> > T2:standard_name = "air_temperature" ;
> > T2:level = "2 m" ;
> > T2:_FillValue = 1.e+20f ;
> > float PSFC(Time, south_north, west_east) ;
> > PSFC:FieldType = 104 ;
> > PSFC:MemoryOrder = "XY " ;
> > PSFC:description = "SFC PRESSURE" ;
> > PSFC:units = "Pa" ;
> > PSFC:stagger = "" ;
> > PSFC:coordinates = "lon lat" ;
> > PSFC:long_name = "Air Pressure at the Surface" ;
> > PSFC:standard_name = "surface_pressure" ;
> > PSFC:level = "surface" ;
> > PSFC:_FillValue = 1.e+20f ;
> > float U10(Time, south_north, west_east) ;
> > U10:FieldType = 104 ;
> > U10:MemoryOrder = "XY " ;
> > U10:description = "U at 10 M" ;
> > U10:units = "m/s" ;
> > U10:stagger = "" ;
> > U10:coordinates = "lon lat" ;
> > U10:long_name = "u-Component of wind at 10 m
(Earth)" ;
> > U10:standard_name = "eastward_wind" ;
> > U10:level = "10 m" ;
> > U10:_FillValue = 1.e+20f ;
> > float V10(Time, south_north, west_east) ;
> > V10:FieldType = 104 ;
> > V10:MemoryOrder = "XY " ;
> > V10:description = "V at 10 M" ;
> > V10:units = "m/s" ;
> > V10:stagger = "" ;
> > V10:coordinates = "lon lat" ;
> > V10:long_name = "v-Component of wind at 10 m
(Earth)" ;
> > V10:standard_name = "northward_wind" ;
> > V10:level = "10 m" ;
> > V10:_FillValue = 1.e+20f ;
> > double time(Time) ;
> > time:standard_name = "time" ;
> > time:long_name = "Time" ;
> > time:calendar = "standard" ;
> > time:units = "hours since 1970-01-01 00:00:00" ;
> > float lat(south_north, west_east) ;
> > lat:long_name = "Latitude" ;
> > lat:standard_name = "latitude" ;
> > lat:units = "degrees_north" ;
> > float lon(south_north, west_east) ;
> > lon:long_name = "Longitude" ;
> > lon:standard_name = "longitude" ;
> > lon:units = "degrees_east" ;
> >
> > // global attributes:
> > :TITLE = " OUTPUT FROM WRF V4.1.1 MODEL" ;
> > :START_DATE = "2021-03-01_00:00:00" ;
> > :WEST-EAST_GRID_DIMENSION = 397 ;
> > :SOUTH-NORTH_GRID_DIMENSION = 481 ;
> > :BOTTOM-TOP_GRID_DIMENSION = 51 ;
> > :DX = 2000.f ;
> > :DY = 2000.f ;
> > :AERCU_OPT = 0 ;
> > :AERCU_FCT = 1.f ;
> > :IDEAL_CASE = 0 ;
> > :DIFF_6TH_SLOPEOPT = 0 ;
> > :AUTO_LEVELS_OPT = 2 ;
> > :DIFF_6TH_THRESH = 0.1f ;
> > :DZBOT = 50.f ;
> > :DZSTRETCH_S = 1.3f ;
> > :DZSTRETCH_U = 1.1f ;
> > :GRIDTYPE = "C" ;
> > :DIFF_OPT = 2 ;
> > :KM_OPT = 4 ;
> > :DAMP_OPT = 3 ;
> > :DAMPCOEF = 0.2f ;
> > :KHDIF = 0.f ;
> > :KVDIF = 0.f ;
> > :MP_PHYSICS = 8 ;
> > :RA_LW_PHYSICS = 4 ;
> > :RA_SW_PHYSICS = 4 ;
> > :SF_SFCLAY_PHYSICS = 5 ;
> > :SF_SURFACE_PHYSICS = 4 ;
> > :BL_PBL_PHYSICS = 6 ;
> > :CU_PHYSICS = 0 ;
> > :SF_LAKE_PHYSICS = 0 ;
> > :SURFACE_INPUT_SOURCE = 3 ;
> > :SST_UPDATE = 0 ;
> > :GRID_FDDA = 0 ;
> > :GFDDA_INTERVAL_M = 0 ;
> > :GFDDA_END_H = 0 ;
> > :GRID_SFDDA = 0 ;
> > :SGFDDA_INTERVAL_M = 0 ;
> > :SGFDDA_END_H = 0 ;
> > :HYPSOMETRIC_OPT = 2 ;
> > :USE_THETA_M = 1 ;
> > :GWD_OPT = 0 ;
> > :SF_URBAN_PHYSICS = 0 ;
> > :SF_SURFACE_MOSAIC = 0 ;
> > :SF_OCEAN_PHYSICS = 0 ;
> > :WEST-EAST_PATCH_START_UNSTAG = 1 ;
> > :WEST-EAST_PATCH_END_UNSTAG = 396 ;
> > :WEST-EAST_PATCH_START_STAG = 1 ;
> > :WEST-EAST_PATCH_END_STAG = 397 ;
> > :SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
> > :SOUTH-NORTH_PATCH_END_UNSTAG = 480 ;
> > :SOUTH-NORTH_PATCH_START_STAG = 1 ;
> > :SOUTH-NORTH_PATCH_END_STAG = 481 ;
> > :BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
> > :BOTTOM-TOP_PATCH_END_UNSTAG = 50 ;
> > :BOTTOM-TOP_PATCH_START_STAG = 1 ;
> > :BOTTOM-TOP_PATCH_END_STAG = 51 ;
> > :GRID_ID = 3 ;
> > :PARENT_ID = 2 ;
> > :I_PARENT_START = 93 ;
> > :J_PARENT_START = 75 ;
> > :PARENT_GRID_RATIO = 3 ;
> > :DT = 8.f ;
> > :CEN_LAT = 38.00961f ;
> > :CEN_LON = -121.0455f ;
> > :TRUELAT1 = 31.f ;
> > :TRUELAT2 = 45.38f ;
> > :MOAD_CEN_LAT = 38.19999f ;
> > :STAND_LON = -121.08f ;
> > :POLE_LAT = 90.f ;
> > :POLE_LON = 0.f ;
> > :GMT = 0.f ;
> > :JULYR = 2021 ;
> > :JULDAY = 60 ;
> > :MAP_PROJ = 1 ;
> > :MAP_PROJ_CHAR = "Lambert Conformal" ;
> > :MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
> > :NUM_LAND_CAT = 21 ;
> > :ISWATER = 17 ;
> > :ISLAKE = 21 ;
> > :ISICE = 15 ;
> > :ISURBAN = 13 ;
> > :ISOILWATER = 14 ;
> > :HYBRID_OPT = 2 ;
> > :ETAC = 0.2f ;
> > :file_type = "NETCDF_NCCF" ;
> > :title = "wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc" ;
> > :Conventions = "CF 1.6, Standard Name Table v19" ;
> > :source = "wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc" ;
> > :notes = "" ;
> > :created_by = "AJEB" ;
> > :institution = "PG&E" ;
> > :creation_date = "Mon Jun 7 17:59:23 UTC 2021" ;
> > :history = "Mon Jun 7 17:59:26 2021: ncap2 -O -s
> > T2F=float(T2F) wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
> > wrfsfc_d03_20210301_0000_f00000-cfcomp.nc\n",
> > "Mon Jun 7 17:59:26 2021: ncap2 -O -s
> > T2F=((T2-273.15)*(9.0/5.0))+32.0
> wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
> > wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
> > :NCO = "netCDF Operators version 4.9.8 (Homepage =
> >
>
https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fnco.sf.net%2F&data=04%7C01%7CAJEB%40pge.com%7Cae1ed37b041e47bfa3a708d929fdd547%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C637586992454384982%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=7Ltw5nmecCqII5PAhzrftQ51SLq6r0wmSDIud9Yb%2Fas%3D&reserved=0,
> Code =
>
https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fgithub.com%2Fnco%2Fnco&data=04%7C01%7CAJEB%40pge.com%7Cae1ed37b041e47bfa3a708d929fdd547%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C637586992454384982%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=ej%2BZ8pSY5CBB3UjOMNyzXgxxhACXPCuizxKShx%2BVDp4%3D&reserved=0)"
> ;
> >
> >
> >
>
>
>
>
------------------------------------------------
Subject: unexpected results from point_stat with NEAREST GRIDPOINT method
From: Eiserloh Jr., A.J.
Time: Tue Jun 08 18:26:31 2021
Hi John,
I got that to work this time, and the image (attached) looks correct
and the map looks good with the data, but not zoomed into my point is
hard for me to tell:
Seems to me like it is interpolating still in some direction, or it is
changing my map projection maybe because I didn't provide enough data?
Here is the full ncdump of my netdf file below if that helps. I did
notice that with a grib2 file of similar model type, it gives me the
expected value at the point after running point_stat. And that grib2
file has all the projection info associated with the lat and lon
variables. But with my netcdf file that I created manually to makeCF
compliant (because that is what is needed for met tool), it is not
giving me the expected value at the nearest location gridpoint as I
described, and my netcdf does not have the projection info as
attributes tied to the lat and lon variables. Rather it is in the
global attributes and I am wondering if Met Tool can't read that info
about the projection?
Thanks,
AJ
netcdf wrfsfc_d03_20210301_0000_f00000-cfcomp {
dimensions:
Time = UNLIMITED ; // (1 currently)
south_north = 480 ;
west_east = 396 ;
DateStrLen = 19 ;
soil_layers_stag = 4 ;
variables:
float T2F(Time, south_north, west_east) ;
T2F:_FillValue = 1.e+20f ;
T2F:FieldType = 104 ;
T2F:MemoryOrder = "XY " ;
T2F:coordinates = "lon lat" ;
T2F:description = "TEMP at 2 M" ;
T2F:level = "2 m" ;
T2F:long_name = "Air Temperature at 2 m" ;
T2F:stagger = "" ;
T2F:standard_name = "air_temperature" ;
T2F:units = "K" ;
float WSPD_10M(Time, south_north, west_east) ;
WSPD_10M:_FillValue = 1.e+20f ;
WSPD_10M:FieldType = 104 ;
WSPD_10M:MemoryOrder = "XY " ;
WSPD_10M:coordinates = "lon lat" ;
WSPD_10M:description = "U at 10 M" ;
WSPD_10M:level = "10 m" ;
WSPD_10M:long_name = "Wind Speed at 10 m" ;
WSPD_10M:stagger = "" ;
WSPD_10M:standard_name = "wind_speed" ;
WSPD_10M:units = "m/s" ;
char Times(Time, DateStrLen) ;
float XLAT(Time, south_north, west_east) ;
XLAT:FieldType = 104 ;
XLAT:MemoryOrder = "XY " ;
XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
XLAT:units = "degrees_north" ;
XLAT:stagger = "" ;
XLAT:coordinates = "XLONG XLAT" ;
XLAT:long_name = "Latitude" ;
XLAT:standard_name = "latitude" ;
float XLONG(Time, south_north, west_east) ;
XLONG:FieldType = 104 ;
XLONG:MemoryOrder = "XY " ;
XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;
XLONG:units = "degrees_east" ;
XLONG:stagger = "" ;
XLONG:coordinates = "XLONG XLAT" ;
XLONG:long_name = "Longitude" ;
XLONG:standard_name = "longitude" ;
float Q2(Time, south_north, west_east) ;
Q2:FieldType = 104 ;
Q2:MemoryOrder = "XY " ;
Q2:description = "QV at 2 M" ;
Q2:units = "kg kg-1" ;
Q2:stagger = "" ;
Q2:coordinates = "lon lat" ;
Q2:long_name = "Specific Humidity at 2 m" ;
Q2:standard_name = "specific_humidity" ;
Q2:level = "2 m" ;
Q2:_FillValue = 1.e+20f ;
float T2(Time, south_north, west_east) ;
T2:FieldType = 104 ;
T2:MemoryOrder = "XY " ;
T2:description = "TEMP at 2 M" ;
T2:units = "K" ;
T2:stagger = "" ;
T2:coordinates = "lon lat" ;
T2:long_name = "Air Temperature at 2 m" ;
T2:standard_name = "air_temperature" ;
T2:level = "2 m" ;
T2:_FillValue = 1.e+20f ;
float PSFC(Time, south_north, west_east) ;
PSFC:FieldType = 104 ;
PSFC:MemoryOrder = "XY " ;
PSFC:description = "SFC PRESSURE" ;
PSFC:units = "Pa" ;
PSFC:stagger = "" ;
PSFC:coordinates = "lon lat" ;
PSFC:long_name = "Air Pressure at the Surface" ;
PSFC:standard_name = "surface_pressure" ;
PSFC:level = "surface" ;
PSFC:_FillValue = 1.e+20f ;
float U10(Time, south_north, west_east) ;
U10:FieldType = 104 ;
U10:MemoryOrder = "XY " ;
U10:description = "U at 10 M" ;
U10:units = "m/s" ;
U10:stagger = "" ;
U10:coordinates = "lon lat" ;
U10:long_name = "u-Component of wind at 10 m (Earth)"
;
U10:standard_name = "eastward_wind" ;
U10:level = "10 m" ;
U10:_FillValue = 1.e+20f ;
float V10(Time, south_north, west_east) ;
V10:FieldType = 104 ;
V10:MemoryOrder = "XY " ;
V10:description = "V at 10 M" ;
V10:units = "m/s" ;
V10:stagger = "" ;
V10:coordinates = "lon lat" ;
V10:long_name = "v-Component of wind at 10 m (Earth)"
;
V10:standard_name = "northward_wind" ;
V10:level = "10 m" ;
V10:_FillValue = 1.e+20f ;
float XTIME(Time) ;
XTIME:FieldType = 104 ;
XTIME:MemoryOrder = "0 " ;
XTIME:description = "minutes since 2021-02-28
21:00:00" ;
XTIME:units = "minutes since 2021-02-28 21:00:00" ;
XTIME:stagger = "" ;
float TSLB(Time, soil_layers_stag, south_north, west_east) ;
TSLB:FieldType = 104 ;
TSLB:MemoryOrder = "XYZ" ;
TSLB:description = "SOIL TEMPERATURE" ;
TSLB:units = "K" ;
TSLB:stagger = "Z" ;
TSLB:coordinates = "lon lat" ;
TSLB:long_name = "Soil Temperature" ;
TSLB:standard_name = "soil_temperature" ;
TSLB:_FillValue = 1.e+20f ;
float SMOIS(Time, soil_layers_stag, south_north, west_east) ;
SMOIS:FieldType = 104 ;
SMOIS:MemoryOrder = "XYZ" ;
SMOIS:description = "SOIL MOISTURE" ;
SMOIS:units = "m3 m-3" ;
SMOIS:stagger = "Z" ;
SMOIS:coordinates = "lon lat" ;
SMOIS:long_name = "Soil Moisture" ;
SMOIS:standard_name = "soil_moisture" ;
SMOIS:_FillValue = 1.e+20f ;
float ACSNOM(Time, south_north, west_east) ;
ACSNOM:FieldType = 104 ;
ACSNOM:MemoryOrder = "XY " ;
ACSNOM:description = "ACCUMULATED MELTED SNOW" ;
ACSNOM:units = "kg m-2" ;
ACSNOM:stagger = "" ;
ACSNOM:coordinates = "lon lat" ;
ACSNOM:long_name = "DOWNWARD SHORT WAVE FLUX AT GROUND
SURFACE" ;
ACSNOM:standard_name = "incoming shortwave raditation
at surface" ;
ACSNOM:_FillValue = 1.e+20f ;
float SNOWH(Time, south_north, west_east) ;
SNOWH:FieldType = 104 ;
SNOWH:MemoryOrder = "XY " ;
SNOWH:description = "PHYSICAL SNOW DEPTH" ;
SNOWH:units = "m" ;
SNOWH:stagger = "" ;
SNOWH:coordinates = "lon lat" ;
SNOWH:long_name = "Physical Snow Depth" ;
SNOWH:standard_name = "snow_depth" ;
SNOWH:_FillValue = 1.e+20f ;
float SWDOWN(Time, south_north, west_east) ;
SWDOWN:FieldType = 104 ;
SWDOWN:MemoryOrder = "XY " ;
SWDOWN:description = "DOWNWARD SHORT WAVE FLUX AT
GROUND SURFACE" ;
SWDOWN:units = "W m-2" ;
SWDOWN:stagger = "" ;
SWDOWN:coordinates = "lon lat" ;
SWDOWN:long_name = "DOWNWARD SHORT WAVE FLUX AT GROUND
SURFACE" ;
SWDOWN:standard_name = "incoming shortwave raditation
at surface" ;
SWDOWN:_FillValue = 1.e+20f ;
float ZNT(Time, south_north, west_east) ;
ZNT:FieldType = 104 ;
ZNT:MemoryOrder = "XY " ;
ZNT:description = "TIME-VARYING ROUGHNESS LENGTH" ;
ZNT:units = "m" ;
ZNT:stagger = "" ;
ZNT:coordinates = "lon lat" ;
ZNT:long_name = "Time Varying Roughness Length" ;
ZNT:standard_name = "roughness_length" ;
ZNT:_FillValue = 1.e+20f ;
float UST(Time, south_north, west_east) ;
UST:FieldType = 104 ;
UST:MemoryOrder = "XY " ;
UST:description = "U* IN SIMILARITY THEORY" ;
UST:units = "m s-1" ;
UST:stagger = "" ;
UST:coordinates = "lon lat" ;
UST:long_name = "U* in Similatiry Theory / Friction
Velocity" ;
UST:standard_name = "friction_velocity" ;
UST:_FillValue = 1.e+20f ;
float PREC_ACC_C(Time, south_north, west_east) ;
PREC_ACC_C:FieldType = 104 ;
PREC_ACC_C:MemoryOrder = "XY " ;
PREC_ACC_C:description = "ACCUMULATED CUMULUS
PRECIPITATION OVER prec_acc_dt PERIODS OF TIME" ;
PREC_ACC_C:units = "mm" ;
PREC_ACC_C:stagger = "" ;
PREC_ACC_C:coordinates = "lon lat" ;
PREC_ACC_C:long_name = "ACCUMULATED CUMULUS
PRECIPITATION" ;
PREC_ACC_C:standard_name = "" ;
PREC_ACC_C:_FillValue = 1.e+20f ;
float PREC_ACC_NC(Time, south_north, west_east) ;
PREC_ACC_NC:FieldType = 104 ;
PREC_ACC_NC:MemoryOrder = "XY " ;
PREC_ACC_NC:description = "ACCUMULATED GRID SCALE
PRECIPITATION OVER prec_acc_dt PERIODS OF TIME" ;
PREC_ACC_NC:units = "mm" ;
PREC_ACC_NC:stagger = "" ;
PREC_ACC_NC:coordinates = "lon lat" ;
PREC_ACC_NC:long_name = "ACCUMULATED GRID SCALE
PRECIPITATION" ;
PREC_ACC_NC:standard_name = "" ;
PREC_ACC_NC:level = "surface" ;
PREC_ACC_NC:_FillValue = 1.e+20f ;
float SNOW_ACC_NC(Time, south_north, west_east) ;
SNOW_ACC_NC:FieldType = 104 ;
SNOW_ACC_NC:MemoryOrder = "XY " ;
SNOW_ACC_NC:description = "ACCUMULATED SNOW WATER
EQUIVALENT OVER prec_acc_dt PERIODS OF TIME" ;
SNOW_ACC_NC:units = "mm" ;
SNOW_ACC_NC:stagger = "" ;
SNOW_ACC_NC:coordinates = "lon lat" ;
SNOW_ACC_NC:long_name = "ACCUMULATED SNOW WATER
EQUIVALENT" ;
SNOW_ACC_NC:standard_name = "" ;
SNOW_ACC_NC:level = "surface" ;
SNOW_ACC_NC:_FillValue = 1.e+20f ;
double time(Time) ;
time:standard_name = "time" ;
time:long_name = "Time" ;
time:calendar = "standard" ;
time:units = "hours since 1970-01-01 00:00:00" ;
float lat(south_north, west_east) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float lon(south_north, west_east) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
// global attributes:
:TITLE = " OUTPUT FROM WRF V4.1.1 MODEL" ;
:START_DATE = "2021-03-01_00:00:00" ;
:BOTTOM-TOP_GRID_DIMENSION = 51 ;
:DX = 2000.f ;
:DY = 2000.f ;
:AERCU_OPT = 0 ;
:AERCU_FCT = 1.f ;
:IDEAL_CASE = 0 ;
:DIFF_6TH_SLOPEOPT = 0 ;
:AUTO_LEVELS_OPT = 2 ;
:DIFF_6TH_THRESH = 0.1f ;
:DZBOT = 50.f ;
:DZSTRETCH_S = 1.3f ;
:DZSTRETCH_U = 1.1f ;
:GRIDTYPE = "C" ;
:DIFF_OPT = 2 ;
:KM_OPT = 4 ;
:DAMP_OPT = 3 ;
:DAMPCOEF = 0.2f ;
:KHDIF = 0.f ;
:KVDIF = 0.f ;
:MP_PHYSICS = 8 ;
:RA_LW_PHYSICS = 4 ;
:RA_SW_PHYSICS = 4 ;
:SF_SFCLAY_PHYSICS = 5 ;
:SF_SURFACE_PHYSICS = 4 ;
:BL_PBL_PHYSICS = 6 ;
:CU_PHYSICS = 0 ;
:SF_LAKE_PHYSICS = 0 ;
:SURFACE_INPUT_SOURCE = 3 ;
:SST_UPDATE = 0 ;
:GRID_FDDA = 0 ;
:GFDDA_INTERVAL_M = 0 ;
:GFDDA_END_H = 0 ;
:GRID_SFDDA = 0 ;
:SGFDDA_INTERVAL_M = 0 ;
:SGFDDA_END_H = 0 ;
:HYPSOMETRIC_OPT = 2 ;
:USE_THETA_M = 1 ;
:GWD_OPT = 0 ;
:SF_URBAN_PHYSICS = 0 ;
:SF_SURFACE_MOSAIC = 0 ;
:SF_OCEAN_PHYSICS = 0 ;
:WEST-EAST_PATCH_START_UNSTAG = 1 ;
:WEST-EAST_PATCH_END_UNSTAG = 396 ;
:WEST-EAST_PATCH_START_STAG = 1 ;
:WEST-EAST_PATCH_END_STAG = 397 ;
:SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
:SOUTH-NORTH_PATCH_END_UNSTAG = 480 ;
:SOUTH-NORTH_PATCH_START_STAG = 1 ;
:SOUTH-NORTH_PATCH_END_STAG = 481 ;
:BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
:BOTTOM-TOP_PATCH_END_UNSTAG = 50 ;
:BOTTOM-TOP_PATCH_START_STAG = 1 ;
:BOTTOM-TOP_PATCH_END_STAG = 51 ;
:GRID_ID = 3 ;
:PARENT_ID = 2 ;
:I_PARENT_START = 93 ;
:J_PARENT_START = 75 ;
:PARENT_GRID_RATIO = 3 ;
:DT = 8.f ;
:CEN_LAT = 38.00961f ;
:CEN_LON = -121.0455f ;
:TRUELAT1 = 31.f ;
:TRUELAT2 = 45.38f ;
:MOAD_CEN_LAT = 38.19999f ;
:STAND_LON = -121.08f ;
:POLE_LAT = 90.f ;
:POLE_LON = 0.f ;
:GMT = 0.f ;
:JULYR = 2021 ;
:JULDAY = 60 ;
:MAP_PROJ = 1 ;
:MAP_PROJ_CHAR = "Lambert Conformal" ;
:MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
:NUM_LAND_CAT = 21 ;
:ISWATER = 17 ;
:ISLAKE = 21 ;
:ISICE = 15 ;
:ISURBAN = 13 ;
:ISOILWATER = 14 ;
:HYBRID_OPT = 2 ;
:ETAC = 0.2f ;
:file_type = "NETCDF_NCCF" ;
:title = "wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
:Conventions = "CF 1.6, Standard Name Table v19" ;
:source = "wrfsfc_d03_20210301_0000_f00000-cfcomp.nc"
;
:notes = "" ;
:created_by = "AJEB" ;
:institution = "PG&E" ;
:creation_date = "Tue Jun 8 23:15:57 UTC 2021" ;
:history = "Tue Jun 8 23:16:00 2021: ncatted -O -a
SOUTH-NORTH_GRID_DIMENSION,global,d,, wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc\n",
"Tue Jun 8 23:16:00 2021: ncatted -O -a WEST-
EAST_GRID_DIMENSION,global,d,, wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc\n",
"Tue Jun 8 23:16:00 2021: ncap2 -O -s
T2F=float(T2F) wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
wrfsfc_d03_20210301_0000_f00000-cfcomp.nc\n",
"Tue Jun 8 23:16:00 2021: ncap2 -O -s
T2F=((T2-273.15)*(9.0/5.0))+32.0 wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
:NCO = "netCDF Operators version 4.9.8 (Homepage =
http://nco.sf.net, Code = http://github.com/nco/nco)" ;
-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Tuesday, June 8, 2021 4:30 PM
To: Eiserloh Jr., A.J. <AJEB at pge.com>
Subject: Re: [rt.rap.ucar.edu #100132] unexpected results from
point_stat with NEAREST GRIDPOINT method
*****CAUTION: This email was sent from an EXTERNAL source. Think
before clicking links or opening attachments.*****
AJ,
Please try this command:
plot_data_plane wrfsfc_d03_20210301_0000_f00000-cfcomp.nc out.ps
'name="WSPD_10M"; level="(0,*,*)";'
That last option is like a mini-inline configuration file that tells
the tool what data to read from that input file.
John
On Tue, Jun 8, 2021 at 5:04 PM Eiserloh Jr., A.J. via RT
<met_help at ucar.edu>
wrote:
>
> <URL:
>
https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Frt.r
>
ap.ucar.edu%2Frt%2FTicket%2FDisplay.html%3Fid%3D100132&data=04%7C0
>
1%7CAJEB%40pge.com%7C0b860e74eccf45e2508508d92ad54e1c%7C44ae661aece641
>
aabc967c2c85a08941%7C0%7C0%7C637587917903778911%7CUnknown%7CTWFpbGZsb3
>
d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7
>
C1000&sdata=eDLjuvbSEv8cosgHtciLwoyf9TaR%2BSTWkLRNejmd9MQ%3D&r
> eserved=0 >
>
> Hi John,
>
> I am unable to ftp to send data to the ftp address unless I get a
> firewall exception request from my company's IT, and that will take
> some time. I tried plotting with the plot_data_plane command but get
> another error that is not helpful.
>
> Thanks,
> -AJ
>
> [root at 698df1a69b34 nccompli_test]# plot_data_plane
> wrfsfc_d03_20210301_0000_f00000-cfcomp.nc out.ps WSPD_10M ERROR :
> ERROR : yyerror() -> syntax error in file "/tmp/met_config_273_0"
> ERROR :
> ERROR : line = 2
> ERROR :
> ERROR : column = 8
> ERROR :
> ERROR : text = "(nul)"
> ERROR :
> ERROR :
> ERROR : (nul)
> ERROR : (nul)
> ERROR :
>
>
>
> -----Original Message-----
> From: John Halley Gotway via RT <met_help at ucar.edu>
> Sent: Monday, June 7, 2021 2:47 PM
> To: Eiserloh Jr., A.J. <AJEB at pge.com>
> Subject: Re: [rt.rap.ucar.edu #100132] unexpected results from
> point_stat with NEAREST GRIDPOINT method
>
> *****CAUTION: This email was sent from an EXTERNAL source. Think
> before clicking links or opening attachments.*****
>
> AJ,
>
> I see that you're not getting the expected results from the nearest
> neighbor interpolation method in point-stat. There's nothing I can
see
> in the Point-Stat config file you sent that is an obvious problem.
>
> If I had the actual model/obs data available to me, I would start by
> running the plot_data_plane tool to plot the 10-m wind data from
your
> NetCDF file.
>
> *plot_data_plane wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
> <
>
https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwrfsf
> c_d03_20210301_0000_f00000-
cfcomp.nc%2F&data=04%7C01%7CAJEB%40pge.
>
com%7C0b860e74eccf45e2508508d92ad54e1c%7C44ae661aece641aabc967c2c85a08
>
941%7C0%7C0%7C637587917903778911%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wL
>
jAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata
> =8H2KGfAsMk%2BQocaxz75QIkeAgE%2Fctl9i0w99nQs0HuA%3D&reserved=0>
> WSPD_10M.ps 'name="WSPD_10M"; level="(0,*,*)";'*
>
> And I'd check to make sure that the data isn't upside-down or
backwards.
> If that looks good, I might also run the plot_point_obs tool. You
can
> actually plot WSPD_10M field from above as a base image and then
> overlay your point.
> Check to make sure that point is really in the right spot on the
grid...
> confirm that sign of the longitude is correct... etc.
>
> If that all checks out, next I'd run point_stat through a debugger
to
> see line-by-line what's going on.
>
> But I'd need sample data to debug in this way. If you'd like to
> provide that, you could post it to our anonymous FTP site following
> these
> instructions:
>
>
https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdtce
> nter.org%2Fcommunity-code%2Fmodel-evaluation-tools-met%2Fmet-help-
desk
>
%23ftp&data=04%7C01%7CAJEB%40pge.com%7C0b860e74eccf45e2508508d92ad
>
54e1c%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C637587917903778911%
>
7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik
>
1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=MqE2mOm9CNXhqn%2F1zAUcsLrLBRWkmc
> LLY%2FCp3Tcu34g%3D&reserved=0
>
> Thanks,
> John Halley Gotway
>
> On Mon, Jun 7, 2021 at 12:45 PM Eiserloh Jr., A.J. via RT <
> met_help at ucar.edu>
> wrote:
>
> >
> > Mon Jun 07 12:45:45 2021: Request 100132 was acted upon.
> > Transaction: Ticket created by AJEB at pge.com
> > Queue: met_help
> > Subject: unexpected results from point_stat with NEAREST
> > GRIDPOINT method
> > Owner: Nobody
> > Requestors: AJEB at pge.com
> > Status: new
> > Ticket <URL:
> >
https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Frt
> > .r
> >
ap.ucar.edu%2Frt%2FTicket%2FDisplay.html%3Fid%3D100132&data=04%7
> > C0
> >
1%7CAJEB%40pge.com%7Cae1ed37b041e47bfa3a708d929fdd547%7C44ae661aece6
> > 41
> >
aabc967c2c85a08941%7C0%7C0%7C637586992454384982%7CUnknown%7CTWFpbGZs
> > b3
> >
d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D
> > %7
> >
C1000&sdata=LuRAF64SY4mRnzYWzWw5YN8%2FKHVDJTkwzWj%2BolIhDmA%3D&a
> > mp
> > ;reserved=0 >
> >
> >
> > Dear Met-Help,
> >
> > I am testing out the point_stat tool (V10.0.0) for 1 station pair,
> > and trying to see if it is grabbing the correct value of windspeed
> > from the nearest gridpoint in my input forecast netcdf grid. In
my
> > ascii obs text I purposefully put the station lat/lon as one of
the
> > gridpoints in my forecast field for lat lon indexes (starting with
0
> > in the ncl fashion) of 218,140. The lat/lon that this corresponds
to
> > is (37.61275, -122.3618). The 10-m windspeed at this point is
> > 3.157048 m/s, however the point_stat*_mpr.txt results of the
matched
> > pair has the windspeed as
> > 3.57623 m/s, which is not correct. In fact, none of the
surrounding
> > points have that value either. Why is point stat giving me a
> > different value? I don't think that it would be doing any kind of
> > vertical interpolation, but not sure. My expectation is that it
should not be.
> > I attached my PointStatConfig file, and the other info I provided
> > below
> that may help.
> >
> > Thanks,
> >
> > A.J. Eiserloh - Data Scientist
> > PG&E - Meteorology Systems and Analytics Applied Technology
Services
> > 3400 Crow Canyon Rd., San Ramon, CA 94583-1393
> > 925-307-4492
> >
> >
> >
> >
> > * Here is my pre ascii2nc ascii input columns for the point:
> >
> > ADPSFC KSFO 20210301_000000 37.61275 -122.3618 0.0 32 0 10.0 NA
> > 11.318018760716
> >
> >
> > * The *mpr.txt results for the one matched pair:
> >
> > OBS_LAT OBS_LON OBS_LVL OBS_ELV FCST
> > OBS OBS_QC CLIMO_MEAN CLIMO_STDEV CLIMO_CDF
> > 37.61275 -122.3618 NA 10
> > 3.57623 11.31802 NA NA
> > NA NA
> >
> >
> > * ncdump of my netcdf file:
> >
> > netcdf wrfsfc_d03_20210301_0000_f00000-cfcomp {
> > dimensions:
> > Time = UNLIMITED ; // (1 currently)
> > south_north = 480 ;
> > west_east = 396 ;
> > DateStrLen = 19 ;
> > soil_layers_stag = 4 ;
> > variables:
> > float WSPD_10M(Time, south_north, west_east) ;
> > WSPD_10M:_FillValue = 1.e+20f ;
> > WSPD_10M:FieldType = 104 ;
> > WSPD_10M:MemoryOrder = "XY " ;
> > WSPD_10M:coordinates = "lon lat" ;
> > WSPD_10M:description = "U at 10 M" ;
> > WSPD_10M:level = "10 m" ;
> > WSPD_10M:long_name = "Wind Speed at 10 m" ;
> > WSPD_10M:stagger = "" ;
> > WSPD_10M:standard_name = "wind_speed" ;
> > WSPD_10M:units = "m/s" ;
> > char Times(Time, DateStrLen) ;
> > float XLAT(Time, south_north, west_east) ;
> > XLAT:FieldType = 104 ;
> > XLAT:MemoryOrder = "XY " ;
> > XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
> > XLAT:units = "degrees_north" ;
> > XLAT:stagger = "" ;
> > XLAT:coordinates = "XLONG XLAT" ;
> > XLAT:long_name = "Latitude" ;
> > XLAT:standard_name = "latitude" ;
> > float XLONG(Time, south_north, west_east) ;
> > XLONG:FieldType = 104 ;
> > XLONG:MemoryOrder = "XY " ;
> > XLONG:description = "LONGITUDE, WEST IS NEGATIVE"
;
> > XLONG:units = "degrees_east" ;
> > XLONG:stagger = "" ;
> > XLONG:coordinates = "XLONG XLAT" ;
> > XLONG:long_name = "Longitude" ;
> > XLONG:standard_name = "longitude" ;
> > float Q2(Time, south_north, west_east) ;
> > Q2:FieldType = 104 ;
> > Q2:MemoryOrder = "XY " ;
> > Q2:description = "QV at 2 M" ;
> > Q2:units = "kg kg-1" ;
> > Q2:stagger = "" ;
> > Q2:coordinates = "lon lat" ;
> > Q2:long_name = "Specific Humidity at 2 m" ;
> > Q2:standard_name = "specific_humidity" ;
> > Q2:level = "2 m" ;
> > Q2:_FillValue = 1.e+20f ;
> > float T2(Time, south_north, west_east) ;
> > T2:FieldType = 104 ;
> > T2:MemoryOrder = "XY " ;
> > T2:description = "TEMP at 2 M" ;
> > T2:units = "K" ;
> > T2:stagger = "" ;
> > T2:coordinates = "lon lat" ;
> > T2:long_name = "Air Temperature at 2 m" ;
> > T2:standard_name = "air_temperature" ;
> > T2:level = "2 m" ;
> > T2:_FillValue = 1.e+20f ;
> > float PSFC(Time, south_north, west_east) ;
> > PSFC:FieldType = 104 ;
> > PSFC:MemoryOrder = "XY " ;
> > PSFC:description = "SFC PRESSURE" ;
> > PSFC:units = "Pa" ;
> > PSFC:stagger = "" ;
> > PSFC:coordinates = "lon lat" ;
> > PSFC:long_name = "Air Pressure at the Surface" ;
> > PSFC:standard_name = "surface_pressure" ;
> > PSFC:level = "surface" ;
> > PSFC:_FillValue = 1.e+20f ;
> > float U10(Time, south_north, west_east) ;
> > U10:FieldType = 104 ;
> > U10:MemoryOrder = "XY " ;
> > U10:description = "U at 10 M" ;
> > U10:units = "m/s" ;
> > U10:stagger = "" ;
> > U10:coordinates = "lon lat" ;
> > U10:long_name = "u-Component of wind at 10 m
(Earth)" ;
> > U10:standard_name = "eastward_wind" ;
> > U10:level = "10 m" ;
> > U10:_FillValue = 1.e+20f ;
> > float V10(Time, south_north, west_east) ;
> > V10:FieldType = 104 ;
> > V10:MemoryOrder = "XY " ;
> > V10:description = "V at 10 M" ;
> > V10:units = "m/s" ;
> > V10:stagger = "" ;
> > V10:coordinates = "lon lat" ;
> > V10:long_name = "v-Component of wind at 10 m
(Earth)" ;
> > V10:standard_name = "northward_wind" ;
> > V10:level = "10 m" ;
> > V10:_FillValue = 1.e+20f ;
> > double time(Time) ;
> > time:standard_name = "time" ;
> > time:long_name = "Time" ;
> > time:calendar = "standard" ;
> > time:units = "hours since 1970-01-01 00:00:00" ;
> > float lat(south_north, west_east) ;
> > lat:long_name = "Latitude" ;
> > lat:standard_name = "latitude" ;
> > lat:units = "degrees_north" ;
> > float lon(south_north, west_east) ;
> > lon:long_name = "Longitude" ;
> > lon:standard_name = "longitude" ;
> > lon:units = "degrees_east" ;
> >
> > // global attributes:
> > :TITLE = " OUTPUT FROM WRF V4.1.1 MODEL" ;
> > :START_DATE = "2021-03-01_00:00:00" ;
> > :WEST-EAST_GRID_DIMENSION = 397 ;
> > :SOUTH-NORTH_GRID_DIMENSION = 481 ;
> > :BOTTOM-TOP_GRID_DIMENSION = 51 ;
> > :DX = 2000.f ;
> > :DY = 2000.f ;
> > :AERCU_OPT = 0 ;
> > :AERCU_FCT = 1.f ;
> > :IDEAL_CASE = 0 ;
> > :DIFF_6TH_SLOPEOPT = 0 ;
> > :AUTO_LEVELS_OPT = 2 ;
> > :DIFF_6TH_THRESH = 0.1f ;
> > :DZBOT = 50.f ;
> > :DZSTRETCH_S = 1.3f ;
> > :DZSTRETCH_U = 1.1f ;
> > :GRIDTYPE = "C" ;
> > :DIFF_OPT = 2 ;
> > :KM_OPT = 4 ;
> > :DAMP_OPT = 3 ;
> > :DAMPCOEF = 0.2f ;
> > :KHDIF = 0.f ;
> > :KVDIF = 0.f ;
> > :MP_PHYSICS = 8 ;
> > :RA_LW_PHYSICS = 4 ;
> > :RA_SW_PHYSICS = 4 ;
> > :SF_SFCLAY_PHYSICS = 5 ;
> > :SF_SURFACE_PHYSICS = 4 ;
> > :BL_PBL_PHYSICS = 6 ;
> > :CU_PHYSICS = 0 ;
> > :SF_LAKE_PHYSICS = 0 ;
> > :SURFACE_INPUT_SOURCE = 3 ;
> > :SST_UPDATE = 0 ;
> > :GRID_FDDA = 0 ;
> > :GFDDA_INTERVAL_M = 0 ;
> > :GFDDA_END_H = 0 ;
> > :GRID_SFDDA = 0 ;
> > :SGFDDA_INTERVAL_M = 0 ;
> > :SGFDDA_END_H = 0 ;
> > :HYPSOMETRIC_OPT = 2 ;
> > :USE_THETA_M = 1 ;
> > :GWD_OPT = 0 ;
> > :SF_URBAN_PHYSICS = 0 ;
> > :SF_SURFACE_MOSAIC = 0 ;
> > :SF_OCEAN_PHYSICS = 0 ;
> > :WEST-EAST_PATCH_START_UNSTAG = 1 ;
> > :WEST-EAST_PATCH_END_UNSTAG = 396 ;
> > :WEST-EAST_PATCH_START_STAG = 1 ;
> > :WEST-EAST_PATCH_END_STAG = 397 ;
> > :SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
> > :SOUTH-NORTH_PATCH_END_UNSTAG = 480 ;
> > :SOUTH-NORTH_PATCH_START_STAG = 1 ;
> > :SOUTH-NORTH_PATCH_END_STAG = 481 ;
> > :BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
> > :BOTTOM-TOP_PATCH_END_UNSTAG = 50 ;
> > :BOTTOM-TOP_PATCH_START_STAG = 1 ;
> > :BOTTOM-TOP_PATCH_END_STAG = 51 ;
> > :GRID_ID = 3 ;
> > :PARENT_ID = 2 ;
> > :I_PARENT_START = 93 ;
> > :J_PARENT_START = 75 ;
> > :PARENT_GRID_RATIO = 3 ;
> > :DT = 8.f ;
> > :CEN_LAT = 38.00961f ;
> > :CEN_LON = -121.0455f ;
> > :TRUELAT1 = 31.f ;
> > :TRUELAT2 = 45.38f ;
> > :MOAD_CEN_LAT = 38.19999f ;
> > :STAND_LON = -121.08f ;
> > :POLE_LAT = 90.f ;
> > :POLE_LON = 0.f ;
> > :GMT = 0.f ;
> > :JULYR = 2021 ;
> > :JULDAY = 60 ;
> > :MAP_PROJ = 1 ;
> > :MAP_PROJ_CHAR = "Lambert Conformal" ;
> > :MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
> > :NUM_LAND_CAT = 21 ;
> > :ISWATER = 17 ;
> > :ISLAKE = 21 ;
> > :ISICE = 15 ;
> > :ISURBAN = 13 ;
> > :ISOILWATER = 14 ;
> > :HYBRID_OPT = 2 ;
> > :ETAC = 0.2f ;
> > :file_type = "NETCDF_NCCF" ;
> > :title = "wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc" ;
> > :Conventions = "CF 1.6, Standard Name Table v19" ;
> > :source = "wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc" ;
> > :notes = "" ;
> > :created_by = "AJEB" ;
> > :institution = "PG&E" ;
> > :creation_date = "Mon Jun 7 17:59:23 UTC 2021" ;
> > :history = "Mon Jun 7 17:59:26 2021: ncap2 -O -s
> > T2F=float(T2F) wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
> > wrfsfc_d03_20210301_0000_f00000-cfcomp.nc\n",
> > "Mon Jun 7 17:59:26 2021: ncap2 -O -s
> > T2F=((T2-273.15)*(9.0/5.0))+32.0
> wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
> > wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
> > :NCO = "netCDF Operators version 4.9.8 (Homepage =
> >
>
https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fnco.s
>
f.net%2F&data=04%7C01%7CAJEB%40pge.com%7C0b860e74eccf45e2508508d92
>
ad54e1c%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C63758791790377891
>
1%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6
>
Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=siex3pLGh5m0T3wXYmbV%2BI4u42Jl
> QcyAfb07ZrRgN2c%3D&reserved=0,
> Code =
>
https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fgithub.com%2Fnco%2Fnco&data=04%7C01%7CAJEB%40pge.com%7C0b860e74eccf45e2508508d92ad54e1c%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C637587917903788907%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=Wob2rFlFVUI%2F9mJWPz7G2nq%2BcNcEoAUL642EH54JoVA%3D&reserved=0)"
> ;
> >
> >
> >
>
>
>
>
------------------------------------------------
Subject: unexpected results from point_stat with NEAREST GRIDPOINT method
From: Eiserloh Jr., A.J.
Time: Tue Jun 08 18:26:31 2021
netcdf wrfsfc_d03_20210301_0000_f00000-cfcomp {
dimensions:
Time = UNLIMITED ; // (1 currently)
south_north = 480 ;
west_east = 396 ;
DateStrLen = 19 ;
soil_layers_stag = 4 ;
variables:
float T2F(Time, south_north, west_east) ;
T2F:_FillValue = 1.e+20f ;
T2F:FieldType = 104 ;
T2F:MemoryOrder = "XY " ;
T2F:coordinates = "lon lat" ;
T2F:description = "TEMP at 2 M" ;
T2F:level = "2 m" ;
T2F:long_name = "Air Temperature at 2 m" ;
T2F:stagger = "" ;
T2F:standard_name = "air_temperature" ;
T2F:units = "K" ;
float WSPD_10M(Time, south_north, west_east) ;
WSPD_10M:_FillValue = 1.e+20f ;
WSPD_10M:FieldType = 104 ;
WSPD_10M:MemoryOrder = "XY " ;
WSPD_10M:coordinates = "lon lat" ;
WSPD_10M:description = "U at 10 M" ;
WSPD_10M:level = "10 m" ;
WSPD_10M:long_name = "Wind Speed at 10 m" ;
WSPD_10M:stagger = "" ;
WSPD_10M:standard_name = "wind_speed" ;
WSPD_10M:units = "m/s" ;
char Times(Time, DateStrLen) ;
float XLAT(Time, south_north, west_east) ;
XLAT:FieldType = 104 ;
XLAT:MemoryOrder = "XY " ;
XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
XLAT:units = "degrees_north" ;
XLAT:stagger = "" ;
XLAT:coordinates = "XLONG XLAT" ;
XLAT:long_name = "Latitude" ;
XLAT:standard_name = "latitude" ;
float XLONG(Time, south_north, west_east) ;
XLONG:FieldType = 104 ;
XLONG:MemoryOrder = "XY " ;
XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;
XLONG:units = "degrees_east" ;
XLONG:stagger = "" ;
XLONG:coordinates = "XLONG XLAT" ;
XLONG:long_name = "Longitude" ;
XLONG:standard_name = "longitude" ;
float Q2(Time, south_north, west_east) ;
Q2:FieldType = 104 ;
Q2:MemoryOrder = "XY " ;
Q2:description = "QV at 2 M" ;
Q2:units = "kg kg-1" ;
Q2:stagger = "" ;
Q2:coordinates = "lon lat" ;
Q2:long_name = "Specific Humidity at 2 m" ;
Q2:standard_name = "specific_humidity" ;
Q2:level = "2 m" ;
Q2:_FillValue = 1.e+20f ;
float T2(Time, south_north, west_east) ;
T2:FieldType = 104 ;
T2:MemoryOrder = "XY " ;
T2:description = "TEMP at 2 M" ;
T2:units = "K" ;
T2:stagger = "" ;
T2:coordinates = "lon lat" ;
T2:long_name = "Air Temperature at 2 m" ;
T2:standard_name = "air_temperature" ;
T2:level = "2 m" ;
T2:_FillValue = 1.e+20f ;
float PSFC(Time, south_north, west_east) ;
PSFC:FieldType = 104 ;
PSFC:MemoryOrder = "XY " ;
PSFC:description = "SFC PRESSURE" ;
PSFC:units = "Pa" ;
PSFC:stagger = "" ;
PSFC:coordinates = "lon lat" ;
PSFC:long_name = "Air Pressure at the Surface" ;
PSFC:standard_name = "surface_pressure" ;
PSFC:level = "surface" ;
PSFC:_FillValue = 1.e+20f ;
float U10(Time, south_north, west_east) ;
U10:FieldType = 104 ;
U10:MemoryOrder = "XY " ;
U10:description = "U at 10 M" ;
U10:units = "m/s" ;
U10:stagger = "" ;
U10:coordinates = "lon lat" ;
U10:long_name = "u-Component of wind at 10 m (Earth)"
;
U10:standard_name = "eastward_wind" ;
U10:level = "10 m" ;
U10:_FillValue = 1.e+20f ;
float V10(Time, south_north, west_east) ;
V10:FieldType = 104 ;
V10:MemoryOrder = "XY " ;
V10:description = "V at 10 M" ;
V10:units = "m/s" ;
V10:stagger = "" ;
V10:coordinates = "lon lat" ;
V10:long_name = "v-Component of wind at 10 m (Earth)"
;
V10:standard_name = "northward_wind" ;
V10:level = "10 m" ;
V10:_FillValue = 1.e+20f ;
float XTIME(Time) ;
XTIME:FieldType = 104 ;
XTIME:MemoryOrder = "0 " ;
XTIME:description = "minutes since 2021-02-28
21:00:00" ;
XTIME:units = "minutes since 2021-02-28 21:00:00" ;
XTIME:stagger = "" ;
float TSLB(Time, soil_layers_stag, south_north, west_east) ;
TSLB:FieldType = 104 ;
TSLB:MemoryOrder = "XYZ" ;
TSLB:description = "SOIL TEMPERATURE" ;
TSLB:units = "K" ;
TSLB:stagger = "Z" ;
TSLB:coordinates = "lon lat" ;
TSLB:long_name = "Soil Temperature" ;
TSLB:standard_name = "soil_temperature" ;
TSLB:_FillValue = 1.e+20f ;
float SMOIS(Time, soil_layers_stag, south_north, west_east) ;
SMOIS:FieldType = 104 ;
SMOIS:MemoryOrder = "XYZ" ;
SMOIS:description = "SOIL MOISTURE" ;
SMOIS:units = "m3 m-3" ;
SMOIS:stagger = "Z" ;
SMOIS:coordinates = "lon lat" ;
SMOIS:long_name = "Soil Moisture" ;
SMOIS:standard_name = "soil_moisture" ;
SMOIS:_FillValue = 1.e+20f ;
float ACSNOM(Time, south_north, west_east) ;
ACSNOM:FieldType = 104 ;
ACSNOM:MemoryOrder = "XY " ;
ACSNOM:description = "ACCUMULATED MELTED SNOW" ;
ACSNOM:units = "kg m-2" ;
ACSNOM:stagger = "" ;
ACSNOM:coordinates = "lon lat" ;
ACSNOM:long_name = "DOWNWARD SHORT WAVE FLUX AT GROUND
SURFACE" ;
ACSNOM:standard_name = "incoming shortwave raditation
at surface" ;
ACSNOM:_FillValue = 1.e+20f ;
float SNOWH(Time, south_north, west_east) ;
SNOWH:FieldType = 104 ;
SNOWH:MemoryOrder = "XY " ;
SNOWH:description = "PHYSICAL SNOW DEPTH" ;
SNOWH:units = "m" ;
SNOWH:stagger = "" ;
SNOWH:coordinates = "lon lat" ;
SNOWH:long_name = "Physical Snow Depth" ;
SNOWH:standard_name = "snow_depth" ;
SNOWH:_FillValue = 1.e+20f ;
float SWDOWN(Time, south_north, west_east) ;
SWDOWN:FieldType = 104 ;
SWDOWN:MemoryOrder = "XY " ;
SWDOWN:description = "DOWNWARD SHORT WAVE FLUX AT
GROUND SURFACE" ;
SWDOWN:units = "W m-2" ;
SWDOWN:stagger = "" ;
SWDOWN:coordinates = "lon lat" ;
SWDOWN:long_name = "DOWNWARD SHORT WAVE FLUX AT GROUND
SURFACE" ;
SWDOWN:standard_name = "incoming shortwave raditation
at surface" ;
SWDOWN:_FillValue = 1.e+20f ;
float ZNT(Time, south_north, west_east) ;
ZNT:FieldType = 104 ;
ZNT:MemoryOrder = "XY " ;
ZNT:description = "TIME-VARYING ROUGHNESS LENGTH" ;
ZNT:units = "m" ;
ZNT:stagger = "" ;
ZNT:coordinates = "lon lat" ;
ZNT:long_name = "Time Varying Roughness Length" ;
ZNT:standard_name = "roughness_length" ;
ZNT:_FillValue = 1.e+20f ;
float UST(Time, south_north, west_east) ;
UST:FieldType = 104 ;
UST:MemoryOrder = "XY " ;
UST:description = "U* IN SIMILARITY THEORY" ;
UST:units = "m s-1" ;
UST:stagger = "" ;
UST:coordinates = "lon lat" ;
UST:long_name = "U* in Similatiry Theory / Friction
Velocity" ;
UST:standard_name = "friction_velocity" ;
UST:_FillValue = 1.e+20f ;
float PREC_ACC_C(Time, south_north, west_east) ;
PREC_ACC_C:FieldType = 104 ;
PREC_ACC_C:MemoryOrder = "XY " ;
PREC_ACC_C:description = "ACCUMULATED CUMULUS
PRECIPITATION OVER prec_acc_dt PERIODS OF TIME" ;
PREC_ACC_C:units = "mm" ;
PREC_ACC_C:stagger = "" ;
PREC_ACC_C:coordinates = "lon lat" ;
PREC_ACC_C:long_name = "ACCUMULATED CUMULUS
PRECIPITATION" ;
PREC_ACC_C:standard_name = "" ;
PREC_ACC_C:_FillValue = 1.e+20f ;
float PREC_ACC_NC(Time, south_north, west_east) ;
PREC_ACC_NC:FieldType = 104 ;
PREC_ACC_NC:MemoryOrder = "XY " ;
PREC_ACC_NC:description = "ACCUMULATED GRID SCALE
PRECIPITATION OVER prec_acc_dt PERIODS OF TIME" ;
PREC_ACC_NC:units = "mm" ;
PREC_ACC_NC:stagger = "" ;
PREC_ACC_NC:coordinates = "lon lat" ;
PREC_ACC_NC:long_name = "ACCUMULATED GRID SCALE
PRECIPITATION" ;
PREC_ACC_NC:standard_name = "" ;
PREC_ACC_NC:level = "surface" ;
PREC_ACC_NC:_FillValue = 1.e+20f ;
float SNOW_ACC_NC(Time, south_north, west_east) ;
SNOW_ACC_NC:FieldType = 104 ;
SNOW_ACC_NC:MemoryOrder = "XY " ;
SNOW_ACC_NC:description = "ACCUMULATED SNOW WATER
EQUIVALENT OVER prec_acc_dt PERIODS OF TIME" ;
SNOW_ACC_NC:units = "mm" ;
SNOW_ACC_NC:stagger = "" ;
SNOW_ACC_NC:coordinates = "lon lat" ;
SNOW_ACC_NC:long_name = "ACCUMULATED SNOW WATER
EQUIVALENT" ;
SNOW_ACC_NC:standard_name = "" ;
SNOW_ACC_NC:level = "surface" ;
SNOW_ACC_NC:_FillValue = 1.e+20f ;
double time(Time) ;
time:standard_name = "time" ;
time:long_name = "Time" ;
time:calendar = "standard" ;
time:units = "hours since 1970-01-01 00:00:00" ;
float lat(south_north, west_east) ;
lat:long_name = "Latitude" ;
lat:standard_name = "latitude" ;
lat:units = "degrees_north" ;
float lon(south_north, west_east) ;
lon:long_name = "Longitude" ;
lon:standard_name = "longitude" ;
lon:units = "degrees_east" ;
// global attributes:
:TITLE = " OUTPUT FROM WRF V4.1.1 MODEL" ;
:START_DATE = "2021-03-01_00:00:00" ;
:BOTTOM-TOP_GRID_DIMENSION = 51 ;
:DX = 2000.f ;
:DY = 2000.f ;
:AERCU_OPT = 0 ;
:AERCU_FCT = 1.f ;
:IDEAL_CASE = 0 ;
:DIFF_6TH_SLOPEOPT = 0 ;
:AUTO_LEVELS_OPT = 2 ;
:DIFF_6TH_THRESH = 0.1f ;
:DZBOT = 50.f ;
:DZSTRETCH_S = 1.3f ;
:DZSTRETCH_U = 1.1f ;
:GRIDTYPE = "C" ;
:DIFF_OPT = 2 ;
:KM_OPT = 4 ;
:DAMP_OPT = 3 ;
:DAMPCOEF = 0.2f ;
:KHDIF = 0.f ;
:KVDIF = 0.f ;
:MP_PHYSICS = 8 ;
:RA_LW_PHYSICS = 4 ;
:RA_SW_PHYSICS = 4 ;
:SF_SFCLAY_PHYSICS = 5 ;
:SF_SURFACE_PHYSICS = 4 ;
:BL_PBL_PHYSICS = 6 ;
:CU_PHYSICS = 0 ;
:SF_LAKE_PHYSICS = 0 ;
:SURFACE_INPUT_SOURCE = 3 ;
:SST_UPDATE = 0 ;
:GRID_FDDA = 0 ;
:GFDDA_INTERVAL_M = 0 ;
:GFDDA_END_H = 0 ;
:GRID_SFDDA = 0 ;
:SGFDDA_INTERVAL_M = 0 ;
:SGFDDA_END_H = 0 ;
:HYPSOMETRIC_OPT = 2 ;
:USE_THETA_M = 1 ;
:GWD_OPT = 0 ;
:SF_URBAN_PHYSICS = 0 ;
:SF_SURFACE_MOSAIC = 0 ;
:SF_OCEAN_PHYSICS = 0 ;
:WEST-EAST_PATCH_START_UNSTAG = 1 ;
:WEST-EAST_PATCH_END_UNSTAG = 396 ;
:WEST-EAST_PATCH_START_STAG = 1 ;
:WEST-EAST_PATCH_END_STAG = 397 ;
:SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
:SOUTH-NORTH_PATCH_END_UNSTAG = 480 ;
:SOUTH-NORTH_PATCH_START_STAG = 1 ;
:SOUTH-NORTH_PATCH_END_STAG = 481 ;
:BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
:BOTTOM-TOP_PATCH_END_UNSTAG = 50 ;
:BOTTOM-TOP_PATCH_START_STAG = 1 ;
:BOTTOM-TOP_PATCH_END_STAG = 51 ;
:GRID_ID = 3 ;
:PARENT_ID = 2 ;
:I_PARENT_START = 93 ;
:J_PARENT_START = 75 ;
:PARENT_GRID_RATIO = 3 ;
:DT = 8.f ;
:CEN_LAT = 38.00961f ;
:CEN_LON = -121.0455f ;
:TRUELAT1 = 31.f ;
:TRUELAT2 = 45.38f ;
:MOAD_CEN_LAT = 38.19999f ;
:STAND_LON = -121.08f ;
:POLE_LAT = 90.f ;
:POLE_LON = 0.f ;
:GMT = 0.f ;
:JULYR = 2021 ;
:JULDAY = 60 ;
:MAP_PROJ = 1 ;
:MAP_PROJ_CHAR = "Lambert Conformal" ;
:MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
:NUM_LAND_CAT = 21 ;
:ISWATER = 17 ;
:ISLAKE = 21 ;
:ISICE = 15 ;
:ISURBAN = 13 ;
:ISOILWATER = 14 ;
:HYBRID_OPT = 2 ;
:ETAC = 0.2f ;
:file_type = "NETCDF_NCCF" ;
:title = "wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
:Conventions = "CF 1.6, Standard Name Table v19" ;
:source = "wrfsfc_d03_20210301_0000_f00000-cfcomp.nc"
;
:notes = "" ;
:created_by = "AJEB" ;
:institution = "PG&E" ;
:creation_date = "Tue Jun 8 23:15:57 UTC 2021" ;
:history = "Tue Jun 8 23:16:00 2021: ncatted -O -a
SOUTH-NORTH_GRID_DIMENSION,global,d,, wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc\n",
"Tue Jun 8 23:16:00 2021: ncatted -O -a WEST-
EAST_GRID_DIMENSION,global,d,, wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc\n",
"Tue Jun 8 23:16:00 2021: ncap2 -O -s
T2F=float(T2F) wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
wrfsfc_d03_20210301_0000_f00000-cfcomp.nc\n",
"Tue Jun 8 23:16:00 2021: ncap2 -O -s
T2F=((T2-273.15)*(9.0/5.0))+32.0 wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
:NCO = "netCDF Operators version 4.9.8 (Homepage =
http://nco.sf.net, Code = http://github.com/nco/nco)" ;
------------------------------------------------
Subject: unexpected results from point_stat with NEAREST GRIDPOINT method
From: John Halley Gotway
Time: Fri Jun 11 17:15:13 2021
AJ,
Sorry for the delay in responding. Ah yes, I think you're onto
something
there... that the issue is in the definition of the projection.
So when MET sees the following global attribute, that tells it to
interpret
this data as being CF-compliant:
* :Conventions = "CF 1.6, Standard Name Table v19" ;*
However, we have a hodge-podge of conventions going on in this file.
This
is a NetCDF WRFOUT in which the projection information is written in a
non-CF-compliant way to the global attributes. Those definitely aren't
read
by MET's library code that processes CF-compliant data. Here's an
example
of how a Lambert Conformal projection is defined in a CF-compliant
way:
https://cfconventions.org/Data/cf-conventions/cf-conventions-1.8/cf-
conventions.html#lambert-conformal-projection
Where the grid_mapping variable attribute defines the name of the
variable
that defines the projection for each variable.
The issue is that you've added the CF conventions attribute without
actually making the file fully CF-compliant.
So what's going on? If there is no projection info defined in a
CF-compliant way, then MET assumes the data is on a lat/lon grid and
attempts to determine that grid by reading the lat/lon variables. I'm
surprised that that didn't produce any warnings or errors... but if
your
data is on a Lambert Conformal grid and not lat/lon then that would
explain
the discrepancy in the interpolated values.
So what should you do? Honestly, I find specifying projection info for
the
CF-convention to be difficult. Instead, I'd recommend that you post-
process
your WRFOUT file using either the Unified Post Processor or the
wrf_interp
utility:
https://www2.mmm.ucar.edu/wrf/users/utilities/util.htm
That'd avoid you having to struggle through all these details on your
own.
UPP writes GRIB2 output which MET reads well. The wrf_interp utility
writes
a slightly modified version of NetCDF that MET can read... but is
limited
in that it can't read staggered variables, like winds.
If all you actually care about it 10M windspeed, you could try the
following command on the raw WRFOUT file:
plot_data_plane wrfsfc_d03_20210301_0000_f00000-cfcomp.nc out.ps
'name="WSPD_10M";
level="(0,*,*)"; file_type=NETCDF_PINT;'
The file_type=NETCDF_PINT; option tells MET to interpret the data as
being
output from the wrf_interp utility (previously called pinterp, thus
the
name 'pint'). And that might work.
If possible, I'd definitely recommend post-processing with UPP. If
not,
using wrf_interp would be my second choice. If not, then I'd recommend
writing a python script to read/format the data to be passed to MET in
memory. And only then, my last choice would be to reformat the NetCDF
file
manually.
Hope that helps.
Thanks,
John
On Tue, Jun 8, 2021 at 6:27 PM Eiserloh Jr., A.J. via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=100132 >
>
> Hi John,
>
> I got that to work this time, and the image (attached) looks correct
and
> the map looks good with the data, but not zoomed into my point is
hard for
> me to tell:
>
> Seems to me like it is interpolating still in some direction, or it
is
> changing my map projection maybe because I didn't provide enough
data?
>
> Here is the full ncdump of my netdf file below if that helps. I did
notice
> that with a grib2 file of similar model type, it gives me the
expected
> value at the point after running point_stat. And that grib2 file has
all
> the projection info associated with the lat and lon variables. But
with my
> netcdf file that I created manually to makeCF compliant (because
that is
> what is needed for met tool), it is not giving me the expected value
at the
> nearest location gridpoint as I described, and my netcdf does not
have the
> projection info as attributes tied to the lat and lon variables.
Rather it
> is in the global attributes and I am wondering if Met Tool can't
read that
> info about the projection?
>
> Thanks,
> AJ
>
> netcdf wrfsfc_d03_20210301_0000_f00000-cfcomp {
> dimensions:
> Time = UNLIMITED ; // (1 currently)
> south_north = 480 ;
> west_east = 396 ;
> DateStrLen = 19 ;
> soil_layers_stag = 4 ;
> variables:
> float T2F(Time, south_north, west_east) ;
> T2F:_FillValue = 1.e+20f ;
> T2F:FieldType = 104 ;
> T2F:MemoryOrder = "XY " ;
> T2F:coordinates = "lon lat" ;
> T2F:description = "TEMP at 2 M" ;
> T2F:level = "2 m" ;
> T2F:long_name = "Air Temperature at 2 m" ;
> T2F:stagger = "" ;
> T2F:standard_name = "air_temperature" ;
> T2F:units = "K" ;
> float WSPD_10M(Time, south_north, west_east) ;
> WSPD_10M:_FillValue = 1.e+20f ;
> WSPD_10M:FieldType = 104 ;
> WSPD_10M:MemoryOrder = "XY " ;
> WSPD_10M:coordinates = "lon lat" ;
> WSPD_10M:description = "U at 10 M" ;
> WSPD_10M:level = "10 m" ;
> WSPD_10M:long_name = "Wind Speed at 10 m" ;
> WSPD_10M:stagger = "" ;
> WSPD_10M:standard_name = "wind_speed" ;
> WSPD_10M:units = "m/s" ;
> char Times(Time, DateStrLen) ;
> float XLAT(Time, south_north, west_east) ;
> XLAT:FieldType = 104 ;
> XLAT:MemoryOrder = "XY " ;
> XLAT:description = "LATITUDE, SOUTH IS NEGATIVE" ;
> XLAT:units = "degrees_north" ;
> XLAT:stagger = "" ;
> XLAT:coordinates = "XLONG XLAT" ;
> XLAT:long_name = "Latitude" ;
> XLAT:standard_name = "latitude" ;
> float XLONG(Time, south_north, west_east) ;
> XLONG:FieldType = 104 ;
> XLONG:MemoryOrder = "XY " ;
> XLONG:description = "LONGITUDE, WEST IS NEGATIVE" ;
> XLONG:units = "degrees_east" ;
> XLONG:stagger = "" ;
> XLONG:coordinates = "XLONG XLAT" ;
> XLONG:long_name = "Longitude" ;
> XLONG:standard_name = "longitude" ;
> float Q2(Time, south_north, west_east) ;
> Q2:FieldType = 104 ;
> Q2:MemoryOrder = "XY " ;
> Q2:description = "QV at 2 M" ;
> Q2:units = "kg kg-1" ;
> Q2:stagger = "" ;
> Q2:coordinates = "lon lat" ;
> Q2:long_name = "Specific Humidity at 2 m" ;
> Q2:standard_name = "specific_humidity" ;
> Q2:level = "2 m" ;
> Q2:_FillValue = 1.e+20f ;
> float T2(Time, south_north, west_east) ;
> T2:FieldType = 104 ;
> T2:MemoryOrder = "XY " ;
> T2:description = "TEMP at 2 M" ;
> T2:units = "K" ;
> T2:stagger = "" ;
> T2:coordinates = "lon lat" ;
> T2:long_name = "Air Temperature at 2 m" ;
> T2:standard_name = "air_temperature" ;
> T2:level = "2 m" ;
> T2:_FillValue = 1.e+20f ;
> float PSFC(Time, south_north, west_east) ;
> PSFC:FieldType = 104 ;
> PSFC:MemoryOrder = "XY " ;
> PSFC:description = "SFC PRESSURE" ;
> PSFC:units = "Pa" ;
> PSFC:stagger = "" ;
> PSFC:coordinates = "lon lat" ;
> PSFC:long_name = "Air Pressure at the Surface" ;
> PSFC:standard_name = "surface_pressure" ;
> PSFC:level = "surface" ;
> PSFC:_FillValue = 1.e+20f ;
> float U10(Time, south_north, west_east) ;
> U10:FieldType = 104 ;
> U10:MemoryOrder = "XY " ;
> U10:description = "U at 10 M" ;
> U10:units = "m/s" ;
> U10:stagger = "" ;
> U10:coordinates = "lon lat" ;
> U10:long_name = "u-Component of wind at 10 m
(Earth)" ;
> U10:standard_name = "eastward_wind" ;
> U10:level = "10 m" ;
> U10:_FillValue = 1.e+20f ;
> float V10(Time, south_north, west_east) ;
> V10:FieldType = 104 ;
> V10:MemoryOrder = "XY " ;
> V10:description = "V at 10 M" ;
> V10:units = "m/s" ;
> V10:stagger = "" ;
> V10:coordinates = "lon lat" ;
> V10:long_name = "v-Component of wind at 10 m
(Earth)" ;
> V10:standard_name = "northward_wind" ;
> V10:level = "10 m" ;
> V10:_FillValue = 1.e+20f ;
> float XTIME(Time) ;
> XTIME:FieldType = 104 ;
> XTIME:MemoryOrder = "0 " ;
> XTIME:description = "minutes since 2021-02-28
21:00:00" ;
> XTIME:units = "minutes since 2021-02-28 21:00:00" ;
> XTIME:stagger = "" ;
> float TSLB(Time, soil_layers_stag, south_north, west_east) ;
> TSLB:FieldType = 104 ;
> TSLB:MemoryOrder = "XYZ" ;
> TSLB:description = "SOIL TEMPERATURE" ;
> TSLB:units = "K" ;
> TSLB:stagger = "Z" ;
> TSLB:coordinates = "lon lat" ;
> TSLB:long_name = "Soil Temperature" ;
> TSLB:standard_name = "soil_temperature" ;
> TSLB:_FillValue = 1.e+20f ;
> float SMOIS(Time, soil_layers_stag, south_north, west_east)
;
> SMOIS:FieldType = 104 ;
> SMOIS:MemoryOrder = "XYZ" ;
> SMOIS:description = "SOIL MOISTURE" ;
> SMOIS:units = "m3 m-3" ;
> SMOIS:stagger = "Z" ;
> SMOIS:coordinates = "lon lat" ;
> SMOIS:long_name = "Soil Moisture" ;
> SMOIS:standard_name = "soil_moisture" ;
> SMOIS:_FillValue = 1.e+20f ;
> float ACSNOM(Time, south_north, west_east) ;
> ACSNOM:FieldType = 104 ;
> ACSNOM:MemoryOrder = "XY " ;
> ACSNOM:description = "ACCUMULATED MELTED SNOW" ;
> ACSNOM:units = "kg m-2" ;
> ACSNOM:stagger = "" ;
> ACSNOM:coordinates = "lon lat" ;
> ACSNOM:long_name = "DOWNWARD SHORT WAVE FLUX AT
GROUND
> SURFACE" ;
> ACSNOM:standard_name = "incoming shortwave
raditation at
> surface" ;
> ACSNOM:_FillValue = 1.e+20f ;
> float SNOWH(Time, south_north, west_east) ;
> SNOWH:FieldType = 104 ;
> SNOWH:MemoryOrder = "XY " ;
> SNOWH:description = "PHYSICAL SNOW DEPTH" ;
> SNOWH:units = "m" ;
> SNOWH:stagger = "" ;
> SNOWH:coordinates = "lon lat" ;
> SNOWH:long_name = "Physical Snow Depth" ;
> SNOWH:standard_name = "snow_depth" ;
> SNOWH:_FillValue = 1.e+20f ;
> float SWDOWN(Time, south_north, west_east) ;
> SWDOWN:FieldType = 104 ;
> SWDOWN:MemoryOrder = "XY " ;
> SWDOWN:description = "DOWNWARD SHORT WAVE FLUX AT
GROUND
> SURFACE" ;
> SWDOWN:units = "W m-2" ;
> SWDOWN:stagger = "" ;
> SWDOWN:coordinates = "lon lat" ;
> SWDOWN:long_name = "DOWNWARD SHORT WAVE FLUX AT
GROUND
> SURFACE" ;
> SWDOWN:standard_name = "incoming shortwave
raditation at
> surface" ;
> SWDOWN:_FillValue = 1.e+20f ;
> float ZNT(Time, south_north, west_east) ;
> ZNT:FieldType = 104 ;
> ZNT:MemoryOrder = "XY " ;
> ZNT:description = "TIME-VARYING ROUGHNESS LENGTH" ;
> ZNT:units = "m" ;
> ZNT:stagger = "" ;
> ZNT:coordinates = "lon lat" ;
> ZNT:long_name = "Time Varying Roughness Length" ;
> ZNT:standard_name = "roughness_length" ;
> ZNT:_FillValue = 1.e+20f ;
> float UST(Time, south_north, west_east) ;
> UST:FieldType = 104 ;
> UST:MemoryOrder = "XY " ;
> UST:description = "U* IN SIMILARITY THEORY" ;
> UST:units = "m s-1" ;
> UST:stagger = "" ;
> UST:coordinates = "lon lat" ;
> UST:long_name = "U* in Similatiry Theory / Friction
> Velocity" ;
> UST:standard_name = "friction_velocity" ;
> UST:_FillValue = 1.e+20f ;
> float PREC_ACC_C(Time, south_north, west_east) ;
> PREC_ACC_C:FieldType = 104 ;
> PREC_ACC_C:MemoryOrder = "XY " ;
> PREC_ACC_C:description = "ACCUMULATED CUMULUS
> PRECIPITATION OVER prec_acc_dt PERIODS OF TIME" ;
> PREC_ACC_C:units = "mm" ;
> PREC_ACC_C:stagger = "" ;
> PREC_ACC_C:coordinates = "lon lat" ;
> PREC_ACC_C:long_name = "ACCUMULATED CUMULUS
PRECIPITATION"
> ;
> PREC_ACC_C:standard_name = "" ;
> PREC_ACC_C:_FillValue = 1.e+20f ;
> float PREC_ACC_NC(Time, south_north, west_east) ;
> PREC_ACC_NC:FieldType = 104 ;
> PREC_ACC_NC:MemoryOrder = "XY " ;
> PREC_ACC_NC:description = "ACCUMULATED GRID SCALE
> PRECIPITATION OVER prec_acc_dt PERIODS OF TIME" ;
> PREC_ACC_NC:units = "mm" ;
> PREC_ACC_NC:stagger = "" ;
> PREC_ACC_NC:coordinates = "lon lat" ;
> PREC_ACC_NC:long_name = "ACCUMULATED GRID SCALE
> PRECIPITATION" ;
> PREC_ACC_NC:standard_name = "" ;
> PREC_ACC_NC:level = "surface" ;
> PREC_ACC_NC:_FillValue = 1.e+20f ;
> float SNOW_ACC_NC(Time, south_north, west_east) ;
> SNOW_ACC_NC:FieldType = 104 ;
> SNOW_ACC_NC:MemoryOrder = "XY " ;
> SNOW_ACC_NC:description = "ACCUMULATED SNOW WATER
> EQUIVALENT OVER prec_acc_dt PERIODS OF TIME" ;
> SNOW_ACC_NC:units = "mm" ;
> SNOW_ACC_NC:stagger = "" ;
> SNOW_ACC_NC:coordinates = "lon lat" ;
> SNOW_ACC_NC:long_name = "ACCUMULATED SNOW WATER
> EQUIVALENT" ;
> SNOW_ACC_NC:standard_name = "" ;
> SNOW_ACC_NC:level = "surface" ;
> SNOW_ACC_NC:_FillValue = 1.e+20f ;
> double time(Time) ;
> time:standard_name = "time" ;
> time:long_name = "Time" ;
> time:calendar = "standard" ;
> time:units = "hours since 1970-01-01 00:00:00" ;
> float lat(south_north, west_east) ;
> lat:long_name = "Latitude" ;
> lat:standard_name = "latitude" ;
> lat:units = "degrees_north" ;
> float lon(south_north, west_east) ;
> lon:long_name = "Longitude" ;
> lon:standard_name = "longitude" ;
> lon:units = "degrees_east" ;
>
> // global attributes:
> :TITLE = " OUTPUT FROM WRF V4.1.1 MODEL" ;
> :START_DATE = "2021-03-01_00:00:00" ;
> :BOTTOM-TOP_GRID_DIMENSION = 51 ;
> :DX = 2000.f ;
> :DY = 2000.f ;
> :AERCU_OPT = 0 ;
> :AERCU_FCT = 1.f ;
> :IDEAL_CASE = 0 ;
> :DIFF_6TH_SLOPEOPT = 0 ;
> :AUTO_LEVELS_OPT = 2 ;
> :DIFF_6TH_THRESH = 0.1f ;
> :DZBOT = 50.f ;
> :DZSTRETCH_S = 1.3f ;
> :DZSTRETCH_U = 1.1f ;
> :GRIDTYPE = "C" ;
> :DIFF_OPT = 2 ;
> :KM_OPT = 4 ;
> :DAMP_OPT = 3 ;
> :DAMPCOEF = 0.2f ;
> :KHDIF = 0.f ;
> :KVDIF = 0.f ;
> :MP_PHYSICS = 8 ;
> :RA_LW_PHYSICS = 4 ;
> :RA_SW_PHYSICS = 4 ;
> :SF_SFCLAY_PHYSICS = 5 ;
> :SF_SURFACE_PHYSICS = 4 ;
> :BL_PBL_PHYSICS = 6 ;
> :CU_PHYSICS = 0 ;
> :SF_LAKE_PHYSICS = 0 ;
> :SURFACE_INPUT_SOURCE = 3 ;
> :SST_UPDATE = 0 ;
> :GRID_FDDA = 0 ;
> :GFDDA_INTERVAL_M = 0 ;
> :GFDDA_END_H = 0 ;
> :GRID_SFDDA = 0 ;
> :SGFDDA_INTERVAL_M = 0 ;
> :SGFDDA_END_H = 0 ;
> :HYPSOMETRIC_OPT = 2 ;
> :USE_THETA_M = 1 ;
> :GWD_OPT = 0 ;
> :SF_URBAN_PHYSICS = 0 ;
> :SF_SURFACE_MOSAIC = 0 ;
> :SF_OCEAN_PHYSICS = 0 ;
> :WEST-EAST_PATCH_START_UNSTAG = 1 ;
> :WEST-EAST_PATCH_END_UNSTAG = 396 ;
> :WEST-EAST_PATCH_START_STAG = 1 ;
> :WEST-EAST_PATCH_END_STAG = 397 ;
> :SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
> :SOUTH-NORTH_PATCH_END_UNSTAG = 480 ;
> :SOUTH-NORTH_PATCH_START_STAG = 1 ;
> :SOUTH-NORTH_PATCH_END_STAG = 481 ;
> :BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
> :BOTTOM-TOP_PATCH_END_UNSTAG = 50 ;
> :BOTTOM-TOP_PATCH_START_STAG = 1 ;
> :BOTTOM-TOP_PATCH_END_STAG = 51 ;
> :GRID_ID = 3 ;
> :PARENT_ID = 2 ;
> :I_PARENT_START = 93 ;
> :J_PARENT_START = 75 ;
> :PARENT_GRID_RATIO = 3 ;
> :DT = 8.f ;
> :CEN_LAT = 38.00961f ;
> :CEN_LON = -121.0455f ;
> :TRUELAT1 = 31.f ;
> :TRUELAT2 = 45.38f ;
> :MOAD_CEN_LAT = 38.19999f ;
> :STAND_LON = -121.08f ;
> :POLE_LAT = 90.f ;
> :POLE_LON = 0.f ;
> :GMT = 0.f ;
> :JULYR = 2021 ;
> :JULDAY = 60 ;
> :MAP_PROJ = 1 ;
> :MAP_PROJ_CHAR = "Lambert Conformal" ;
> :MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
> :NUM_LAND_CAT = 21 ;
> :ISWATER = 17 ;
> :ISLAKE = 21 ;
> :ISICE = 15 ;
> :ISURBAN = 13 ;
> :ISOILWATER = 14 ;
> :HYBRID_OPT = 2 ;
> :ETAC = 0.2f ;
> :file_type = "NETCDF_NCCF" ;
> :title = "wrfsfc_d03_20210301_0000_f00000-cfcomp.nc"
;
> :Conventions = "CF 1.6, Standard Name Table v19" ;
> :source = "wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc" ;
> :notes = "" ;
> :created_by = "AJEB" ;
> :institution = "PG&E" ;
> :creation_date = "Tue Jun 8 23:15:57 UTC 2021" ;
> :history = "Tue Jun 8 23:16:00 2021: ncatted -O -a
> SOUTH-NORTH_GRID_DIMENSION,global,d,,
> wrfsfc_d03_20210301_0000_f00000-cfcomp.nc\n",
> "Tue Jun 8 23:16:00 2021: ncatted -O -a
> WEST-EAST_GRID_DIMENSION,global,d,,
> wrfsfc_d03_20210301_0000_f00000-cfcomp.nc\n",
> "Tue Jun 8 23:16:00 2021: ncap2 -O -s
> T2F=float(T2F) wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
> wrfsfc_d03_20210301_0000_f00000-cfcomp.nc\n",
> "Tue Jun 8 23:16:00 2021: ncap2 -O -s
> T2F=((T2-273.15)*(9.0/5.0))+32.0 wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc
> wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
> :NCO = "netCDF Operators version 4.9.8 (Homepage =
> http://nco.sf.net, Code = http://github.com/nco/nco)" ;
>
> -----Original Message-----
> From: John Halley Gotway via RT <met_help at ucar.edu>
> Sent: Tuesday, June 8, 2021 4:30 PM
> To: Eiserloh Jr., A.J. <AJEB at pge.com>
> Subject: Re: [rt.rap.ucar.edu #100132] unexpected results from
point_stat
> with NEAREST GRIDPOINT method
>
> *****CAUTION: This email was sent from an EXTERNAL source. Think
before
> clicking links or opening attachments.*****
>
> AJ,
>
> Please try this command:
>
> plot_data_plane wrfsfc_d03_20210301_0000_f00000-cfcomp.nc out.ps
> 'name="WSPD_10M"; level="(0,*,*)";'
>
> That last option is like a mini-inline configuration file that tells
the
> tool what data to read from that input file.
>
> John
>
> On Tue, Jun 8, 2021 at 5:04 PM Eiserloh Jr., A.J. via RT <
> met_help at ucar.edu>
> wrote:
>
> >
> > <URL:
> >
https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Frt.r
> >
ap.ucar.edu%2Frt%2FTicket%2FDisplay.html%3Fid%3D100132&data=04%7C0
> >
1%7CAJEB%40pge.com%7C0b860e74eccf45e2508508d92ad54e1c%7C44ae661aece641
> >
aabc967c2c85a08941%7C0%7C0%7C637587917903778911%7CUnknown%7CTWFpbGZsb3
> >
d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7
> >
C1000&sdata=eDLjuvbSEv8cosgHtciLwoyf9TaR%2BSTWkLRNejmd9MQ%3D&r
> > eserved=0 >
> >
> > Hi John,
> >
> > I am unable to ftp to send data to the ftp address unless I get a
> > firewall exception request from my company's IT, and that will
take
> > some time. I tried plotting with the plot_data_plane command but
get
> > another error that is not helpful.
> >
> > Thanks,
> > -AJ
> >
> > [root at 698df1a69b34 nccompli_test]# plot_data_plane
> > wrfsfc_d03_20210301_0000_f00000-cfcomp.nc out.ps WSPD_10M ERROR :
> > ERROR : yyerror() -> syntax error in file "/tmp/met_config_273_0"
> > ERROR :
> > ERROR : line = 2
> > ERROR :
> > ERROR : column = 8
> > ERROR :
> > ERROR : text = "(nul)"
> > ERROR :
> > ERROR :
> > ERROR : (nul)
> > ERROR : (nul)
> > ERROR :
> >
> >
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT <met_help at ucar.edu>
> > Sent: Monday, June 7, 2021 2:47 PM
> > To: Eiserloh Jr., A.J. <AJEB at pge.com>
> > Subject: Re: [rt.rap.ucar.edu #100132] unexpected results from
> > point_stat with NEAREST GRIDPOINT method
> >
> > *****CAUTION: This email was sent from an EXTERNAL source. Think
> > before clicking links or opening attachments.*****
> >
> > AJ,
> >
> > I see that you're not getting the expected results from the
nearest
> > neighbor interpolation method in point-stat. There's nothing I can
see
> > in the Point-Stat config file you sent that is an obvious problem.
> >
> > If I had the actual model/obs data available to me, I would start
by
> > running the plot_data_plane tool to plot the 10-m wind data from
your
> > NetCDF file.
> >
> > *plot_data_plane wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
> > <
> >
https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwrfsf
> > c_d03_20210301_0000_f00000-
cfcomp.nc%2F&data=04%7C01%7CAJEB%40pge.
> >
com%7C0b860e74eccf45e2508508d92ad54e1c%7C44ae661aece641aabc967c2c85a08
> >
941%7C0%7C0%7C637587917903778911%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wL
> >
jAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata
> >
=8H2KGfAsMk%2BQocaxz75QIkeAgE%2Fctl9i0w99nQs0HuA%3D&reserved=0>
> > WSPD_10M.ps 'name="WSPD_10M"; level="(0,*,*)";'*
> >
> > And I'd check to make sure that the data isn't upside-down or
backwards.
> > If that looks good, I might also run the plot_point_obs tool. You
can
> > actually plot WSPD_10M field from above as a base image and then
> > overlay your point.
> > Check to make sure that point is really in the right spot on the
grid...
> > confirm that sign of the longitude is correct... etc.
> >
> > If that all checks out, next I'd run point_stat through a debugger
to
> > see line-by-line what's going on.
> >
> > But I'd need sample data to debug in this way. If you'd like to
> > provide that, you could post it to our anonymous FTP site
following
> > these
> > instructions:
> >
> >
https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fdtce
> > nter.org%2Fcommunity-code%2Fmodel-evaluation-tools-met%2Fmet-help-
desk
> >
%23ftp&data=04%7C01%7CAJEB%40pge.com%7C0b860e74eccf45e2508508d92ad
> >
54e1c%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C637587917903778911%
> >
7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik
> >
1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=MqE2mOm9CNXhqn%2F1zAUcsLrLBRWkmc
> > LLY%2FCp3Tcu34g%3D&reserved=0
> >
> > Thanks,
> > John Halley Gotway
> >
> > On Mon, Jun 7, 2021 at 12:45 PM Eiserloh Jr., A.J. via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > >
> > > Mon Jun 07 12:45:45 2021: Request 100132 was acted upon.
> > > Transaction: Ticket created by AJEB at pge.com
> > > Queue: met_help
> > > Subject: unexpected results from point_stat with NEAREST
> > > GRIDPOINT method
> > > Owner: Nobody
> > > Requestors: AJEB at pge.com
> > > Status: new
> > > Ticket <URL:
> > >
https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Frt
> > > .r
> > >
ap.ucar.edu%2Frt%2FTicket%2FDisplay.html%3Fid%3D100132&data=04%7
> > > C0
> > >
1%7CAJEB%40pge.com%7Cae1ed37b041e47bfa3a708d929fdd547%7C44ae661aece6
> > > 41
> > >
aabc967c2c85a08941%7C0%7C0%7C637586992454384982%7CUnknown%7CTWFpbGZs
> > > b3
> > >
d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D
> > > %7
> > >
C1000&sdata=LuRAF64SY4mRnzYWzWw5YN8%2FKHVDJTkwzWj%2BolIhDmA%3D&a
> > > mp
> > > ;reserved=0 >
> > >
> > >
> > > Dear Met-Help,
> > >
> > > I am testing out the point_stat tool (V10.0.0) for 1 station
pair,
> > > and trying to see if it is grabbing the correct value of
windspeed
> > > from the nearest gridpoint in my input forecast netcdf grid. In
my
> > > ascii obs text I purposefully put the station lat/lon as one of
the
> > > gridpoints in my forecast field for lat lon indexes (starting
with 0
> > > in the ncl fashion) of 218,140. The lat/lon that this
corresponds to
> > > is (37.61275, -122.3618). The 10-m windspeed at this point is
> > > 3.157048 m/s, however the point_stat*_mpr.txt results of the
matched
> > > pair has the windspeed as
> > > 3.57623 m/s, which is not correct. In fact, none of the
surrounding
> > > points have that value either. Why is point stat giving me a
> > > different value? I don't think that it would be doing any kind
of
> > > vertical interpolation, but not sure. My expectation is that it
should
> not be.
> > > I attached my PointStatConfig file, and the other info I
provided
> > > below
> > that may help.
> > >
> > > Thanks,
> > >
> > > A.J. Eiserloh - Data Scientist
> > > PG&E - Meteorology Systems and Analytics Applied Technology
Services
> > > 3400 Crow Canyon Rd., San Ramon, CA 94583-1393
> > > 925-307-4492
> > >
> > >
> > >
> > >
> > > * Here is my pre ascii2nc ascii input columns for the point:
> > >
> > > ADPSFC KSFO 20210301_000000 37.61275 -122.3618 0.0 32 0 10.0 NA
> > > 11.318018760716
> > >
> > >
> > > * The *mpr.txt results for the one matched pair:
> > >
> > > OBS_LAT OBS_LON OBS_LVL OBS_ELV FCST
> > > OBS OBS_QC CLIMO_MEAN CLIMO_STDEV CLIMO_CDF
> > > 37.61275 -122.3618 NA 10
> > > 3.57623 11.31802 NA NA
> > > NA NA
> > >
> > >
> > > * ncdump of my netcdf file:
> > >
> > > netcdf wrfsfc_d03_20210301_0000_f00000-cfcomp {
> > > dimensions:
> > > Time = UNLIMITED ; // (1 currently)
> > > south_north = 480 ;
> > > west_east = 396 ;
> > > DateStrLen = 19 ;
> > > soil_layers_stag = 4 ;
> > > variables:
> > > float WSPD_10M(Time, south_north, west_east) ;
> > > WSPD_10M:_FillValue = 1.e+20f ;
> > > WSPD_10M:FieldType = 104 ;
> > > WSPD_10M:MemoryOrder = "XY " ;
> > > WSPD_10M:coordinates = "lon lat" ;
> > > WSPD_10M:description = "U at 10 M" ;
> > > WSPD_10M:level = "10 m" ;
> > > WSPD_10M:long_name = "Wind Speed at 10 m" ;
> > > WSPD_10M:stagger = "" ;
> > > WSPD_10M:standard_name = "wind_speed" ;
> > > WSPD_10M:units = "m/s" ;
> > > char Times(Time, DateStrLen) ;
> > > float XLAT(Time, south_north, west_east) ;
> > > XLAT:FieldType = 104 ;
> > > XLAT:MemoryOrder = "XY " ;
> > > XLAT:description = "LATITUDE, SOUTH IS NEGATIVE"
;
> > > XLAT:units = "degrees_north" ;
> > > XLAT:stagger = "" ;
> > > XLAT:coordinates = "XLONG XLAT" ;
> > > XLAT:long_name = "Latitude" ;
> > > XLAT:standard_name = "latitude" ;
> > > float XLONG(Time, south_north, west_east) ;
> > > XLONG:FieldType = 104 ;
> > > XLONG:MemoryOrder = "XY " ;
> > > XLONG:description = "LONGITUDE, WEST IS
NEGATIVE" ;
> > > XLONG:units = "degrees_east" ;
> > > XLONG:stagger = "" ;
> > > XLONG:coordinates = "XLONG XLAT" ;
> > > XLONG:long_name = "Longitude" ;
> > > XLONG:standard_name = "longitude" ;
> > > float Q2(Time, south_north, west_east) ;
> > > Q2:FieldType = 104 ;
> > > Q2:MemoryOrder = "XY " ;
> > > Q2:description = "QV at 2 M" ;
> > > Q2:units = "kg kg-1" ;
> > > Q2:stagger = "" ;
> > > Q2:coordinates = "lon lat" ;
> > > Q2:long_name = "Specific Humidity at 2 m" ;
> > > Q2:standard_name = "specific_humidity" ;
> > > Q2:level = "2 m" ;
> > > Q2:_FillValue = 1.e+20f ;
> > > float T2(Time, south_north, west_east) ;
> > > T2:FieldType = 104 ;
> > > T2:MemoryOrder = "XY " ;
> > > T2:description = "TEMP at 2 M" ;
> > > T2:units = "K" ;
> > > T2:stagger = "" ;
> > > T2:coordinates = "lon lat" ;
> > > T2:long_name = "Air Temperature at 2 m" ;
> > > T2:standard_name = "air_temperature" ;
> > > T2:level = "2 m" ;
> > > T2:_FillValue = 1.e+20f ;
> > > float PSFC(Time, south_north, west_east) ;
> > > PSFC:FieldType = 104 ;
> > > PSFC:MemoryOrder = "XY " ;
> > > PSFC:description = "SFC PRESSURE" ;
> > > PSFC:units = "Pa" ;
> > > PSFC:stagger = "" ;
> > > PSFC:coordinates = "lon lat" ;
> > > PSFC:long_name = "Air Pressure at the Surface" ;
> > > PSFC:standard_name = "surface_pressure" ;
> > > PSFC:level = "surface" ;
> > > PSFC:_FillValue = 1.e+20f ;
> > > float U10(Time, south_north, west_east) ;
> > > U10:FieldType = 104 ;
> > > U10:MemoryOrder = "XY " ;
> > > U10:description = "U at 10 M" ;
> > > U10:units = "m/s" ;
> > > U10:stagger = "" ;
> > > U10:coordinates = "lon lat" ;
> > > U10:long_name = "u-Component of wind at 10 m
(Earth)" ;
> > > U10:standard_name = "eastward_wind" ;
> > > U10:level = "10 m" ;
> > > U10:_FillValue = 1.e+20f ;
> > > float V10(Time, south_north, west_east) ;
> > > V10:FieldType = 104 ;
> > > V10:MemoryOrder = "XY " ;
> > > V10:description = "V at 10 M" ;
> > > V10:units = "m/s" ;
> > > V10:stagger = "" ;
> > > V10:coordinates = "lon lat" ;
> > > V10:long_name = "v-Component of wind at 10 m
(Earth)" ;
> > > V10:standard_name = "northward_wind" ;
> > > V10:level = "10 m" ;
> > > V10:_FillValue = 1.e+20f ;
> > > double time(Time) ;
> > > time:standard_name = "time" ;
> > > time:long_name = "Time" ;
> > > time:calendar = "standard" ;
> > > time:units = "hours since 1970-01-01 00:00:00" ;
> > > float lat(south_north, west_east) ;
> > > lat:long_name = "Latitude" ;
> > > lat:standard_name = "latitude" ;
> > > lat:units = "degrees_north" ;
> > > float lon(south_north, west_east) ;
> > > lon:long_name = "Longitude" ;
> > > lon:standard_name = "longitude" ;
> > > lon:units = "degrees_east" ;
> > >
> > > // global attributes:
> > > :TITLE = " OUTPUT FROM WRF V4.1.1 MODEL" ;
> > > :START_DATE = "2021-03-01_00:00:00" ;
> > > :WEST-EAST_GRID_DIMENSION = 397 ;
> > > :SOUTH-NORTH_GRID_DIMENSION = 481 ;
> > > :BOTTOM-TOP_GRID_DIMENSION = 51 ;
> > > :DX = 2000.f ;
> > > :DY = 2000.f ;
> > > :AERCU_OPT = 0 ;
> > > :AERCU_FCT = 1.f ;
> > > :IDEAL_CASE = 0 ;
> > > :DIFF_6TH_SLOPEOPT = 0 ;
> > > :AUTO_LEVELS_OPT = 2 ;
> > > :DIFF_6TH_THRESH = 0.1f ;
> > > :DZBOT = 50.f ;
> > > :DZSTRETCH_S = 1.3f ;
> > > :DZSTRETCH_U = 1.1f ;
> > > :GRIDTYPE = "C" ;
> > > :DIFF_OPT = 2 ;
> > > :KM_OPT = 4 ;
> > > :DAMP_OPT = 3 ;
> > > :DAMPCOEF = 0.2f ;
> > > :KHDIF = 0.f ;
> > > :KVDIF = 0.f ;
> > > :MP_PHYSICS = 8 ;
> > > :RA_LW_PHYSICS = 4 ;
> > > :RA_SW_PHYSICS = 4 ;
> > > :SF_SFCLAY_PHYSICS = 5 ;
> > > :SF_SURFACE_PHYSICS = 4 ;
> > > :BL_PBL_PHYSICS = 6 ;
> > > :CU_PHYSICS = 0 ;
> > > :SF_LAKE_PHYSICS = 0 ;
> > > :SURFACE_INPUT_SOURCE = 3 ;
> > > :SST_UPDATE = 0 ;
> > > :GRID_FDDA = 0 ;
> > > :GFDDA_INTERVAL_M = 0 ;
> > > :GFDDA_END_H = 0 ;
> > > :GRID_SFDDA = 0 ;
> > > :SGFDDA_INTERVAL_M = 0 ;
> > > :SGFDDA_END_H = 0 ;
> > > :HYPSOMETRIC_OPT = 2 ;
> > > :USE_THETA_M = 1 ;
> > > :GWD_OPT = 0 ;
> > > :SF_URBAN_PHYSICS = 0 ;
> > > :SF_SURFACE_MOSAIC = 0 ;
> > > :SF_OCEAN_PHYSICS = 0 ;
> > > :WEST-EAST_PATCH_START_UNSTAG = 1 ;
> > > :WEST-EAST_PATCH_END_UNSTAG = 396 ;
> > > :WEST-EAST_PATCH_START_STAG = 1 ;
> > > :WEST-EAST_PATCH_END_STAG = 397 ;
> > > :SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
> > > :SOUTH-NORTH_PATCH_END_UNSTAG = 480 ;
> > > :SOUTH-NORTH_PATCH_START_STAG = 1 ;
> > > :SOUTH-NORTH_PATCH_END_STAG = 481 ;
> > > :BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
> > > :BOTTOM-TOP_PATCH_END_UNSTAG = 50 ;
> > > :BOTTOM-TOP_PATCH_START_STAG = 1 ;
> > > :BOTTOM-TOP_PATCH_END_STAG = 51 ;
> > > :GRID_ID = 3 ;
> > > :PARENT_ID = 2 ;
> > > :I_PARENT_START = 93 ;
> > > :J_PARENT_START = 75 ;
> > > :PARENT_GRID_RATIO = 3 ;
> > > :DT = 8.f ;
> > > :CEN_LAT = 38.00961f ;
> > > :CEN_LON = -121.0455f ;
> > > :TRUELAT1 = 31.f ;
> > > :TRUELAT2 = 45.38f ;
> > > :MOAD_CEN_LAT = 38.19999f ;
> > > :STAND_LON = -121.08f ;
> > > :POLE_LAT = 90.f ;
> > > :POLE_LON = 0.f ;
> > > :GMT = 0.f ;
> > > :JULYR = 2021 ;
> > > :JULDAY = 60 ;
> > > :MAP_PROJ = 1 ;
> > > :MAP_PROJ_CHAR = "Lambert Conformal" ;
> > > :MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
> > > :NUM_LAND_CAT = 21 ;
> > > :ISWATER = 17 ;
> > > :ISLAKE = 21 ;
> > > :ISICE = 15 ;
> > > :ISURBAN = 13 ;
> > > :ISOILWATER = 14 ;
> > > :HYBRID_OPT = 2 ;
> > > :ETAC = 0.2f ;
> > > :file_type = "NETCDF_NCCF" ;
> > > :title = "wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc" ;
> > > :Conventions = "CF 1.6, Standard Name Table v19"
;
> > > :source = "wrfsfc_d03_20210301_0000_f00000-
cfcomp.nc"
> ;
> > > :notes = "" ;
> > > :created_by = "AJEB" ;
> > > :institution = "PG&E" ;
> > > :creation_date = "Mon Jun 7 17:59:23 UTC 2021"
;
> > > :history = "Mon Jun 7 17:59:26 2021: ncap2 -O
-s
> > > T2F=float(T2F) wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
> > > wrfsfc_d03_20210301_0000_f00000-cfcomp.nc\n",
> > > "Mon Jun 7 17:59:26 2021: ncap2 -O -s
> > > T2F=((T2-273.15)*(9.0/5.0))+32.0
> > wrfsfc_d03_20210301_0000_f00000-cfcomp.nc
> > > wrfsfc_d03_20210301_0000_f00000-cfcomp.nc" ;
> > > :NCO = "netCDF Operators version 4.9.8 (Homepage
=
> > >
> >
https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fnco.s
> >
f.net%2F&data=04%7C01%7CAJEB%40pge.com%7C0b860e74eccf45e2508508d92
> >
ad54e1c%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C63758791790377891
> >
1%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6
> >
Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=siex3pLGh5m0T3wXYmbV%2BI4u42Jl
> > QcyAfb07ZrRgN2c%3D&reserved=0,
> > Code =
> >
>
https://nam10.safelinks.protection.outlook.com/?url=http%3A%2F%2Fgithub.com%2Fnco%2Fnco&data=04%7C01%7CAJEB%40pge.com%7C0b860e74eccf45e2508508d92ad54e1c%7C44ae661aece641aabc967c2c85a08941%7C0%7C0%7C637587917903788907%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=Wob2rFlFVUI%2F9mJWPz7G2nq%2BcNcEoAUL642EH54JoVA%3D&reserved=0
> )"
> > ;
> > >
> > >
> > >
> >
> >
> >
> >
>
>
>
------------------------------------------------
More information about the Met_help
mailing list