Decoders error

NOTE: The decoders mailing list is no longer active. The list archives are made available for historical reasons.

Dan,

I found the source of the decoders error message that you sent me awhile
back.  The sao2nc decoder was using the wrong literal. Attached is the
fixed version of sao2nc.

Robb...

==============================================================================
Robb Kambic                                Unidata Program Center
Software Engineer III                      Univ. Corp for Atmospheric Research
rkambic@xxxxxxxxxxxxxxxx                   WWW: http://www.unidata.ucar.edu/
==============================================================================
#! /usr/local/bin/perl
#
#  usage: sao2nc [cdlfile] netcdf
#
# This sao to netCDF decoder was created by using a sao decoder provided by
# Mark Albright as a starting point.  Robb Kambic added a netCDF wrapper to
# change the output from ascii to netCDF format.
#
use NetCDF ;

# process input parameters
if( $#ARGV == 0 ) {
        $ncfile = $ARGV[ 0 ] ;
} elsif( $#ARGV == 1 ) {
        $cdlfile = $ARGV[ 0 ] ;
        $ncfile = $ARGV[ 1 ] ;
} else {
        die "Wrong number of parameters " ;
}

# set interrupt handler
$SIG{ 'INT' }  = 'cleanup' ;
$SIG{ 'TERM' }  = 'cleanup' ;
$SIG{ 'KILL' }  = 'cleanup' ;
$SIG{ 'QUIT' }  = 'cleanup' ;

# default value
$F = -9999 ;

# change to home
chdir() ;
chdir( "/home/rkambic/code/decoders/perldec/marka" ) ;

# open or create ncfiles
if( -e $ncfile ) {
        $ncid = NetCDF::open( "$ncfile", WRITE ) ;
        $record_id = NetCDF::dimid( $ncid, "record" ) ;
        $name_id =  "xxxxxx"  ;
        $recnum =  -1  ;
        NetCDF::diminq( $ncid, $record_id, $name_id, $recnum ) ;

} else {
        die "Wrong or missing cdlfile parameter" unless  -e $cdlfile ;
        if( index( $ncfile, '/' ) != -1 ) {
                $dir = substr( $ncfile, 0, rindex( $ncfile, '/' ) ) ;
                if( ! -e $dir ) {
                        system( "mkdir -p $dir" ) ;
                }
        }

        if( -e "util/ncgen" ) {
                $ncgen = "util/ncgen" ;
        } elsif( -e "/usr/local/ldm/util/ncgen" ) {
                $ncgen = "/usr/local/ldm/util/ncgen" ;
        } elsif( -e "/upc/netcdf/bin/ncgen" ) {
                $ncgen = "/upc/netcdf/bin/ncgen" ;
        } elsif( -e "./ncgen" ) {
                $ncgen = "./ncgen" ;
        } else {
                open( NCGEN, "which ncgen |" ) ;
                $ncgen = <NCGEN> ;
                close( NCGEN ) ;

                if( $ncgen =~ /no ncgen/ ) {
                        die "Can't find NetCDF utility 'ncgen' in PATH, 
util/ncgen
        /usr/local/ldm/util/ncgen, /upc/netcdf/bin/ncgen, or ./ncgen : $!\n" ;
                } else {
                        $ncgen = "ncgen" ;
                }
        }

        system( "$ncgen -o $ncfile $cdlfile" ) ;
        $ncid = NetCDF::open( "$ncfile", WRITE ) ;
        # NetCDF record counter
        $recnum = 0 ;
}

# set fill mode
NetCDF::setfill($ncid, NetCDF::NOFILL) == 0 || die "Couldn't set fill mode\n";

select( STDOUT ) ;

# read in station data
if( -e "etc/city_sfc.dat" ) {
        $sfile = "etc/city_sfc.dat" ;
} elsif( -e "./city_sfc.dat" ) {
        $sfile = "./city_sfc.dat" ;
} else {
        die "Can't file city_sfc.dat station file" ;
}
open( STATION, "$sfile" ) || die "could not open $sfile: $!\n" ;

while( <STATION> ) {
        @station = split ;
        @station = reverse( @station ) ;
        $station[ 4 ] =~ /(\d)/ ;
        
        # set these vars ( $lat, $lon, $elev, $priority, $wmo_id ) 
        $STATIONS{ "$station[ 5 ]" } = 
                "$station[ 3 ] $station[ 2 ] $station[ 1 ] $1 $station[ 0 ]" ;
}
close STATION ;

