Inconsistent parameters between projections extracted by wrf-python and xwrf

[thaolinhtran]

Hi there,

I wanted to utilize the projection extracted by xwrf for further calculations. However, I found an inconsistency in the value of "+lon_0" in the wrf_projection created by xwrf compared to the one from get_cartopy in wrf-python. Please see the details in the code snippets below.

I would appreciate any suggestions to address this issue.

Thank you in advance.

Best regards,
Linh

from netCDF4 import Dataset
import xwrf
from wrf import (get_cartopy, getvar)

ds        = xr.open_mfdataset(path+"wrfout_d01*", engine='netcdf4', concat_dim="Time", combine="nested" ); ds=ds.xwrf.postprocess(); ds=ds.xwrf.destagger()
rootgroup = Dataset(infile)

proj1     = ds['wrf_projection'].values.item()
proj2     = get_cartopy(getvar(rootgroup, "pw", timeidx=0))

print(proj1) ## printed values: +proj=merc +x_0=0 +y_0=0 +a=6370000 +b=6370000 +lon_0=172.90992736816406 +lat_ts=-15.0 +type=crs
print(proj2) ## printed values: +proj=merc +a=6370000.0 +b=6370000.0 +nadgrids=@null +lon_0=152.0 +lat_ts=-15.0 +k=1 +units=m +no_defs +type=crs

[lpilz]

Hi Linh,

thanks for reporting this. We are checking our projections against salem at them moment. Would you be able to generate a projection using salem as well and compare the results with that, too?

Thanks!