search

COSIMA / mom6-panan

Pan-Antarctic regional configuration of MOM6
MIT License
6 stars 6 forks source link

Regridding the 1/20th Boundary forcing to 1/40th #33

Closed willaguiar closed 1 year ago

willaguiar commented 1 year ago

I am trying to regrid the panan 1/20th boundary forcing to the 1/40th grid, but xesmf keeps returning this error:

can't currently regrid a grid that contains a DE of width less than 2

Which seems to be caused by xesmf requiring Lon,Lat,z coordinates, while the boundary forcing has only Lon,Z, and Lat is a single value. (It seems like this is a problem people found before)

Has anyone here found this problem before and have a way to fix it? Or perhaps another suggested method for interpolation. @angus-g @adele157 @fabiobdias @micaeljtoliveira

1/20th boundary forcing: /g/data/x77/ahg157/inputs/mom6/panan/forcing_obc_shifted.nc Issue relative to: https://github.com/COSIMA/mom6-panan/issues/29#issuecomment-1469051318

angus-g commented 1 year ago

You might be able to use xarray's interp (although it might not handle the periodic boundary properly). Otherwise, you could just do 50/50 linear interpolation between the existing points since it's a double-resolution grid.

willaguiar commented 1 year ago

Thanks. For the 50/50 case, can it br a simple interpolation, or we need to be conserving ?(Interpolating T,S,U,V).