[Met_help] [rt.rap.ucar.edu #62455] History for Problem with TC_PAIRS
John Halley Gotway via RT
met_help at ucar.edu
Tue Aug 27 10:49:53 MDT 2013
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Here we go again with a new issue. I'm testing the tropical cyclone tracker for an upcoming project we've been asked to do. So I gathered some sample files from NHC's FTP site and tried to run the code. I got the TC_DLAND to run just fine, but something appears to prevent TC_PAIRS from running, and I can't figure out what it is.
Here's the log file as I see it. My main concerns are why Model Suffix is (nul), and why the ADECK file claims to be read but none of the records are used. I am using the wwpts file included in 4.1; I can't find a newer one on the NHC FTP site-I looked in many directories there but didn't see this file. The storm I'm looking at has no entries in the wwpts file, presumably because the wwpts file predates the aal and bal files.
DEBUG 1: [Source 1 of 1] ADECK Source: /home/qcteam/METv4.1/data/aal142012.dat, Model Suffix: (nul)
DEBUG 1: [Source 1 of 1] BDECK Source: /home/qcteam/METv4.1/data/bal142012.dat, Model Suffix: (nul)
DEBUG 1: Config File Default: /home/qcteam/METv4.1/data/config/TCPairsConfig_default
DEBUG 1: Config File User: /home/qcteam/METv4.1/data/config/TCPairs.cfg
DEBUG 1: Distance to land file: /home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc
DEBUG 4: Latitude/Longitude Grid Data:
DEBUG 4: lat_ll = 0
DEBUG 4: lon_ll = -180
DEBUG 4: delta_lat = 0.1
DEBUG 4: delta_lon = 0.1
DEBUG 4: Nlat = 601
DEBUG 4: Nlon = 1801
DEBUG 2: Processing 1 ADECK file(s).
DEBUG 4: [File 1 of 1] Used 0 of 202726 lines read from file "/home/qcteam/METv4.1/data/aal142012.dat"
DEBUG 3: Used 0 of 202726 lines read from 1 file(s).
DEBUG 3: Identified 0 track(s).
DEBUG 2: Found 0 ADECK track(s).
DEBUG 2: Filtering 0 ADECK tracks based on config file settings.
DEBUG 3: Total tracks read = 0
DEBUG 3: Total tracks kept = 0
DEBUG 3: Rejected for model = 0
DEBUG 3: Rejected for storm id = 0
DEBUG 3: Rejected for basin = 0
DEBUG 3: Rejected for cyclone = 0
DEBUG 3: Rejected for storm name = 0
DEBUG 3: Rejected for init time = 0
DEBUG 3: Rejected for init hour = 0
DEBUG 3: Rejected for valid time = 0
DEBUG 3: Rejected for init mask = 0
DEBUG 3: Rejected for valid mask = 0
DEBUG 2: Processing 1 BDECK file(s).
DEBUG 4: Updating storm name from "INVEST" to "FOURTEEN" for AL142012.
DEBUG 4: Updating storm name from "FOURTEEN" to "NADINE" for AL142012.
DEBUG 4: [File 1 of 1] Used 184 of 184 lines read from file "/home/qcteam/METv4.1/data/bal142012.dat"
DEBUG 3: Used 184 of 184 lines read from 1 file(s).
DEBUG 3: Identified 1 track(s).
DEBUG 4: [Track 1 of 1] TrackInfo: StormId = "AL142012", Basin = "AL", Cyclone = "14", StormName = "NADINE", TechniqueNumber = -9999, Technique = "BEST", Initials = "(nul)", InitTime = NA, MinValidTime = 20120910_120000, MaxValidTime = 20121004_060000, NPoints = 96, NAlloc = 100, NTrackLines = 0
DEBUG 2: Found 1 BDECK track(s).
DEBUG 2: Deriving 12-hour interpolated tracks.
DEBUG 2: Finished deriving 0 12-hour interpolated tracks.
DEBUG 2: Deriving 0 ADECK consensus model(s).
DEBUG 2: Added 0 ADECK consensus tracks(s).
DEBUG 2: Deriving 0 ADECK lag model(s).
DEBUG 2: Added 0 ADECK lag tracks(s).
DEBUG 2: Deriving 0 CLIPER/SHIFOR baseline model(s).
DEBUG 2: Added 0 CLIPER/SHIFOR baseline track(s).
DEBUG 2: Matching 0 ADECK tracks to 1 BDECK tracks.
DEBUG 1: Watch/Warning file: /home/qcteam/METv4.1/data/tc_data/wwpts_us.txt
DEBUG 1: Output file: ./out_tcmpr.tcst
My config file contents:
model = ["AEMN"];
storm_id = ["AL142012"];
basin = ["AL"];
cyclone = ["14"];
storm_name = ["NADINE"];
init_beg = "";
init_end = "";
init_inc = [];
init_exc = [];
valid_beg = "";
valid_end = "";
init_hour = [];
init_mask = "";
valid_mask = "";
check_dup = FALSE;
interp12 = TRUE;
consensus = [];
lag_time = [];
best_baseline = [];
oper_baseline = [];
match_points = FALSE;
dland_file = "/home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc";
watch_warn = {
file_name = "/home/qcteam/METv4.1/data/tc_data/wwpts_us.txt";
time_offset = -14400;
};
version = "V4.1";
Hopefully you can suggest something to get this working!
Thanks,
Matt
// SIGNED //
Matthew C. Sittel
Associate Scientist
University Corporation for Atmospheric Research
16WS/WXN, Offutt AFB, NE
(402) 294-3473 DSN: 271-3473
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
From: John Halley Gotway
Time: Thu Aug 01 09:37:10 2013
Hello Matt,
Try removing the storm name from the config file, like this:
storm_name = [];
When you rerun, you should get matching tracks.
However, I think there may be a bug in the config file parsing code,
and I'll let you know what I find.
Thanks,
John
On 07/31/2013 03:40 PM, matthew.sittel.ctr at offutt.af.mil via RT wrote:
>
> Wed Jul 31 15:40:07 2013: Request 62455 was acted upon.
> Transaction: Ticket created by matthew.sittel.ctr at offutt.af.mil
> Queue: met_help
> Subject: Problem with TC_PAIRS
> Owner: Nobody
> Requestors: matthew.sittel.ctr at offutt.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>
>
> Here we go again with a new issue. I'm testing the tropical cyclone
tracker for an upcoming project we've been asked to do. So I gathered
some sample files from NHC's FTP site and tried to run the code. I
got the TC_DLAND to run just fine, but something appears to prevent
TC_PAIRS from running, and I can't figure out what it is.
>
> Here's the log file as I see it. My main concerns are why Model
Suffix is (nul), and why the ADECK file claims to be read but none of
the records are used. I am using the wwpts file included in 4.1; I
can't find a newer one on the NHC FTP site-I looked in many
directories there but didn't see this file. The storm I'm looking at
has no entries in the wwpts file, presumably because the wwpts file
predates the aal and bal files.
>
> DEBUG 1: [Source 1 of 1] ADECK Source:
/home/qcteam/METv4.1/data/aal142012.dat, Model Suffix: (nul)
> DEBUG 1: [Source 1 of 1] BDECK Source:
/home/qcteam/METv4.1/data/bal142012.dat, Model Suffix: (nul)
> DEBUG 1: Config File Default:
/home/qcteam/METv4.1/data/config/TCPairsConfig_default
> DEBUG 1: Config File User:
/home/qcteam/METv4.1/data/config/TCPairs.cfg
> DEBUG 1: Distance to land file:
/home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc
> DEBUG 4: Latitude/Longitude Grid Data:
> DEBUG 4: lat_ll = 0
> DEBUG 4: lon_ll = -180
> DEBUG 4: delta_lat = 0.1
> DEBUG 4: delta_lon = 0.1
> DEBUG 4: Nlat = 601
> DEBUG 4: Nlon = 1801
> DEBUG 2: Processing 1 ADECK file(s).
> DEBUG 4: [File 1 of 1] Used 0 of 202726 lines read from file
"/home/qcteam/METv4.1/data/aal142012.dat"
> DEBUG 3: Used 0 of 202726 lines read from 1 file(s).
> DEBUG 3: Identified 0 track(s).
> DEBUG 2: Found 0 ADECK track(s).
> DEBUG 2: Filtering 0 ADECK tracks based on config file settings.
> DEBUG 3: Total tracks read = 0
> DEBUG 3: Total tracks kept = 0
> DEBUG 3: Rejected for model = 0
> DEBUG 3: Rejected for storm id = 0
> DEBUG 3: Rejected for basin = 0
> DEBUG 3: Rejected for cyclone = 0
> DEBUG 3: Rejected for storm name = 0
> DEBUG 3: Rejected for init time = 0
> DEBUG 3: Rejected for init hour = 0
> DEBUG 3: Rejected for valid time = 0
> DEBUG 3: Rejected for init mask = 0
> DEBUG 3: Rejected for valid mask = 0
> DEBUG 2: Processing 1 BDECK file(s).
> DEBUG 4: Updating storm name from "INVEST" to "FOURTEEN" for
AL142012.
> DEBUG 4: Updating storm name from "FOURTEEN" to "NADINE" for
AL142012.
> DEBUG 4: [File 1 of 1] Used 184 of 184 lines read from file
"/home/qcteam/METv4.1/data/bal142012.dat"
> DEBUG 3: Used 184 of 184 lines read from 1 file(s).
> DEBUG 3: Identified 1 track(s).
> DEBUG 4: [Track 1 of 1] TrackInfo: StormId = "AL142012", Basin =
"AL", Cyclone = "14", StormName = "NADINE", TechniqueNumber = -9999,
Technique = "BEST", Initials = "(nul)", InitTime = NA, MinValidTime =
20120910_120000, MaxValidTime = 20121004_060000, NPoints = 96, NAlloc
= 100, NTrackLines = 0
> DEBUG 2: Found 1 BDECK track(s).
> DEBUG 2: Deriving 12-hour interpolated tracks.
> DEBUG 2: Finished deriving 0 12-hour interpolated tracks.
> DEBUG 2: Deriving 0 ADECK consensus model(s).
> DEBUG 2: Added 0 ADECK consensus tracks(s).
> DEBUG 2: Deriving 0 ADECK lag model(s).
> DEBUG 2: Added 0 ADECK lag tracks(s).
> DEBUG 2: Deriving 0 CLIPER/SHIFOR baseline model(s).
> DEBUG 2: Added 0 CLIPER/SHIFOR baseline track(s).
> DEBUG 2: Matching 0 ADECK tracks to 1 BDECK tracks.
> DEBUG 1: Watch/Warning file:
/home/qcteam/METv4.1/data/tc_data/wwpts_us.txt
> DEBUG 1: Output file: ./out_tcmpr.tcst
>
> My config file contents:
>
> model = ["AEMN"];
>
> storm_id = ["AL142012"];
>
> basin = ["AL"];
>
> cyclone = ["14"];
>
> storm_name = ["NADINE"];
>
> init_beg = "";
> init_end = "";
> init_inc = [];
> init_exc = [];
>
> valid_beg = "";
> valid_end = "";
>
> init_hour = [];
>
> init_mask = "";
> valid_mask = "";
>
> check_dup = FALSE;
>
> interp12 = TRUE;
>
> consensus = [];
>
> lag_time = [];
>
> best_baseline = [];
> oper_baseline = [];
>
> match_points = FALSE;
>
> dland_file =
"/home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc";
>
> watch_warn = {
> file_name = "/home/qcteam/METv4.1/data/tc_data/wwpts_us.txt";
> time_offset = -14400;
> };
>
> version = "V4.1";
>
> Hopefully you can suggest something to get this working!
>
> Thanks,
> Matt
>
> // SIGNED //
> Matthew C. Sittel
> Associate Scientist
> University Corporation for Atmospheric Research
> 16WS/WXN, Offutt AFB, NE
> (402) 294-3473 DSN: 271-3473
>
------------------------------------------------
Subject: Problem with TC_PAIRS
From: matthew.sittel.ctr at offutt.af.mil
Time: Thu Aug 01 09:43:08 2013
Hey John,
Yeah, that change took care of it.
Our ensembles guys are going to make ATCF-formatted tracks for their
WRF members, the ensemble mean, and then compare them to a
deterministic WRF model... all to be compared to best tracks for
storms. TS Gil is in the east Pacific so they've got something to
work with. This will be our chance to utilize MET (which is in my
statement of work, to be the "subject matter expert" in house).
Thanks,
Matt
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Thursday, August 01, 2013 10:37 AM
To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
Hello Matt,
Try removing the storm name from the config file, like this:
storm_name = [];
When you rerun, you should get matching tracks.
However, I think there may be a bug in the config file parsing code,
and I'll let you know what I find.
Thanks,
John
On 07/31/2013 03:40 PM, matthew.sittel.ctr at offutt.af.mil via RT wrote:
>
> Wed Jul 31 15:40:07 2013: Request 62455 was acted upon.
> Transaction: Ticket created by matthew.sittel.ctr at offutt.af.mil
> Queue: met_help
> Subject: Problem with TC_PAIRS
> Owner: Nobody
> Requestors: matthew.sittel.ctr at offutt.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455
> >
>
>
> Here we go again with a new issue. I'm testing the tropical cyclone
tracker for an upcoming project we've been asked to do. So I gathered
some sample files from NHC's FTP site and tried to run the code. I
got the TC_DLAND to run just fine, but something appears to prevent
TC_PAIRS from running, and I can't figure out what it is.
>
> Here's the log file as I see it. My main concerns are why Model
Suffix is (nul), and why the ADECK file claims to be read but none of
the records are used. I am using the wwpts file included in 4.1; I
can't find a newer one on the NHC FTP site-I looked in many
directories there but didn't see this file. The storm I'm looking at
has no entries in the wwpts file, presumably because the wwpts file
predates the aal and bal files.
>
> DEBUG 1: [Source 1 of 1] ADECK Source:
> /home/qcteam/METv4.1/data/aal142012.dat, Model Suffix: (nul) DEBUG
1:
> [Source 1 of 1] BDECK Source:
/home/qcteam/METv4.1/data/bal142012.dat,
> Model Suffix: (nul) DEBUG 1: Config File Default:
> /home/qcteam/METv4.1/data/config/TCPairsConfig_default
> DEBUG 1: Config File User:
> /home/qcteam/METv4.1/data/config/TCPairs.cfg
> DEBUG 1: Distance to land file:
> /home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc
> DEBUG 4: Latitude/Longitude Grid Data:
> DEBUG 4: lat_ll = 0
> DEBUG 4: lon_ll = -180
> DEBUG 4: delta_lat = 0.1
> DEBUG 4: delta_lon = 0.1
> DEBUG 4: Nlat = 601
> DEBUG 4: Nlon = 1801
> DEBUG 2: Processing 1 ADECK file(s).
> DEBUG 4: [File 1 of 1] Used 0 of 202726 lines read from file
"/home/qcteam/METv4.1/data/aal142012.dat"
> DEBUG 3: Used 0 of 202726 lines read from 1 file(s).
> DEBUG 3: Identified 0 track(s).
> DEBUG 2: Found 0 ADECK track(s).
> DEBUG 2: Filtering 0 ADECK tracks based on config file settings.
> DEBUG 3: Total tracks read = 0
> DEBUG 3: Total tracks kept = 0
> DEBUG 3: Rejected for model = 0
> DEBUG 3: Rejected for storm id = 0
> DEBUG 3: Rejected for basin = 0
> DEBUG 3: Rejected for cyclone = 0
> DEBUG 3: Rejected for storm name = 0
> DEBUG 3: Rejected for init time = 0
> DEBUG 3: Rejected for init hour = 0
> DEBUG 3: Rejected for valid time = 0
> DEBUG 3: Rejected for init mask = 0
> DEBUG 3: Rejected for valid mask = 0
> DEBUG 2: Processing 1 BDECK file(s).
> DEBUG 4: Updating storm name from "INVEST" to "FOURTEEN" for
AL142012.
> DEBUG 4: Updating storm name from "FOURTEEN" to "NADINE" for
AL142012.
> DEBUG 4: [File 1 of 1] Used 184 of 184 lines read from file
"/home/qcteam/METv4.1/data/bal142012.dat"
> DEBUG 3: Used 184 of 184 lines read from 1 file(s).
> DEBUG 3: Identified 1 track(s).
> DEBUG 4: [Track 1 of 1] TrackInfo: StormId = "AL142012", Basin =
"AL",
> Cyclone = "14", StormName = "NADINE", TechniqueNumber = -9999,
Technique = "BEST", Initials = "(nul)", InitTime = NA, MinValidTime =
20120910_120000, MaxValidTime = 20121004_060000, NPoints = 96, NAlloc
= 100, NTrackLines = 0 DEBUG 2: Found 1 BDECK track(s).
> DEBUG 2: Deriving 12-hour interpolated tracks.
> DEBUG 2: Finished deriving 0 12-hour interpolated tracks.
> DEBUG 2: Deriving 0 ADECK consensus model(s).
> DEBUG 2: Added 0 ADECK consensus tracks(s).
> DEBUG 2: Deriving 0 ADECK lag model(s).
> DEBUG 2: Added 0 ADECK lag tracks(s).
> DEBUG 2: Deriving 0 CLIPER/SHIFOR baseline model(s).
> DEBUG 2: Added 0 CLIPER/SHIFOR baseline track(s).
> DEBUG 2: Matching 0 ADECK tracks to 1 BDECK tracks.
> DEBUG 1: Watch/Warning file:
> /home/qcteam/METv4.1/data/tc_data/wwpts_us.txt
> DEBUG 1: Output file: ./out_tcmpr.tcst
>
> My config file contents:
>
> model = ["AEMN"];
>
> storm_id = ["AL142012"];
>
> basin = ["AL"];
>
> cyclone = ["14"];
>
> storm_name = ["NADINE"];
>
> init_beg = "";
> init_end = "";
> init_inc = [];
> init_exc = [];
>
> valid_beg = "";
> valid_end = "";
>
> init_hour = [];
>
> init_mask = "";
> valid_mask = "";
>
> check_dup = FALSE;
>
> interp12 = TRUE;
>
> consensus = [];
>
> lag_time = [];
>
> best_baseline = [];
> oper_baseline = [];
>
> match_points = FALSE;
>
> dland_file =
> "/home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc";
>
> watch_warn = {
> file_name = "/home/qcteam/METv4.1/data/tc_data/wwpts_us.txt";
> time_offset = -14400;
> };
>
> version = "V4.1";
>
> Hopefully you can suggest something to get this working!
>
> Thanks,
> Matt
>
> // SIGNED //
> Matthew C. Sittel
> Associate Scientist
> University Corporation for Atmospheric Research 16WS/WXN, Offutt
AFB,
> NE
> (402) 294-3473 DSN: 271-3473
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
From: John Halley Gotway
Time: Thu Aug 01 09:52:59 2013
Matt,
Sounds interesting. I've actually switched gears recently too. I'm
working on the real-time evaluation for the HFIP products. Currently,
we're mostly just plotting track data, but the plan is to
start running MET-TC in near-realtime mode soon.
I did find a copy-and-paste bug in the tc_pairs tool. I'll post a
bugfix for it shortly and write you back. You should still continue
to omit the storm name though. While the ATCF file format has a
place for the storm name, it usually is not included in the ADECK
files. If you list a storm name in the config file, it'll only keep
tracks that include it. Since ADECK files often don't have it,
you won't keep many of them.
Also, don't worry about the model name suffix being "(nul)". We put
in this option of appending a string to the model names found in the
track files. At one point, we had multiple models to evaluate
that all made use of the same ATCF ID. To distinguish them in the
output, we added a suffix to their ATCF ID's to make them unique. If
you're not in that situation, you don't need to use the suffix.
Thanks,
John
On 08/01/2013 09:43 AM, matthew.sittel.ctr at offutt.af.mil via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>
> Hey John,
>
> Yeah, that change took care of it.
>
> Our ensembles guys are going to make ATCF-formatted tracks for their
WRF members, the ensemble mean, and then compare them to a
deterministic WRF model... all to be compared to best tracks for
storms. TS Gil is in the east Pacific so they've got something to
work with. This will be our chance to utilize MET (which is in my
statement of work, to be the "subject matter expert" in house).
>
> Thanks,
> Matt
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, August 01, 2013 10:37 AM
> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>
> Hello Matt,
>
> Try removing the storm name from the config file, like this:
> storm_name = [];
>
> When you rerun, you should get matching tracks.
>
> However, I think there may be a bug in the config file parsing code,
and I'll let you know what I find.
>
> Thanks,
> John
>
>
> On 07/31/2013 03:40 PM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>
>> Wed Jul 31 15:40:07 2013: Request 62455 was acted upon.
>> Transaction: Ticket created by matthew.sittel.ctr at offutt.af.mil
>> Queue: met_help
>> Subject: Problem with TC_PAIRS
>> Owner: Nobody
>> Requestors: matthew.sittel.ctr at offutt.af.mil
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455
>>>
>>
>>
>> Here we go again with a new issue. I'm testing the tropical
cyclone tracker for an upcoming project we've been asked to do. So I
gathered some sample files from NHC's FTP site and tried to run the
code. I got the TC_DLAND to run just fine, but something appears to
prevent TC_PAIRS from running, and I can't figure out what it is.
>>
>> Here's the log file as I see it. My main concerns are why Model
Suffix is (nul), and why the ADECK file claims to be read but none of
the records are used. I am using the wwpts file included in 4.1; I
can't find a newer one on the NHC FTP site-I looked in many
directories there but didn't see this file. The storm I'm looking at
has no entries in the wwpts file, presumably because the wwpts file
predates the aal and bal files.
>>
>> DEBUG 1: [Source 1 of 1] ADECK Source:
>> /home/qcteam/METv4.1/data/aal142012.dat, Model Suffix: (nul) DEBUG
1:
>> [Source 1 of 1] BDECK Source:
/home/qcteam/METv4.1/data/bal142012.dat,
>> Model Suffix: (nul) DEBUG 1: Config File Default:
>> /home/qcteam/METv4.1/data/config/TCPairsConfig_default
>> DEBUG 1: Config File User:
>> /home/qcteam/METv4.1/data/config/TCPairs.cfg
>> DEBUG 1: Distance to land file:
>> /home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc
>> DEBUG 4: Latitude/Longitude Grid Data:
>> DEBUG 4: lat_ll = 0
>> DEBUG 4: lon_ll = -180
>> DEBUG 4: delta_lat = 0.1
>> DEBUG 4: delta_lon = 0.1
>> DEBUG 4: Nlat = 601
>> DEBUG 4: Nlon = 1801
>> DEBUG 2: Processing 1 ADECK file(s).
>> DEBUG 4: [File 1 of 1] Used 0 of 202726 lines read from file
"/home/qcteam/METv4.1/data/aal142012.dat"
>> DEBUG 3: Used 0 of 202726 lines read from 1 file(s).
>> DEBUG 3: Identified 0 track(s).
>> DEBUG 2: Found 0 ADECK track(s).
>> DEBUG 2: Filtering 0 ADECK tracks based on config file settings.
>> DEBUG 3: Total tracks read = 0
>> DEBUG 3: Total tracks kept = 0
>> DEBUG 3: Rejected for model = 0
>> DEBUG 3: Rejected for storm id = 0
>> DEBUG 3: Rejected for basin = 0
>> DEBUG 3: Rejected for cyclone = 0
>> DEBUG 3: Rejected for storm name = 0
>> DEBUG 3: Rejected for init time = 0
>> DEBUG 3: Rejected for init hour = 0
>> DEBUG 3: Rejected for valid time = 0
>> DEBUG 3: Rejected for init mask = 0
>> DEBUG 3: Rejected for valid mask = 0
>> DEBUG 2: Processing 1 BDECK file(s).
>> DEBUG 4: Updating storm name from "INVEST" to "FOURTEEN" for
AL142012.
>> DEBUG 4: Updating storm name from "FOURTEEN" to "NADINE" for
AL142012.
>> DEBUG 4: [File 1 of 1] Used 184 of 184 lines read from file
"/home/qcteam/METv4.1/data/bal142012.dat"
>> DEBUG 3: Used 184 of 184 lines read from 1 file(s).
>> DEBUG 3: Identified 1 track(s).
>> DEBUG 4: [Track 1 of 1] TrackInfo: StormId = "AL142012", Basin =
"AL",
>> Cyclone = "14", StormName = "NADINE", TechniqueNumber = -9999,
Technique = "BEST", Initials = "(nul)", InitTime = NA, MinValidTime =
20120910_120000, MaxValidTime = 20121004_060000, NPoints = 96, NAlloc
= 100, NTrackLines = 0 DEBUG 2: Found 1 BDECK track(s).
>> DEBUG 2: Deriving 12-hour interpolated tracks.
>> DEBUG 2: Finished deriving 0 12-hour interpolated tracks.
>> DEBUG 2: Deriving 0 ADECK consensus model(s).
>> DEBUG 2: Added 0 ADECK consensus tracks(s).
>> DEBUG 2: Deriving 0 ADECK lag model(s).
>> DEBUG 2: Added 0 ADECK lag tracks(s).
>> DEBUG 2: Deriving 0 CLIPER/SHIFOR baseline model(s).
>> DEBUG 2: Added 0 CLIPER/SHIFOR baseline track(s).
>> DEBUG 2: Matching 0 ADECK tracks to 1 BDECK tracks.
>> DEBUG 1: Watch/Warning file:
>> /home/qcteam/METv4.1/data/tc_data/wwpts_us.txt
>> DEBUG 1: Output file: ./out_tcmpr.tcst
>>
>> My config file contents:
>>
>> model = ["AEMN"];
>>
>> storm_id = ["AL142012"];
>>
>> basin = ["AL"];
>>
>> cyclone = ["14"];
>>
>> storm_name = ["NADINE"];
>>
>> init_beg = "";
>> init_end = "";
>> init_inc = [];
>> init_exc = [];
>>
>> valid_beg = "";
>> valid_end = "";
>>
>> init_hour = [];
>>
>> init_mask = "";
>> valid_mask = "";
>>
>> check_dup = FALSE;
>>
>> interp12 = TRUE;
>>
>> consensus = [];
>>
>> lag_time = [];
>>
>> best_baseline = [];
>> oper_baseline = [];
>>
>> match_points = FALSE;
>>
>> dland_file =
>> "/home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc";
>>
>> watch_warn = {
>> file_name =
"/home/qcteam/METv4.1/data/tc_data/wwpts_us.txt";
>> time_offset = -14400;
>> };
>>
>> version = "V4.1";
>>
>> Hopefully you can suggest something to get this working!
>>
>> Thanks,
>> Matt
>>
>> // SIGNED //
>> Matthew C. Sittel
>> Associate Scientist
>> University Corporation for Atmospheric Research 16WS/WXN, Offutt
AFB,
>> NE
>> (402) 294-3473 DSN: 271-3473
>>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
From: John Halley Gotway
Time: Thu Aug 01 10:06:35 2013
Matt,
OK, I've posted the bugfix for the cut-and-paste error I found:
http://www.dtcenter.org/met/users/support/known_issues/METv4.1/index.php
Just let us know if more questions or problems arise.
Thanks,
John
On 08/01/2013 09:43 AM, matthew.sittel.ctr at offutt.af.mil via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>
> Hey John,
>
> Yeah, that change took care of it.
>
> Our ensembles guys are going to make ATCF-formatted tracks for their
WRF members, the ensemble mean, and then compare them to a
deterministic WRF model... all to be compared to best tracks for
storms. TS Gil is in the east Pacific so they've got something to
work with. This will be our chance to utilize MET (which is in my
statement of work, to be the "subject matter expert" in house).
>
> Thanks,
> Matt
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, August 01, 2013 10:37 AM
> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>
> Hello Matt,
>
> Try removing the storm name from the config file, like this:
> storm_name = [];
>
> When you rerun, you should get matching tracks.
>
> However, I think there may be a bug in the config file parsing code,
and I'll let you know what I find.
>
> Thanks,
> John
>
>
> On 07/31/2013 03:40 PM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>
>> Wed Jul 31 15:40:07 2013: Request 62455 was acted upon.
>> Transaction: Ticket created by matthew.sittel.ctr at offutt.af.mil
>> Queue: met_help
>> Subject: Problem with TC_PAIRS
>> Owner: Nobody
>> Requestors: matthew.sittel.ctr at offutt.af.mil
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455
>>>
>>
>>
>> Here we go again with a new issue. I'm testing the tropical
cyclone tracker for an upcoming project we've been asked to do. So I
gathered some sample files from NHC's FTP site and tried to run the
code. I got the TC_DLAND to run just fine, but something appears to
prevent TC_PAIRS from running, and I can't figure out what it is.
>>
>> Here's the log file as I see it. My main concerns are why Model
Suffix is (nul), and why the ADECK file claims to be read but none of
the records are used. I am using the wwpts file included in 4.1; I
can't find a newer one on the NHC FTP site-I looked in many
directories there but didn't see this file. The storm I'm looking at
has no entries in the wwpts file, presumably because the wwpts file
predates the aal and bal files.
>>
>> DEBUG 1: [Source 1 of 1] ADECK Source:
>> /home/qcteam/METv4.1/data/aal142012.dat, Model Suffix: (nul) DEBUG
1:
>> [Source 1 of 1] BDECK Source:
/home/qcteam/METv4.1/data/bal142012.dat,
>> Model Suffix: (nul) DEBUG 1: Config File Default:
>> /home/qcteam/METv4.1/data/config/TCPairsConfig_default
>> DEBUG 1: Config File User:
>> /home/qcteam/METv4.1/data/config/TCPairs.cfg
>> DEBUG 1: Distance to land file:
>> /home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc
>> DEBUG 4: Latitude/Longitude Grid Data:
>> DEBUG 4: lat_ll = 0
>> DEBUG 4: lon_ll = -180
>> DEBUG 4: delta_lat = 0.1
>> DEBUG 4: delta_lon = 0.1
>> DEBUG 4: Nlat = 601
>> DEBUG 4: Nlon = 1801
>> DEBUG 2: Processing 1 ADECK file(s).
>> DEBUG 4: [File 1 of 1] Used 0 of 202726 lines read from file
"/home/qcteam/METv4.1/data/aal142012.dat"
>> DEBUG 3: Used 0 of 202726 lines read from 1 file(s).
>> DEBUG 3: Identified 0 track(s).
>> DEBUG 2: Found 0 ADECK track(s).
>> DEBUG 2: Filtering 0 ADECK tracks based on config file settings.
>> DEBUG 3: Total tracks read = 0
>> DEBUG 3: Total tracks kept = 0
>> DEBUG 3: Rejected for model = 0
>> DEBUG 3: Rejected for storm id = 0
>> DEBUG 3: Rejected for basin = 0
>> DEBUG 3: Rejected for cyclone = 0
>> DEBUG 3: Rejected for storm name = 0
>> DEBUG 3: Rejected for init time = 0
>> DEBUG 3: Rejected for init hour = 0
>> DEBUG 3: Rejected for valid time = 0
>> DEBUG 3: Rejected for init mask = 0
>> DEBUG 3: Rejected for valid mask = 0
>> DEBUG 2: Processing 1 BDECK file(s).
>> DEBUG 4: Updating storm name from "INVEST" to "FOURTEEN" for
AL142012.
>> DEBUG 4: Updating storm name from "FOURTEEN" to "NADINE" for
AL142012.
>> DEBUG 4: [File 1 of 1] Used 184 of 184 lines read from file
"/home/qcteam/METv4.1/data/bal142012.dat"
>> DEBUG 3: Used 184 of 184 lines read from 1 file(s).
>> DEBUG 3: Identified 1 track(s).
>> DEBUG 4: [Track 1 of 1] TrackInfo: StormId = "AL142012", Basin =
"AL",
>> Cyclone = "14", StormName = "NADINE", TechniqueNumber = -9999,
Technique = "BEST", Initials = "(nul)", InitTime = NA, MinValidTime =
20120910_120000, MaxValidTime = 20121004_060000, NPoints = 96, NAlloc
= 100, NTrackLines = 0 DEBUG 2: Found 1 BDECK track(s).
>> DEBUG 2: Deriving 12-hour interpolated tracks.
>> DEBUG 2: Finished deriving 0 12-hour interpolated tracks.
>> DEBUG 2: Deriving 0 ADECK consensus model(s).
>> DEBUG 2: Added 0 ADECK consensus tracks(s).
>> DEBUG 2: Deriving 0 ADECK lag model(s).
>> DEBUG 2: Added 0 ADECK lag tracks(s).
>> DEBUG 2: Deriving 0 CLIPER/SHIFOR baseline model(s).
>> DEBUG 2: Added 0 CLIPER/SHIFOR baseline track(s).
>> DEBUG 2: Matching 0 ADECK tracks to 1 BDECK tracks.
>> DEBUG 1: Watch/Warning file:
>> /home/qcteam/METv4.1/data/tc_data/wwpts_us.txt
>> DEBUG 1: Output file: ./out_tcmpr.tcst
>>
>> My config file contents:
>>
>> model = ["AEMN"];
>>
>> storm_id = ["AL142012"];
>>
>> basin = ["AL"];
>>
>> cyclone = ["14"];
>>
>> storm_name = ["NADINE"];
>>
>> init_beg = "";
>> init_end = "";
>> init_inc = [];
>> init_exc = [];
>>
>> valid_beg = "";
>> valid_end = "";
>>
>> init_hour = [];
>>
>> init_mask = "";
>> valid_mask = "";
>>
>> check_dup = FALSE;
>>
>> interp12 = TRUE;
>>
>> consensus = [];
>>
>> lag_time = [];
>>
>> best_baseline = [];
>> oper_baseline = [];
>>
>> match_points = FALSE;
>>
>> dland_file =
>> "/home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc";
>>
>> watch_warn = {
>> file_name =
"/home/qcteam/METv4.1/data/tc_data/wwpts_us.txt";
>> time_offset = -14400;
>> };
>>
>> version = "V4.1";
>>
>> Hopefully you can suggest something to get this working!
>>
>> Thanks,
>> Matt
>>
>> // SIGNED //
>> Matthew C. Sittel
>> Associate Scientist
>> University Corporation for Atmospheric Research 16WS/WXN, Offutt
AFB,
>> NE
>> (402) 294-3473 DSN: 271-3473
>>
>
>
------------------------------------------------
Subject: Problem with TC_PAIRS
From: matthew.sittel.ctr at offutt.af.mil
Time: Thu Aug 01 10:12:30 2013
It's nice to be doing some other stuff, just for variety's sake. I'm
working on a low-pressure tracking routine as well that will evaluate
model skill. It will have similarities to MET, and ultimately I may
use MET to do the aggregation of statistics... not sure just yet, it's
still in the 'proof of concept' stage.
I can't get tc_stat to run... well, it runs but produces an empty
output file. Here's what I am running:
tc_stat -lookin /home/qcteam/METv4.1/data/tc_data -config
/home/qcteam/METv4.1/data/config/TCStat.cfg -v 4
The "-lookin" directory has the .tcst output file I produced with
TC_PAIRS. All I get in response is this:
DEBUG 1: [Source 1] Lookin: /home/qcteam/METv4.1/data/tc_data
DEBUG 1: Config File Default:
/home/qcteam/METv4.1/data/config/TCStatConfig_default
DEBUG 1: Config File User: /home/qcteam/METv4.1/data/config/TCStat.cfg
DEBUG 1:
DEBUG 1: Creating output file:
/home/qcteam/METv4.1/data/tc_data/tc_stat_output.txt
Here's what is in the config file I created:
amodel = ["AEMN"];
bmodel = ["BEST"];
storm_id = ["AL142012"];
basin = ["AL"];
cyclone = ["14"];
storm_name = [];
init_beg = "";
init_end = "";
init_inc = [];
init_exc = [];
valid_beg = "";
valid_end = "";
valid_inc = [];
valid_exc = [];
lead = [];
init_mask = [];
valid_mask = [];
line_type = [];
track_watch_warn = [];
column_thresh_name = [];
column_thresh_val = [];
column_str_name = [];
column_str_val = [];
init_thresh_name = [];
init_thresh_val = [];
init_str_name = [];
init_str_val = [];
water_only = FALSE;
rapid_inten = FALSE;
rapid_inten_thresh = >=30.0;
landfall = FALSE;
landfall_beg = -86400;
landfall_end = 0;
match_points = FALSE;
event_equal = FALSE;
out_init_mask = "";
out_valid_mask = "";
jobs = [];
version = "V4.1";
Do you see anything obvious that I missed?
Note that the TC 4.1 documentation says "lookin" rather than "-lookin"
when listing the arguments for running TC_STAT.
Thanks,
Matt
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Thursday, August 01, 2013 10:53 AM
To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
Matt,
Sounds interesting. I've actually switched gears recently too. I'm
working on the real-time evaluation for the HFIP products. Currently,
we're mostly just plotting track data, but the plan is to start
running MET-TC in near-realtime mode soon.
I did find a copy-and-paste bug in the tc_pairs tool. I'll post a
bugfix for it shortly and write you back. You should still continue
to omit the storm name though. While the ATCF file format has a place
for the storm name, it usually is not included in the ADECK files. If
you list a storm name in the config file, it'll only keep tracks that
include it. Since ADECK files often don't have it, you won't keep
many of them.
Also, don't worry about the model name suffix being "(nul)". We put
in this option of appending a string to the model names found in the
track files. At one point, we had multiple models to evaluate that
all made use of the same ATCF ID. To distinguish them in the output,
we added a suffix to their ATCF ID's to make them unique. If you're
not in that situation, you don't need to use the suffix.
Thanks,
John
On 08/01/2013 09:43 AM, matthew.sittel.ctr at offutt.af.mil via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>
> Hey John,
>
> Yeah, that change took care of it.
>
> Our ensembles guys are going to make ATCF-formatted tracks for their
WRF members, the ensemble mean, and then compare them to a
deterministic WRF model... all to be compared to best tracks for
storms. TS Gil is in the east Pacific so they've got something to
work with. This will be our chance to utilize MET (which is in my
statement of work, to be the "subject matter expert" in house).
>
> Thanks,
> Matt
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, August 01, 2013 10:37 AM
> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>
> Hello Matt,
>
> Try removing the storm name from the config file, like this:
> storm_name = [];
>
> When you rerun, you should get matching tracks.
>
> However, I think there may be a bug in the config file parsing code,
and I'll let you know what I find.
>
> Thanks,
> John
>
>
> On 07/31/2013 03:40 PM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>
>> Wed Jul 31 15:40:07 2013: Request 62455 was acted upon.
>> Transaction: Ticket created by matthew.sittel.ctr at offutt.af.mil
>> Queue: met_help
>> Subject: Problem with TC_PAIRS
>> Owner: Nobody
>> Requestors: matthew.sittel.ctr at offutt.af.mil
>> Status: new
>> Ticket <URL:
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455
>>>
>>
>>
>> Here we go again with a new issue. I'm testing the tropical
cyclone tracker for an upcoming project we've been asked to do. So I
gathered some sample files from NHC's FTP site and tried to run the
code. I got the TC_DLAND to run just fine, but something appears to
prevent TC_PAIRS from running, and I can't figure out what it is.
>>
>> Here's the log file as I see it. My main concerns are why Model
Suffix is (nul), and why the ADECK file claims to be read but none of
the records are used. I am using the wwpts file included in 4.1; I
can't find a newer one on the NHC FTP site-I looked in many
directories there but didn't see this file. The storm I'm looking at
has no entries in the wwpts file, presumably because the wwpts file
predates the aal and bal files.
>>
>> DEBUG 1: [Source 1 of 1] ADECK Source:
>> /home/qcteam/METv4.1/data/aal142012.dat, Model Suffix: (nul) DEBUG
1:
>> [Source 1 of 1] BDECK Source:
>> /home/qcteam/METv4.1/data/bal142012.dat,
>> Model Suffix: (nul) DEBUG 1: Config File Default:
>> /home/qcteam/METv4.1/data/config/TCPairsConfig_default
>> DEBUG 1: Config File User:
>> /home/qcteam/METv4.1/data/config/TCPairs.cfg
>> DEBUG 1: Distance to land file:
>> /home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc
>> DEBUG 4: Latitude/Longitude Grid Data:
>> DEBUG 4: lat_ll = 0
>> DEBUG 4: lon_ll = -180
>> DEBUG 4: delta_lat = 0.1
>> DEBUG 4: delta_lon = 0.1
>> DEBUG 4: Nlat = 601
>> DEBUG 4: Nlon = 1801
>> DEBUG 2: Processing 1 ADECK file(s).
>> DEBUG 4: [File 1 of 1] Used 0 of 202726 lines read from file
"/home/qcteam/METv4.1/data/aal142012.dat"
>> DEBUG 3: Used 0 of 202726 lines read from 1 file(s).
>> DEBUG 3: Identified 0 track(s).
>> DEBUG 2: Found 0 ADECK track(s).
>> DEBUG 2: Filtering 0 ADECK tracks based on config file settings.
>> DEBUG 3: Total tracks read = 0
>> DEBUG 3: Total tracks kept = 0
>> DEBUG 3: Rejected for model = 0
>> DEBUG 3: Rejected for storm id = 0
>> DEBUG 3: Rejected for basin = 0
>> DEBUG 3: Rejected for cyclone = 0
>> DEBUG 3: Rejected for storm name = 0
>> DEBUG 3: Rejected for init time = 0
>> DEBUG 3: Rejected for init hour = 0
>> DEBUG 3: Rejected for valid time = 0
>> DEBUG 3: Rejected for init mask = 0
>> DEBUG 3: Rejected for valid mask = 0
>> DEBUG 2: Processing 1 BDECK file(s).
>> DEBUG 4: Updating storm name from "INVEST" to "FOURTEEN" for
AL142012.
>> DEBUG 4: Updating storm name from "FOURTEEN" to "NADINE" for
AL142012.
>> DEBUG 4: [File 1 of 1] Used 184 of 184 lines read from file
"/home/qcteam/METv4.1/data/bal142012.dat"
>> DEBUG 3: Used 184 of 184 lines read from 1 file(s).
>> DEBUG 3: Identified 1 track(s).
>> DEBUG 4: [Track 1 of 1] TrackInfo: StormId = "AL142012", Basin =
>> "AL", Cyclone = "14", StormName = "NADINE", TechniqueNumber =
-9999, Technique = "BEST", Initials = "(nul)", InitTime = NA,
MinValidTime = 20120910_120000, MaxValidTime = 20121004_060000,
NPoints = 96, NAlloc = 100, NTrackLines = 0 DEBUG 2: Found 1 BDECK
track(s).
>> DEBUG 2: Deriving 12-hour interpolated tracks.
>> DEBUG 2: Finished deriving 0 12-hour interpolated tracks.
>> DEBUG 2: Deriving 0 ADECK consensus model(s).
>> DEBUG 2: Added 0 ADECK consensus tracks(s).
>> DEBUG 2: Deriving 0 ADECK lag model(s).
>> DEBUG 2: Added 0 ADECK lag tracks(s).
>> DEBUG 2: Deriving 0 CLIPER/SHIFOR baseline model(s).
>> DEBUG 2: Added 0 CLIPER/SHIFOR baseline track(s).
>> DEBUG 2: Matching 0 ADECK tracks to 1 BDECK tracks.
>> DEBUG 1: Watch/Warning file:
>> /home/qcteam/METv4.1/data/tc_data/wwpts_us.txt
>> DEBUG 1: Output file: ./out_tcmpr.tcst
>>
>> My config file contents:
>>
>> model = ["AEMN"];
>>
>> storm_id = ["AL142012"];
>>
>> basin = ["AL"];
>>
>> cyclone = ["14"];
>>
>> storm_name = ["NADINE"];
>>
>> init_beg = "";
>> init_end = "";
>> init_inc = [];
>> init_exc = [];
>>
>> valid_beg = "";
>> valid_end = "";
>>
>> init_hour = [];
>>
>> init_mask = "";
>> valid_mask = "";
>>
>> check_dup = FALSE;
>>
>> interp12 = TRUE;
>>
>> consensus = [];
>>
>> lag_time = [];
>>
>> best_baseline = [];
>> oper_baseline = [];
>>
>> match_points = FALSE;
>>
>> dland_file =
>> "/home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc";
>>
>> watch_warn = {
>> file_name =
"/home/qcteam/METv4.1/data/tc_data/wwpts_us.txt";
>> time_offset = -14400;
>> };
>>
>> version = "V4.1";
>>
>> Hopefully you can suggest something to get this working!
>>
>> Thanks,
>> Matt
>>
>> // SIGNED //
>> Matthew C. Sittel
>> Associate Scientist
>> University Corporation for Atmospheric Research 16WS/WXN, Offutt
AFB,
>> NE
>> (402) 294-3473 DSN: 271-3473
>>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
From: John Halley Gotway
Time: Thu Aug 01 10:26:53 2013
Matt,
Yes, you're missing the "jobs" setting. You need to define what type
of statistics you'd like computed. This is just like the
stat_analysis tool - you can run a single job on the command line or
multiple jobs via a config file. Clearly though, I should add an
error message when a user passes in a config file with no jobs
defined. I'll just add that to the development version.
While there are a few different job types defined for stat_analysis,
currently the only supported job types for tc_stat are "filter" and
"summary". Here's an example of how you might run a summary
job on the command line without a config file:
METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
summary -column TK_ERR,ALTK_ERR,CRTK_ERR
That'll extract the 3 columns (TK_ERR, ALTK_ERR, and CRTK_ERR) and
print summary information (min/max/mean/median, and so on).
Next, try this one:
METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
summary -column TK_ERR -by LEAD
This will summarize the track errors but will give you output for each
unique entry in the "LEAD" column. So you'll see track errors by lead
time. As you'd expect, the track errors grow by lead time.
Hope that helps.
Thanks,
John
On 08/01/2013 10:12 AM, matthew.sittel.ctr at offutt.af.mil via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>
> It's nice to be doing some other stuff, just for variety's sake.
I'm working on a low-pressure tracking routine as well that will
evaluate model skill. It will have similarities to MET, and
ultimately I may use MET to do the aggregation of statistics... not
sure just yet, it's still in the 'proof of concept' stage.
>
> I can't get tc_stat to run... well, it runs but produces an empty
output file. Here's what I am running:
>
> tc_stat -lookin /home/qcteam/METv4.1/data/tc_data -config
/home/qcteam/METv4.1/data/config/TCStat.cfg -v 4
>
> The "-lookin" directory has the .tcst output file I produced with
TC_PAIRS. All I get in response is this:
>
> DEBUG 1: [Source 1] Lookin: /home/qcteam/METv4.1/data/tc_data
> DEBUG 1: Config File Default:
/home/qcteam/METv4.1/data/config/TCStatConfig_default
> DEBUG 1: Config File User:
/home/qcteam/METv4.1/data/config/TCStat.cfg
> DEBUG 1:
> DEBUG 1: Creating output file:
/home/qcteam/METv4.1/data/tc_data/tc_stat_output.txt
>
> Here's what is in the config file I created:
>
> amodel = ["AEMN"];
> bmodel = ["BEST"];
> storm_id = ["AL142012"];
> basin = ["AL"];
> cyclone = ["14"];
> storm_name = [];
> init_beg = "";
> init_end = "";
> init_inc = [];
> init_exc = [];
> valid_beg = "";
> valid_end = "";
> valid_inc = [];
> valid_exc = [];
> lead = [];
> init_mask = [];
> valid_mask = [];
> line_type = [];
> track_watch_warn = [];
> column_thresh_name = [];
> column_thresh_val = [];
> column_str_name = [];
> column_str_val = [];
> init_thresh_name = [];
> init_thresh_val = [];
> init_str_name = [];
> init_str_val = [];
> water_only = FALSE;
> rapid_inten = FALSE;
> rapid_inten_thresh = >=30.0;
> landfall = FALSE;
> landfall_beg = -86400;
> landfall_end = 0;
> match_points = FALSE;
> event_equal = FALSE;
> out_init_mask = "";
> out_valid_mask = "";
> jobs = [];
> version = "V4.1";
>
> Do you see anything obvious that I missed?
>
> Note that the TC 4.1 documentation says "lookin" rather than "-
lookin" when listing the arguments for running TC_STAT.
>
> Thanks,
> Matt
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, August 01, 2013 10:53 AM
> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>
> Matt,
>
> Sounds interesting. I've actually switched gears recently too. I'm
working on the real-time evaluation for the HFIP products. Currently,
we're mostly just plotting track data, but the plan is to start
running MET-TC in near-realtime mode soon.
>
> I did find a copy-and-paste bug in the tc_pairs tool. I'll post a
bugfix for it shortly and write you back. You should still continue
to omit the storm name though. While the ATCF file format has a place
for the storm name, it usually is not included in the ADECK files. If
you list a storm name in the config file, it'll only keep tracks that
include it. Since ADECK files often don't have it, you won't keep
many of them.
>
> Also, don't worry about the model name suffix being "(nul)". We put
in this option of appending a string to the model names found in the
track files. At one point, we had multiple models to evaluate that
all made use of the same ATCF ID. To distinguish them in the output,
we added a suffix to their ATCF ID's to make them unique. If you're
not in that situation, you don't need to use the suffix.
>
> Thanks,
> John
>
> On 08/01/2013 09:43 AM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>>
>> Hey John,
>>
>> Yeah, that change took care of it.
>>
>> Our ensembles guys are going to make ATCF-formatted tracks for
their WRF members, the ensemble mean, and then compare them to a
deterministic WRF model... all to be compared to best tracks for
storms. TS Gil is in the east Pacific so they've got something to
work with. This will be our chance to utilize MET (which is in my
statement of work, to be the "subject matter expert" in house).
>>
>> Thanks,
>> Matt
>>
>> -----Original Message-----
>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> Sent: Thursday, August 01, 2013 10:37 AM
>> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>>
>> Hello Matt,
>>
>> Try removing the storm name from the config file, like this:
>> storm_name = [];
>>
>> When you rerun, you should get matching tracks.
>>
>> However, I think there may be a bug in the config file parsing
code, and I'll let you know what I find.
>>
>> Thanks,
>> John
>>
>>
>> On 07/31/2013 03:40 PM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>>
>>> Wed Jul 31 15:40:07 2013: Request 62455 was acted upon.
>>> Transaction: Ticket created by matthew.sittel.ctr at offutt.af.mil
>>> Queue: met_help
>>> Subject: Problem with TC_PAIRS
>>> Owner: Nobody
>>> Requestors: matthew.sittel.ctr at offutt.af.mil
>>> Status: new
>>> Ticket <URL:
>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455
>>>>
>>>
>>>
>>> Here we go again with a new issue. I'm testing the tropical
cyclone tracker for an upcoming project we've been asked to do. So I
gathered some sample files from NHC's FTP site and tried to run the
code. I got the TC_DLAND to run just fine, but something appears to
prevent TC_PAIRS from running, and I can't figure out what it is.
>>>
>>> Here's the log file as I see it. My main concerns are why Model
Suffix is (nul), and why the ADECK file claims to be read but none of
the records are used. I am using the wwpts file included in 4.1; I
can't find a newer one on the NHC FTP site-I looked in many
directories there but didn't see this file. The storm I'm looking at
has no entries in the wwpts file, presumably because the wwpts file
predates the aal and bal files.
>>>
>>> DEBUG 1: [Source 1 of 1] ADECK Source:
>>> /home/qcteam/METv4.1/data/aal142012.dat, Model Suffix: (nul) DEBUG
1:
>>> [Source 1 of 1] BDECK Source:
>>> /home/qcteam/METv4.1/data/bal142012.dat,
>>> Model Suffix: (nul) DEBUG 1: Config File Default:
>>> /home/qcteam/METv4.1/data/config/TCPairsConfig_default
>>> DEBUG 1: Config File User:
>>> /home/qcteam/METv4.1/data/config/TCPairs.cfg
>>> DEBUG 1: Distance to land file:
>>> /home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc
>>> DEBUG 4: Latitude/Longitude Grid Data:
>>> DEBUG 4: lat_ll = 0
>>> DEBUG 4: lon_ll = -180
>>> DEBUG 4: delta_lat = 0.1
>>> DEBUG 4: delta_lon = 0.1
>>> DEBUG 4: Nlat = 601
>>> DEBUG 4: Nlon = 1801
>>> DEBUG 2: Processing 1 ADECK file(s).
>>> DEBUG 4: [File 1 of 1] Used 0 of 202726 lines read from file
"/home/qcteam/METv4.1/data/aal142012.dat"
>>> DEBUG 3: Used 0 of 202726 lines read from 1 file(s).
>>> DEBUG 3: Identified 0 track(s).
>>> DEBUG 2: Found 0 ADECK track(s).
>>> DEBUG 2: Filtering 0 ADECK tracks based on config file settings.
>>> DEBUG 3: Total tracks read = 0
>>> DEBUG 3: Total tracks kept = 0
>>> DEBUG 3: Rejected for model = 0
>>> DEBUG 3: Rejected for storm id = 0
>>> DEBUG 3: Rejected for basin = 0
>>> DEBUG 3: Rejected for cyclone = 0
>>> DEBUG 3: Rejected for storm name = 0
>>> DEBUG 3: Rejected for init time = 0
>>> DEBUG 3: Rejected for init hour = 0
>>> DEBUG 3: Rejected for valid time = 0
>>> DEBUG 3: Rejected for init mask = 0
>>> DEBUG 3: Rejected for valid mask = 0
>>> DEBUG 2: Processing 1 BDECK file(s).
>>> DEBUG 4: Updating storm name from "INVEST" to "FOURTEEN" for
AL142012.
>>> DEBUG 4: Updating storm name from "FOURTEEN" to "NADINE" for
AL142012.
>>> DEBUG 4: [File 1 of 1] Used 184 of 184 lines read from file
"/home/qcteam/METv4.1/data/bal142012.dat"
>>> DEBUG 3: Used 184 of 184 lines read from 1 file(s).
>>> DEBUG 3: Identified 1 track(s).
>>> DEBUG 4: [Track 1 of 1] TrackInfo: StormId = "AL142012", Basin =
>>> "AL", Cyclone = "14", StormName = "NADINE", TechniqueNumber =
-9999, Technique = "BEST", Initials = "(nul)", InitTime = NA,
MinValidTime = 20120910_120000, MaxValidTime = 20121004_060000,
NPoints = 96, NAlloc = 100, NTrackLines = 0 DEBUG 2: Found 1 BDECK
track(s).
>>> DEBUG 2: Deriving 12-hour interpolated tracks.
>>> DEBUG 2: Finished deriving 0 12-hour interpolated tracks.
>>> DEBUG 2: Deriving 0 ADECK consensus model(s).
>>> DEBUG 2: Added 0 ADECK consensus tracks(s).
>>> DEBUG 2: Deriving 0 ADECK lag model(s).
>>> DEBUG 2: Added 0 ADECK lag tracks(s).
>>> DEBUG 2: Deriving 0 CLIPER/SHIFOR baseline model(s).
>>> DEBUG 2: Added 0 CLIPER/SHIFOR baseline track(s).
>>> DEBUG 2: Matching 0 ADECK tracks to 1 BDECK tracks.
>>> DEBUG 1: Watch/Warning file:
>>> /home/qcteam/METv4.1/data/tc_data/wwpts_us.txt
>>> DEBUG 1: Output file: ./out_tcmpr.tcst
>>>
>>> My config file contents:
>>>
>>> model = ["AEMN"];
>>>
>>> storm_id = ["AL142012"];
>>>
>>> basin = ["AL"];
>>>
>>> cyclone = ["14"];
>>>
>>> storm_name = ["NADINE"];
>>>
>>> init_beg = "";
>>> init_end = "";
>>> init_inc = [];
>>> init_exc = [];
>>>
>>> valid_beg = "";
>>> valid_end = "";
>>>
>>> init_hour = [];
>>>
>>> init_mask = "";
>>> valid_mask = "";
>>>
>>> check_dup = FALSE;
>>>
>>> interp12 = TRUE;
>>>
>>> consensus = [];
>>>
>>> lag_time = [];
>>>
>>> best_baseline = [];
>>> oper_baseline = [];
>>>
>>> match_points = FALSE;
>>>
>>> dland_file =
>>> "/home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc";
>>>
>>> watch_warn = {
>>> file_name =
"/home/qcteam/METv4.1/data/tc_data/wwpts_us.txt";
>>> time_offset = -14400;
>>> };
>>>
>>> version = "V4.1";
>>>
>>> Hopefully you can suggest something to get this working!
>>>
>>> Thanks,
>>> Matt
>>>
>>> // SIGNED //
>>> Matthew C. Sittel
>>> Associate Scientist
>>> University Corporation for Atmospheric Research 16WS/WXN, Offutt
AFB,
>>> NE
>>> (402) 294-3473 DSN: 271-3473
>>>
>>
>>
>
>
------------------------------------------------
Subject: Problem with TC_PAIRS
From: matthew.sittel.ctr at offutt.af.mil
Time: Thu Aug 01 11:12:10 2013
That took care of it. Thanks.
Is there a field for wind error? That's the other one I need to
build. Basically the wind difference "forecast minus best track
speed". I was thinking I could get at it, indirectly, by summarizing
AMAX_WIND and BMAX_WIND, but that would not give me pairwise
differences or any sort of stratification of errors as a function of
wind speed.
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Thursday, August 01, 2013 11:27 AM
To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
Matt,
Yes, you're missing the "jobs" setting. You need to define what type
of statistics you'd like computed. This is just like the
stat_analysis tool - you can run a single job on the command line or
multiple jobs via a config file. Clearly though, I should add an
error message when a user passes in a config file with no jobs
defined. I'll just add that to the development version.
While there are a few different job types defined for stat_analysis,
currently the only supported job types for tc_stat are "filter" and
"summary". Here's an example of how you might run a summary job on
the command line without a config file:
METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
summary -column TK_ERR,ALTK_ERR,CRTK_ERR
That'll extract the 3 columns (TK_ERR, ALTK_ERR, and CRTK_ERR) and
print summary information (min/max/mean/median, and so on).
Next, try this one:
METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
summary -column TK_ERR -by LEAD
This will summarize the track errors but will give you output for each
unique entry in the "LEAD" column. So you'll see track errors by lead
time. As you'd expect, the track errors grow by lead time.
Hope that helps.
Thanks,
John
On 08/01/2013 10:12 AM, matthew.sittel.ctr at offutt.af.mil via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>
> It's nice to be doing some other stuff, just for variety's sake.
I'm working on a low-pressure tracking routine as well that will
evaluate model skill. It will have similarities to MET, and
ultimately I may use MET to do the aggregation of statistics... not
sure just yet, it's still in the 'proof of concept' stage.
>
> I can't get tc_stat to run... well, it runs but produces an empty
output file. Here's what I am running:
>
> tc_stat -lookin /home/qcteam/METv4.1/data/tc_data -config
> /home/qcteam/METv4.1/data/config/TCStat.cfg -v 4
>
> The "-lookin" directory has the .tcst output file I produced with
TC_PAIRS. All I get in response is this:
>
> DEBUG 1: [Source 1] Lookin: /home/qcteam/METv4.1/data/tc_data DEBUG
1:
> Config File Default:
> /home/qcteam/METv4.1/data/config/TCStatConfig_default
> DEBUG 1: Config File User:
/home/qcteam/METv4.1/data/config/TCStat.cfg
> DEBUG 1:
> DEBUG 1: Creating output file:
> /home/qcteam/METv4.1/data/tc_data/tc_stat_output.txt
>
> Here's what is in the config file I created:
>
> amodel = ["AEMN"];
> bmodel = ["BEST"];
> storm_id = ["AL142012"];
> basin = ["AL"];
> cyclone = ["14"];
> storm_name = [];
> init_beg = "";
> init_end = "";
> init_inc = [];
> init_exc = [];
> valid_beg = "";
> valid_end = "";
> valid_inc = [];
> valid_exc = [];
> lead = [];
> init_mask = [];
> valid_mask = [];
> line_type = [];
> track_watch_warn = [];
> column_thresh_name = [];
> column_thresh_val = [];
> column_str_name = [];
> column_str_val = [];
> init_thresh_name = [];
> init_thresh_val = [];
> init_str_name = [];
> init_str_val = [];
> water_only = FALSE;
> rapid_inten = FALSE;
> rapid_inten_thresh = >=30.0;
> landfall = FALSE;
> landfall_beg = -86400;
> landfall_end = 0;
> match_points = FALSE;
> event_equal = FALSE;
> out_init_mask = "";
> out_valid_mask = "";
> jobs = [];
> version = "V4.1";
>
> Do you see anything obvious that I missed?
>
> Note that the TC 4.1 documentation says "lookin" rather than "-
lookin" when listing the arguments for running TC_STAT.
>
> Thanks,
> Matt
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, August 01, 2013 10:53 AM
> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>
> Matt,
>
> Sounds interesting. I've actually switched gears recently too. I'm
working on the real-time evaluation for the HFIP products. Currently,
we're mostly just plotting track data, but the plan is to start
running MET-TC in near-realtime mode soon.
>
> I did find a copy-and-paste bug in the tc_pairs tool. I'll post a
bugfix for it shortly and write you back. You should still continue
to omit the storm name though. While the ATCF file format has a place
for the storm name, it usually is not included in the ADECK files. If
you list a storm name in the config file, it'll only keep tracks that
include it. Since ADECK files often don't have it, you won't keep
many of them.
>
> Also, don't worry about the model name suffix being "(nul)". We put
in this option of appending a string to the model names found in the
track files. At one point, we had multiple models to evaluate that
all made use of the same ATCF ID. To distinguish them in the output,
we added a suffix to their ATCF ID's to make them unique. If you're
not in that situation, you don't need to use the suffix.
>
> Thanks,
> John
>
> On 08/01/2013 09:43 AM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>>
>> Hey John,
>>
>> Yeah, that change took care of it.
>>
>> Our ensembles guys are going to make ATCF-formatted tracks for
their WRF members, the ensemble mean, and then compare them to a
deterministic WRF model... all to be compared to best tracks for
storms. TS Gil is in the east Pacific so they've got something to
work with. This will be our chance to utilize MET (which is in my
statement of work, to be the "subject matter expert" in house).
>>
>> Thanks,
>> Matt
>>
>> -----Original Message-----
>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> Sent: Thursday, August 01, 2013 10:37 AM
>> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>>
>> Hello Matt,
>>
>> Try removing the storm name from the config file, like this:
>> storm_name = [];
>>
>> When you rerun, you should get matching tracks.
>>
>> However, I think there may be a bug in the config file parsing
code, and I'll let you know what I find.
>>
>> Thanks,
>> John
>>
>>
>> On 07/31/2013 03:40 PM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>>
>>> Wed Jul 31 15:40:07 2013: Request 62455 was acted upon.
>>> Transaction: Ticket created by matthew.sittel.ctr at offutt.af.mil
>>> Queue: met_help
>>> Subject: Problem with TC_PAIRS
>>> Owner: Nobody
>>> Requestors: matthew.sittel.ctr at offutt.af.mil
>>> Status: new
>>> Ticket <URL:
>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455
>>>>
>>>
>>>
>>> Here we go again with a new issue. I'm testing the tropical
cyclone tracker for an upcoming project we've been asked to do. So I
gathered some sample files from NHC's FTP site and tried to run the
code. I got the TC_DLAND to run just fine, but something appears to
prevent TC_PAIRS from running, and I can't figure out what it is.
>>>
>>> Here's the log file as I see it. My main concerns are why Model
Suffix is (nul), and why the ADECK file claims to be read but none of
the records are used. I am using the wwpts file included in 4.1; I
can't find a newer one on the NHC FTP site-I looked in many
directories there but didn't see this file. The storm I'm looking at
has no entries in the wwpts file, presumably because the wwpts file
predates the aal and bal files.
>>>
>>> DEBUG 1: [Source 1 of 1] ADECK Source:
>>> /home/qcteam/METv4.1/data/aal142012.dat, Model Suffix: (nul) DEBUG
1:
>>> [Source 1 of 1] BDECK Source:
>>> /home/qcteam/METv4.1/data/bal142012.dat,
>>> Model Suffix: (nul) DEBUG 1: Config File Default:
>>> /home/qcteam/METv4.1/data/config/TCPairsConfig_default
>>> DEBUG 1: Config File User:
>>> /home/qcteam/METv4.1/data/config/TCPairs.cfg
>>> DEBUG 1: Distance to land file:
>>> /home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc
>>> DEBUG 4: Latitude/Longitude Grid Data:
>>> DEBUG 4: lat_ll = 0
>>> DEBUG 4: lon_ll = -180
>>> DEBUG 4: delta_lat = 0.1
>>> DEBUG 4: delta_lon = 0.1
>>> DEBUG 4: Nlat = 601
>>> DEBUG 4: Nlon = 1801
>>> DEBUG 2: Processing 1 ADECK file(s).
>>> DEBUG 4: [File 1 of 1] Used 0 of 202726 lines read from file
"/home/qcteam/METv4.1/data/aal142012.dat"
>>> DEBUG 3: Used 0 of 202726 lines read from 1 file(s).
>>> DEBUG 3: Identified 0 track(s).
>>> DEBUG 2: Found 0 ADECK track(s).
>>> DEBUG 2: Filtering 0 ADECK tracks based on config file settings.
>>> DEBUG 3: Total tracks read = 0
>>> DEBUG 3: Total tracks kept = 0
>>> DEBUG 3: Rejected for model = 0
>>> DEBUG 3: Rejected for storm id = 0
>>> DEBUG 3: Rejected for basin = 0
>>> DEBUG 3: Rejected for cyclone = 0
>>> DEBUG 3: Rejected for storm name = 0 DEBUG 3: Rejected for init
time
>>> = 0 DEBUG 3: Rejected for init hour = 0 DEBUG 3: Rejected for
valid
>>> time = 0 DEBUG 3: Rejected for init mask = 0 DEBUG 3: Rejected
for
>>> valid mask = 0 DEBUG 2: Processing 1 BDECK file(s).
>>> DEBUG 4: Updating storm name from "INVEST" to "FOURTEEN" for
AL142012.
>>> DEBUG 4: Updating storm name from "FOURTEEN" to "NADINE" for
AL142012.
>>> DEBUG 4: [File 1 of 1] Used 184 of 184 lines read from file
"/home/qcteam/METv4.1/data/bal142012.dat"
>>> DEBUG 3: Used 184 of 184 lines read from 1 file(s).
>>> DEBUG 3: Identified 1 track(s).
>>> DEBUG 4: [Track 1 of 1] TrackInfo: StormId = "AL142012", Basin =
>>> "AL", Cyclone = "14", StormName = "NADINE", TechniqueNumber =
-9999, Technique = "BEST", Initials = "(nul)", InitTime = NA,
MinValidTime = 20120910_120000, MaxValidTime = 20121004_060000,
NPoints = 96, NAlloc = 100, NTrackLines = 0 DEBUG 2: Found 1 BDECK
track(s).
>>> DEBUG 2: Deriving 12-hour interpolated tracks.
>>> DEBUG 2: Finished deriving 0 12-hour interpolated tracks.
>>> DEBUG 2: Deriving 0 ADECK consensus model(s).
>>> DEBUG 2: Added 0 ADECK consensus tracks(s).
>>> DEBUG 2: Deriving 0 ADECK lag model(s).
>>> DEBUG 2: Added 0 ADECK lag tracks(s).
>>> DEBUG 2: Deriving 0 CLIPER/SHIFOR baseline model(s).
>>> DEBUG 2: Added 0 CLIPER/SHIFOR baseline track(s).
>>> DEBUG 2: Matching 0 ADECK tracks to 1 BDECK tracks.
>>> DEBUG 1: Watch/Warning file:
>>> /home/qcteam/METv4.1/data/tc_data/wwpts_us.txt
>>> DEBUG 1: Output file: ./out_tcmpr.tcst
>>>
>>> My config file contents:
>>>
>>> model = ["AEMN"];
>>>
>>> storm_id = ["AL142012"];
>>>
>>> basin = ["AL"];
>>>
>>> cyclone = ["14"];
>>>
>>> storm_name = ["NADINE"];
>>>
>>> init_beg = "";
>>> init_end = "";
>>> init_inc = [];
>>> init_exc = [];
>>>
>>> valid_beg = "";
>>> valid_end = "";
>>>
>>> init_hour = [];
>>>
>>> init_mask = "";
>>> valid_mask = "";
>>>
>>> check_dup = FALSE;
>>>
>>> interp12 = TRUE;
>>>
>>> consensus = [];
>>>
>>> lag_time = [];
>>>
>>> best_baseline = [];
>>> oper_baseline = [];
>>>
>>> match_points = FALSE;
>>>
>>> dland_file =
>>> "/home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc";
>>>
>>> watch_warn = {
>>> file_name =
"/home/qcteam/METv4.1/data/tc_data/wwpts_us.txt";
>>> time_offset = -14400;
>>> };
>>>
>>> version = "V4.1";
>>>
>>> Hopefully you can suggest something to get this working!
>>>
>>> Thanks,
>>> Matt
>>>
>>> // SIGNED //
>>> Matthew C. Sittel
>>> Associate Scientist
>>> University Corporation for Atmospheric Research 16WS/WXN, Offutt
>>> AFB, NE
>>> (402) 294-3473 DSN: 271-3473
>>>
>>
>>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
From: John Halley Gotway
Time: Thu Aug 01 11:53:50 2013
Matt,
Yes, the "-column" option supports differences. For example...
METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
summary -column AMAX_WIND-BMAX_WIND -by LEAD
However when I ran this job, I just got all NA's in the output. Upon
closer inspection, I see that the AMAX_WIND column is all NA's in the
.tcst output of tc_pairs. Taking a step back, I see that
the AEMN track data in the ATCF file does not contain valid wind
speed. Is that perhaps a track only ensemble mean forecast?
Thanks,
John
On 08/01/2013 11:12 AM, matthew.sittel.ctr at offutt.af.mil via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>
> That took care of it. Thanks.
>
> Is there a field for wind error? That's the other one I need to
build. Basically the wind difference "forecast minus best track
speed". I was thinking I could get at it, indirectly, by summarizing
AMAX_WIND and BMAX_WIND, but that would not give me pairwise
differences or any sort of stratification of errors as a function of
wind speed.
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, August 01, 2013 11:27 AM
> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>
> Matt,
>
> Yes, you're missing the "jobs" setting. You need to define what
type of statistics you'd like computed. This is just like the
stat_analysis tool - you can run a single job on the command line or
multiple jobs via a config file. Clearly though, I should add an
error message when a user passes in a config file with no jobs
defined. I'll just add that to the development version.
>
> While there are a few different job types defined for stat_analysis,
currently the only supported job types for tc_stat are "filter" and
"summary". Here's an example of how you might run a summary job on
the command line without a config file:
>
> METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
summary -column TK_ERR,ALTK_ERR,CRTK_ERR
>
> That'll extract the 3 columns (TK_ERR, ALTK_ERR, and CRTK_ERR) and
print summary information (min/max/mean/median, and so on).
>
> Next, try this one:
> METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
summary -column TK_ERR -by LEAD
>
> This will summarize the track errors but will give you output for
each unique entry in the "LEAD" column. So you'll see track errors by
lead time. As you'd expect, the track errors grow by lead time.
>
> Hope that helps.
>
> Thanks,
> John
>
>
> On 08/01/2013 10:12 AM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>>
>> It's nice to be doing some other stuff, just for variety's sake.
I'm working on a low-pressure tracking routine as well that will
evaluate model skill. It will have similarities to MET, and
ultimately I may use MET to do the aggregation of statistics... not
sure just yet, it's still in the 'proof of concept' stage.
>>
>> I can't get tc_stat to run... well, it runs but produces an empty
output file. Here's what I am running:
>>
>> tc_stat -lookin /home/qcteam/METv4.1/data/tc_data -config
>> /home/qcteam/METv4.1/data/config/TCStat.cfg -v 4
>>
>> The "-lookin" directory has the .tcst output file I produced with
TC_PAIRS. All I get in response is this:
>>
>> DEBUG 1: [Source 1] Lookin: /home/qcteam/METv4.1/data/tc_data DEBUG
1:
>> Config File Default:
>> /home/qcteam/METv4.1/data/config/TCStatConfig_default
>> DEBUG 1: Config File User:
/home/qcteam/METv4.1/data/config/TCStat.cfg
>> DEBUG 1:
>> DEBUG 1: Creating output file:
>> /home/qcteam/METv4.1/data/tc_data/tc_stat_output.txt
>>
>> Here's what is in the config file I created:
>>
>> amodel = ["AEMN"];
>> bmodel = ["BEST"];
>> storm_id = ["AL142012"];
>> basin = ["AL"];
>> cyclone = ["14"];
>> storm_name = [];
>> init_beg = "";
>> init_end = "";
>> init_inc = [];
>> init_exc = [];
>> valid_beg = "";
>> valid_end = "";
>> valid_inc = [];
>> valid_exc = [];
>> lead = [];
>> init_mask = [];
>> valid_mask = [];
>> line_type = [];
>> track_watch_warn = [];
>> column_thresh_name = [];
>> column_thresh_val = [];
>> column_str_name = [];
>> column_str_val = [];
>> init_thresh_name = [];
>> init_thresh_val = [];
>> init_str_name = [];
>> init_str_val = [];
>> water_only = FALSE;
>> rapid_inten = FALSE;
>> rapid_inten_thresh = >=30.0;
>> landfall = FALSE;
>> landfall_beg = -86400;
>> landfall_end = 0;
>> match_points = FALSE;
>> event_equal = FALSE;
>> out_init_mask = "";
>> out_valid_mask = "";
>> jobs = [];
>> version = "V4.1";
>>
>> Do you see anything obvious that I missed?
>>
>> Note that the TC 4.1 documentation says "lookin" rather than "-
lookin" when listing the arguments for running TC_STAT.
>>
>> Thanks,
>> Matt
>>
>> -----Original Message-----
>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> Sent: Thursday, August 01, 2013 10:53 AM
>> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>>
>> Matt,
>>
>> Sounds interesting. I've actually switched gears recently too.
I'm working on the real-time evaluation for the HFIP products.
Currently, we're mostly just plotting track data, but the plan is to
start running MET-TC in near-realtime mode soon.
>>
>> I did find a copy-and-paste bug in the tc_pairs tool. I'll post a
bugfix for it shortly and write you back. You should still continue
to omit the storm name though. While the ATCF file format has a place
for the storm name, it usually is not included in the ADECK files. If
you list a storm name in the config file, it'll only keep tracks that
include it. Since ADECK files often don't have it, you won't keep
many of them.
>>
>> Also, don't worry about the model name suffix being "(nul)". We
put in this option of appending a string to the model names found in
the track files. At one point, we had multiple models to evaluate
that all made use of the same ATCF ID. To distinguish them in the
output, we added a suffix to their ATCF ID's to make them unique. If
you're not in that situation, you don't need to use the suffix.
>>
>> Thanks,
>> John
>>
>> On 08/01/2013 09:43 AM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>>>
>>> Hey John,
>>>
>>> Yeah, that change took care of it.
>>>
>>> Our ensembles guys are going to make ATCF-formatted tracks for
their WRF members, the ensemble mean, and then compare them to a
deterministic WRF model... all to be compared to best tracks for
storms. TS Gil is in the east Pacific so they've got something to
work with. This will be our chance to utilize MET (which is in my
statement of work, to be the "subject matter expert" in house).
>>>
>>> Thanks,
>>> Matt
>>>
>>> -----Original Message-----
>>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>>> Sent: Thursday, August 01, 2013 10:37 AM
>>> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
>>> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>>>
>>> Hello Matt,
>>>
>>> Try removing the storm name from the config file, like this:
>>> storm_name = [];
>>>
>>> When you rerun, you should get matching tracks.
>>>
>>> However, I think there may be a bug in the config file parsing
code, and I'll let you know what I find.
>>>
>>> Thanks,
>>> John
>>>
>>>
>>> On 07/31/2013 03:40 PM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>>>
>>>> Wed Jul 31 15:40:07 2013: Request 62455 was acted upon.
>>>> Transaction: Ticket created by matthew.sittel.ctr at offutt.af.mil
>>>> Queue: met_help
>>>> Subject: Problem with TC_PAIRS
>>>> Owner: Nobody
>>>> Requestors: matthew.sittel.ctr at offutt.af.mil
>>>> Status: new
>>>> Ticket <URL:
>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455
>>>>>
>>>>
>>>>
>>>> Here we go again with a new issue. I'm testing the tropical
cyclone tracker for an upcoming project we've been asked to do. So I
gathered some sample files from NHC's FTP site and tried to run the
code. I got the TC_DLAND to run just fine, but something appears to
prevent TC_PAIRS from running, and I can't figure out what it is.
>>>>
>>>> Here's the log file as I see it. My main concerns are why Model
Suffix is (nul), and why the ADECK file claims to be read but none of
the records are used. I am using the wwpts file included in 4.1; I
can't find a newer one on the NHC FTP site-I looked in many
directories there but didn't see this file. The storm I'm looking at
has no entries in the wwpts file, presumably because the wwpts file
predates the aal and bal files.
>>>>
>>>> DEBUG 1: [Source 1 of 1] ADECK Source:
>>>> /home/qcteam/METv4.1/data/aal142012.dat, Model Suffix: (nul)
DEBUG 1:
>>>> [Source 1 of 1] BDECK Source:
>>>> /home/qcteam/METv4.1/data/bal142012.dat,
>>>> Model Suffix: (nul) DEBUG 1: Config File Default:
>>>> /home/qcteam/METv4.1/data/config/TCPairsConfig_default
>>>> DEBUG 1: Config File User:
>>>> /home/qcteam/METv4.1/data/config/TCPairs.cfg
>>>> DEBUG 1: Distance to land file:
>>>> /home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc
>>>> DEBUG 4: Latitude/Longitude Grid Data:
>>>> DEBUG 4: lat_ll = 0
>>>> DEBUG 4: lon_ll = -180
>>>> DEBUG 4: delta_lat = 0.1
>>>> DEBUG 4: delta_lon = 0.1
>>>> DEBUG 4: Nlat = 601
>>>> DEBUG 4: Nlon = 1801
>>>> DEBUG 2: Processing 1 ADECK file(s).
>>>> DEBUG 4: [File 1 of 1] Used 0 of 202726 lines read from file
"/home/qcteam/METv4.1/data/aal142012.dat"
>>>> DEBUG 3: Used 0 of 202726 lines read from 1 file(s).
>>>> DEBUG 3: Identified 0 track(s).
>>>> DEBUG 2: Found 0 ADECK track(s).
>>>> DEBUG 2: Filtering 0 ADECK tracks based on config file settings.
>>>> DEBUG 3: Total tracks read = 0
>>>> DEBUG 3: Total tracks kept = 0
>>>> DEBUG 3: Rejected for model = 0
>>>> DEBUG 3: Rejected for storm id = 0
>>>> DEBUG 3: Rejected for basin = 0
>>>> DEBUG 3: Rejected for cyclone = 0
>>>> DEBUG 3: Rejected for storm name = 0 DEBUG 3: Rejected for init
time
>>>> = 0 DEBUG 3: Rejected for init hour = 0 DEBUG 3: Rejected for
valid
>>>> time = 0 DEBUG 3: Rejected for init mask = 0 DEBUG 3: Rejected
for
>>>> valid mask = 0 DEBUG 2: Processing 1 BDECK file(s).
>>>> DEBUG 4: Updating storm name from "INVEST" to "FOURTEEN" for
AL142012.
>>>> DEBUG 4: Updating storm name from "FOURTEEN" to "NADINE" for
AL142012.
>>>> DEBUG 4: [File 1 of 1] Used 184 of 184 lines read from file
"/home/qcteam/METv4.1/data/bal142012.dat"
>>>> DEBUG 3: Used 184 of 184 lines read from 1 file(s).
>>>> DEBUG 3: Identified 1 track(s).
>>>> DEBUG 4: [Track 1 of 1] TrackInfo: StormId = "AL142012", Basin =
>>>> "AL", Cyclone = "14", StormName = "NADINE", TechniqueNumber =
-9999, Technique = "BEST", Initials = "(nul)", InitTime = NA,
MinValidTime = 20120910_120000, MaxValidTime = 20121004_060000,
NPoints = 96, NAlloc = 100, NTrackLines = 0 DEBUG 2: Found 1 BDECK
track(s).
>>>> DEBUG 2: Deriving 12-hour interpolated tracks.
>>>> DEBUG 2: Finished deriving 0 12-hour interpolated tracks.
>>>> DEBUG 2: Deriving 0 ADECK consensus model(s).
>>>> DEBUG 2: Added 0 ADECK consensus tracks(s).
>>>> DEBUG 2: Deriving 0 ADECK lag model(s).
>>>> DEBUG 2: Added 0 ADECK lag tracks(s).
>>>> DEBUG 2: Deriving 0 CLIPER/SHIFOR baseline model(s).
>>>> DEBUG 2: Added 0 CLIPER/SHIFOR baseline track(s).
>>>> DEBUG 2: Matching 0 ADECK tracks to 1 BDECK tracks.
>>>> DEBUG 1: Watch/Warning file:
>>>> /home/qcteam/METv4.1/data/tc_data/wwpts_us.txt
>>>> DEBUG 1: Output file: ./out_tcmpr.tcst
>>>>
>>>> My config file contents:
>>>>
>>>> model = ["AEMN"];
>>>>
>>>> storm_id = ["AL142012"];
>>>>
>>>> basin = ["AL"];
>>>>
>>>> cyclone = ["14"];
>>>>
>>>> storm_name = ["NADINE"];
>>>>
>>>> init_beg = "";
>>>> init_end = "";
>>>> init_inc = [];
>>>> init_exc = [];
>>>>
>>>> valid_beg = "";
>>>> valid_end = "";
>>>>
>>>> init_hour = [];
>>>>
>>>> init_mask = "";
>>>> valid_mask = "";
>>>>
>>>> check_dup = FALSE;
>>>>
>>>> interp12 = TRUE;
>>>>
>>>> consensus = [];
>>>>
>>>> lag_time = [];
>>>>
>>>> best_baseline = [];
>>>> oper_baseline = [];
>>>>
>>>> match_points = FALSE;
>>>>
>>>> dland_file =
>>>> "/home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc";
>>>>
>>>> watch_warn = {
>>>> file_name =
"/home/qcteam/METv4.1/data/tc_data/wwpts_us.txt";
>>>> time_offset = -14400;
>>>> };
>>>>
>>>> version = "V4.1";
>>>>
>>>> Hopefully you can suggest something to get this working!
>>>>
>>>> Thanks,
>>>> Matt
>>>>
>>>> // SIGNED //
>>>> Matthew C. Sittel
>>>> Associate Scientist
>>>> University Corporation for Atmospheric Research 16WS/WXN, Offutt
>>>> AFB, NE
>>>> (402) 294-3473 DSN: 271-3473
>>>>
>>>
>>>
>>
>>
>
>
------------------------------------------------
Subject: Problem with TC_PAIRS
From: matthew.sittel.ctr at offutt.af.mil
Time: Thu Aug 01 12:57:59 2013
Sure looks that way. I'm testing using a single ensemble member
instead as those do appear to have winds... I think our cases will.
I didn't realize I could do differences... I didn't gather that from
the PDF... maybe I should re-read it!
-----Original Message-----
From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Sent: Thursday, August 01, 2013 12:54 PM
To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
Matt,
Yes, the "-column" option supports differences. For example...
METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
summary -column AMAX_WIND-BMAX_WIND -by LEAD
However when I ran this job, I just got all NA's in the output. Upon
closer inspection, I see that the AMAX_WIND column is all NA's in the
.tcst output of tc_pairs. Taking a step back, I see that the AEMN
track data in the ATCF file does not contain valid wind speed. Is
that perhaps a track only ensemble mean forecast?
Thanks,
John
On 08/01/2013 11:12 AM, matthew.sittel.ctr at offutt.af.mil via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>
> That took care of it. Thanks.
>
> Is there a field for wind error? That's the other one I need to
build. Basically the wind difference "forecast minus best track
speed". I was thinking I could get at it, indirectly, by summarizing
AMAX_WIND and BMAX_WIND, but that would not give me pairwise
differences or any sort of stratification of errors as a function of
wind speed.
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, August 01, 2013 11:27 AM
> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>
> Matt,
>
> Yes, you're missing the "jobs" setting. You need to define what
type of statistics you'd like computed. This is just like the
stat_analysis tool - you can run a single job on the command line or
multiple jobs via a config file. Clearly though, I should add an
error message when a user passes in a config file with no jobs
defined. I'll just add that to the development version.
>
> While there are a few different job types defined for stat_analysis,
currently the only supported job types for tc_stat are "filter" and
"summary". Here's an example of how you might run a summary job on
the command line without a config file:
>
> METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
> summary -column TK_ERR,ALTK_ERR,CRTK_ERR
>
> That'll extract the 3 columns (TK_ERR, ALTK_ERR, and CRTK_ERR) and
print summary information (min/max/mean/median, and so on).
>
> Next, try this one:
> METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
> summary -column TK_ERR -by LEAD
>
> This will summarize the track errors but will give you output for
each unique entry in the "LEAD" column. So you'll see track errors by
lead time. As you'd expect, the track errors grow by lead time.
>
> Hope that helps.
>
> Thanks,
> John
>
>
> On 08/01/2013 10:12 AM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>>
>> It's nice to be doing some other stuff, just for variety's sake.
I'm working on a low-pressure tracking routine as well that will
evaluate model skill. It will have similarities to MET, and
ultimately I may use MET to do the aggregation of statistics... not
sure just yet, it's still in the 'proof of concept' stage.
>>
>> I can't get tc_stat to run... well, it runs but produces an empty
output file. Here's what I am running:
>>
>> tc_stat -lookin /home/qcteam/METv4.1/data/tc_data -config
>> /home/qcteam/METv4.1/data/config/TCStat.cfg -v 4
>>
>> The "-lookin" directory has the .tcst output file I produced with
TC_PAIRS. All I get in response is this:
>>
>> DEBUG 1: [Source 1] Lookin: /home/qcteam/METv4.1/data/tc_data DEBUG
1:
>> Config File Default:
>> /home/qcteam/METv4.1/data/config/TCStatConfig_default
>> DEBUG 1: Config File User:
>> /home/qcteam/METv4.1/data/config/TCStat.cfg
>> DEBUG 1:
>> DEBUG 1: Creating output file:
>> /home/qcteam/METv4.1/data/tc_data/tc_stat_output.txt
>>
>> Here's what is in the config file I created:
>>
>> amodel = ["AEMN"];
>> bmodel = ["BEST"];
>> storm_id = ["AL142012"];
>> basin = ["AL"];
>> cyclone = ["14"];
>> storm_name = [];
>> init_beg = "";
>> init_end = "";
>> init_inc = [];
>> init_exc = [];
>> valid_beg = "";
>> valid_end = "";
>> valid_inc = [];
>> valid_exc = [];
>> lead = [];
>> init_mask = [];
>> valid_mask = [];
>> line_type = [];
>> track_watch_warn = [];
>> column_thresh_name = [];
>> column_thresh_val = [];
>> column_str_name = [];
>> column_str_val = [];
>> init_thresh_name = [];
>> init_thresh_val = [];
>> init_str_name = [];
>> init_str_val = [];
>> water_only = FALSE;
>> rapid_inten = FALSE;
>> rapid_inten_thresh = >=30.0;
>> landfall = FALSE;
>> landfall_beg = -86400;
>> landfall_end = 0;
>> match_points = FALSE;
>> event_equal = FALSE;
>> out_init_mask = "";
>> out_valid_mask = "";
>> jobs = [];
>> version = "V4.1";
>>
>> Do you see anything obvious that I missed?
>>
>> Note that the TC 4.1 documentation says "lookin" rather than "-
lookin" when listing the arguments for running TC_STAT.
>>
>> Thanks,
>> Matt
>>
>> -----Original Message-----
>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> Sent: Thursday, August 01, 2013 10:53 AM
>> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>>
>> Matt,
>>
>> Sounds interesting. I've actually switched gears recently too.
I'm working on the real-time evaluation for the HFIP products.
Currently, we're mostly just plotting track data, but the plan is to
start running MET-TC in near-realtime mode soon.
>>
>> I did find a copy-and-paste bug in the tc_pairs tool. I'll post a
bugfix for it shortly and write you back. You should still continue
to omit the storm name though. While the ATCF file format has a place
for the storm name, it usually is not included in the ADECK files. If
you list a storm name in the config file, it'll only keep tracks that
include it. Since ADECK files often don't have it, you won't keep
many of them.
>>
>> Also, don't worry about the model name suffix being "(nul)". We
put in this option of appending a string to the model names found in
the track files. At one point, we had multiple models to evaluate
that all made use of the same ATCF ID. To distinguish them in the
output, we added a suffix to their ATCF ID's to make them unique. If
you're not in that situation, you don't need to use the suffix.
>>
>> Thanks,
>> John
>>
>> On 08/01/2013 09:43 AM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>>>
>>> Hey John,
>>>
>>> Yeah, that change took care of it.
>>>
>>> Our ensembles guys are going to make ATCF-formatted tracks for
their WRF members, the ensemble mean, and then compare them to a
deterministic WRF model... all to be compared to best tracks for
storms. TS Gil is in the east Pacific so they've got something to
work with. This will be our chance to utilize MET (which is in my
statement of work, to be the "subject matter expert" in house).
>>>
>>> Thanks,
>>> Matt
>>>
>>> -----Original Message-----
>>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>>> Sent: Thursday, August 01, 2013 10:37 AM
>>> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
>>> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>>>
>>> Hello Matt,
>>>
>>> Try removing the storm name from the config file, like this:
>>> storm_name = [];
>>>
>>> When you rerun, you should get matching tracks.
>>>
>>> However, I think there may be a bug in the config file parsing
code, and I'll let you know what I find.
>>>
>>> Thanks,
>>> John
>>>
>>>
>>> On 07/31/2013 03:40 PM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>>>
>>>> Wed Jul 31 15:40:07 2013: Request 62455 was acted upon.
>>>> Transaction: Ticket created by matthew.sittel.ctr at offutt.af.mil
>>>> Queue: met_help
>>>> Subject: Problem with TC_PAIRS
>>>> Owner: Nobody
>>>> Requestors: matthew.sittel.ctr at offutt.af.mil
>>>> Status: new
>>>> Ticket <URL:
>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455
>>>>>
>>>>
>>>>
>>>> Here we go again with a new issue. I'm testing the tropical
cyclone tracker for an upcoming project we've been asked to do. So I
gathered some sample files from NHC's FTP site and tried to run the
code. I got the TC_DLAND to run just fine, but something appears to
prevent TC_PAIRS from running, and I can't figure out what it is.
>>>>
>>>> Here's the log file as I see it. My main concerns are why Model
Suffix is (nul), and why the ADECK file claims to be read but none of
the records are used. I am using the wwpts file included in 4.1; I
can't find a newer one on the NHC FTP site-I looked in many
directories there but didn't see this file. The storm I'm looking at
has no entries in the wwpts file, presumably because the wwpts file
predates the aal and bal files.
>>>>
>>>> DEBUG 1: [Source 1 of 1] ADECK Source:
>>>> /home/qcteam/METv4.1/data/aal142012.dat, Model Suffix: (nul)
DEBUG 1:
>>>> [Source 1 of 1] BDECK Source:
>>>> /home/qcteam/METv4.1/data/bal142012.dat,
>>>> Model Suffix: (nul) DEBUG 1: Config File Default:
>>>> /home/qcteam/METv4.1/data/config/TCPairsConfig_default
>>>> DEBUG 1: Config File User:
>>>> /home/qcteam/METv4.1/data/config/TCPairs.cfg
>>>> DEBUG 1: Distance to land file:
>>>> /home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc
>>>> DEBUG 4: Latitude/Longitude Grid Data:
>>>> DEBUG 4: lat_ll = 0
>>>> DEBUG 4: lon_ll = -180
>>>> DEBUG 4: delta_lat = 0.1
>>>> DEBUG 4: delta_lon = 0.1
>>>> DEBUG 4: Nlat = 601
>>>> DEBUG 4: Nlon = 1801
>>>> DEBUG 2: Processing 1 ADECK file(s).
>>>> DEBUG 4: [File 1 of 1] Used 0 of 202726 lines read from file
"/home/qcteam/METv4.1/data/aal142012.dat"
>>>> DEBUG 3: Used 0 of 202726 lines read from 1 file(s).
>>>> DEBUG 3: Identified 0 track(s).
>>>> DEBUG 2: Found 0 ADECK track(s).
>>>> DEBUG 2: Filtering 0 ADECK tracks based on config file settings.
>>>> DEBUG 3: Total tracks read = 0
>>>> DEBUG 3: Total tracks kept = 0
>>>> DEBUG 3: Rejected for model = 0
>>>> DEBUG 3: Rejected for storm id = 0
>>>> DEBUG 3: Rejected for basin = 0
>>>> DEBUG 3: Rejected for cyclone = 0
>>>> DEBUG 3: Rejected for storm name = 0 DEBUG 3: Rejected for init
>>>> time = 0 DEBUG 3: Rejected for init hour = 0 DEBUG 3: Rejected
for
>>>> valid time = 0 DEBUG 3: Rejected for init mask = 0 DEBUG 3:
>>>> Rejected for valid mask = 0 DEBUG 2: Processing 1 BDECK file(s).
>>>> DEBUG 4: Updating storm name from "INVEST" to "FOURTEEN" for
AL142012.
>>>> DEBUG 4: Updating storm name from "FOURTEEN" to "NADINE" for
AL142012.
>>>> DEBUG 4: [File 1 of 1] Used 184 of 184 lines read from file
"/home/qcteam/METv4.1/data/bal142012.dat"
>>>> DEBUG 3: Used 184 of 184 lines read from 1 file(s).
>>>> DEBUG 3: Identified 1 track(s).
>>>> DEBUG 4: [Track 1 of 1] TrackInfo: StormId = "AL142012", Basin =
>>>> "AL", Cyclone = "14", StormName = "NADINE", TechniqueNumber =
-9999, Technique = "BEST", Initials = "(nul)", InitTime = NA,
MinValidTime = 20120910_120000, MaxValidTime = 20121004_060000,
NPoints = 96, NAlloc = 100, NTrackLines = 0 DEBUG 2: Found 1 BDECK
track(s).
>>>> DEBUG 2: Deriving 12-hour interpolated tracks.
>>>> DEBUG 2: Finished deriving 0 12-hour interpolated tracks.
>>>> DEBUG 2: Deriving 0 ADECK consensus model(s).
>>>> DEBUG 2: Added 0 ADECK consensus tracks(s).
>>>> DEBUG 2: Deriving 0 ADECK lag model(s).
>>>> DEBUG 2: Added 0 ADECK lag tracks(s).
>>>> DEBUG 2: Deriving 0 CLIPER/SHIFOR baseline model(s).
>>>> DEBUG 2: Added 0 CLIPER/SHIFOR baseline track(s).
>>>> DEBUG 2: Matching 0 ADECK tracks to 1 BDECK tracks.
>>>> DEBUG 1: Watch/Warning file:
>>>> /home/qcteam/METv4.1/data/tc_data/wwpts_us.txt
>>>> DEBUG 1: Output file: ./out_tcmpr.tcst
>>>>
>>>> My config file contents:
>>>>
>>>> model = ["AEMN"];
>>>>
>>>> storm_id = ["AL142012"];
>>>>
>>>> basin = ["AL"];
>>>>
>>>> cyclone = ["14"];
>>>>
>>>> storm_name = ["NADINE"];
>>>>
>>>> init_beg = "";
>>>> init_end = "";
>>>> init_inc = [];
>>>> init_exc = [];
>>>>
>>>> valid_beg = "";
>>>> valid_end = "";
>>>>
>>>> init_hour = [];
>>>>
>>>> init_mask = "";
>>>> valid_mask = "";
>>>>
>>>> check_dup = FALSE;
>>>>
>>>> interp12 = TRUE;
>>>>
>>>> consensus = [];
>>>>
>>>> lag_time = [];
>>>>
>>>> best_baseline = [];
>>>> oper_baseline = [];
>>>>
>>>> match_points = FALSE;
>>>>
>>>> dland_file =
>>>> "/home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc";
>>>>
>>>> watch_warn = {
>>>> file_name =
"/home/qcteam/METv4.1/data/tc_data/wwpts_us.txt";
>>>> time_offset = -14400;
>>>> };
>>>>
>>>> version = "V4.1";
>>>>
>>>> Hopefully you can suggest something to get this working!
>>>>
>>>> Thanks,
>>>> Matt
>>>>
>>>> // SIGNED //
>>>> Matthew C. Sittel
>>>> Associate Scientist
>>>> University Corporation for Atmospheric Research 16WS/WXN, Offutt
>>>> AFB, NE
>>>> (402) 294-3473 DSN: 271-3473
>>>>
>>>
>>>
>>
>>
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
From: John Halley Gotway
Time: Thu Aug 01 13:34:43 2013
Matt,
I checked and you are right, there is no mention of differences in the
MET-TC documentation. This our first release of the
code/documentation, so I'm not surprised we overlooked. Please let us
know
of other instances of incomplete/unclear documentation. I'll add a
task for us to improve that for the next release.
The other option not documented is the absolute value. You can put
"abs()" around any of the columns or differences of columns:
METv4.1/bin/tc_stat -lookin tcmpr_al142012.tcst -job summary
-column ALAT-BLAT -column abs\(ALAT-BLAT\)
Thanks,
John
On 08/01/2013 12:57 PM, matthew.sittel.ctr at offutt.af.mil via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>
> Sure looks that way. I'm testing using a single ensemble member
instead as those do appear to have winds... I think our cases will.
>
> I didn't realize I could do differences... I didn't gather that from
the PDF... maybe I should re-read it!
>
> -----Original Message-----
> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Sent: Thursday, August 01, 2013 12:54 PM
> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>
> Matt,
>
> Yes, the "-column" option supports differences. For example...
>
> METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
summary -column AMAX_WIND-BMAX_WIND -by LEAD
>
> However when I ran this job, I just got all NA's in the output.
Upon closer inspection, I see that the AMAX_WIND column is all NA's in
the .tcst output of tc_pairs. Taking a step back, I see that the AEMN
track data in the ATCF file does not contain valid wind speed. Is
that perhaps a track only ensemble mean forecast?
>
> Thanks,
> John
>
>
> On 08/01/2013 11:12 AM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>>
>> That took care of it. Thanks.
>>
>> Is there a field for wind error? That's the other one I need to
build. Basically the wind difference "forecast minus best track
speed". I was thinking I could get at it, indirectly, by summarizing
AMAX_WIND and BMAX_WIND, but that would not give me pairwise
differences or any sort of stratification of errors as a function of
wind speed.
>>
>> -----Original Message-----
>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>> Sent: Thursday, August 01, 2013 11:27 AM
>> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
>> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>>
>> Matt,
>>
>> Yes, you're missing the "jobs" setting. You need to define what
type of statistics you'd like computed. This is just like the
stat_analysis tool - you can run a single job on the command line or
multiple jobs via a config file. Clearly though, I should add an
error message when a user passes in a config file with no jobs
defined. I'll just add that to the development version.
>>
>> While there are a few different job types defined for
stat_analysis, currently the only supported job types for tc_stat are
"filter" and "summary". Here's an example of how you might run a
summary job on the command line without a config file:
>>
>> METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
>> summary -column TK_ERR,ALTK_ERR,CRTK_ERR
>>
>> That'll extract the 3 columns (TK_ERR, ALTK_ERR, and CRTK_ERR) and
print summary information (min/max/mean/median, and so on).
>>
>> Next, try this one:
>> METv4.1_bugfix/bin/tc_stat -lookin tcmpr_al142012.tcst -job
>> summary -column TK_ERR -by LEAD
>>
>> This will summarize the track errors but will give you output for
each unique entry in the "LEAD" column. So you'll see track errors by
lead time. As you'd expect, the track errors grow by lead time.
>>
>> Hope that helps.
>>
>> Thanks,
>> John
>>
>>
>> On 08/01/2013 10:12 AM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>>
>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>>>
>>> It's nice to be doing some other stuff, just for variety's sake.
I'm working on a low-pressure tracking routine as well that will
evaluate model skill. It will have similarities to MET, and
ultimately I may use MET to do the aggregation of statistics... not
sure just yet, it's still in the 'proof of concept' stage.
>>>
>>> I can't get tc_stat to run... well, it runs but produces an empty
output file. Here's what I am running:
>>>
>>> tc_stat -lookin /home/qcteam/METv4.1/data/tc_data -config
>>> /home/qcteam/METv4.1/data/config/TCStat.cfg -v 4
>>>
>>> The "-lookin" directory has the .tcst output file I produced with
TC_PAIRS. All I get in response is this:
>>>
>>> DEBUG 1: [Source 1] Lookin: /home/qcteam/METv4.1/data/tc_data
DEBUG 1:
>>> Config File Default:
>>> /home/qcteam/METv4.1/data/config/TCStatConfig_default
>>> DEBUG 1: Config File User:
>>> /home/qcteam/METv4.1/data/config/TCStat.cfg
>>> DEBUG 1:
>>> DEBUG 1: Creating output file:
>>> /home/qcteam/METv4.1/data/tc_data/tc_stat_output.txt
>>>
>>> Here's what is in the config file I created:
>>>
>>> amodel = ["AEMN"];
>>> bmodel = ["BEST"];
>>> storm_id = ["AL142012"];
>>> basin = ["AL"];
>>> cyclone = ["14"];
>>> storm_name = [];
>>> init_beg = "";
>>> init_end = "";
>>> init_inc = [];
>>> init_exc = [];
>>> valid_beg = "";
>>> valid_end = "";
>>> valid_inc = [];
>>> valid_exc = [];
>>> lead = [];
>>> init_mask = [];
>>> valid_mask = [];
>>> line_type = [];
>>> track_watch_warn = [];
>>> column_thresh_name = [];
>>> column_thresh_val = [];
>>> column_str_name = [];
>>> column_str_val = [];
>>> init_thresh_name = [];
>>> init_thresh_val = [];
>>> init_str_name = [];
>>> init_str_val = [];
>>> water_only = FALSE;
>>> rapid_inten = FALSE;
>>> rapid_inten_thresh = >=30.0;
>>> landfall = FALSE;
>>> landfall_beg = -86400;
>>> landfall_end = 0;
>>> match_points = FALSE;
>>> event_equal = FALSE;
>>> out_init_mask = "";
>>> out_valid_mask = "";
>>> jobs = [];
>>> version = "V4.1";
>>>
>>> Do you see anything obvious that I missed?
>>>
>>> Note that the TC 4.1 documentation says "lookin" rather than "-
lookin" when listing the arguments for running TC_STAT.
>>>
>>> Thanks,
>>> Matt
>>>
>>> -----Original Message-----
>>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>>> Sent: Thursday, August 01, 2013 10:53 AM
>>> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
>>> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>>>
>>> Matt,
>>>
>>> Sounds interesting. I've actually switched gears recently too.
I'm working on the real-time evaluation for the HFIP products.
Currently, we're mostly just plotting track data, but the plan is to
start running MET-TC in near-realtime mode soon.
>>>
>>> I did find a copy-and-paste bug in the tc_pairs tool. I'll post a
bugfix for it shortly and write you back. You should still continue
to omit the storm name though. While the ATCF file format has a place
for the storm name, it usually is not included in the ADECK files. If
you list a storm name in the config file, it'll only keep tracks that
include it. Since ADECK files often don't have it, you won't keep
many of them.
>>>
>>> Also, don't worry about the model name suffix being "(nul)". We
put in this option of appending a string to the model names found in
the track files. At one point, we had multiple models to evaluate
that all made use of the same ATCF ID. To distinguish them in the
output, we added a suffix to their ATCF ID's to make them unique. If
you're not in that situation, you don't need to use the suffix.
>>>
>>> Thanks,
>>> John
>>>
>>> On 08/01/2013 09:43 AM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>>>
>>>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455 >
>>>>
>>>> Hey John,
>>>>
>>>> Yeah, that change took care of it.
>>>>
>>>> Our ensembles guys are going to make ATCF-formatted tracks for
their WRF members, the ensemble mean, and then compare them to a
deterministic WRF model... all to be compared to best tracks for
storms. TS Gil is in the east Pacific so they've got something to
work with. This will be our chance to utilize MET (which is in my
statement of work, to be the "subject matter expert" in house).
>>>>
>>>> Thanks,
>>>> Matt
>>>>
>>>> -----Original Message-----
>>>> From: John Halley Gotway via RT [mailto:met_help at ucar.edu]
>>>> Sent: Thursday, August 01, 2013 10:37 AM
>>>> To: Sittel, Matthew C CTR USAF AFWA 16 WS/WXN
>>>> Subject: Re: [rt.rap.ucar.edu #62455] Problem with TC_PAIRS
>>>>
>>>> Hello Matt,
>>>>
>>>> Try removing the storm name from the config file, like this:
>>>> storm_name = [];
>>>>
>>>> When you rerun, you should get matching tracks.
>>>>
>>>> However, I think there may be a bug in the config file parsing
code, and I'll let you know what I find.
>>>>
>>>> Thanks,
>>>> John
>>>>
>>>>
>>>> On 07/31/2013 03:40 PM, matthew.sittel.ctr at offutt.af.mil via RT
wrote:
>>>>>
>>>>> Wed Jul 31 15:40:07 2013: Request 62455 was acted upon.
>>>>> Transaction: Ticket created by matthew.sittel.ctr at offutt.af.mil
>>>>> Queue: met_help
>>>>> Subject: Problem with TC_PAIRS
>>>>> Owner: Nobody
>>>>> Requestors: matthew.sittel.ctr at offutt.af.mil
>>>>> Status: new
>>>>> Ticket <URL:
>>>>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62455
>>>>>>
>>>>>
>>>>>
>>>>> Here we go again with a new issue. I'm testing the tropical
cyclone tracker for an upcoming project we've been asked to do. So I
gathered some sample files from NHC's FTP site and tried to run the
code. I got the TC_DLAND to run just fine, but something appears to
prevent TC_PAIRS from running, and I can't figure out what it is.
>>>>>
>>>>> Here's the log file as I see it. My main concerns are why Model
Suffix is (nul), and why the ADECK file claims to be read but none of
the records are used. I am using the wwpts file included in 4.1; I
can't find a newer one on the NHC FTP site-I looked in many
directories there but didn't see this file. The storm I'm looking at
has no entries in the wwpts file, presumably because the wwpts file
predates the aal and bal files.
>>>>>
>>>>> DEBUG 1: [Source 1 of 1] ADECK Source:
>>>>> /home/qcteam/METv4.1/data/aal142012.dat, Model Suffix: (nul)
DEBUG 1:
>>>>> [Source 1 of 1] BDECK Source:
>>>>> /home/qcteam/METv4.1/data/bal142012.dat,
>>>>> Model Suffix: (nul) DEBUG 1: Config File Default:
>>>>> /home/qcteam/METv4.1/data/config/TCPairsConfig_default
>>>>> DEBUG 1: Config File User:
>>>>> /home/qcteam/METv4.1/data/config/TCPairs.cfg
>>>>> DEBUG 1: Distance to land file:
>>>>> /home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc
>>>>> DEBUG 4: Latitude/Longitude Grid Data:
>>>>> DEBUG 4: lat_ll = 0
>>>>> DEBUG 4: lon_ll = -180
>>>>> DEBUG 4: delta_lat = 0.1
>>>>> DEBUG 4: delta_lon = 0.1
>>>>> DEBUG 4: Nlat = 601
>>>>> DEBUG 4: Nlon = 1801
>>>>> DEBUG 2: Processing 1 ADECK file(s).
>>>>> DEBUG 4: [File 1 of 1] Used 0 of 202726 lines read from file
"/home/qcteam/METv4.1/data/aal142012.dat"
>>>>> DEBUG 3: Used 0 of 202726 lines read from 1 file(s).
>>>>> DEBUG 3: Identified 0 track(s).
>>>>> DEBUG 2: Found 0 ADECK track(s).
>>>>> DEBUG 2: Filtering 0 ADECK tracks based on config file settings.
>>>>> DEBUG 3: Total tracks read = 0
>>>>> DEBUG 3: Total tracks kept = 0
>>>>> DEBUG 3: Rejected for model = 0
>>>>> DEBUG 3: Rejected for storm id = 0
>>>>> DEBUG 3: Rejected for basin = 0
>>>>> DEBUG 3: Rejected for cyclone = 0
>>>>> DEBUG 3: Rejected for storm name = 0 DEBUG 3: Rejected for init
>>>>> time = 0 DEBUG 3: Rejected for init hour = 0 DEBUG 3: Rejected
for
>>>>> valid time = 0 DEBUG 3: Rejected for init mask = 0 DEBUG 3:
>>>>> Rejected for valid mask = 0 DEBUG 2: Processing 1 BDECK file(s).
>>>>> DEBUG 4: Updating storm name from "INVEST" to "FOURTEEN" for
AL142012.
>>>>> DEBUG 4: Updating storm name from "FOURTEEN" to "NADINE" for
AL142012.
>>>>> DEBUG 4: [File 1 of 1] Used 184 of 184 lines read from file
"/home/qcteam/METv4.1/data/bal142012.dat"
>>>>> DEBUG 3: Used 184 of 184 lines read from 1 file(s).
>>>>> DEBUG 3: Identified 1 track(s).
>>>>> DEBUG 4: [Track 1 of 1] TrackInfo: StormId = "AL142012", Basin =
>>>>> "AL", Cyclone = "14", StormName = "NADINE", TechniqueNumber =
-9999, Technique = "BEST", Initials = "(nul)", InitTime = NA,
MinValidTime = 20120910_120000, MaxValidTime = 20121004_060000,
NPoints = 96, NAlloc = 100, NTrackLines = 0 DEBUG 2: Found 1 BDECK
track(s).
>>>>> DEBUG 2: Deriving 12-hour interpolated tracks.
>>>>> DEBUG 2: Finished deriving 0 12-hour interpolated tracks.
>>>>> DEBUG 2: Deriving 0 ADECK consensus model(s).
>>>>> DEBUG 2: Added 0 ADECK consensus tracks(s).
>>>>> DEBUG 2: Deriving 0 ADECK lag model(s).
>>>>> DEBUG 2: Added 0 ADECK lag tracks(s).
>>>>> DEBUG 2: Deriving 0 CLIPER/SHIFOR baseline model(s).
>>>>> DEBUG 2: Added 0 CLIPER/SHIFOR baseline track(s).
>>>>> DEBUG 2: Matching 0 ADECK tracks to 1 BDECK tracks.
>>>>> DEBUG 1: Watch/Warning file:
>>>>> /home/qcteam/METv4.1/data/tc_data/wwpts_us.txt
>>>>> DEBUG 1: Output file: ./out_tcmpr.tcst
>>>>>
>>>>> My config file contents:
>>>>>
>>>>> model = ["AEMN"];
>>>>>
>>>>> storm_id = ["AL142012"];
>>>>>
>>>>> basin = ["AL"];
>>>>>
>>>>> cyclone = ["14"];
>>>>>
>>>>> storm_name = ["NADINE"];
>>>>>
>>>>> init_beg = "";
>>>>> init_end = "";
>>>>> init_inc = [];
>>>>> init_exc = [];
>>>>>
>>>>> valid_beg = "";
>>>>> valid_end = "";
>>>>>
>>>>> init_hour = [];
>>>>>
>>>>> init_mask = "";
>>>>> valid_mask = "";
>>>>>
>>>>> check_dup = FALSE;
>>>>>
>>>>> interp12 = TRUE;
>>>>>
>>>>> consensus = [];
>>>>>
>>>>> lag_time = [];
>>>>>
>>>>> best_baseline = [];
>>>>> oper_baseline = [];
>>>>>
>>>>> match_points = FALSE;
>>>>>
>>>>> dland_file =
>>>>>
"/home/qcteam/METv4.1/data/tc_data/dland_nw_hem_tenth_degree.nc";
>>>>>
>>>>> watch_warn = {
>>>>> file_name =
"/home/qcteam/METv4.1/data/tc_data/wwpts_us.txt";
>>>>> time_offset = -14400;
>>>>> };
>>>>>
>>>>> version = "V4.1";
>>>>>
>>>>> Hopefully you can suggest something to get this working!
>>>>>
>>>>> Thanks,
>>>>> Matt
>>>>>
>>>>> // SIGNED //
>>>>> Matthew C. Sittel
>>>>> Associate Scientist
>>>>> University Corporation for Atmospheric Research 16WS/WXN, Offutt
>>>>> AFB, NE
>>>>> (402) 294-3473 DSN: 271-3473
>>>>>
>>>>
>>>>
>>>
>>>
>>
>>
>
>
------------------------------------------------
More information about the Met_help
mailing list