Skip to content

Tempterc #613

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
remicousin merged 2 commits into from
Apr 20, 2026
Merged

Tempterc #613

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
ceef48f
config for temp
remicousin Mar 31, 2026
5d468b9
temperature tercile fcst mr
remicousin Apr 16, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
6 changes: 6 additions & 0 deletions ccsr/config-defaults.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -30,4 +30,10 @@ maprooms:
variable: Precipitation
var: precip
issue_format: "%b %Y"
- title: terc_fcst_temp
core_path: terc_fcst_temp
forecast_path: /Data/data21/iri/FD/NMME_Seasonal_Forecast/Temperature_ELR
variable: Temperature
var: tref
issue_format: "%b %Y"

6 changes: 6 additions & 0 deletions ccsr/config-iri.yaml
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -43,3 +43,9 @@ maprooms:
variable: Precipitation
var: precip
issue_format: "%b %Y"
- title: terc_fcst_temp
core_path: terc_fcst_temp
forecast_path: /Data/data21/iri/FD/NMME_Seasonal_Forecast/Temperature_ELR
variable: Temperature
var: tref
issue_format: "%b %Y"
20 changes: 10 additions & 10 deletions ccsr/terc_fcst/layout.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -8,7 +8,7 @@
IRI_GRAY = "rgb(113,112,116)"
LIGHT_GRAY = "#eeeeee"

def app_layout():
def app_layout(config):

return lou.app_1(

Expand All @@ -19,18 +19,18 @@ def app_layout():
),

