Flight plan - KA-20240827a

Contents

Flight plan - KA-20240827a#

in_cloud spiral meteor

Crew#

The flight is planned to take off at 2024-08-27 15:30:00+00:00.

Job

Name

PI

Robert Oscar David

Instrument operator

Sorin Ghemulet

Instrument operator

Marius Corbu

Flight plan#

Arrive @FL200 at 16:09 UTC at WP1 - holding pattern until 16:29 UTC.

At 16:29 UTC: WP1 -> WP2 at cloud top (max FL280, ~20min). HALO passes over us at 16:39 UTC at FL450 at HALO_meet (13.0 N, -24.04 W). Holding pattern at WP2 for 5 min to allow for westward scan direction of SEA-POL radar on RV Meteor.

16:55 UTC: WP2 -> WP3 -> WP2 -> …: Descent with 1000 ft/min from cloud top down to FL010.

Arrive at WP4 (RV Meteor) at FL010 (safe altitude from SEA-POL radar).

Spiral ascent from FL010 to cloud top (max. FL280). Ascent as race-track pattern with one straight leg between WP5 and WP6: - WP5: -23.4785 W, 13.0 N - WP6: -23.5215 W, 13.0 N

Return home to RAI: descent through potential clouds on the way home.

Hide code cell source 
from orcestra.flightplan import sal, bco, LatLon, IntoCircle, path_preview, plot_cwv
from datetime import datetime
import intake
import easygems.healpix as egh

cat = intake.open_catalog("https://tcodata.mpimet.mpg.de/internal.yaml")

# Define dates for forecast initialization and flight

issued_time = datetime(2024, 8, 27, 0, 0, 0)
issued_time_str = issued_time.strftime('%Y-%m-%d')

flight_time = datetime(2024, 8, 27, 12, 0, 0)
flight_time_str = flight_time.strftime('%Y-%m-%d')
flight_index = f"KA-{flight_time.strftime('%Y%m%d')}a"

print("Initalization date of IFS forecast: " + issued_time_str + "\nFlight date: " + flight_time_str + "\nFlight index: " + flight_index)

airport = LatLon(lat=14.945, lon=-23.4863889, label='RAI')
wp1 = LatLon(lat=13.0, lon=-24.47, label='wp1')
halo_meet = LatLon(lat=13.0, lon=-24.04, label='HALO_meet')
wp2 = LatLon(lat=13.0, lon=-23.61, label='wp2')
wp3 = LatLon(lat=13.0, lon=-23.8135, label='wp3')
meteor = LatLon(lat=13.0, lon=-23.5, label='meteor')
wp5 = LatLon(lat=13.0, lon=-23.4785, label='wp5')
wp6 = LatLon(lat=13.0, lon=-23.5215, label='wp6')


spiral = LatLon(lat=13.0, lon=-23.5, label='spiral')

leg_out = [
     airport,
     wp1
]

leg_halo = [
      wp1,
      halo_meet,
      wp2
]

leg_seapol = [
      wp2,
      wp3,
      wp2
]

leg_meteor = [
      wp2,
      meteor,
      spiral,
      wp5,
      wp6
]

leg_home = [
     wp6,
     airport
]

path = leg_out + leg_halo + leg_seapol + leg_meteor + leg_home 

cat = intake.open_catalog("https://tcodata.mpimet.mpg.de/internal.yaml")
ds = cat.HIFS(datetime=issued_time).to_dask().pipe(egh.attach_coords)
cwv_flight_time = ds["tcwv"].sel(time=flight_time, method = "nearest")

ax = path_preview(path)

#ax.scatter(-23.07972222, 14.33555556, c='red', marker='o', s=200, ec = 'r', alpha = 0.3)

plot_cwv(cwv_flight_time)

Initalization date of IFS forecast: 2024-08-27
Flight date: 2024-08-27
Flight index: KA-20240827a
../_images/7068bc45f5cb47b6e96f0c9f5173c61ae5f39c461eb9e4f5d2ae8483880bf4ab.png
Hide code cell source 
import pandas as pd
from dataclasses import asdict

pd.DataFrame.from_records(map(asdict, [airport, wp1, halo_meet, wp2, wp3, meteor, wp5, wp6, airport])).set_index("label")
lat lon fl time note
label
RAI 14.945 -23.486389 None None None
wp1 13.000 -24.470000 None None None
HALO_meet 13.000 -24.040000 None None None
wp2 13.000 -23.610000 None None None
wp3 13.000 -23.813500 None None None
meteor 13.000 -23.500000 None None None
wp5 13.000 -23.478500 None None None
wp6 13.000 -23.521500 None None None
RAI 14.945 -23.486389 None None None
Contents