Environment Class Usage#

Here we explore the different environment models available by default in RocketPy. Let’s start by importing the rocketpy module.

[1]:

from rocketpy import Environment, SolidMotor, Rocket, Flight

If you are using Jupyter Notebooks, it is recommended to run the following line to make matplotlib plots which will be shown later interactive and higher quality.

[2]:

%matplotlib inline

Datetime#

RocketPy requires date and time information specifying year, month, day and hour to properly forecast the weather. Alongside this data, a timezone may also be specified as follows.

[3]:

env = Environment(date=(2022, 2, 16, 18), timezone="America/New_York")

It is also possible to omit the timezone, in which case RocketPy will assume the datetime is given in standard UTC time, as done in the examples that follow.

Elevation#

The elevation of the launch site of a rocket is particularly important for determining its interaction with the atmosphere and Earth’s gravitational field.

Sometimes, together with the latitude and longitude of the launch site, the elevation which is measured in meters above sea level is known and can be used to initialize an Environment class instance as follows.

[4]:

env = Environment(
    date=(2019, 2, 10, 18),
    latitude=-21.960641,
    longitude=-47.482122,
    elevation=110,
)

However, this is not always the case. Lucky for us, there are alternatives to find an approximate value of the elevation of the launch site if the latitude and longitude values are know.

One option is to use the value supplied by some atmospheric models. Since elevation data is crucial for Numerical Weather Programing, some weather models make elevation data available together with other variables. But this will be covered later.

Another very useful and handy option is to use Open-Elevation, a free and open-source elevation API. It is built right into RocketPy and one can use it as follows.

First, we initialize a new Environment.

[5]:

env = Environment(date=(2019, 2, 10, 18), latitude=-21.960641, longitude=-47.482122)

Then, we use the Environment.setElevation method with the input parameter as ‘Open-Elevation’. Note, however, that Open-Elevation servers are not always available.

[6]:

env.set_elevation("Open-Elevation")

Fetching elevation from open-elevation.com...
Elevation received: 651.0

To get information from the Environment, just use the Environment.info() or the Environment.all_info() methods.

[7]:

env.info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.785541145259767 m/s²


Launch Site Details

Launch Date: 2019-02-10 18:00:00 UTC
Launch Site Latitude: -21.96064°
Launch Site Longitude: -47.48212°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 243664.48 W    7569452.20 S
Launch Site UTM zone: 23K
Launch Site Surface Elevation: 651.0 m


Atmospheric Model Details

Atmospheric Model Type: standard_atmosphere
standard_atmosphere Maximum Height: 80.000 km


Surface Atmospheric Conditions

Surface Wind Speed: 0.00 m/s
Surface Wind Direction: 0.00°
Surface Wind Heading: 0.00°
Surface Pressure: 937.46 hPa
Surface Temperature: 283.93 K
Surface Air Density: 1.150 kg/m³
Surface Speed of Sound: 337.62 m/s


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_15_1.png

Using the Standard Atmosphere#

By default, when initializing an environment class, the International Standard Atmosphere as defined by ISO 2533 is initialized.

Note that the International Standard Atmosphere only has temperature and pressure profiles properly specified. Other profiles can be derived from it, however, winds are automatically set to 0 m/s.

[8]:

EnvISA = Environment()

[9]:

EnvISA.info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.7803253359 m/s²


Launch Site Details

Launch Site Latitude: 0.00000°
Launch Site Longitude: 0.00000°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 166021.44 W|E    0.00 N
Launch Site UTM zone: 31N
Launch Site Surface Elevation: 0.0 m


Atmospheric Model Details

Atmospheric Model Type: standard_atmosphere
standard_atmosphere Maximum Height: 80.000 km


Surface Atmospheric Conditions

Surface Wind Speed: 0.00 m/s
Surface Wind Direction: 0.00°
Surface Wind Heading: 0.00°
Surface Pressure: 1013.25 hPa
Surface Temperature: 288.15 K
Surface Air Density: 1.225 kg/m³
Surface Speed of Sound: 340.29 m/s


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_19_1.png

The International Standard Atmosphere can also be reset at any time by using the following method.

[10]:

EnvISA.set_atmospheric_model(type="standard_atmosphere")

Using a Custom Atmosphere#

The user can also set a completely customized atmosphere. This is particularly useful for setting custom values of wind in both directions.

First, we initialize a new Environment.

