[Met_help] [rt.rap.ucar.edu #93267] History for Fcst time in NetCDF Forecast File
John Halley Gotway via RT
met_help at ucar.edu
Mon Dec 2 14:07:24 MST 2019
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
John, hope the snow isn't too bad. Question on NETCDF Forecast files. I am verifying NETCDF precip forecast files. The header for the files is listed below:
netcdf PS.557WW_SC.U_DI.C_GP.GALWEM-GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.012_DF {
dimensions:
latitude = 721 ;
longitude = 1440 ;
time = 1 ;
variables:
float latitude(latitude) ;
latitude:units = "degrees_north" ;
latitude:long_name = "latitude" ;
float longitude(longitude) ;
longitude:units = "degrees_east" ;
longitude:long_name = "longitude" ;
int time(time) ;
time:units = "hours since 1900-01-01 00:00:0.0" ;
time:long_name = "time" ;
time:calendar = "gregorian" ;
double tprec(time, latitude, longitude) ;
tprec:_FillValue = -32767. ;
tprec:missing_value = -32767s ;
tprec:long_name = "Total precipitation" ;
// global attributes:
:Conventions = "CF-1.0" ;
:history = "Wed Nov 27 00:19:41 2019: /h/ENSEMBLE/bin/ncbo --op_typ=- PS.557WW_SC.U_DI.C_GP.GALWEM-GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.012_DF.nc PS.557WW_SC.U
_DI.C_GP.GALWEM-GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.006_DF.nc ../6hr_netcdf/PS.557WW_SC.U_DI.C_GP.GALWEM-GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.012_DF.nc\n",
"2019-11-26 21:46:00 GMT by grib_to_netcdf-1.13.0: /h/3P/bin/grib_to_netcdf PS.557WW_SC.U_DI.C_GP.GALWEM-GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.012_DF.GR
2 -o netcdf/PS.557WW_SC.U_DI.C_GP.GALWEM-GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.012_DF.GR2" ;
:NCO = "\"4.6.1\"" ;
data:
latitude = -90, -89.75, -89.5, -89.25, -89, -88.75, -88.5, -88.25, -88,
-87.75, -87.5, -87.25, -87, -86.75, -86.5, -86.25, -86, -85.75, -85.5,
When I verify these files met thinks the forecast file is a 0hr forecast. I get the following:
DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file: /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file: /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
DEBUG 1: Default Config File: /h/MET/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: /h/data/global/WXQC/data/met/config/met_config/PointStatConfig_TP_NC_controlrun_updated
WARNING:
WARNING: InterpInfo::validate() -> Resetting interpolation width 0 from 2 to 1 for interpolation method "NEAREST".
WARNING:
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=18446744072342855261
DEBUG 1: Forecast File: /home/franksc/cs_proj/April2019_controlRun/6hr_netcdf/PS.557WW_SC.U_DI.C_GP.GALWEM-GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.012_DF.nc
DEBUG 1: Observation File: /h/data/global/WXQC/data/met/tmp/c72cfd4c-5c6c-4ffc-9496-6976b7baeafc
WARNING:
WARNING: MetNcCFDataFile::data_plane(VarInfo &, DataPlane &) -> for "tprec" variable, the lead time does not match the requested lead time: (000000 != 120000)
WARNING:
DEBUG 3: Applying conversion function.
What should the time field look like for this to work? I will send the forecast NETCDF by AMRDEC.
Thanks
Bob Craig
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Fcst time in NetCDF Forecast File
From: John Halley Gotway
Time: Wed Nov 27 11:21:38 2019
Bob,
In CF-compliant NetCDF files, the forecast reference time variable
defines
the model initialization time. MET reads that and the time variable
and
computes the lead time as their difference.
So adding the forecast reference time variable should do the trick.
I’m out of the office the rest of this week, but you can look at the
CF
documentation and search for “forecast_reference_time”.
http://cfconventions.org/cf-conventions/cf-conventions.html
Thanks
John
On Wed, Nov 27, 2019 at 7:58 AM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> Wed Nov 27 07:58:11 2019: Request 93267 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
> Queue: met_help
> Subject: Fcst time in NetCDF Forecast File
> Owner: Nobody
> Requestors: robert.craig.2 at us.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=93267 >
>
>
> John, hope the snow isn't too bad. Question on NETCDF Forecast
files. I
> am verifying NETCDF precip forecast files. The header for the files
is
> listed below:
>
> netcdf
> PS.557WW_SC.U_DI.C_GP.GALWEM-
GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.012_DF
> {
> dimensions:
> latitude = 721 ;
> longitude = 1440 ;
> time = 1 ;
> variables:
> float latitude(latitude) ;
> latitude:units = "degrees_north" ;
> latitude:long_name = "latitude" ;
> float longitude(longitude) ;
> longitude:units = "degrees_east" ;
> longitude:long_name = "longitude" ;
> int time(time) ;
> time:units = "hours since 1900-01-01 00:00:0.0" ;
> time:long_name = "time" ;
> time:calendar = "gregorian" ;
> double tprec(time, latitude, longitude) ;
> tprec:_FillValue = -32767. ;
> tprec:missing_value = -32767s ;
> tprec:long_name = "Total precipitation" ;
>
> // global attributes:
> :Conventions = "CF-1.0" ;
> :history = "Wed Nov 27 00:19:41 2019:
/h/ENSEMBLE/bin/ncbo
> --op_typ=-
> PS.557WW_SC.U_DI.C_GP.GALWEM-
GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.012_DF.nc
> PS.557WW_SC.U
> _DI.C_GP.GALWEM-
GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.006_DF.nc
> ../6hr_netcdf/
> PS.557WW_SC.U_DI.C_GP.GALWEM-
GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.012_DF.nc
> \n",
> "2019-11-26 21:46:00 GMT by grib_to_netcdf-
1.13.0:
> /h/3P/bin/grib_to_netcdf
> PS.557WW_SC.U_DI.C_GP.GALWEM-
GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.012_DF.GR
> 2 -o
> netcdf/PS.557WW_SC.U_DI.C_GP.GALWEM-
GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.012_DF.GR2"
> ;
> :NCO = "\"4.6.1\"" ;
> data:
>
> latitude = -90, -89.75, -89.5, -89.25, -89, -88.75, -88.5, -88.25,
-88,
> -87.75, -87.5, -87.25, -87, -86.75, -86.5, -86.25, -86, -85.75,
-85.5,
>
>
> When I verify these files met thinks the forecast file is a 0hr
forecast.
> I get the following:
>
> DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file:
> /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
> DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file:
> /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
> DEBUG 1: Default Config File:
> /h/MET/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File:
>
/h/data/global/WXQC/data/met/config/met_config/PointStatConfig_TP_NC_controlrun_updated
> WARNING:
> WARNING: InterpInfo::validate() -> Resetting interpolation width 0
from 2
> to 1 for interpolation method "NEAREST".
> WARNING:
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=18446744072342855261
> DEBUG 1: Forecast File:
> /home/franksc/cs_proj/April2019_controlRun/6hr_netcdf/
> PS.557WW_SC.U_DI.C_GP.GALWEM-
GD_GR.C0P25DEG_AR.GLOBAL_DD.20190401_CY.00_FH.012_DF.nc
> DEBUG 1: Observation File:
> /h/data/global/WXQC/data/met/tmp/c72cfd4c-5c6c-4ffc-9496-
6976b7baeafc
> WARNING:
> WARNING: MetNcCFDataFile::data_plane(VarInfo &, DataPlane &) -> for
> "tprec" variable, the lead time does not match the requested lead
time:
> (000000 != 120000)
> WARNING:
> DEBUG 3: Applying conversion function.
>
> What should the time field look like for this to work? I will send
the
> forecast NETCDF by AMRDEC.
>
> Thanks
> Bob Craig
>
>
------------------------------------------------
Subject: Fcst time in NetCDF Forecast File
From: John Halley Gotway
Time: Mon Dec 02 11:00:32 2019
Bob,
Just wanted to follow up on this ticket. I did pull the file you sent
via
ARL SAFE but haven't done anything with it. Just let me know if you
have
issues or questions.
Thanks,
John
On Wed, Nov 27, 2019 at 11:21 AM The RT System itself via RT <
met_help at ucar.edu> wrote:
>
> Wed Nov 27 11:21:38 2019: Request 93267 was acted upon.
> Transaction: Given to johnhg (John Halley Gotway) by RT_System
> Queue: met_help
> Subject: Fcst time in NetCDF Forecast File
> Owner: johnhg
> Requestors: robert.craig.2 at us.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=93267 >
>
>
> This transaction appears to have no content
>
------------------------------------------------
Subject: RE: [Non-DoD Source] Re: [rt.rap.ucar.edu #93267] Fcst time in NetCDF Forecast File
From: robert.craig.2 at us.af.mil
Time: Mon Dec 02 11:28:30 2019
Thanks, you answered the question already.
-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Monday, December 2, 2019 12:01 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Cc: WICKARD, BRAEDI J CTR USAF AFMC AFLCMC/HBAW-OL
<braedi.wickard.ctr at us.af.mil>; FRANKS, CHRISTOPHER C CTR USAF AFMC
AFLCMC/HBAW-OL <christopher.franks.2.ctr at us.af.mil>
Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #93267] Fcst time in
NetCDF Forecast File
Bob,
Just wanted to follow up on this ticket. I did pull the file you sent
via ARL SAFE but haven't done anything with it. Just let me know if
you have issues or questions.
Thanks,
John
On Wed, Nov 27, 2019 at 11:21 AM The RT System itself via RT <
met_help at ucar.edu> wrote:
>
> Wed Nov 27 11:21:38 2019: Request 93267 was acted upon.
> Transaction: Given to johnhg (John Halley Gotway) by RT_System
> Queue: met_help
> Subject: Fcst time in NetCDF Forecast File
> Owner: johnhg
> Requestors: robert.craig.2 at us.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=93267
> >
>
>
> This transaction appears to have no content
>
------------------------------------------------
More information about the Met_help
mailing list