lou.description(
"Precipitation Terciles Seasonal Forecast",
"""
The seasonal forecast for above-, below- and near-normal precipitation
f"{config['variable']} Terciles Seasonal Forecast",
f"""
The seasonal forecast for above-, below- and near-normal {config['variable']}
from the IRI.
""",
html.P(
"""
The default map shows globally the seasonal precipitation forecast
tercile probability. The historical climatology used is 1982-2010 up
to the forecast issued in August 2021, and is 1991-2020 from the
forecast issued in September 2021. The forecast shown is the latest
forecast made (e.g. Dec 2017) for the next season to come
f"""
The default map shows globally the seasonal {config['variable']}
forecast tercile probability. The historical climatology used is
1982-2010 up to the forecast issued in August 2021, and is 1991-2020
from the forecast issued in September 2021. The forecast shown is
the latest forecast made (e.g. Dec 2017) for the next season to come
(e.g. Jan-Mar 2018). Four different seasons are forecasted and it is
also possible to consult forecasts made previously. The forecasts
are directly computed from the extended logistic regression model as
Expand Down
35 changes: 27 additions & 8 deletions ccsr/terc_fcst/maproom.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -35,7 +35,7 @@ def register(FLASK, config):
)
APP.title = "Tercile Forecast"

APP.layout = layout.app_layout()
APP.layout = layout.app_layout(config)

@APP.callback(
Output("start_div", "children"),
Expand All @@ -60,13 +60,18 @@ def initialize(path):
options=[ld["value"] for ld in target_dict],
labels=[ld["label"] for ld in target_dict],
)
if config["variable"] == "Precipitation":
cs_below = "prcp_terciles_below"
cs_above = "prcp_terciles_above"
else:
cs_below = "temp_terciles_below"
cs_above = "temp_terciles_above"
return (
issue_select, lead_select,
# while the data is classified and mapped according to classification
# colormap, I draw the colorscale from 2 colormaps (for below and above)
# corresponding to those classifications and with the probability values
CMAPS["prcp_terciles_below"].to_dash_leaflet(),
CMAPS["prcp_terciles_above"].to_dash_leaflet(),
CMAPS[cs_below].to_dash_leaflet(), CMAPS[cs_above].to_dash_leaflet(),
) + mru.initialize_map(fcst_ds)


Expand Down Expand Up @@ -141,8 +146,18 @@ def local_plots(marker_pos, start_date, lead_time, targets):
lng_units = "˚E" if (fcst_ds['X'] >= 0) else "˚W"
lat_units = "˚N" if (fcst_ds['Y'] >= 0) else "˚S"
#Get colors from the colormap to assign each bar
if config["variable"] == "Precipitation":
cs_below = "prcp_terciles_below"
cs_above = "prcp_terciles_above"
low_below = "lightyellow"
low_above = "rgb(220, 255, 220)"
else:
cs_below = "temp_terciles_below"
cs_above = "temp_terciles_above"
low_below = "azure"
low_above = "lightpink"
if fcst_ds["prob"].isel(cat=0) >= 37.5:
below_color = CMAPS["prcp_terciles_below"].to_rgba_array()[
below_color = CMAPS[cs_below].to_rgba_array()[
int(
0.5
+ 255
Expand All @@ -155,9 +170,9 @@ def local_plots(marker_pos, start_date, lead_time, targets):
f"rgb({below_color[0]}, {below_color[1]}, {below_color[2]})"
)
else:
below_color = "lightyellow"
below_color = low_below
if fcst_ds["prob"].isel(cat=2) >= 37.5:
above_color = CMAPS["prcp_terciles_above"].to_rgba_array()[
above_color = CMAPS[cs_above].to_rgba_array()[
int(
0.5
+ 255
Expand All @@ -170,7 +185,7 @@ def local_plots(marker_pos, start_date, lead_time, targets):
f"rgb({above_color[0]}, {above_color[1]}, {above_color[2]})"
)
else:
above_color = "rgb(220, 255, 220)"
above_color = low_above
local_graph = px.bar(
(
fcst_ds["prob"]
Expand Down Expand Up @@ -254,7 +269,11 @@ def fcst_tiles(tz, tx, ty, start_date, lead_time):
],
[lambda x : x, lambda x : x - 4, 6],
)
dominant_fcst_class.attrs["colormap"] = CMAPS["prcp_terciles"]
if config["variable"] == "Precipitation":
cs_terc = "prcp_terciles"
else:
cs_terc = "temp_terciles"
dominant_fcst_class.attrs["colormap"] = CMAPS[cs_terc]
dominant_fcst_class.attrs["scale_min"] = 0
dominant_fcst_class.attrs["scale_max"] = 11
return pingrid.tile(
Expand Down
48 changes: 48 additions & 0 deletions pingrid/impl.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -880,6 +880,51 @@ def produce_shape_tile(
[37.5, 42.5, 42.5, 47.5, 47.7, 57.5, 57.5, 67.5, 67.5, 100],
)


_TEMP_TERCILES_CS = ColorScale(
"temp_terciles",
[
WHITE, WHITE,
Color(200, 240, 250), Color(200, 240, 250),
Color(177, 202, 250), Color(177, 202, 250),
Color(153, 163, 250), Color(153, 163, 250),
Color(129, 124, 250), Color(129, 124, 250),
Color(85, 85, 200), Color(85, 85, 200),
GRAY, GRAY,
Color(255, 210, 210), Color(255, 210, 210),
Color(242, 170, 170), Color(242, 170, 170),
Color(250, 120, 120), Color(250, 120, 120),
Color(255, 40, 40), Color(255, 40, 40),
Color(150, 19, 19), Color(150, 19, 18),
],
[0, 0.5, 0.5, 1.5, 1.5, 2.5, 2.5, 3.5, 3.5, 4.5, 4.5, 5.5,
5.5, 6.5, 6.5, 7.5, 7.5, 8.5, 8.5, 9.5, 9.5, 10.5, 10.5, 11],
)

_TEMP_TERCILES_BELOW_CS = ColorScale(
"temp_terciles_below",
[
Color(200, 240, 250), Color(200, 240, 250),
Color(177, 202, 250), Color(177, 202, 250),
Color(153, 163, 250), Color(153, 163, 250),
Color(129, 124, 250), Color(129, 124, 250),
Color(85, 85, 200), Color(85, 85, 200),
],
[37.5, 42.5, 42.5, 47.5, 47.7, 57.5, 57.5, 67.5, 67.5, 100],
)

_TEMP_TERCILES_ABOVE_CS = ColorScale(
"temp_terciles_above",
[
Color(255, 210, 210), Color(255, 210, 210),
Color(242, 170, 170), Color(242, 170, 170),
Color(250, 120, 120), Color(250, 120, 120),
Color(255, 40, 40), Color(255, 40, 40),
Color(150, 19, 19), Color(150, 19, 18),
],
[37.5, 42.5, 42.5, 47.5, 47.7, 57.5, 57.5, 67.5, 67.5, 100],
)

_RAIN_POE_CS = ColorScale(
"rain_poe",
[BLACK, BROWN, ORANGE, YELLOW, MOCCASIN, MOCCASIN, LIMEGREEN, TURQUOISE, BLUE, PURPLE],
Expand Down Expand Up @@ -1034,6 +1079,9 @@ def produce_shape_tile(
_STD_ANOMALY_CS,
_TEMP_CS,
_TEMP_ANOMALY_CS,
_TEMP_TERCILES_CS,
_TEMP_TERCILES_ABOVE_CS,
_TEMP_TERCILES_BELOW_CS,
_VULN_CS,
]}

Expand Down
Loading