I am running an MPAS Limited Area simulation over CONUS with ERA5 boundary conditions and find that the lateral boundary condition files produced by MPAS (lbc*nc) contain wrong values that are not in line with the underlying ERA5 data they were created from. Specifically, I can see that all 55 levels contain the same values for lbc_qv which should vary with altitude. Strangely, I do not see the same values for all 55 levels for other variables (e.g. lbc_u).
The result is also evident in the simulation result. There is no visible inflow of moisture in the Pacific Northwest that is clearly visible in ERA5:
The log-file for the model run shows that the lateral boundary conditions are correctly updated, but I am confused about the log.init_atmosphere file used for creating the boundary conditions because it shows that it assumes model levels and interpolates values onto pressure levels, although my intermediate files are structured in pressure levels.
I used this python script to create the intermediate files: GitHub - NCAR/era5_to_int: A simple Python script for converting ERA5 model-level netCDF files to the WPS intermediate format and those files are located in /glade/derecho/scratch/kukulies/era5/met_data/.
The lbc*nc files are located in /glade/campaign/mmm/c3we/juliak/MPAS_C404/limited_area_3km/ and the run directory for the simulation (and creation of lateral boundary conditions) is /glade/work/kukulies/MPAS_v8.2.3/MPAS-Model.
Many thanks in advance!
Julia
The result is also evident in the simulation result. There is no visible inflow of moisture in the Pacific Northwest that is clearly visible in ERA5:
The log-file for the model run shows that the lateral boundary conditions are correctly updated, but I am confused about the log.init_atmosphere file used for creating the boundary conditions because it shows that it assumes model levels and interpolates values onto pressure levels, although my intermediate files are structured in pressure levels.
I used this python script to create the intermediate files: GitHub - NCAR/era5_to_int: A simple Python script for converting ERA5 model-level netCDF files to the WPS intermediate format and those files are located in /glade/derecho/scratch/kukulies/era5/met_data/.
The lbc*nc files are located in /glade/campaign/mmm/c3we/juliak/MPAS_C404/limited_area_3km/ and the run directory for the simulation (and creation of lateral boundary conditions) is /glade/work/kukulies/MPAS_v8.2.3/MPAS-Model.
Many thanks in advance!
Julia