#!/usr/bin/perl

if($#ARGV != 0)
{
   $prgName = `basename $0`;
   chomp($prgName);
   die <<EOHELP;

Generates a BCL load script based on the family file.

Usage: $prgName [all|fam]

EOHELP
}

$type=$ARGV[0];

@files = `cat family.$type`;
#@files = `ls dimrmd.struct???.iupac.pdb`;
#@files = (
#          "dimrmd.struct002.iupac.pdb",
#          "dimrmd.struct003.iupac.pdb",
#          "dimrmd.struct005.iupac.pdb",
#          "dimrmd.struct006.iupac.pdb",
#          "dimrmd.struct008.iupac.pdb",
#          "dimrmd.struct015.iupac.pdb",
#         );

$firstTime = 1;
foreach $file (@files)
{
   chomp($file);
   $name = $file;
   $name =~ s/.*f36g//;
   $name =~ s/\..*//;
   $head = "r";
   $head = "m" if $firstTime;
   print <<EOSECTION;
$head:Get Molecule PDB User $file f36g_$name -Heteroatom -Reference_Object
EOSECTION
   $firstTime = 0;
}

$firstTime = 1;
$firstHead = 1;
foreach $file (@files)
{
   chomp($file);
   $name = $file;
   $name =~ s/.*f36g//;
   $name =~ s/\..*//;
   if($firstTime)
   {
      $firstName = $name;
      $firstTime = 0;
   }
   else
   {
      $head = "r";
      $head = "m" if($firstHead);
      print <<EOSECTION;
$head:Superimpose -End_Definition Backbone "Label Mode" f36g_$name:3-15,25-35,46-53,63-73 f36g_$firstName:3-15,25-35,46-53,63-73
r:Superimpose End_Definition
EOSECTION
      $firstHead = 0;
   }
}

print <<EOSECTION;
m:Display Molecule Off Atoms Specified *
r:Display Molecule On Atoms Backbone * 
m:Color Molecule Atoms *:1-15 Specified Specification 255,0,0
r:Color Molecule Atoms *:16-24 Specified Specification 0,255,255
r:Color Molecule Atoms *:25-35 Specified Specification 0,255,0
r:Color Molecule Atoms *:36-45 Specified Specification 255,0,255
r:Color Molecule Atoms *:46-54 Specified Specification 0,0,255 
r:Color Molecule Atoms *:55-62 Specified Specification 255,255,0
r:Color Molecule Atoms *:63-75 Specified Specification 255,255,255
EOSECTION

$firstTime = 1;
$firstHead = 1;
LINE:
foreach $file (@files)
{
   chomp($file);
   $name = $file;
   $name =~ s/.*f36g//;
   $name =~ s/\..*//;
   if($firstTime)
   {
      $firstName = $name;
      $firstTime = 0;
   }
   else
   {
      print <<EOSECTION;
m:Associate Assembly f36g_$firstName f36g_$name f36g
EOSECTION
      last LINE;
   }
}

$firstHead = 1;
foreach $file (@files)
{
   next if(++$count < 3);
   chomp($file);
   $name = $file;
   $name =~ s/.*f36g//;
   $name =~ s/\..*//;
   $head = "r";
   $head = "m" if($firstHead);
   print <<EOSECTION;
$head:Add Assembly f36g_$name f36g
EOSECTION
   $firstHead = 0;
}

print "m:DepthCue Object * -DepthCueing\n";
print "m:Clip Auto Molecule_Objects Contour_Objects Grid_Objects ";
print "Graph_Objects User_Objects * 10 -Perspective INSIGHTII\n";