/********************************************************************** /********************************************************************** /* Name: genhrap.aml /* Purpose: Generate polygon coverage(s) from user specified input /* file(s) (i.e. inputgc.COD)generated by genhrap.f, project the /* polygon coverage into chosen projection. Create an INFO file with /* HRAP-IDs (given hrap.COD.dat), and join this INFO file to the PAT of /* the projected polygon coverage. /* /* Notes: (1) This program should be run following the execution of /* a compiled version of genhrap.f has been executed. This program /* prompts the user for the three character suffix used in genhrap.f /* and the name of a projection file. An example projection /* (albdd.prj) is listed in the comments at the end /* of this file. Any other output projection parameters could be /* substituted for those listed. The input projection needs to be /* geographic but the choice of an input datum is a moot point because /* NEXRAD data are generated in a spherical coordinate system -- see /* our report for further explanation. /********************************************************************** /********************************************************************** &sv suff = [response 'Enter the 3 character suffix used to ID hrap ~ files:'] &sv prjfile [response 'Enter the name of your projection file'] &sv covgc = %suff%geocc &sv inputgc = inputgc.%suff% &if [exists %covgc% -cover] &then kill %covgc% all generate %covgc% &if [exists %inputgc% -file] &then input %inputgc% &else &type Can't find input file. polys /* must quit out of the GENERATE sub-program quit clean %covgc% &sv covgcprj = %covgc%prj &if [exists %covgcprj% -cover] &then kill %covgcprj% all project cover %covgc% %covgcprj% %prjfile% clean %covgcprj% tables &if [exists hrapxy2.dat -info] &then &sv delvar = [delete hrapxy2.dat -info] /* Add data to the INFO file hrapxy2.dat from the file hrap.***.dat /* created by the FORTRAN program create.f &sv addfile = hrap.%suff%.dat /*add from %addfile% define hrapxy2.dat %covgcprj%-id 5 5 i hrapx 4 4 i hrapy 4 4 i ~ add from %addfile% quit /* Join the newly created INFO file to the PAT, creating two new colums /* in the HRAP polygon coverage joinitem %covgcprj%.pat hrapxy2.dat %covgcprj%.pat %covgcprj%-id ~ %covgcprj%-id ordered &return /* Listing of albdd.prj /*input /*projection geographic /*units dd /*datum wgs72 /*parameters /*output /*projection albers /*units meters /*datum wgs72 /*parameters /*29 30 00 /*45 30 00 /*-96 00 00 /*23 00 00 /*0.0 /*0.0 /*end