$cover{'SCT'} = 3 ;
$cover{'BKN'} = 6 ;
$cover{'OVC'} = 8 ;
$cover{'0VC'} = 8 ;     # Occasionally observers misspell OVC

$pre1995 = "" ;         # Hard wired to process sao obs from 1995

# Now begin parsing file and decoding observations breaking on cntrl C
$/ = "\cC" ;

# set select processing here
while( 1 ) {
        open( STDIN, '-' ) ;
        vec($rin,fileno(STDIN),1) = 1;
        $timeout = 600 ;
        $nfound = 0 ;
        $nfound = select( $rout = $rin, undef, undef, $timeout );
        # timed out
        if( ! $nfound ) {
                #print "nfound = 0 \n" ;
                &cleanup() ;
        }

        &cleanup() if( eof( STDIN ) ) ;

   $_ = <STDIN> ;     # Process each bulletin
   s#\cM##g ;
   s#\cC##g ;

   if( s/^\s*\001\s*\d\d\d\s+S\w(\w\w)\d{1,2} [CK]\w\w\w \d{2}(\d{2})(\d{2})( 
RTD\d\d)?\s+//s ) {
      # Ignore bulletins of Canadian and Mexican hourlies
      #print "$&\n" ;
      $header_hr  = $2 ;
      $header_min = $3 ;
   }
   else {
      #print STDERR "Couldn't match header.\n" ;
      #print STDERR "$_\n\n" ;
      next ;      # Advance to the next bulletin
   }

   # Usually reports are separated by a RS character
   if( /\036/ ) {
      s/[ \n]RSNK /\036RSNK / ;         # Separate RSNK from the HMS report   
      @reports = split ( /\036/ ) ;
   }
   elsif( /=\n/ ) {
#      s/[ \n]RSNK /=RSNK / ;            # Separate RSNK from the HMS report
      @reports = split( /=\n/ ) ;
   }    
   else {
#      s/[ \n]RSNK /\nRSNK / ;           # Separate RSNK from the HMS report
      @reports = split ( /\n/ ) ;
   }

   for ( @reports ) {                   # Process each report in the bulletin  
      $raw = $lat = $lon = $elev = $priority = "" ;
      $id = $type = $hr = $min = $ceiling = $skyc = $vis = $wea = $pres =
      $t = $td = $wdir = $wspd = $wgst = $alt = $ptend = $deltap =
      $pcp1 = $pcp3_6 = $xt = $pcp24 = $snow = $tmax = $tmin = $comment = "" ;
      $ptend_found = $read_cloud = $snowfall = $snowh2o = "" ;
      $asos = $ramos = $amos = $awos = $correction = $dashx = "" ;
      $one_inch = $two_inch = $three_inch = "" ;

      # save original undecoded report
      $raw = $_ ;

              # Remove multiple spaces, newlines and trailing = from report
      s#[\s=\$]+$## ;    s#/\n\s*#/ # ;   s#\n# #g ;   s#\s{2,}# #g ;
#      s#\s*\n\s*/#/#g ; 

      if( /FINO/ ) { next ; }  
      $hold_report = $_ ;

      if( s/^(\w{3,4}) (SA|RS|SP|USP) (COR |RTD )?([0-2]\d)([0-5]\d) (COR |RTD 
)?// ) {
         $id   = $1 ;
#      if( !($id eq "CNK" || $id eq "ALS" || $id eq "DDC" || $id eq "GLD" )) {
#      next ; }           # These 2 lines used to check ASOS decoding
         $type = $2 ;
         $correction = $3 ;
         $hr   = $4 ;
         $min  = $5 ;
      }
      elsif( s/^RSNK // ) {   # Treat Rattlesnake Mtn. as special case
         $id = "RSN" ;
         $type = "SA" ;
         $hr = $header_hr ;
         $min = $header_min ;
         if( s#^([-\d]+)/(\d\d)(\d\d)(G(\d+))?## ) {
            $t = $1 ;
            $wdir = $2 ;
            $wspd = $3 ;
            $wgst = $5 ;
         }
      }
      else { next ; }

      if( $id eq "LAX" ) {  # Handles city temperature appended at end of report
         if( s#/ ?(CITY\s*[\d/]{2,3})## ) {
            $comment = " $1" ;
         }
      }

#       Extract station's latitude, longitude, elevation, priority.

        ( $lat, $lon, $elev, $priority, $wmo_id ) = 
                        split( ' ', $STATIONS{ "$id" } ) ;


      if( $lat eq "" ) {        # Station not found in directory
         #print STDERR "Couldn't find $id in station directory.\n" ;
         next ;
      }

      if( s/^AUTOB // ) {  # for now skip AUTOB reports for SDB, DRT, and INW
         next ;
      }

      if( s/^AMOS // ) {  # Stations such as SMP
         $amos = "TRUE" ;
      }

      if( s#^RAMOS /## ) { # Process RAMOS reports somewhat differently
         $ramos = "TRUE" ;
      }

      if( s/^A[O0]2A? // ) { # Process ASOS reports somewhat differently
        $asos = "TRUE" ;
        if( s/CLR BLO 120 // ) {
           $ceiling = "888" ;
        }
        s/^MM // ;
      }

      if( s/^AWOS // ) {   # Process AWOS reports somewhat differently
         $awos = "TRUE" ;
         if( s/CLR BLO 120 // ) {    # Sky cover is unknown for this case
#            $ceiling = ">120" ;
            $ceiling = "888" ;         # Above value confusing to fortran io 
         }
         if( s/^M // ) {        # sky cover info missing
           $ceiling = "M" ;
         }
         if( s/SLP (\d\d\d)// ) {      # OTH reports in this manner
            $pres = $1 ;
         }
      }

      if( s/^-X // ) {
         $dashx = "TRUE" ;
         $skyc = 0 ;
      }

      if( s/^W ?(\d{1,2}) ?X // ) {        # Sky obscured
         $skyc = 9 ;
         $ceiling = $1 ;
      }
      elsif( s/^CLR // ) {       # Clear
         $skyc = 0 ;
         $ceiling = 999 ;
      }
      else {
         # Capture each cloud layer
         while ( s/^[ME]?(\d{1,3})V? ?(-?)(SCT|BKN|[0O]VC) // ) {
            $thin = $2 ;
            if( ( $3 eq "BKN" || $3 eq "OVC"  || $3 eq "0VC" ) && !$ceiling ) {
               if( !$thin ) { $ceiling = $1 ; }
            }
            $skyc = $cover{$3} ;
            $read_cloud = "TRUE" ;
         }
         if( !$ceiling && ( $read_cloud || $dashx ) ) {
            $ceiling = 999 ;
         }
      }

   #This section assumes sky info accompanied by vis and pres weather
      if( $ceiling ne "" ) {
         if( $ceiling eq "M" ) { $ceiling = "" ; } 
         # Now decode visibility and present weather
         if( s#^((<?\d{1,2}/\d{1,2})|(\d{1,3}\+?)) ?([A-Z+-]*) ## ) {
            $vis = $1 ;
            $wea = $4 ;
         } 
         if( s#^(\d)/(\d{1,2})## ) {  # Check whether more visibility follows
            if( $2 > $1 ) {            # Otherwise must be T and Td 
               $vis .= $1 . "/" . $2 ;
               if( s#^([A-Z]+[A-Z+-]*) ## ) {
                  $wea = $1 ;
               }
               else {
                  s#^ ## ;     # If no present weather then remove leading space
               }
            }
            else {           # T and Td are both single digits
               $t = $1 ;
               $td= $2 ;
               if( s#^/E?(\d\d)(\d\d)(G(\d+))?/## ) {
                  $wdir = $1 ;
                  $wspd = $2 ;
                  $wgst = $4 ;
               }
            }
         }
      }
      $wea =~ s/^V// ;  # Remove reference to variable visibility from $wea
         
      if( $vis =~ m#^<1/4# ) {        # AWOS stations report in this way
         $vis = "1/8" ;
      }
#  Use this code to change visibility to decimal notation if needed
      if( $vis =~ m#^(\d?)(\d)/(\d{1,2})# ) {       # Fractional visibility
         if( $3 != 0 ) { $vis = $1 + $2/$3 ; }
         else { $vis = "" ; }
      }

      # Regular(SA or RS) obs on the hour, or amos/asos report
      if( $type ne "SP" || $amos  || $asos || $awos ) { 
         if( s#^E?(\w+)/ ?([\w-]+)/([\w-]+)/E?(\d\d)(\d\d)(G(\d+)| M)?/## ){ 
            $pres = $1 ;
            $t    = $2 ;
            $td   = $3 ;
            $wdir = $4 ;
            $wspd = $5 ;
            $wgst = $7 ;
         }
         elsif( s#^ ?([\w-]+)/([\w-]+)/E?(\d\d)(\d\d)(G(\d+)| M)?/## ){ 
            $t    = $1 ;
            $td   = $2 ;
            $wdir = $3 ;
            $wspd = $4 ;
            $wgst = $6 ;
            $td =~ s/M//g ;     # Check for missing value
         }
         else {
            $skyc = $ceiling = $vis = $wea = "" ;
         }
         if( $pres =~ /[a-zA-Z]/ || length( $pres ) != 3 ) { $pres = "" ; }
         if( $t    =~ /[a-zA-Z]/ ) { $t  = "" ; }
         if( $td   =~ /[a-zA-Z]/ ) { $td = "" ; }         
      }
      else {                     # Special report(SP)
         if( s#^E?(\d\d)(\d\d)(G(\d+))?/## ) {
            $wdir = $1 ;
            $wspd = $2 ;
            $wgst = $4 ;
         }
      }

      if( $wdir ne "" ) {
         $wdir *= 10 ;
         if( $wdir == 0 ) { $wdir = "00" ; }
      }

      if( $ramos && s#PK WND (\d{2,3})## ) {
         $wgst = $1 ;
         $wgst =~ s/^0(\d)$/$1/ ;
      }

      if( $amos && s/PK WND (\d{2,3}) \d\d\d// ) {
         unless ( $wgst ) { 
            $wgst = $1 ;
#            $wgst =~ s/^0(\d\d)/$1/ ;
         }
      }

      if( $wspd ne "" ) { $wspd *= 1 ; }
      if( $wgst ne "" ) { $wgst *= 1 ; }

      if( $wspd ne "" ) { # If wind was found then search for altimeter setting
         if( s#^(\d\d\d|M)/?# # ) {    # Altimeter setting
            $alt = $1 ;   $alt =~ s#M##g ;
         }
      }

      if( $asos ) {
         if( s/ PCPN (\d\d\d\d|M)( |$)/ / ) {
            $pcp1 = $1 ;
            if( $pcp1 eq "0000" ) { $pcp1 = " T" ; }
            if( $pcp1 =~ /^\d+$/ ) {
               $pcp1 *= 1 ;   # Remove leading zeroes
               if( $pcp1 < 10 ) { $pcp1 = "0" . $pcp1 ; }
            }
         }
         else { $pcp1 = " 0" ; }
      }
         
      if( s/ NOSPL// ) {      # no specials(SP) between hours
         $comment .= $& ;
      }

      if( $hr == 7 && $min > 45 ) {
         if( s# 98(\d\d\d)([ /]|$)# # ) {
            $sunshine = $1/60. + .05 ;
            $sunshine =~ s/(\.\d)\d*/$1/ ;
            $comment = " sunshine=$sunshine hrs" . $comment ;
         }
      }
      if( $hr >= 4 && $hr <= 11 ) {   # Midnight Local Standard Time
         if( s# 4(\d)(\d{3})(\d)(\d{3})([ /]|$)# # ) {
            $tmax = $1 ? -$2*1 : $2*1 ;
            $tmin = $3 ? -$4*1 : $4*1 ;
            $comment .= " Mx=$tmax Mn=$tmin" ;
            if( $id eq "RSN" ) { $comment = "Hanford:" . $comment ; }
         }
      }

      $hourly_3 = $hourly_6 = $hourly_12 = "" ;
      if( ($hr+1)%3 == 0  &&  $min>20  &&  $type ne "SP" ) {
         $hourly_3 = "TRUE" ;
         if( ( $hr+1 )% 6 == 0 ) { $hourly_6  = "TRUE" ; }
         if( ( $hr+1 )%12 == 0 ) { $hourly_12 = "TRUE" ; }
      }
      elsif( $hr%3 == 0  &&  $min<30  &&  $type ne "SP" ) {
         $hourly_3 = "TRUE" ;
         if( $hr% 6 == 0 ) { $hourly_6  = "TRUE" ; }
         if( $hr%12 == 0 ) { $hourly_12 = "TRUE" ; }
      }

  # Process only 3-hourly SA or RS reports in the next block
      if( $hourly_3 && !$ramos ) {

         if( $hourly_12 ) {   # Process only 12-hourly SA/RS reports here
            if( s/(( [A-Z]{3})? RADAT \w+)// ){
               $comment .= $1 ;
            }
         }

         if( $asos || ! $pre1995 ) {
            if( s# 5(\d)(\d\d\d)([ /]|$)# # ) {
               $ptend = $1 ;
               $deltap = $2 ;
               $deltap =~ s/0(\d\d)/$1/ ;
               $ptend_found = "TRUE" ;
            }

            if( s# 6(\d\d\d)/([ /]|$)# # ) {
               $pcp3_6 = $1 ;
               if( $pcp3_6 eq "000" ) { $pcp3_6 = "T" ; }
               $pcp3_6 =~ s/0(\d\d)/$1/ ;
            }
            elsif( s# 6////( |$)# # ) { $pcp3_6 = "" ; }
            elsif( $ptend_found ) { $pcp3_6 = "0" ; }

            if( $hr == 11 ) {   # 12 GMT only
               if( s# 7(\d{4})([ /]|$)# # ) {
                  $pcp24 = $1 ;
                  $pcp24 *= 1 ;
                  if( $pcp24 < 10 ) { $pcp24 = "0" . $pcp24 ; }
               }
               elsif( s# 7////([ /]|$)# # ) { $pcp24 = "" ; }
               elsif( $ptend_found ) { $pcp24 = "0" ; }
            }

            if( s# 1(\d)(\d\d\d)([ /]|$)# # ) {
               $tmax = $1 ? -$2*1 : $2*1 ;
            }
            else { s# 1////([ /]|$)# # ; }

            if( s# 2(\d)(\d\d\d)([ /]|$)# # ) {
               $tmin = $1 ? -$2*1 : $2*1 ;
            }
            else { s# 2////([ /]|$)# # ; }

            if( $hr <= 6 || $hr == 23 ) {  # 00 and 06 GMT
               $xt = $tmax ;
            }
            else { $xt = $tmin ; }      # 12 and 18 GMT
            
            if( s# 8/([\d/]{3})([ /]|$)# # ) {
               $comment .= " C$1" ;
            }
            
#            if( s/ 901(\d\d)( |$)/ / ) { # New snowfall in inches
#               $snowfall = $1 ;
#            }

            if( s# 933(\d\d\d)([ /]|$)# # ) { # Water equivalent snow depth in
               $snowh2o = $1/10. ;         # tenths of inch
            }
            if( $snowh2o ) {
               $snowh2o += .05 ;
               $snowh2o =~ s/(\.\d)\d*/$1/ ;
            }

            if( s# 4/(\d\d\d)([ /]|$)# # ) {   # Snow depth in inches
               $snow = $1 ;
               $snow =~ s#^0+## ;
            }
         }              # end of ASOS processing
         else {         # non ASOS stations
            if( s/ 901(\d\d)( |$)/ / ) { # New snowfall in inches
               $snowfall = $1 ;
            }
            if( s/ 902(\d\d)( |$)/ / ) { # Water equivalent snow depth in
               $snowh2o = $1/10. ;         # tenths of inch
            }
            if( s/ 903(\d\d)( |$)/ / ) { # Water equivalent snow depth in
               $snowh2o += $1 ;            # whole inches
            }
            if( $snowh2o ) {
               $snowh2o += .05 ;
               $snowh2o =~ s/(\.\d)\d*/$1/ ;
            }
            while ( s/ 904(\d\d)( |$)/ / ) {   # Snow depth in inches
               if( $1 == 99 ) { $snow = $snow + 100 ; }
               else { $snow = $snow + $1 ; }
            }   
                # probably Max/Min temp and 24 hr precipitation
            if( s/ 4(\d\d)(\d\d)( 2(\d\d\d\d))?( |$)/ / ) {
               if( $1 eq "00" ) {            # found ptend/precip group instead
                  $ptend = "4" ;
                  $deltap = $1 ;
                  $pcp3_6 = $2 ;
                  $ptend_found = 1 ;
               }
               else {                           # found Max and Min temperature
                  $tmax = $1 ;
                  $tmin = $2 ;
                  if( $pcp24 = $4 ) {
                     $pcp24 *= 1 ;               # Remove leading zeroes ;
                     if( $pcp24 < 10 ) { $pcp24 = "0" . $pcp24 ; }
                  }
               }
            }

            if( s# [A-Z]{1,2}R[\d/]{2}[ \w]*$## ) { # Runway condition
               $comment .= $& ;
            }

#            s#[\s/]+$## ; # Remove trailing slashes such as found on OTH report

            if( $hourly_12 ) {
               if( s/ 2(\d\d\d\d)\s*$// ) {    # 24 hr precipitation
                  $pcp24 = $1 ;
                  $pcp24 *= 1 ;               # Remove leading zeroes ;
                  if( $pcp24 < 10 ) { $pcp24 = "0" . $pcp24 ; }
               }
            }

            if( s# 1\d[\d/][\d/]( |$)# # ) {  # Cloud types
               $comment = $& . " " . $comment ;
            }

            if( s/ (\d{1,2})\s*$/ / ) {
                $xt = $1 ;
            }

            if( s/ 99(\d\d\d)\s*$/ / ) { # Pressure change exceeds 9.9 mb
               $deltap = $1 ;
            }
 
            if( s/ ONE\s*$/ / )   { $one_inch = "TRUE" ; }
            elsif( s/ TWO\s*$/ / )   { $two_inch = "TRUE" ; }
            elsif( s/ THREE\s*$/ / ) { $three_inch = "TRUE" ; }

            if( !$ptend_found && s#[ /](\d)(\d\d)(\d\d)?\s*$# # ) {
               $ptend  = $1 ;
               if( !$deltap ) { $deltap = $2 ; }
               $pcp3_6 = $3 ;
               if( $one_inch )   { $pcp3_6 += 100 ; }
               elsif( $two_inch )   { $pcp3_6 += 200 ; }
               elsif( $three_inch ) { $pcp3_6 += 300 ; }
               if( !$pcp3_6 ) {
                  $pcp3_6 = "0" ;
               }
            # Military stations do not report precip at 03, 09, 15, and 21GMT.
               if( !$hourly_6 ) {
                  if( $id eq "GRF"||$id eq "TCM"||$id eq "NUW"||$id eq "SKA" ) {
                     $pcp3_6 = "" ;
                  }
               }
            }
            if( $pcp3_6 eq "00" ) { $pcp3_6 = "-1" ; } # Trace precip coded as 
-1
         }
      }                  # end of 3-hourly processing

      # Clean up in preparation for adding to comment field
      s#^[/ ]*## ;  s#[/ ]*$## ; s# / # #g ;  s#\s+# # ;
      s# #_#g ; s#,#;#g ;  # Substitute ; for , and underscore for spaces

      $comment = $_ . $comment ;

      if( $snowh2o ) {
         $comment = "snowH2O=$snowh2o " . $comment ;  
      }
  
      $comment = $correction . $comment ;
      $comment =~ s#^\s+## ;
      $comment = $pcp1 . " " . $comment ;
#      if( $asos ) {
#         $pcp1str = sprintf ( "%3.3s ", $pcp1 ) ;
#         $comment = $pcp1str . $comment ;
#      }

      # Needed for using with list directed fortran io
      #$id = "\"" . $id . "\"" ;

      $t =~ s/(-?)(0+)([1-9]+)/$1$3/ ;  
      $td =~ s/(-?)(0+)([1-9]+)/$1$3/ ;  

      # expand abreviations and partial numbers
      if( $alt && $alt < 500 ) {
        $alt = "3" . "$alt" ;
        $alt /= 100 ;
      } elsif( $alt && $alt > 500 ) {
        $alt = "2" . "$alt" ;
        $alt /= 100 ;
      }
      if( $pres && $pres < 600 ) {
        $pres = "10" . "$pres" ;
        $pres /= 10 ;
      } elsif( $pres && $pres > 600 ) {
        $pres = "9" . "$pres" ;
        $pres /= 10 ;
      }

if( 0 ) {
      $report_dec = join(',', $id, $type, $priority, $hr, $min, $lat, $lon,
                               $elev, $ceiling, $skyc, $vis, $wea, $pres, $t,
                               $td, $wdir, $wspd, $wgst, $alt, $ptend, $deltap,
                               $pcp3_6, $xt, $pcp24, $snow, $tmax, $tmin,
                               $comment ) ;

      print "\n---$hold_report---\n" ;
      print OUT $report_dec, "\n" ;
}

        # set defaults for NetCDF
        $hr = $F unless $hr ;
        $t = $F unless $t ;
        $td = $F unless $td ;
        $pres = $F unless $pres ;
        $wdir = $F unless $wdir ;
        $wspd = $F unless $wspd ;
        $wgst = $F unless $wgst ;
        $alt = $F unless $alt ;
        $vis = $F unless $vis ;
        $wea = $F unless $wea ;
        $ptend = $F unless $ptend ;
        $deltap = $F unless $deltap ;
        $pcp1 = $F unless $pcp1 ;
        $pcp3 = $F unless $pcp3 ;
        $pcp6 = $F unless $pcp6 ;
        $pcp24 = $F unless $pcp24 ;
        $tmax = $F unless $tmax ;
        $tmin = $F unless $tmin ;


        # output the NetCDF data here
        $datap[ 0 ] = \$wmo_id ;
        #&charord( \@id, $id, 4 ) ;
        &padstr( \$id, 4 ) ;
        $datap[ 1 ] = \$id ;
        $datap[ 2 ] = \$lat ;
        $datap[ 3 ] = \$lon ;
        $datap[ 4 ] = \$elev ;
        $datap[ 5 ] = \$hr ;
        $datap[ 6 ] = \$hr ;
        $datap[ 7 ] = \$t ;
        $datap[ 8 ] = \$td ;
        $datap[ 9 ] = \$pres ;
        $datap[ 10 ] = \$wdir ;
        $datap[ 11 ] = \$wspd ;
        $datap[ 12 ] = \$wgst ;
        $datap[ 13 ] = \$alt ;
        $datap[ 14 ] = \$vis ;
#print "$wea\n" ;
        @wea = ( $wea, 
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0) ;
        $datap[ 15 ] = \@wea ;
        &padstr( \$stn_type, 1 ) ;
        $datap[ 16 ] = \$stn_type ;
        &padstr( \$type, 2 ) ;
        $datap[ 17 ] = \$type ;
        $datap[ 18 ] = \$ptend ;
        $datap[ 19 ] = \$deltap ;
        $datap[ 20 ] = \$pcp1 ;
        $datap[ 21 ] = \$pcp3 ;
        $datap[ 22 ] = \$pcp6 ;
        $datap[ 23 ] = \$pcp24 ;
        $datap[ 24 ] = \$tmax ;
        $datap[ 25 ] = \$tmin ;
        &padarr( \@cc, 5, 8 ) ;
        $datap[ 26 ] = \@cc ;
        $datap[ 27 ] = \[ 0, 0, 0, 0, 0 ] ;
        &padarr( \@cloudtype, 5, 1 ) ;
        $datap[ 28 ] = \@cloudtype ;
        #&padstr( \$comment, 256 ) ;
        #$datap[ 29 ] = \$comment ;
        &padstr( \$raw, 256 ) ;
        $datap[ 29 ] = \$raw ;

        $result = NetCDF::recput( $ncid, $recnum, [ @datap ] ) ;
        #$result = NetCDF::recput( $ncid, $recnum, [ \$wmo_id ] ) ;
        #print STDOUT "NetCDF::recput result = $result\n" ;

        $recnum++ ;
        #last ;

   } # end of for ( @reports ) loop

   &cleanup() if( eof( STDIN ) ) ;

} # end while( 1 )
&cleanup();
exit( 0 ) ; #should never get here

sub cleanup
{

local( $sig ) = @_ ;

#print "Caught a SIG$sig --shutting down\n" if( $sig ) ;
$result = NetCDF::close( $ncid ) ;
#print STDOUT "NetCDF::close result = $result\n" ;

#close OUT ;      
exit( 0 ) ;

}

# pad str to correct length
sub padstr
{

( $str, $len ) = @_ ;

local( $size, $i ) ;

$size = length( $$str ) ;

for( $i = $size; $i < $len; $i++ ) {
        $$str .= "\0" ;
        #print "$$str,\n" ;
}

}

# pad arr to correct length
sub padarr
{

( $arr, $x, $y ) = @_ ;

local( $size, $i, $j ) ;

for( $i = 0; $i < $x; $i++ ) {
        $size = length( $$arr[ $i ] ) ;
        for( $j = $size; $j < $y; $j++ ) {
                $$arr[ $i ] .= "\0" ;
                #print STDOUT ",$$arr[ $i ],\n" ;
        }
}
}

  • 1999 messages navigation, sorted by:
    1. Thread
    2. Subject
    3. Author
    4. Date
    5. ↑ Table Of Contents
  • Search the decoders archives: