#!/usr/local/bin/perl

 $pi= 3.141592654;
 $co= $pi/180.;
 $erad= 6370.8;    # earth radius in kilometers
 $arad= 6378.139;  # earth radius at equator in kilometers
 $brad= 6356.75;   # earth radius an north pole in kilometers

 $narg= @ARGV;

 open(FILE2, ">gcen");
  while (@ARGV) {
     $dat1= .0; $dat2= .0;
     $dat8= .0; $dat9= .0;
 $maximum1= .0;
 $minimum1= 1000000.e3;
    $file= shift @ARGV;
   $file1= substr($file, 8, 3);
   $file1= lc ($file1); open(FILE, $file);
    open(FILE1, ">in/$file1");
    @lines= <FILE>;
   $nlines= @lines;
     @word= split(' ', $lines[0]);
      $lon= - ( $word[1] );
     @word= split(' ', $lines[1]);
     $lat1= $word[1];
      $lat= atan2( $brad**2/$arad**2 *sin($lat1 *$co)/ cos($lat1 *$co), 1. )/$co;
       $ix= 0;
      $sum= 0;
    for ($i= 2; $i < $nlines - 1;  $i++) {
       $line= $lines[$i];
      $month= substr($line, 0, 2);
       $year= substr($line, 3, 4);
       $time= $year + ($month + .5) /12.;
       @word= split(' ', $line);
       $dat1= $word[1] *1.2335;
        $dat= $word[1] *1.2335 /1.e6;
#      if ($dat > .0) {
           printf FILE1 "%9.4f %23.15e %23.15e  %3s_  lat  96jan01  \n", $time, $dat, $dat, $file1;
#          printf FILE1 "%9.4f %23.15e %23.15e  %3s_  lat  96jan01  \n", $time, $dat,    .1, $file1;
           $ix++;
           $sum= $sum + $dat1;
#         }
       if ($line =~ /^10\/2011/) {
          $dat8= $dat1;
          }
       if ($line =~ /^10\/2015/) {
          $dat9= $dat1;
          }
       if ($dat1 > $maximum1) {
          $maximum1= $dat1;
          }
       if ($dat1 < $minimum1) {
          if ($line =~ /^10\/2017/) { }
          else {
             $minimum1= $dat1;
             }
          }
       }
   if ($ix > 0) {
      $mean= $sum /$ix;
      }
#  printf "%s  %9.3f %9.3f  %4d  %9d\n", $file1, $lat1, $lon, $ix, $mean ;
#  if ($dat8 > 0 ) { if ($dat9 > 0) {
   printf " %s %8.3f %9.3f  %4d  %5d %5d %5d %6d %6d\n", $file1, $lat1, $lon, $ix, $maximum1 /1.e3, $mean /1.e3, $minimum1 /1.e3, ($minimum1 - $maximum1) /1.e3, ($dat9 - $dat8) /1.e3 ;
   $sum8= $sum8 + $dat8;
   $sum9= $sum9 + $dat9;
#  } } 
   printf FILE2 "%s__res %15.9f %15.9f\n", $file1, $lat, $lon;
   }
   printf "total                                           %6d\n", ($sum9 - $sum8) /1.e3 ;