wrf-model/WPS

Wrong NLATS in processing NCAR RDA's ECMWF "ds113.1" dataset grib

ben108472 opened this issue · 1 comments

I have processed grib files from RDA's ECMWF "ds113.1" dataset by using ungrib and metgrid.
But I found the wrong information with intermediate files.
FIELD = GHT
UNITS = m DESCRIPTION = Height
DATE = 2019-10-05_00:00:00 FCST = 0.000000
SOURCE = ECMWF
LEVEL = 10000.000000
I,J DIMS = 1849, 783
IPROJ = 4 PROJECTION = GAUSSIAN
REF_X, REF_Y = 1.000000, 1.000000
REF_LAT, REF_LON = 60.000004, 80.016006
NLATS, DLON = 1280.000000, 0.070000
EARTH_RADIUS = 6367.470215
DATA(1,1)=16054.108398

NLATS is 1280, not -0.07.
So the metgrid shows the error.
Oops, something is not right with the Gaussian latitude computation.
The input data gave the starting latitude as 60.000.
This routine computed the starting latitude as +- 89.946.
The difference is larger than 0.01 degrees, which is not expected.
ERROR: Gaussian_latitude_computation

Is there any soultion?

Use the global region can fix this problem.