[11]:

EnvCA = Environment()

Then, we set the atmospheric model as ‘custom_atmosphere’.

Leaving the pressure field as None means we want the International Standard Atmosphere’s pressure profile to be used. We could have done the same with temperature, but to showcase how floats can be used, we set the temperature field as a constant 300 K profile.

For the wind, we need to specify its value in both U (east) and V (north) direction. In this case, we used arrays to specify points. Consider a wind U profile of 5 m/s at 0 m and 10 m/s at 1000 m. For the wind V, we used -2 m/s at 0 m, 3 m/s at 500 m and 2 m/s at 1000 m.

[12]:

EnvCA.set_atmospheric_model(
    type="custom_atmosphere",
    pressure=None,
    temperature=300,
    wind_u=[(0, 5), (1000, 10)],
    wind_v=[(0, -2), (500, 3), (1600, 2)],
)

[13]:

EnvCA.all_info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.7803253359 m/s²


Launch Site Details

Launch Site Latitude: 0.00000°
Launch Site Longitude: 0.00000°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 166021.44 W|E    0.00 N
Launch Site UTM zone: 31N
Launch Site Surface Elevation: 0.0 m


Atmospheric Model Details

Atmospheric Model Type: custom_atmosphere
custom_atmosphere Maximum Height: 1.600 km


Surface Atmospheric Conditions

Surface Wind Speed: 5.39 m/s
Surface Wind Direction: 291.80°
Surface Wind Heading: 111.80°
Surface Pressure: 1013.25 hPa
Surface Temperature: 300.00 K
Surface Air Density: 1.177 kg/m³
Surface Speed of Sound: 347.22 m/s



Gravity Model Plots
../../_images/notebooks_environment_environment_class_usage_27_1.png 


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_27_3.png

Keep in mind that importing .csv files for Custom Atmospheres is also possible. See the documentation on Environment.set_atmospheric_model for more information.

Using Wyoming Upper Air Soundings#

Another option to define atmospheric conditions is to use upper air soundings. The University of Wyoming - College of Engineering - Department of Atmospheric Sciences has a great compilation of atmospheric soundings on their website which can be found here.

For this example, we will use the sounding from 83779 SBMT Marte Civ Observations at 00Z 04 Feb 2019 which can be accessed using the URL generated by the above website: http://weather.uwyo.edu/cgi-bin/sounding?region=samer&TYPE=TEXT%3ALIST&YEAR=2019&MONTH=02&FROM=0500&TO=0512&STNM=83779.

First, we initialize a new Environment. There is no need to specify dates, latitude and longitude here, so we will skip them for know.

[14]:

EnvWUAS = Environment()

Now, we set use the URL above as the input to our atmospheric model.

[15]:

URL = "http://weather.uwyo.edu/cgi-bin/sounding?region=samer&TYPE=TEXT%3ALIST&YEAR=2019&MONTH=02&FROM=0500&TO=0512&STNM=83779"

EnvWUAS.set_atmospheric_model(type="wyoming_sounding", file=URL)

Next, notice how even though we have not specified a surface elevation for our Environment, it has one set! The reason is that Wyoming Upper Air Soundings have elevation data for the sounding station, which is used to set the elevation of the Environment. If you do not want this elevation, you can change using Environment.setElevation(desired elevation in meters).

[16]:

EnvWUAS.all_info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.778096408189615 m/s²


Launch Site Details

Launch Site Latitude: 0.00000°
Launch Site Longitude: 0.00000°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 166021.44 W|E    0.00 N
Launch Site UTM zone: 31N
Launch Site Surface Elevation: 722.0 m


Atmospheric Model Details

Atmospheric Model Type: wyoming_sounding
wyoming_sounding Maximum Height: 22.077 km


Surface Atmospheric Conditions

Surface Wind Speed: 3.09 m/s
Surface Wind Direction: 110.00°
Surface Wind Heading: 290.00°
Surface Pressure: 936.00 hPa
Surface Temperature: 291.75 K
Surface Air Density: 1.118 kg/m³
Surface Speed of Sound: 342.41 m/s



Gravity Model Plots
../../_images/notebooks_environment_environment_class_usage_34_1.png 


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_34_3.png

Using NOAA’s Ruc Soundings#

Another option for upper air soundings is NOAA’s Ruc Soundings servers. It allows users to download virtual soundings from numerical weather prediction models such as GFS, RAP and NAM, and also allow for the download of real soundings from the Integrated Global Radiosonde Archive (IGRA).

