#! /usr/bin/perl

use WXP;

WXP::Prog::readResrcFile;

$mspec = new WXP::ModelSpec( "nam" );
$file = new WXP::Filename( $mspec->getName(0), "la" );

$grib = new WXP::GribFile( $file );
if( !$grib->isOpen ){
    printf( "File not found: %s\n", $file->getName );
    exit 1;
}
$grids = new WXP::Grids;
$values = new WXP::FloatArray;
$ept = new WXP::EarthPoint;
$city = new WXP::CityFile;
$cdata = new WXP::CityData;
$date = new WXP::Date;
$date->set();
$date->print;

$grids->addSpec( "+24hr-12acc::prec:in" );
#$grids->addSpec( "+24hr-6acc::prec:in" );
#$grids->addSpec( "+24hr:sl:prese:mb" );
$grids->addSpec( "+24hr:sl:slpres:mb" );
$grids->addSpec( "+24hr:2m_ag:temp:C" );
$grids->addSpec( "+24hr:1000mb:hght:m" );
$grids->addSpec( "+24hr:500mb:hght:m" );

#$grib->read( $grids );
$gridread = new WXP::GridRead();

$gridread->read( "nam", $date, $grids );

$city->open( "/home/wxp/etc/sao.cty" );

while( $city->read( $cdata )){
   next if $cdata->getPrior > 3 || $cdata->getRegion ne "US";
   next if $cdata->getPrior > 1;
   $ept->set( $cdata->getLat, $cdata->getLon );
   $grids->getValue( $ept, $values );

   printf "%4s %5.2f %6.1f %5.1f %4.0f %4.0f\n", $cdata->getId, $values->get(0),
      $values->get(1), $values->get(2), $values->get(3), $values->get(4);

}