Flight plan - HALO-20240903a

ec_under ec_track c_north c_mid c_south c_atr

Crew

The flight is planned to take off at 2024-09-03 11:04:13+00:00.

Job	Name
PI	Bjorn Stevens
WALES	Tanya Bodenbach
HAMP	Clara Bayley
Dropsondes	Nina Robbins
Smart/VELOX	Sophie Rosenburg
SpecMACS	Zekican Demiralay
Flight Documentation	Marius Rixen
Ground contact	Karsten Peters

Flight plan

Show code cell source

Hide code cell source

from datetime import datetime
import orcestra.sat
from orcestra.flightplan import bco, sal, point_on_track, LatLon, IntoCircle, FlightPlan

# Some fixed coordinates

lon_min, lon_max, lat_min, lat_max = -65, -5, -5, 25

band    = "east"
airport = sal if band == "east" else bco

radius = 72e3*1.852

atr_radius = 38e3*1.852
c_atr_nw = LatLon(17.433,-23.500, label = "c_atr")
c_atr_se = LatLon(16.080,-21.715, label = "c_atr")

# Define dates for flight

flight_time   = datetime(2024, 9, 3, 12, 0, 0)
flight_index = f"HALO-{flight_time.strftime('%Y%m%d')}a"

# Load satellite tracks 

ec_fcst_time  = "2024-09-02"
ec_track = orcestra.sat.SattrackLoader("EARTHCARE", ec_fcst_time, kind="PRE") \
    .get_track_for_day(f"{flight_time:%Y-%m-%d}")\
    .sel(time=slice(f"{flight_time:%Y-%m-%d} 14:00", None))

# Create elements of track

ec_north = point_on_track(ec_track,lat= 14.75).assign(label = "ec_north") 
ec_south = point_on_track(ec_track,lat=  2.00).assign(label = "ec_south")
c_north  = point_on_track(ec_track,lat= 13.00).assign(label = "c_north")
c_south  = point_on_track(ec_track,lat=  4.00).assign(label = "c_south")
c_mid    = point_on_track(ec_track,lat=c_south.towards(c_north).lat).assign(label = "c_mid")
c_atr    = LatLon(ec_north.lat+atr_radius/1852./60.,-22.1, label = "c_atr")
j_beg    = LatLon(lat=ec_north.lat,lon=-25.).assign(label = "junction")
ec_under = point_on_track(ec_track, lat=c_mid.towards(c_north,fraction=2./5.).lat, with_time=True).assign(label = "ec_under", note = "meet EarthCARE, drop sonde")

# Additional Waypoints

x_cnn = point_on_track(ec_track,lat= 13.67).assign(label="x_cnn", note="Northern alternative center for c_north")
x_cns = point_on_track(ec_track,lat= 11.50).assign(label="x_cns", note="Southern alternative center for c_north")
x_cms = point_on_track(ec_track,lat=  7.75).assign(label="x_cms", note="Southern alternative center for c_mid")
x_xsd = c_south.towards(c_mid).assign(label="x_sonde", note="Position for an extra sonde heading south")
j_atr = c_atr.towards(sal,distance=-atr_radius).assign(label="j_atr", note = "ATR circle start and end")
extra_waypoints = [LatLon(17.5,-24.00).assign(label="x_atr", note="alternative ATR circle")]

# Define Flight Paths

waypoints = [
     airport.assign(fl=0),
     j_beg.assign(fl=390),
     ec_north.assign(fl=410),
     c_mid.assign(fl=410),
     x_xsd.assign(fl=410),
     IntoCircle(c_south.assign(fl=410), radius, -360),
     x_cms.assign(fl=450).assign(note='might be to soon to rise to FL450, if so rise to FL430 and higher later'),
     IntoCircle(c_mid.assign(fl=450), radius, -360),
     ec_under.assign(fl=450),
     x_cns.assign(fl=450),
     IntoCircle(c_north.assign(fl=450), radius, -360,enter=90),
     x_cnn.assign(fl=450),
     ec_north.assign(fl=450),
     j_beg.assign(fl=450),
     j_atr.assign(fl=350),
     IntoCircle(c_atr.assign(fl=350), atr_radius,-360),
     airport.assign(fl=0),
     ]

# Crew

crew = {'Mission PI': 'Bjorn Stevens',
        'DropSondes': 'Nina Robbins',
        'HAMP': 'Clara Bayley',
        'SMART/VELOX': 'Sophie Rosenburg',
        'SpecMACS': 'Zekican Demiralay',
        'WALES' : 'Tanja Bodenbach',
        'Flight Documentation': 'Marius Rixen',
        'Ground Support': 'Karsten Peters',
        }

# Plan

