#======================================================================
# . I S O P Y C N A L _ T S . G A M M A _ N
# doc: Tue Dec 13 21:50:18 2005
# dlm: Mon Dec 19 13:04:15 2005
# (c) 2005 A.M. Thurnherr
# uE-Info: 27 30 NIL 0 0 72 0 2 4 NIL ofnI
#======================================================================
# HISTORY:
# Dec 14, 2005: - created
# Dec 19, 2005: - finalized
# NOTES:
# - requires %lat/%lon PARAMs
require "$ANTS/libgamma.pl"; # load equation of state
unless (defined($P{ITS})) {
&antsInfo("using default %ITS=90");
&antsAddParams(ITS,90);
}
$gamma::temp_scale = $P{ITS};
croak("$0: need %lat/%lon\n")
unless defined($P{lat}) && defined($P{lon});
croak("$0: need %press\n")
unless defined($P{press});
sub density($$)
{
my($S,$T) = @_;
my($gamma) = gamma::gamma_n($S,$T,$P{press},$P{lat},$P{lon});
return $gamma > 0 ? $gamma : undef;
}
1;