From Code to Community: Sponsoring The Perl and Raku Conference 2025 Learn more

PUCKERING / Weather-GHCN-Fetch-v0.0.011 / bin / ghcn_extremes 

            

              
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
—
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156

              
#!/usr/bin/perl
# ghcn_extremes - Report temperature extremes from ghcn_fetch output
use v5.18;  # minimum for Object::Pad
our $VERSION = 'v0.0.011';
use Weather::GHCN::App::Extremes;
Weather::GHCN::App::Extremes->run( \@ARGV );
########################################################################
# Documentation
########################################################################
__END__
=head1 NAME
ghcn_extremes - Report temperature extremes from ghcn_fetch output
=head1 VERSION
version v0.0.011
=head1 SYNOPSIS
    ghcn_extremes [-limit <int>] [-ndays <int>]
                  [-peryear] [-daycounts] [-cold] [-nogaps]
                  [-outclip] [ file... ]
    ghcn_extremes [--help | -? | --usage]
=head1 DESCRIPTION
Report patterns of temperature extremes (heatwaves or coldwaves) by
analyzing daily temperature records and looking for consecutive days
of extreme temperatures.  By default the script reports a count of
heatwaves in each year of the input data, where a heatwave is defined
as 5 or more consecutive days where the maximum temperature is 30
Celsius or higher.
Using the -cold option, it can report coldwaves; by default, where
the temperatures are at or below -20 C for 5 days or more days.
When the -daycounts option is given, the script reports the
year-month-day the heatwave (or coldwave) began along with a count of
the number of consecutive days it lasted.
The input data must be tab-separated, in date order, and contain the
following fields:
    year, month, day, decade,     # decade is 10 * int( year % 100 / 10 )
    s_decade, s_year, s_qtr,      # seasonal period (Jul-Jun years)
    tmax, tmin, tavg, qflags,     # qflags are ignored, so may be empty
    stn_id, station_name
All trailing columns of data are ignored.  This format is the one
generated by ghcn_fetch -report detail; i.e. it is a report
of daily weather data for each station id.
In typical use, you would extract the data of interest with
    ghcn_fetch -report detail <other_options>
and then either pipe it to ghcn_extremes or save the output to a file
and then use that file as input to ghcn_extremes.
=head1 PARAMETERS
Getoptions::Long is used, so either - or -- may be used.  Parameter
names may be abbreviated, so long as they remains unambiguous.  Flag
options may appear after filenames.
=over 4
=item -limit <int>
Sets the limit (in degrees Celsius) at or above which the TMAX temperature
for the day is counted as part of a heatwave.  Defaults to 30 Celsius.
When -cold is given, the limit is the temperature at or below which the
TMIN temperature for the day is counted as part of a coldwave.  Defaults
to -20 Celsius.
=item -ndays <int>
Consider any consecutive run of <int> or more days where the temperature
exceeds that specified by -limit as a heatwave (or coldwave).
=item -peryear
Report number of heatwaves per year.
=item -daycounts
Report the number of days in each heatwave (or coldwave, rather than
the number of heatwaves (or coldwaves) per year.
=item -cold
Report coldwaves instead of heatwaves.  Days with temperatures less
than or equal to the value specified by -limit will be counted.  The
default for -limit is set to -20 Celsius when this option is selected.
=item -nogaps
When a year has no heat (or cold) waves, there will be no data for
that year in the output.  Thus, when charting, those years are missing
from the x-axis, which not only distorts the picture, it can lead
to inaccurate trend analysis.
When -nogaps is specified, the script will fill in these gaps by
adding rows at the end with the missing years.  It will be necessary
to manually sort the table in Excel, or specify an ascending sort
order for the year axis, in order for charts to look correct.
=item -outclip
Send output to the Windows clipboard.
=item -h | -help
Display this documentation.
=back
=head1 INPUT
The input data is expected to be tab-separated, with columns year,
month, day, then four other columns, then TMAX and TMIN columns.  TMAX
values are used for determining heatwaves.  TMIN values are used for
coldwaves.
The first row of data can be a row header, but column 1 must have
the value 'Year' and columns 6 and 7 must begin with 'TMAX' and 'TMIN'
respectively.
=head1 RELATED SCRIPTS
This script was designed to use the output from ghcn_fetch      ,
specifically when the option "-report daily" is used.  It generates
the correct level and format of data for input to this script.
The option "-report detail" can also be used, but if the data contains
multiple stations then this script will miscount.
=head1 AUTHOR
Gary Puckering (jgpuckering@rogers.com)
=head1 LICENSE AND COPYRIGHT
Copyright 2022, Gary Puckering
=cut