anchor.diff:  tc15(:,:) = -9999.
   tc55(:,:) = -9999.
   lst(:,:) = -9999.
  +qsfc = qsfc/1000.
   !call fsin
   call wndo
   
  
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: read(11,rec=1) spfh
  read(12,rec=1) tsfc
  read(13,rec=1) psfc
  read(14,rec=1) qsfc1
  close(10)
  close(11)
  close(12)
  