.isopycnal_TS.sigma4
author A.M. Thurnherr <athurnherr@yahoo.com>
Sat, 24 Jul 2021 09:38:16 -0400
changeset 46 70e566505a12
parent 39 56bdfe65a697
permissions -rw-r--r--
V7.3
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 4 
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: Wed Sep 16 10:19:10 2009
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: 26 34 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
#	Sep 16, 2009: 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},4000);
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;