anchor.diff:  			ny=size(domain%u,3)
   			if (options%external_winds) then
   				call ext_winds_init(domain,bc,options)
  +				call smooth_wind(domain%u,1,3)
  +				call smooth_wind(domain%v,1,3)
   			elseif (options%remove_lowres_linear) then
   				call remove_linear_winds(domain,bc,options,file_list(curfile),curstep)
   			elseif (options%mean_winds) then
  
positive.code: 		use_boundary=.True.
  		if (options%external_winds) then
  			call update_ext_winds(bc,options)
  		elseif (options%remove_lowres_linear) then
  			call remove_linear_winds(bc%next_domain,bc,options,file_list(curfile),curstep)
  		elseif (options%mean_winds) then
  
negative.code: 		else
  			call read_var(bc%next_domain%u,    file_list(curfile),"U",      bc%geolut,curstep,use_interior)
  			call read_var(bc%next_domain%v,    file_list(curfile),"V",      bc%geolut,curstep,use_interior)
  			call smooth_wind(bc%next_domain%u,smoothing_window,3)
  			call smooth_wind(bc%next_domain%v,smoothing_window,3)
  		endif
  		call read_var(bc%next_domain%p,    file_list(curfile),"P",      bc%geolut,curstep,use_interior)
  		call read_var(bc%next_domain%th,   file_list(curfile),"T",      bc%geolut,curstep,use_boundary)
  