TOPCONS default

#!/usr/bin/perl -w

use SOAP::Lite;

my $soap_TopconsService = SOAP::Lite
    ->uri('topcons')
    ->proxy('http://wsdl.sbc.su.se/cgi-bin/topcons.cgi');


my $strseq = ">1a91A\nMENLNMDLLYMAAAVMMGLAAIGAAIGIGILGGKFLEGAARQPDLIPLLRTQFFIVMGLV\nDAIPMIAVGLGLYVMFAVA\n";

my $result = $soap_TopconsService -> runTopconsDefault($strseq);


# check for error
unless ($result->fault)
#unless (0)
{
  # check sanity of output
  my $scampi_msa = 0;
  my $scampi_seq = 0;
  my $prodiv = 0;
  my $pro = 0;
  my $octopus = 0;
  my $topcons = 0;
  my $zpred = 0;
  my $deltag = 0;
  my $rel = 0;

  my $output = $result->result();
  my @data = split(/\n/, $output);
#   open(DEBUG, "/tmp/topcons_tmp_t2gj/query.txt");
#   my @data = <DEBUG>;
#   close(DEBUG);
  for (my $i = 0; $i < $#data; $i++) {
      my $line = $data[$i];
      my $nextline = $data[$i+1];
      chomp($line);
      if ( $line =~ /^SCAMPI-seq/ && $nextline =~ /^[oiM]+/ ) {
	  $scampi_seq = 1;
      }
      elsif ( $line =~ /^SCAMPI-msa/ && $nextline =~ /^[oiM]+/ ) {
	  $scampi_msa = 1;
      }
      elsif ( $line =~ /^PRODIV predicted/ && $nextline =~ /^[oiM]+/ ) {
	  $prodiv = 1;
      }
      elsif ( $line =~ /^PRO predicted/ && $nextline =~ /^[oiM]+/ ) {
	  $pro = 1;
      }
      elsif ( $line =~ /^OCTOPUS/ && $nextline =~ /^[oiM]+/ ) {
	  $octopus = 1;
      }
      elsif ( $line =~ /^TOPCONS/ && $nextline =~ /^[oiM]+/ ) {
	  $topcons = 1;
      }
      elsif ( $line =~ /^Predicted Z-coordinates/ && $nextline =~ /^\d+\s+\d+\.\d+/ ) {
	  $zpred = 1;
      }
      elsif ( $line =~ /^Predicted Delta-G-values/ && $nextline =~ /^\d+\s+\d+\.\d+/ ) {
	  $deltag = 1;
      }
      elsif ( $line =~ /^Predicted TOPCONS reliability/ && $nextline =~ /^\d+\s+\d+\.\d+/ ) {
	  $rel = 1;
      }
  }
  print STDERR "Preds: $scampi_seq, $scampi_msa, $prodiv, $pro, $octopus, $topcons, $zpred, $deltag, $rel\n";
  # return 0 if everything ok
  if ( $scampi_seq && $scampi_msa && $prodiv && $pro && $octopus && $topcons && $zpred && $deltag && $rel ) {
      print "0\n";
  }
  else {
      print "1\n";
  }
} 
else 
{
  # report error code
  print join ', ',
    $result->faultcode,
      $result->faultstring,
	$result->faultdetail;
  print "1\n";
}