This archive contains answers to questions sent to Unidata support through mid-2025. Note that the archive is no longer being updated. We provide the archive for reference; many of the answers presented here remain technically correct, even if somewhat outdated. For the most up-to-date information on the use of NSF Unidata software and data services, please consult the Software Documentation first.
*text was cut-off in my previous reply* Kevin, RE: gridded hail size, you're going to need write hail size data into SFLIST format for importing into a GEMPAK surface file. in short: 1) format hail size into readable format 2) create empty surface file with SFCFIL http://www.unidata.ucar.edu/cgi-bin/gempak/manual/programs_index?sfcfil 3) import hail size into this surface file with SFEDIT http://www.unidata.ucar.edu/cgi-bin/gempak/manual/programs_index?sfedit 4) create an empty grid file with GDCFIL http://www.unidata.ucar.edu/cgi-bin/gempak/manual/programs_index?gdcfil 5) perform an objective analysis on your surface data using OABSFC http://www.unidata.ucar.edu/cgi-bin/gempak/manual/programs_index?oabsfc -------------- more detailed: 1) as of now you have hail size saved in a surface file of the form (making an example of a 2-station file): KSTL 724340 ST._LOUIS/LAMBERT MO US 3875 -9037 172 0 0. KCGI 723489 CAPE_GIRARDEAU_MUNI MO US 3723 -8957 104 0 2.35 what we need is a separate file of the form: PARM = HAIL KSTL 100607/1200 0.00 KCGI 100607/1200 2.35 with date/time specified. the first line specifies the parm names for columns that come after date/time, so if you wanted to plot both hail size and tmpf, for example, your file would look like. PARM = TMPF;HAIL KSTL 100607/1200 77.7 0.00 KCGI 100607/1200 88.8 2.35 *NOTE* - if you have both TMPF and HAIL specified in your file, like this example, and you want to add both TMPF and HAIL to your GEMPAK surface file, you will need to specify both when creating your file with SFCFIL. make sense? (let's assume we only have hail size and not tmpf for the rest of this problem) so now you have your original "rucstns6_max" file and this new file we'll call "rucstns6.fil" -------------- 2) your next step is to create a surface gempak file using SFCFIL. GEMPAK-SFCFIL> SFOUTF = hail.sfc SFPRMF = HAIL STNFIL = rucstns6_max SHIPFL = NO TIMSTN = 1/0 SFFSRC = text GEMPAK-SFCFIL>r the new (empty) surface file is called "hail.sfc" and "rucstns6_max" is specified as the station information file. run this and you should have an empty surface file approx. ~17k -------------- 3) run SFEDIT to add hail size values to hail.sfc with SFEFIL = rucstns6.fil GEMPAK-SFEDIT>l SFEFIL = rucstns6.fil SFFILE = hail.sfc GEMPAK-SFEDIT>r you should see: SFEDIT PARAMETERS: Edit file: rucstns6.fil Output surface file: hail.sfc Parameters to be added to file: HAIL Enter <cr> to accept parameters or type EXIT: -------------- now exit SFEDIT and check your file with SFLIST GEMPAK-SFLIST>l SFFILE = hail.sfc AREA = dset DATTIM = all SFPARM = HAIL OUTPUT = T IDNTYP = STID GEMPAK-SFLIST>r you should see PARM = HAIL STN YYMMDD/HHMM HAIL KSTL 100607/1200 0.00 KCGI 100607/1200 2.35 * note here that 100607/1200 is the tmp date/time I specified when creating rucstns6.fil -------------- 4) create an empty grid file with GDCFIL http://www.unidata.ucar.edu/cgi-bin/gempak/manual/programs_index?gdcfil GEMPAK-GDCFIL>l GDOUTF = hail.gem PROJ = MER GRDAREA = us KXKY = 10;10 MAXGRD = 200 CPYFIL = copy.gem ANLYSS = 4/2;2;2;2 GEMPAK-GDCFIL> where copy.gem is an existing grid/model file specified for navigation. you should see: Creating process: gplt for queue 557058 GDCFIL PARAMETERS: New grid file: hail.gem File copied: copy.gem Grid projection: MER Grid area: us Grid dimension: 10;10 Grid anlysis: 4/2;2;2;2 Maximum number of grids that can be added to file 200 Enter <cr> to accept parameters or type EXIT: -------------- 5) perform an objective analysis on your surface data using OABSFC outlined in detail in the following two links: http://www.unidata.ucar.edu/cgi-bin/gempak/manual/programs_index?oabsfc http://www.unidata.ucar.edu/software/gempak/tutorial/barnes3.html Ticket Details =================== Ticket ID: NZZ-240543 Department: Support GEMPAK Priority: Critical Status: Open