#$=head1 NAME
#$
#$mask1 - simple masking operator
#$
#$=head1 SYNOPSIS
#$
#$call mask1_init(m)
#$
#$C<ierr=mask1_lop(adj,add,model,data)>
#$
#$=head1 PARAMETERS
#$
#$=over 4
#$
#$=item m  - C<logical(:)>  
#$
#$      Masking operator
#$
#$=item adj,add,model,data -
#$
#$       Standard operator interfaces
#$
#$=back
#$
#$=head1 DESCRIPTION
#$
#$Simple masking opertor.  data=mod when mask=.true.
#$
#$=head1 LIBRARY
#$
#$B<geef90>
#$
#$=cut
#$
#$
module mask1 { # masking operator
  logical, dimension( :), pointer :: m
#% _init( m)
#% _lop( x, y)
  if( adj) 
	   where( m) x += y
  else  #
           where( m) y += x
}