All of these options can be retrieved as a text file which follows the GSD format. By generating such a file through the link above, the file’s URL can be used to import the atmospheric data into RocketPy.

We will use the same sounding station as we did for Wyoming Soundings. Selecting ROABs as the initial data source, specifying the station through it’s WMO-ID and opting for the ASCII (GSD format) button, the following page opens up: https://rucsoundings.noaa.gov/get_raobs.cgi?data_source=RAOB&latest=latest&start_year=2019&start_month_name=Feb&start_mday=5&start_hour=12&start_min=0&n_hrs=1.0&fcst_len=shortest&airport=83779&text=Ascii%20text%20%28GSD%20format%29&hydrometeors=false&start=latest

Let’s initialize a new Environment to use this as it’s source for atmospheric conditions. There is no need to specify dates, latitude and longitude here as well, so we will skip them for know.

[17]:

EnvRUC = Environment()

[18]:

URL = "https://rucsoundings.noaa.gov/get_raobs.cgi?data_source=RAOB&latest=latest&start_year=2019&start_month_name=Feb&start_mday=5&start_hour=12&start_min=0&n_hrs=1.0&fcst_len=shortest&airport=83779&text=Ascii%20text%20%28GSD%20format%29&hydrometeors=false&start=latest"

EnvRUC.set_atmospheric_model(type="NOAARucSounding", file=URL)

[19]:

EnvRUC.all_info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.778096408189615 m/s²


Launch Site Details

Launch Site Latitude: 0.00000°
Launch Site Longitude: 0.00000°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 166021.44 W|E    0.00 N
Launch Site UTM zone: 31N
Launch Site Surface Elevation: 722.0 m


Atmospheric Model Details

Atmospheric Model Type: NOAARucSounding
NOAARucSounding Maximum Height: 20.116 km


Surface Atmospheric Conditions

Surface Wind Speed: 1.03 m/s
Surface Wind Direction: 180.00°
Surface Wind Heading: 0.00°
Surface Pressure: 936.00 hPa
Surface Temperature: 294.75 K
Surface Air Density: 1.106 kg/m³
Surface Speed of Sound: 344.17 m/s



Gravity Model Plots
../../_images/notebooks_environment_environment_class_usage_38_1.png 


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_38_3.png

Using Weather Forecasts and Reanalysis#

Weather Forecasts and Reanalysis can also be used to set the atmospheric model. Here we will showcase how to import global forecasts such as GFS and also local forecasts such as NAM and RAP for North America, all of these available through OPeNDAP in the NOAA’s NCEP NOMADS website. Keep in mind that other generic forecasts can also be imported.

Furthermore, reanalysis such as ERA5, which can be easily downloaded from ECMWF’s Climate Data Store, will also be showcased.

Setting up Date, Latitude and Longitude#

All Environments which use weather forecasts or reanalysis must have date, latitude and longitude set up correctly during or after initialization in order for everything to work smoothly.

Let’s first initialize a new Environment for the GFS weather forecast by specifying date, latitude and longitude during initialization. Keep in mind that date is specified using UTC time in year, month, day, UTC hour format.

[20]:

import datetime

tomorrow = datetime.date.today() + datetime.timedelta(days=1)

date_info = (tomorrow.year, tomorrow.month, tomorrow.day, 9)  # Hour given in UTC time

print("Tomorrow's date:", date_info)

Tomorrow's date: (2023, 8, 12, 9)

[21]:

EnvGFS = Environment(date=date_info, latitude=-21.960641, longitude=-47.482122)

Another approach is to specify date, latitude and longitude values after initialization. Let’s do this when creating a new Environment instance that will be using the NAM weather forecast.

[22]:

EnvNAM = Environment()

[23]:

EnvNAM.set_date(date_info)

[24]:

EnvNAM.set_location(-21.960641, -47.482122)

GFS#

Using the latest forecast from GFS is really simple. All we need to do is specify that GFS is the forecast we want. Note that since data is downloaded from NOMADS server, this line of code can take longer than usual.

[25]:

EnvGFS.set_atmospheric_model(type="Forecast", file="GFS")

[26]:

EnvGFS.all_info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.785478852100368 m/s²


Launch Site Details

