#!/usr/local/bin/perl
# PBWIJ3   Perl script to test BWIJ3    99/03/09

  $HSOCON  = '/raid2/C089/A10/HSO.089';   #  units (cm)
  $HSOGHG  = '/raid2/C090/A10/HSO.090';   #  units (cm)
  $MOSICON = '/raid2/C089/A10/MOSI.089';  #  units (cm)
  $MOSIGHG = '/raid2/C090/A10/MOSI.090';  #  units (cm)

  open (BW,"| BWIJ3 $HSOCON $HSOGHG $MOSICON");
  print BW "D = Y5-X5 - Y2+X2\n";  #  D1 = dHSO 199X-196X Ann
  print BW "NY $MOSIGHG\n";        #  Y file is now $MOSIGHG
  print BW "D = Y5-Z5 - Y2+Z2\n";  #  D2 = dMOSI 199X-196X Ann
  print BW "D = D1 - 100.*D2\n";   #  D3 = dTHERMAL EXPANSION
  print BW "S 4 3 23 0\n";
  print BW "C DIAG1388\n";
  print BW "T\n";
  print BW "SEA LEVEL CHANGE due to THERMAL EXPANSIO".
           "N (cm)   GHG-Con   199X-196X DJF        \n";
  print BW "PP D3\n";
  print BW "Q\n";
  close(BW);