.isopycnal_TS.sigma0
author Andreas Thurnherr <ant@ldeo.columbia.edu>
Mon, 13 Apr 2020 11:06:22 -0400
changeset 40 c1803ae2540f
parent 39 56bdfe65a697
permissions -rw-r--r--
.
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
39
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     1
#======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     2
#                    . I S O P Y C N A L _ T S . S I G M A 0 
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     3
#                    doc: Mon Dec 19 12:28:01 2005
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     4
#                    dlm: Mon Dec 19 13:27:24 2005
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     5
#                    (c) 2005 A.M. Thurnherr
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     6
#                    uE-Info: 19 56 NIL 0 0 72 0 2 4 NIL ofnI
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     7
#======================================================================
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     8
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
     9
# HISTORY:
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    10
#	Dec 19, 2005: - created
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    11
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    12
require "$ANTS/libEOS83.pl";				# load equation of state
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    13
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    14
unless (defined($P{ITS})) {
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    15
	&antsInfo("using default %ITS=90");
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    16
	&antsAddParams(ITS,90);
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    17
}
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    18
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    19
&antsInfo("Warning: ignoring non-zero %press = $P{press}")
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    20
	if ($P{press} != 0);
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    21
$P{press} = 0;
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    22
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    23
sub density($$)
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    24
{
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    25
	my($S,$T) = @_;
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    26
	return sigma($S,$T,$P{press},0);
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    27
}
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    28
A.M. Thurnherr <athurnherr@yahoo.com>
parents:
diff changeset
    29
1;