wrf-model / WPS

The official repository for the WRF Preprocessing System (WPS)
202 stars 159 forks source link

Adding support for CMA-GFS in WPS #236

Open andrewsoong opened 12 months ago

andrewsoong commented 12 months ago

The CMA-GFS grib2 datasets is now able to be ingested by the ungrib program. The GRAPES_GFS(here in after referred to as CMA_GFS), the new generation of Global/Regional Assimilation and Prediction Enhanced System which is independently developped by China Meteorological Administration (CMA), has been put into formal operation. So we are committed to solving the problem of WPS correctly handling CMA_GFS forecast data.

one Vtable files(Vtable.CMA_GFS) has been created with the fields required by the WRF model. Modifications to the ungrib source code: Adding support for code 738 in rd_grid2.F because model levels of CMA_GFS are on generalized vertical height coordinates.

The CMA_GFS real time data is here:http://data.wis.cma.cn/DCPC_WMC_BJ/open/nwp/gmf_gra/

I've tested the program with an Intel compiler and it passed without errors. CMA_GFS forecast data were also used to test, the results are normal.

jimbresch commented 4 months ago

A test of the modified code on a dataset pulled from that CMA web site segfaulted. I haven't found the cause, but note that CMA uses several vertical level types that aren't decoded by rd_grib2. Their files also cause problems for g2print. So, it will take debugging and mods to get it working.

I also suggest renaming Vtable.CMA_GFS to something else, since it may be confused with the NCEP GFS by some users. Vtable.CMA_GRAPES perhaps? Otherwise we may need to have a naming convention of the form Vtable.center_model and rename all the NCEP Vtables.

andrewsoong commented 4 months ago

A test of the modified code on a dataset pulled from that CMA web site segfaulted. I haven't found the cause, but note that CMA uses several vertical level types that aren't decoded by rd_grib2. Their files also cause problems for g2print. So, it will take debugging and mods to get it working.

I also suggest renaming Vtable.CMA_GFS to something else, since it may be confused with the NCEP GFS by some users. Vtable.CMA_GRAPES perhaps? Otherwise we may need to have a naming convention of the form Vtable.center_model and rename all the NCEP Vtables.

Hello, Jim. The China Meteorological Administration(CMA) has changed its name from GRAPES_GFS to CMA_GFS. To differentiate from NCEP GFS, it is necessary to remind users. Do you now need to change Vtable.CMA_GFS to Vtable.CMA_GRAPES? Also, I think we need to have a naming convention of the form Vtable.center_model. Thanks !

jimbresch commented 4 months ago

Hi Andrew, OK, we'll keep the name as Vtable.CMA_GFS. The decoding is failing inside an NCEP g2 library routine. Do you have a dataset that has been successfully ungribbed with WPS? Every file that I've tested from the CMA web site is very large and has failed to be decoded. The date tag on your code modification is from 2018. If CMA increased the resolution of their GFS since then, it might explain the failure.

andrewsoong commented 4 months ago

Hi Andrew, OK, we'll keep the name as Vtable.CMA_GFS. The decoding is failing inside an NCEP g2 library routine. Do you have a dataset that has been successfully ungribbed with WPS? Every file that I've tested from the CMA web site is very large and has failed to be decoded. The date tag on your code modification is from 2018. If CMA increased the resolution of their GFS since then, it might explain the failure.

Hi,Jim, I have just tested the latest CMA-GFS forecast data from:http://data.wis.cma.cn/DCPC_WMC_BJ/open/nwp/gmf_gra/t1200/f0_f240_6h/, WPS-V4.5 runs well, generated normal data. See: hsong_2024-05-21_8_09_45.txt for test information. My Linux platform is Rock Linux 9.2, and the Intel compiler version is 2021.10.0. By the way, do you need me to mail you any data files for your tests? What's your e-mail address? Feel free to email me. hsong_2024-05-21_8_09_45.txt

Haiqing Soong (宋海清) Ph.D, Associate Professor, Ecological and Agricultural Meteorology Center of Inner Mongolia & Inner Mongolia Meteorological Satellite Remote Sensing Center, Building A, Inner Mongolia Meteorological Bureau, Hailaer Road, Xincheng District, Hohhot, 010051,China Tel: 0471-3335230 E-Mail: haiqingsong2010@163.com, and haiqingsong@emails.imau.edu.cn

andrewsoong commented 3 months ago

@jimbresch Do We need to resubmit this PR based on WPS-V4.6.0?

weiwangncar commented 3 months ago

@andrewsoong You should be able to continue to amend this PR, unless you are addressing a different issue.