anchor.diff:  c
   c	d_vs_vl=vs(ig_refmax)-vl(ig_refmax)	! Olivier
   c	kvsl=ifix((d_vs_vl/vnyq_el)*0.5)+5	! Olivier
  -c        print *,'    d_VS_VL :',d_vs_vl,' -> KVSL :',kvsl	
  +c        print *,'    d_VS_VL :',d_vs_vl,' -> KVSL :',kvsl
   c	if(kvsl.ge.1.and.kvsl.le.9)then		! Olivier
   c	  vs_depl=vs(ig_refmax)+xms(kvsl)*vnyq_el	! Olivier
   c	  vl_depl=vl(ig_refmax)+xml(kvsl)*vnyq_el	! Olivier
  
positive.code:      &               ,sqrt(rms_var_vinsitu(12)/swadvinsitu_tot)
              else
                print *,' !!!! DVDOP_insitu -> sWGHTs:'
       &               ,swdvinsitu_tot,' !!!!' 
              endif
  c
  c******************************************************************
  
negative.code:             if(sp(i,j).gt.spmin)then
                sz(i,j)=sz(i,j)/sp(i,j)
                nin=nin+1
  	    else
  	      sz(i,j)=-999.
              endif
           enddo
  