Dfs2 - Global Forecasting System.ipynb result was backward

[Junik4]

I just wanna be convert netcdf data from ECMWF to dfs2, see plot in below :

then i running the script Dfs2 - Global Forecasting System.ipynb result was backward, see plot in below :

basically my script is below:

!pip install mikeio
import mikeio

from google.colab import drive
drive.mount('/content/drive') # Mount to the base directory first
import os

Create the target directory if it doesn't exist

os.makedirs('/content/drive/MyDrive/Project_XXX/', exist_ok=True)

Now you can access files within this directory

Example:

my_file_path = '/content/drive/MyDrive/Project_XXX/my_file.txt'

pip install netcdf4 xarray

import numpy as np
import pandas as pd
import xarray as xr
import mikeio

!ls /content/drive/MyDrive/Project_XXX/

from google.colab import drive
drive.mount('/content/drive')

from google.colab import drive
drive.mount('/content/drive', force_remount=True) # Force remount to ensure connection

import xarray as xr

Assuming the file is in your Google Drive under "Project_PBR"

file_path = '/content/drive/MyDrive/Project_XXX/Wave ECMWF XXX.nc'

Print the current working directory for debugging

import os
print(os.getcwd())

List files in the target directory to verify file existence

!ls '/content/drive/MyDrive/Project_XXX/'

Try opening the dataset

ds = xr.open_dataset(file_path)
ds

ds["swh"]

ds.swh

ds.swh.isel(valid_time=0).plot();

time = pd.DatetimeIndex(ds.valid_time)
time

mwd = ds.mwd.values
mwp = ds.mwp.values
swh = ds.swh.values

geometry = mikeio.Grid2D(x=ds.longitude.values, y=np.sort(ds.latitude.values), projection="LONG/LAT")

from mikeio import ItemInfo, EUMType, EUMUnit

mwd_da = mikeio.DataArray(data=mwd,time=time, geometry=geometry, item=ItemInfo("Mean Wave Direction", EUMType.Wave_direction, EUMUnit.degree))
mwp_da = mikeio.DataArray(data=mwp,time=time, geometry=geometry, item=ItemInfo("Mean Wave Period", EUMType.Wave_period, EUMUnit.second))
swh_da = mikeio.DataArray(data=swh,time=time, geometry=geometry, item=ItemInfo("Significant Wave Height", EUMType.Wave_height, EUMUnit.meter))

mds = mikeio.Dataset([mwd_da, mwp_da, swh_da])
mds

mds.to_dfs("wave ecmwf xxx.dfs2")

---

is there any error in the script above?
please give me solution, Thanks

[ecomodeller]

The error is that you forgot to flip the data upside down along the y axis.