Model

The Model class provides access to weather forecast data from Jua's weather prediction models. It allows you to retrieve future weather predictions with flexible spatial and temporal selection options.

Historical weather data (known as hindcasts) are also available directly directly through the Model class. Hindcasts are past model runs that allow you to analyze historical weather patterns, evaluate model performance, or test your pipelines. All model runs that are older that 90 days (i.e. an having an init_time that is more than 90 days in the past) are classified as hindcasts.

Information about model/query pricing is available in the pricing page.

Overview

The Model class serves as your interface to Jua's weather forecasting capabilities. It provides methods to:

• Retrieve forecast data for global coverage or specific locations

• Access forecasts for specific initialization times or the latest available forecast

• Query metadata about available forecasts

• Check forecast availability and readiness

List Available Models

Code

from jua.weather import Models

print("Available models:")
for model in Models:
    print(f"  - {model}")

Output

Request model metadata

Code

import json

from jua import JuaClient
from jua.weather import Models

client = JuaClient()
model = client.weather.get_model(Models.EPT1_5)
metadata = model.get_metadata()
print(json.dumps(metadata.model_dump(), indent=4, default=str))

Output

Getting the Forecasts Available for a Model

Code

from jua import JuaClient
from jua.weather import Models

client = JuaClient()
model = client.weather.get_model(Models.EPT2)

# The available forecasts will be returned in descending order,
#  from most recent to oldest; this returns the 10 most recent forecasts
available_forecasts = model.get_available_forecasts(limit=10)
for forecast in available_forecasts:
    print(f" - {forecast.init_time}")

Output

The get_available_forecasts forecast response is paginated. So to obtain all available forecasts for a model, you need to iterate through the responses.

Code

from jua import JuaClient
from jua.weather import Models

client = JuaClient()
model = client.weather.get_model(Models.EPT2)

# Paginate through the response to get all available init_times
result = model.get_available_forecasts(limit=1000)
available_forecasts = list(result.forecasts)
while result.has_more:
    result = result.next()
    available_forecasts.extend(result.forecasts)

print(f"Available init times for {model}: {len(available_forecasts)}")

Output

You can also specify a range of dates for which you want to get available init_times.

Code

from datetime import datetime

from jua import JuaClient
from jua.weather import Models

client = JuaClient()
model = client.weather.get_model(Models.EPT2)

# get the available hindcasts from EPT-2 for January 2025
available_forecasts = model.get_available_forecasts(
    since=datetime(2025, 1, 1),
    before=datetime(2025, 2, 1),
    limit=200,
)
for forecast in available_forecasts:
    print(f" - {forecast.init_time}")

Output

Check if the forecasted data is available

Code

from jua import JuaClient
from jua.weather import Models

client = JuaClient()
model = client.weather.get_model(Models.EPT1_5)

# init_time might be any valid datetime object or datetime string
if model.is_ready(forecasted_hours=480, init_time="latest"):
    print(f"{model.name} has at least 480 hours of forecast data available")
else:
    print(f"{model.name} does not yet have 480 hours of forecast data available")

Output

Requesting the forecast data

Point-Based

from jua import JuaClient
from jua.types.geo import LatLon
from jua.weather import Models, Variables

# Setup JuaClient
client = JuaClient()

# Choose the model
model = client.weather.get_model(Models.EPT1_5)

# Defining the sites that we are interested in
sites = [
    # label is optional
    LatLon(55.06, 13.00, label="Kriegers Flak"),
    LatLon(54.04, 5.96, label="Gemini Wind Farm"),
    LatLon(51.71, 2.91, label="Borssele III & IV"),
    LatLon(53.89, 1.79, label="Hornsea Two"),
]

# Define the variables we are interested in
# If set to `None` all available variables will be loaded
variables = [
    Variables.WIND_SPEED_AT_HEIGHT_LEVEL_100M,
]

site_forecasts = model.get_forecasts(
    init_time="latest", # or any valid datetime / datetime string
    points=sites,
    # The below parameters are optional
    variables=variables,
    min_lead_time=0, # hours
    max_lead_time=72, # hours
    method="bilinear",
)

# Convert the response to an xarray dataset
print(site_forecasts.to_xarray()) 

<xarray.Dataset> Size: 2kB
Dimensions:                          (points: 4, init_time: 1,
                                      prediction_timedelta: 73)
Coordinates:
  * points                           (points) object 32B 'borssele_iii_&_iv' ...
  * init_time                        (init_time) datetime64[ns] 8B 2025-10-30...
  * prediction_timedelta             (prediction_timedelta) timedelta64[ns] 584B ...
    latitude                         (points) float32 16B 51.71 54.04 ... 55.06
    longitude                        (points) float32 16B 2.91 5.96 1.79 13.0
    requested_lat                    (points) float64 32B 51.71 54.04 ... 55.06
    requested_lon                    (points) float64 32B 2.91 5.96 1.79 13.0
Data variables:
    wind_speed_at_height_level_100m  (points, init_time, prediction_timedelta) float32 1kB ...

⚠️ When querying data with method="nearest" the returned coordinates will be slightly different to the requested ones, as it will return the coordinates of the nearest grid points for the model.

Region-Based

Requesting Hindcasts for Models

Hindcasts can be obtained directly through the same get_forecasts method.

Point-Based

from datetime import datetime
from jua import JuaClient
from jua.types.geo import LatLon
from jua.weather import Models, Variables

# Setup JuaClient
client = JuaClient()

# Choose the model
model = client.weather.get_model(Models.EPT2)

# Defining the sites that we are interested in
sites = [
    # label is optional
    LatLon(55.06, 13.00, label="Kriegers Flak"),
    LatLon(54.04, 5.96, label="Gemini Wind Farm"),
    LatLon(51.71, 2.91, label="Borssele III & IV"),
    LatLon(53.89, 1.79, label="Hornsea Two"),
]

# Define the variables we are interested in
# If set to `None` all available variables will be loaded
variables = [
    Variables.WIND_SPEED_AT_HEIGHT_LEVEL_10M,
]

site_forecasts = model.get_forecasts(
    init_time=slice(
        datetime(2025, month=1, day=1),
        datetime(2025, month=3, day=31, hour=23, minute=59),
    ),
    points=sites,
    # The below parameters are optional
    variables=variables,
    min_lead_time=0, # hours
    max_lead_time=72, # hours
    method="bilinear",
)

# Convert the response to an xarray dataset
print(site_forecasts.to_xarray()) 

<xarray.Dataset> Size: 424kB
Dimensions:                         (points: 4, init_time: 360,
                                     prediction_timedelta: 73)
Coordinates:
  * points                          (points) object 32B 'borssele_iii_&_iv' ....
  * init_time                       (init_time) datetime64[ns] 3kB 2025-01-01...
  * prediction_timedelta            (prediction_timedelta) timedelta64[ns] 584B ...
    latitude                        (points) float32 16B 51.71 54.04 53.89 55.06
    longitude                       (points) float32 16B 2.91 5.96 1.79 13.0
    requested_lat                   (points) float64 32B 51.71 54.04 53.89 55.06
    requested_lon                   (points) float64 32B 2.91 5.96 1.79 13.0
Data variables:
    wind_speed_at_height_level_10m  (points, init_time, prediction_timedelta) float32 420kB ...

⚠️ When querying data with method="nearest" the returned coordinates will be slightly different to the requested ones, as it will return the coordinates of the nearest grid points for the model.

Region-Based