anchor.diff:  read(11,rec=1) spfh
   read(12,rec=1) tsfc
   read(13,rec=1) psfc
  -read(14,rec=1) qsfc/1000.  ! for testing
  +read(14,rec=1) qsfc1
   close(10)
   close(11)
   close(12)
  
positive.code: character(len=400) :: tprof, qprof, tsfcfile, presfile, qsfcfile,tradfile,viewfile
  integer :: tindex
  character(len=255) :: ifile
  
  call getarg(1,tprof)
  call getarg(2,qprof)
  
negative.code: tc15(:,:) = -9999.
  tc55(:,:) = -9999.
  lst(:,:) = -9999.
  qsfc = qsfc/1000.
  !call fsin
  call wndo
  
  