anchor.diff:  		
   		if (options%restart) then
   			call load_restart_file(domain,options%restart_file)
  +			if (options%external_winds) then
  +				call ext_winds_init(domain,bc,options)
  +			endif
   		else
   			boundary_value=.False.
   			nx=size(domain%u,1)
  
positive.code: 		logical :: boundary_value
  		integer::nx,ny,nz,i
  		real::domainsize
  		
  		curfile=1
  		if (options%restart) then
  
negative.code: 		type(options_type),intent(in)::options
  		integer,dimension(io_maxDims)::dims	!note, io_maxDims is included from io_routines.
  		logical :: use_boundary,use_interior
  		! MODULE variables : ext_winds_ curstep, curfile, nfiles, steps_in_file, file_list
  		ext_winds_curstep=ext_winds_curstep+1
  		if (ext_winds_curstep>ext_winds_steps_in_file) then
  			ext_winds_curfile=ext_winds_curfile+1
  