1 #======================================================================

     2 #                    . I S O P Y C N A L _ T S . S I G M A 4 

     3 #                    doc: Mon Dec 19 12:28:01 2005

     4 #                    dlm: Wed Sep 16 10:19:10 2009

     5 #                    (c) 2005 A.M. Thurnherr

     6 #                    uE-Info: 26 34 NIL 0 0 72 0 2 4 NIL ofnI

     7 #======================================================================

     8 

     9 # HISTORY:

    10 #	Sep 16, 2009: created

    11 

    12 require "$ANTS/libEOS83.pl";				# load equation of state

    13 

    14 unless (defined($P{ITS})) {

    15 	&antsInfo("using default %ITS=90");

    16 	&antsAddParams(ITS,90);

    17 }

    18 

    19 &antsInfo("Warning: ignoring non-zero %press = $P{press}")

    20 	if ($P{press} != 0);

    21 $P{press} = 0;

    22 

    23 sub density($$)

    24 {

    25 	my($S,$T) = @_;

    26 	return sigma($S,$T,$P{press},4000);

    27 }

    28 

    29 1;