Metpoint Exec and Conversion of Lat-Lon to Grid Point

[kmsxkuse]

I am currently parsing HYSPLIT archival weather data to pull relevant surface and altitude data for the continental United States using CHK_DATA.F which was available here: https://www.ready.noaa.gov/data/archive ... chk_data.f

After successfully parsing the records of the meteorological data header, I need to know which grid point corresponds to a specific latitude and longitude provided in an input parameter. The exec program Metpoint does this, but I would prefer to not have to bundle a specific executable when everything else is written in python.

I would like to know what the procedure is for converting a latitude and longitude into a grid point using the ARL-HYSPLIT data format.
Is the process identical to this module used for WRF?
https://wrf-python.readthedocs.io/en/la ... _proj.html

[sonny.zinn]

Please refer to https://www.arl.noaa.gov/hysplit/cmapf- ... -routines/. The METPOINT program initializes the mapping routine and then converts latitude and longitude to grid indices by calling GBL2XY() or CLLX2XY() depending on the coordinate system and truncating the resulting X, Y coordinates into integers. The source code for the mapping routines including GBL2XY() and others may be downloaded from the website. A detailed description of coordinate transforms is given in an article under Answer #1.