Launch Date: 2023-08-12 09:00:00 UTC
Launch Site Latitude: -21.96064°
Launch Site Longitude: -47.48212°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 243664.48 W    7569452.20 S
Launch Site UTM zone: 23K
Launch Site Surface Elevation: 671.2 m


Atmospheric Model Details

Atmospheric Model Type: Forecast
Forecast Maximum Height: 78.272 km
Forecast Time Period: From  2023-08-11 12:00:00  to  2023-08-27 12:00:00  UTC
Forecast Hour Interval: 3  hrs
Forecast Latitude Range: From  -90.0 ° To  90.0 °
Forecast Longitude Range: From  0.0 ° To  359.75 °


Surface Atmospheric Conditions

Surface Wind Speed: 1.47 m/s
Surface Wind Direction: 95.72°
Surface Wind Heading: 275.72°
Surface Pressure: 940.78 hPa
Surface Temperature: 295.08 K
Surface Air Density: 1.111 kg/m³
Surface Speed of Sound: 344.36 m/s



Gravity Model Plots
../../_images/notebooks_environment_environment_class_usage_49_1.png 


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_49_3.png

NAM and RAP#

We can also request the latest forecasts from NAM and RAP. Since they are regional models for North America, we need to specify latitude and longitude points in North America. We will use SpacePort America for this, represented by coordinates 32.988528, -106.975056 or 32°59’18.7”N 106°58’30.2”W.

[27]:

EnvNAM = Environment(date=date_info, latitude=32.988517, longitude=-106.975062)

EnvRAP = Environment(date=date_info, latitude=32.988517, longitude=-106.975062)

[28]:

EnvNAM.set_atmospheric_model(type="Forecast", file="NAM")

EnvRAP.set_atmospheric_model(type="Forecast", file="RAP")

[29]:

EnvNAM.info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.791284255443486 m/s²


Launch Site Details

Launch Date: 2023-08-12 09:00:00 UTC
Launch Site Latitude: 32.98852°
Launch Site Longitude: -106.97506°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 315459.04 W    3651746.15 N
Launch Site UTM zone: 13S
Launch Site Surface Elevation: 1415.4 m


Atmospheric Model Details

Atmospheric Model Type: Forecast
Forecast Maximum Height: 31.642 km
Forecast Time Period: From  2023-08-11 18:00:00  to  2023-08-14 06:00:00  UTC
Forecast Hour Interval: 3  hrs
Forecast Latitude Range: From  21.140671 ° To  52.61339827273 °
Forecast Longitude Range: From  -134.096127 ° To  -60.93723597484 °


Surface Atmospheric Conditions

Surface Wind Speed: 2.65 m/s
Surface Wind Direction: 263.06°
Surface Wind Heading: 83.06°
Surface Pressure: 861.50 hPa
Surface Temperature: 294.74 K
Surface Air Density: 1.018 kg/m³
Surface Speed of Sound: 344.16 m/s


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_53_1.png 
[30]:

EnvRAP.info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.791069448766502 m/s²


Launch Site Details

Launch Date: 2023-08-12 09:00:00 UTC
Launch Site Latitude: 32.98852°
Launch Site Longitude: -106.97506°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 315459.04 W    3651746.15 N
Launch Site UTM zone: 13S
Launch Site Surface Elevation: 1485.0 m


Atmospheric Model Details

Atmospheric Model Type: Forecast
Forecast Maximum Height: 16.760 km
Forecast Time Period: From  2023-08-11 21:00:00  to  2023-08-12 18:00:00  UTC
Forecast Hour Interval: 1  hrs
Forecast Latitude Range: From  16.322011 ° To  58.25097463636 °
Forecast Longitude Range: From  -139.856603 ° To  -57.49840949379 °


Surface Atmospheric Conditions

Surface Wind Speed: 3.30 m/s
Surface Wind Direction: 20.02°
Surface Wind Heading: 200.02°
Surface Pressure: 853.48 hPa
Surface Temperature: 296.16 K
Surface Air Density: 1.004 kg/m³
Surface Speed of Sound: 344.99 m/s


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_54_1.png

ERA-5 Reanalysis#

Using data from ERA-5 reanalysis is as simple as specifying the file location and the dictionary to be used. This is much faster since the data has already been downloaded and is stored locally.

[31]:

EnvERA = Environment(
    date=(2018, 2, 17, 18),
    latitude=32.988517,
    longitude=-106.975062,
    elevation=1400,
)

