anchor.diff: +program test
  +call test_missing_values
  +stop
  +end
   subroutine test_missing_values
   !
   !  rm -f missing.fst ; MISSING_VALUE_FLAGS="99.0 127 255" ./a.out
  
positive.code: sa(2)=sm   ; sa(ASIZE-1)=sm ; sa(4)=sm   ; sa(ASIZE-3)=sm
  ba(2)=bm   ; ba(ASIZE-1)=bm ; ba(1)=bm
  uia(3)=uim ; uia(ASIZE-2)=uim
  usa(3)=usm ; usa(ASIZE-2)=usm ; usa(ASIZE)=usm
  uba(3)=ubm ; uba(ASIZE-2)=ubm
  ia(ASIZE/2)=125      ! set to 127 to force error message
  sa(ASIZE/2)=124      ! set to 127 to force error message
  ba(ASIZE/2)=123      ! set to 127 to force error message
  
negative.code: call fstecr(da,work,-64,11,0,0,0,ASIZE,1,1,18,0,0,'XX','YYYY','ETIKET','X',0,0,0,0,69,.false.) ! IEEE64 with missing
  
  call fst_data_length(2)
  call fstecr(usa,work,-16,11,0,0,0,ASIZE,1,1,9,0,0,'XX','YYYY','ETIKET','X',0,0,0,0,2,.false.)  ! unsigned short
  call fst_data_length(2)
  call fstecr(usa,work,-16,11,0,0,0,ASIZE,1,1,10,0,0,'XX','YYYY','ETIKET','X',0,0,0,0,66,.false.)! unsigned short with missing
  
  call fst_data_length(2)
  call fstecr(sa,work,-8,11,0,0,0,ASIZE,1,1,11,0,0,'XX','YYYY','ETIKET','X',0,0,0,0,4,.false.)  ! signed short
  