plan = FlightPlan(waypoints, flight_index, extra_waypoints=extra_waypoints, crew=crew)
print(f"Flight index: {plan.flight_id}")
print(f"Take-off: {plan.takeoff_time:%H:%M %Z}\nLanding:  {plan.landing_time:%H:%M %Z}\n")

/home/runner/miniconda3/envs/orcestra_book/lib/python3.12/site-packages/orcestra/sat.py:195: UserWarning: You are using an old forecast (issued on 2024-09-02) for EARTHCARE on 2024-09-03! The newest forecast issued so far was issued on 2024-09-03. It's a RES forecast.
  warnings.warn(

Flight index: HALO-20240903a
Take-off: 11:04 UTC
Landing:  19:44 UTC

/home/runner/miniconda3/envs/orcestra_book/lib/python3.12/site-packages/orcestra/flightplan.py:214: FutureWarning: The aircraft attribute of this FlightPlan is not set, but it's required for aircraft performance calculations. Currently, the code will proceed using the default aircraft performance, but the code may break in future.
  warn(

Show code cell source

Hide code cell source

import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import easygems.healpix as egh
import intake
from orcestra.flightplan import plot_cwv, plot_path

plt.figure(figsize = (14, 8))
ax = plt.axes(projection=ccrs.PlateCarree())
ax.set_extent([lon_min, lon_max, lat_min, lat_max], crs=ccrs.PlateCarree())
ax.coastlines(alpha=1.0)
ax.gridlines(draw_labels=True, dms=True, x_inline=False, y_inline=False, alpha = 0.25)

ifs_fcst_time = datetime(2024, 9, 1, 0, 0, 0)
ds = intake.open_catalog("https://tcodata.mpimet.mpg.de/internal.yaml").HIFS(datetime = ifs_fcst_time).to_dask().pipe(egh.attach_coords)
cwv_flight_time = ds["tcwv"].sel(time=flight_time, method = "nearest")
plot_cwv(cwv_flight_time, levels = [50.0, 60.0], ax=ax)
plt.title(f"{flight_time}\n(CWV forecast issued on {ifs_fcst_time})")

plt.plot(ec_track.lon, ec_track.lat, c='k', ls='dotted')
plot_path(plan, ax, color="C1")

Show code cell source

Hide code cell source

from orcestra.flightplan import vertical_preview

vertical_preview(waypoints)

Show code cell source

Hide code cell source

plan.show_details()
plan.export()

Detailed Overview:
              SAL          N16 44.07, W022 56.64, FL000, 11:04:13 UTC, 
to            junction     N14 45.00, W025 00.00, FL390, 11:30:52 UTC, 
to            ec_north     N14 45.00, W029 38.68, FL410, 12:06:37 UTC, 
to            c_mid        N08 30.09, W030 50.49, FL410, 12:56:38 UTC, 
to            x_sonde      N06 15.06, W031 15.95, FL410, 13:14:39 UTC, Position for an extra sonde heading south
circle around c_south      N04 00.00, W031 41.20, FL410, 13:23:10 UTC - 14:22:44 UTC, radius: 72 nm, 360° CCW, enter from north
to            x_cms        N07 45.00, W030 59.00, FL450, 14:43:04 UTC, might be to soon to rise to FL450, if so rise to FL430 and higher later
circle around c_mid        N08 30.09, W030 50.49, FL450, 14:46:29 UTC - 15:45:01 UTC, radius: 72 nm, 360° CCW, enter from south
to            ec_under     N10 18.08, W030 30.02, FL450, 16:08:30 UTC, meet EarthCARE, drop sonde
to            x_cns        N11 30.00, W030 16.30, FL450, 16:17:56 UTC, Southern alternative center for c_north
circle around c_north      N13 00.00, W029 59.02, FL450, 16:39:04 UTC - 17:37:36 UTC, radius: 72 nm, 360° CCW, enter from north, fly through circle before entering
to            x_cnn        N13 40.20, W029 51.26, FL450, 17:41:39 UTC, Northern alternative center for c_north
to            ec_north     N14 45.00, W029 38.68, FL450, 17:50:09 UTC, 
to            junction     N14 45.00, W025 00.00, FL450, 18:25:06 UTC, 
to            j_atr        N14 50.28, W021 45.79, FL350, 18:50:04 UTC, ATR circle start and end
circle around c_atr        N15 23.00, W022 06.00, FL350, 18:50:04 UTC - 19:22:42 UTC, radius: 38 nm, 360° CCW, enter from south east
to            SAL          N16 44.07, W022 56.64, FL000, 19:44:03 UTC, 

/home/runner/miniconda3/envs/orcestra_book/lib/python3.12/site-packages/orcestra/flightplan.py:214: FutureWarning: The aircraft attribute of this FlightPlan is not set, but it's required for aircraft performance calculations. Currently, the code will proceed using the default aircraft performance, but the code may break in future.
  warn(