[32]:

EnvERA.set_atmospheric_model(
    type="Reanalysis",
    file="../../../data/weather/SpaceportAmerica_2018_ERA-5.nc",
    dictionary="ECMWF",
)

[33]:

EnvERA.info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.791331652638586 m/s²


Launch Site Details

Launch Date: 2018-02-17 18:00:00 UTC
Launch Site Latitude: 32.98852°
Launch Site Longitude: -106.97506°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 315459.04 W    3651746.15 N
Launch Site UTM zone: 13S
Launch Site Surface Elevation: 1400.0 m


Atmospheric Model Details

Atmospheric Model Type: Reanalysis
Reanalysis Maximum Height: 4.366 km
Reanalysis Time Period: From  2018-01-01 00:00:00  to  2018-12-31 18:00:00  UTC
Reanalysis Hour Interval: 6  hrs
Reanalysis Latitude Range: From  34.5 ° To  31.5 °
Reanalysis Longitude Range: From  252.0 ° To  254.25 °


Surface Atmospheric Conditions

Surface Wind Speed: 2.41 m/s
Surface Wind Direction: 1.56°
Surface Wind Heading: 181.56°
Surface Pressure: 864.76 hPa
Surface Temperature: 282.81 K
Surface Air Density: 1.065 kg/m³
Surface Speed of Sound: 337.13 m/s


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_58_1.png

Other Weather Forecasts or Reanalysis - HIRESW#

We can use other forecasts or reanalysis as well. As an example, we will use the HIRESW CONUS forecast available at NOAA’s NCEP NOMADS website.

[34]:

import datetime

today = datetime.date.today()
date_info = (today.year, today.month, today.day, 12)  # Hour given in UTC time
print("Today's date:", date_info)

EnvHIRESW = Environment(date=date_info, latitude=32.988517, longitude=-106.975062)

Today's date: (2023, 8, 11, 12)

This model can be accessed through the following link: https://nomads.ncep.noaa.gov:9090/dods/hiresw/hiresw20190206/hiresw_conusarw_12z.

Since it is available from NOAA’s NCEP NOMADS, we could use the NOAA dictionary However, we will create one from scratch to showcase the methodology.

The variables from the dataset can also be seen through the link above. We need to determine the name used by the model of the following dimensions and variables:

  • Time

  • Latitude

  • Longitude

  • Pressure Levels

  • Geopotential Height (as a function of Time, Pressure Levels, Latitude and Longitude)

  • Surface Geopotential Height (as a function of Time, Latitude and Longitude)

  • Wind - U Component (as a function of Time, Pressure Levels, Latitude and Longitude)

  • Wind - V Component (as a function of Time, Pressure Levels, Latitude and Longitude)

Looking through the variable list in the link above, we find the following correspondence:

  • Time = ‘time’

  • Latitude = ‘lat’

  • Longitude = ‘lon’

  • Pressure Levels = ‘lev’

  • Geopotential Height = ‘hgtprs’

  • Surface Geopotential Height = ‘hgtsfc’

  • Wind - U Component = ‘ugrdprs’

  • Wind - V Component = ‘vgrdprs’

Therefore, we use the dictionary shown below.

[35]:

HIRESW_dictionary = {
    "time": "time",
    "latitude": "lat",
    "longitude": "lon",
    "level": "lev",
    "temperature": "tmpprs",
    "surface_geopotential_height": "hgtsfc",
    "geopotential_height": "hgtprs",
    "u_wind": "ugrdprs",
    "v_wind": "vgrdprs",
}

This dictionary can then be used in the Environment.set_atmospheric_model method.

[36]:

date_string = f"{date_info[0]}{date_info[1]:02}{date_info[2]:02}"
print("Today's date string:", date_string)

EnvHIRESW.set_atmospheric_model(
    type="Forecast",
    file=f"https://nomads.ncep.noaa.gov/dods/hiresw/hiresw{date_string}/hiresw_conusarw_12z",
    dictionary=HIRESW_dictionary,
)

Today's date string: 20230811

[37]:

EnvHIRESW.info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.79128516880994 m/s²


Launch Site Details

Launch Date: 2023-08-11 12:00:00 UTC
Launch Site Latitude: 32.98852°
Launch Site Longitude: -106.97506°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 315459.04 W    3651746.15 N
Launch Site UTM zone: 13S
Launch Site Surface Elevation: 1415.1 m


