#======================================================================

#                    L I B N O D C . P L 

#                    doc: Mon Aug 28 11:07:47 2000

#                    dlm: Sun Jul  2 00:16:26 2006

#                    (c) 2000 A.M. Thurnherr

#                    uE-Info: 117 0 NIL 0 0 72 2 2 4 NIL ofnI

#======================================================================

# HISTORY:

#	Aug 28, 2000: - created

#	Sep 05, 2000: - allow spaces instead of 0es in lat/lon to accomodate

#					for Talley's OCEANUS 24S files

#	Oct 16, 2000: - added &DD[D]MMSSh2d()

#	Feb 28, 2001: - changed &depth to &obs_depth to remove clash with

#					[libEOS83]

#	Aug  1, 2001: - BUG: obs() could not handle Reid and Mantyla -ve values

#					correctly (such as -80 with precision 3!)

#   Jul  1, 2006: - Version 3.3 [HISTORY]

require "$ANTS/libconv.pl";

#----------------------------------------------------------------------

# Lat/Lon

#----------------------------------------------------------------------

sub DDMMXh2d(@)			# NODC SD2 header info

{

	my($DDMMX,$H) = &antsFunUsage(2,"51","DDMMX H",@_);

	$DDMMX =~ s/ /0/g;

	return &dmh2d(substr($DDMMX,0,2),

				  substr($DDMMX,2,2) . "." . substr($DDMMX,4,1),

				  $H);

}

sub DDDMMXh2d(@)		# NODC SD2 header info

{

	my($DDDMMX,$H) = &antsFunUsage(2,"61","DDDMMX H",@_);

	$DDDMMX =~ s/ /0/g;

	return &dmh2d(substr($DDDMMX,0,3),

				  substr($DDDMMX,3,2) . "." . substr($DDDMMX,5,1),

				  $H);

}

sub DDMMSSh2d(@)		# NODC detailed inventory info

{

	my($DDMMSS,$H) = &antsFunUsage(2,"61","DDMMSS H",@_);

	$DDMMSS =~ s/ /0/g;

	return &dmsh2d(substr($DDMMSS,0,2),

				   substr($DDMMSS,2,2),

				   substr($DDMMSS,4,2),

				   $H);

}

sub DDDMMSSh2d(@)		# NODC detailed inventory info

{

	my($DDDMMSS,$H) = &antsFunUsage(2,"71","DDDMMSS H",@_);

	$DDDMMSS =~ s/ /0/g;

	return &dmsh2d(substr($DDDMMSS,0,3),

				   substr($DDDMMSS,3,2),

				   substr($DDDMMSS,5,2),

				   $H);

}

#----------------------------------------------------------------------

# date/time

#----------------------------------------------------------------------

sub YYMMDD(@)		# 6 digit date

{

	my($YYMMDD) = &antsFunUsage(1,"6","YYMMDD",@_);

	return substr($YYMMDD,2,2) . "/" .

		   substr($YYMMDD,4,2) . "/19" . substr($YYMMDD,0,2);

}

sub HHt(@)			# 3 digits (hours to tenths)

{

	my($HHt) = &antsFunUsage(1,"3","HHt",@_);

	return sprintf("%02d:%02d",substr($HHt,0,2),substr($HHt,2,1)*6);

}

#----------------------------------------------------------------------

# depth

#----------------------------------------------------------------------

sub obs_depth(@)									# good depth only

{

	my($obs,$quality,$t_flag) =

		&antsFunUsage(3,"c..","obs quality t_flag",@_);

	return (isnan($quality) && ($t_flag ne 'T'))

		? $obs : nan;

}

sub wire_out(@)										# wire-out

{

	my($obs,$quality,$t_flag) =

		&antsFunUsage(3,"c..","obs quality t_flag",@_);

	return (($quality == 6) && ($t_flag ne 'T'))

		? $obs : nan;

}

sub t_depth(@)										# good thermometric depth

{

	my($obs,$quality,$t_flag) =

		&antsFunUsage(3,"c..","obs quality t_flag",@_);

	return (isnan($quality) && ($t_flag eq 'T'))

		? $obs : nan;

}

#----------------------------------------------------------------------

# temp, salin, O2, ...

#----------------------------------------------------------------------

sub obs(@)

{

	my($obs,$prec,$qual) =

		&antsFunUsage(3,".1.","obs prec qual",@_);

	return nan if isnan($obs);

	return nan if isnan($qual);						# spc->nan==OK

	my($fac) = 1;									# Reid and Mantyla weird fmt

	if ($obs =~ /^-/) {

		$fac = -1;

		$obs = $';

	}

	$obs = sprintf("%0${prec}d",$obs);				# pre-pad missing 0es

	substr($obs,-$prec,0) = ".";					# PERL is wonderful...

	return $fac * $obs;

}

1;