#!perl #
#
# Perl script to display the colour saturation of the snapshots
#
# Version: 5.50
# Date: 2001, March 27
#
# Usage:
#     perl fdmax.pl SEP DIV ROT PS
# Parameters:
#     SEP...  Name of the history file.
#     DIV...  Filename of the 4-D data cube of the "normalized" wavefield
#             divergence (P waves).
#     ROT...  Filename of the 4-D data cube of the X2 component of the
#             "normalized" wavefield rotation (S waves).
#     PS...   Output PostScript figure to be displayed and edited.
#             Refer to 'print' commands at the end of this file.
#
# ======================================================================
#
  @ARGVsave=@ARGV;
  @ARGV=();
  require 'go.pl';
  &CHK("forms/","sep.pl");
  require 'sep.pl';
  @ARGV=@ARGVsave;
#
  &CHK("fd/","fdmax.pro");
  &CHK(""   ,"$ARGV[0]" );
  &CHK(""   ,"$ARGV[1]" );
  &CHK(""   ,"$ARGV[2]" );
#
  &RSEP1($ARGV[0]);
  &RSEP3("N4"   ,$N4   ,1);
  &RSEP3("D4"   ,$D4   ,1);
  &RSEP3("SAMP" ,$SAMP ,1);
  &RSEP3("VSAT" ,$VSAT ,1);
  &RSEP3("VSAT4",$VSAT4,0);
#
  &ECHO( ">$ARGV[3]","%!");
  &ECHO(">>$ARGV[3]","/N4    $N4    def");
  &ECHO(">>$ARGV[3]","/D4    $D4    def");
  &ECHO(">>$ARGV[3]","/SAMP  $SAMP  def");
  &ECHO(">>$ARGV[3]","/VSAT  $VSAT  def");
  &ECHO(">>$ARGV[3]","/VSAT4 $VSAT4 def");
#
  open(LU,">fdmax.h");
  print LU "N1=111    N3=41     N4=321 \n";
  print LU "D1=1.000  D3=1.000  D4=0.004300 \n";
  print LU "O1=0.000  O3=0.000  O4=0.000 \n";
  print LU "N1NEW=1  N2NEW=1  N3NEW=1  GNORM=999. \n";
  print LU "GRD='$ARGV[1]'  GRDNEW='fdmax.tmp' \n";
  close (LU) || die "Error in writing history file fdmax.h";
  &RUN("grdnorm","'fdmax.h' /");
  &APPEND("$ARGV[3]","fdmax.tmp");
#
  open(LU,">>fdmax.h");
  print LU "GRD='$ARGV[2]'  GRDNEW='fdmax.tmp' \n";
  close (LU) || die "Error in writing history file fdmax.h";
  &RUN("grdnorm","'fdmax.h' /");
  &APPEND("$ARGV[3]","fdmax.tmp");
  &APPEND("$ARGV[3]","fdmax.pro");
#
  print "-----------------------------------------------------------------\n";
  print "Display and edit file '$ARGV[3]' to adjust the colour saturation!\n\n";
  print "Black line in figure '$ARGV[3]' corresponds to the maximum colour\n";
  print "saturation, green dots to P waves and blue dots to S waves.\n";
  print "Change SAMP to amplify S waves with respect to P waves,\n";
  print "VSAT to adjust the initial level of the saturation line and\n";
  print "VSAT4 to adjust its exponential decay.\n\n";
  print "When figure '$ARGV[3]' is OK, modify SAMP, VSAT and VSAT4 in the\n";
  print "history file accordingly!\n";
# ======================================================================
1;                                                               #