Atmospheric Model Details

Atmospheric Model Type: Forecast
Forecast Maximum Height: 12.494 km
Forecast Time Period: From  2023-08-11 12:00:00  to  2023-08-13 12:00:00  UTC
Forecast Hour Interval: 1  hrs
Forecast Latitude Range: From  12.202469 ° To  61.23790536364 °
Forecast Longitude Range: From  -152.852997 ° To  -49.40219446122 °


Surface Atmospheric Conditions

Surface Wind Speed: 2.86 m/s
Surface Wind Direction: 163.98°
Surface Wind Heading: 343.98°
Surface Pressure: 860.83 hPa
Surface Temperature: 293.56 K
Surface Air Density: 1.022 kg/m³
Surface Speed of Sound: 343.47 m/s


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_65_1.png

Using Windy Atmosphere#

Windy.com is a website that lets you see and get weather and atmospheric forecasts in any place in the world. The same atmospheric predictions and data that are available on windy.com can be set as an atmospheric model in RocketPy.

A few different models are accepted:

  • ECMWF-HRES

  • GFS

  • ICON-Global

  • ICON-EU.

Setting up Date, latitude and Longitude#

Just like when using Forecasts and Reanalysis, all Environments which use windy must have date, latitude and longitude set up correctly during or after initialization in order for everything to work smoothly.

Lets first initialize a few different Envioronmets with the same date. Each one for a different atmospheric model.

[38]:

import datetime

tomorrow = datetime.date.today() + datetime.timedelta(days=1)

date_info = (tomorrow.year, tomorrow.month, tomorrow.day, 12)  # Hour given in UTC time

print("Tomorrow's date:", date_info)

Tomorrow's date: (2023, 8, 12, 12)

[39]:

EnvECMWF = Environment(date=date_info, latitude=-21.960641, longitude=-47.482122)

EnvGFS = Environment(date=date_info, latitude=-21.960641, longitude=-47.482122)

EnvICON = Environment(date=date_info, latitude=-21.960641, longitude=-47.482122)

EnvICONEU = Environment(date=date_info, latitude=48.137154, longitude=11.576124)

Setting up the Atmospheric Models#

Now we set up each different atmospheric model

[40]:

EnvECMWF.set_atmospheric_model(type="Windy", file="ECMWF")

EnvGFS.set_atmospheric_model(type="Windy", file="GFS")

EnvICON.set_atmospheric_model(type="Windy", file="ICON")

EnvICONEU.set_atmospheric_model(type="Windy", file="ICONEU")

And printing and plotting everything availabe (just for ECMWF model for simplicity)

[41]:

EnvICONEU.all_info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.807415796147971 m/s²


Launch Site Details

Launch Date: 2023-08-12 12:00:00 UTC
Launch Site Latitude: 48.13715°
Launch Site Longitude: 11.57612°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 691650.37 E    5334754.25 N
Launch Site UTM zone: 32U
Launch Site Surface Elevation: 524.0 m


Atmospheric Model Details

Atmospheric Model Type: Windy
Windy Maximum Height: 13.888 km


Surface Atmospheric Conditions

Surface Wind Speed: 2.19 m/s
Surface Wind Direction: 272.69°
Surface Wind Heading: 92.69°
Surface Pressure: 957.52 hPa
Surface Temperature: 300.78 K
Surface Air Density: 1.109 kg/m³
Surface Speed of Sound: 347.67 m/s



Gravity Model Plots
../../_images/notebooks_environment_environment_class_usage_73_1.png 


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_73_3.png

Using Ensemble Forecasts and Reanalysis#

Ensemble weather forecasts and reanalysis can also be imported. Since ensembles are just a gathering of multiple weather predictions, they are great for dispersion analysis and Monte Carlo simulations. One member forecast of the ensemble can be selected to run each flight simulation.

Two global ensembles are available by default: GEFS and CMC Ensemble, but others can also be used. To showcase another example, an ERA5 Reanalysis ensemble will also be shown.

GEFS#

[42]:

import datetime

tomorrow = datetime.date.today() + datetime.timedelta(days=1)
date_info = (tomorrow.year, tomorrow.month, tomorrow.day, 12)  # Hour given in UTC time

EnvGEFS = Environment(
    date=date_info,
    latitude=-21.960641,
    longitude=-47.482122,
    elevation=640,
)

[43]:

EnvGEFS.set_atmospheric_model(type="Ensemble", file="GEFS")

By default, the 0th forecast member is activated, which means that using EnvGFS will be exactly the same as using the 0th forecast member, and all other forecasts are kept away stored. Changing this is as easy as the following line of code.

[44]:

EnvGEFS.select_ensemble_member(1)

Now, a new forecast is active and EnvGFS behaves as if ensemble member 1 was its only forecast.

A very useful thing is to compare the forecast of all ensemble members, as well as confirming how many are there. All of this can be done by using Environment.all_info().

[45]:

EnvGEFS.all_info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.785575094392566 m/s²


Launch Site Details

Launch Date: 2023-08-12 12:00:00 UTC
Launch Site Latitude: -21.96064°
Launch Site Longitude: -47.48212°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 243664.48 W    7569452.20 S
Launch Site UTM zone: 23K
Launch Site Surface Elevation: 640.0 m


Atmospheric Model Details

Atmospheric Model Type: Ensemble
Ensemble Maximum Height: 31.184 km
Ensemble Time Period: From  2023-08-11 12:00:00  to  2023-08-27 12:00:00  UTC
Ensemble Hour Interval: 6  hrs
Ensemble Latitude Range: From  -90.0 ° To  90.0 °
Ensemble Longitude Range: From  0.0 ° To  359.5 °
Number of Ensemble Members: 21
Selected Ensemble Member: 1  (Starts from 0)


Surface Atmospheric Conditions

Surface Wind Speed: 2.09 m/s
Surface Wind Direction: 30.01°
Surface Wind Heading: 210.01°
Surface Pressure: 947.13 hPa
Surface Temperature: 297.20 K
Surface Air Density: 1.110 kg/m³
Surface Speed of Sound: 345.60 m/s



Gravity Model Plots
../../_images/notebooks_environment_environment_class_usage_81_1.png 


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_81_3.png 


Ensemble Members Comparison
../../_images/notebooks_environment_environment_class_usage_81_5.png

CMC Ensemble#

[46]:

EnvCMC = Environment(
    date=date_info,
    latitude=-21.960641,
    longitude=-47.482122,
    elevation=640,
)

[ ]:

EnvCMC.set_atmospheric_model(type="Ensemble", file="CMC")

[48]:

EnvCMC.all_info()


Gravity Details

Acceleration of Gravity at Lauch Site: 9.785575094392566 m/s²


Launch Site Details

Launch Date: 2023-08-12 12:00:00 UTC
Launch Site Latitude: -21.96064°
Launch Site Longitude: -47.48212°
Reference Datum: SIRGAS2000
Launch Site UTM coordinates: 243664.48 W    7569452.20 S
Launch Site UTM zone: 23K
Launch Site Surface Elevation: 640.0 m


Atmospheric Model Details

Atmospheric Model Type: Ensemble
Ensemble Maximum Height: 31.199 km
Ensemble Time Period: From  2023-08-11 12:00:00  to  2023-08-27 12:00:00  UTC
Ensemble Hour Interval: 6  hrs
Ensemble Latitude Range: From  -90.0 ° To  90.0 °
Ensemble Longitude Range: From  0.0 ° To  359.5 °
Number of Ensemble Members: 21
Selected Ensemble Member: 0  (Starts from 0)


Surface Atmospheric Conditions

Surface Wind Speed: 6.05 m/s
Surface Wind Direction: 31.30°
Surface Wind Heading: 211.30°
Surface Pressure: 944.91 hPa
Surface Temperature: 295.75 K
Surface Air Density: 1.113 kg/m³
Surface Speed of Sound: 344.75 m/s



Gravity Model Plots
../../_images/notebooks_environment_environment_class_usage_85_1.png 


Atmospheric Model Plots
../../_images/notebooks_environment_environment_class_usage_85_3.png 


Ensemble Members Comparison
../../_images/notebooks_environment_environment_class_usage_85_5.png

ERA5 Ensemble#

[49]:

EnvERA5 = Environment(date=(2018, 2, 17, 18), latitude=-21.960641, longitude=-47.482122)

[ ]:

EnvERA5.set_atmospheric_model(
    type="Ensemble",
    file="../../../data/weather/CuritibaRioSaoPauloEnsemble_2018_ERA-5.nc",
    dictionary="ECMWF",
)

[ ]:

EnvERA5.all_info()