initialize.f90 803 Bytes
  ! Initilisation
  zeta(:,:) = 0.0
  psi(:,:)  = 0.0
  v(:,:)    = 0.0
  rhs(:,:)  = 0.0
  jacoV(:,:)= 0.0
  jacoZ(:,:)= 0.0

  call readField(realV,'Vini.txt')
  call readField(realU,'Uini.txt')
  call writeField("realV.bin", realV)
  call writeField("realU.bin",realU)
  
  ! set initialisation and iguess to 0
  intlD     = 0
  iguessD   = 0
  intlN     = 0 
  iguessN   = 0
 
  write(*,*) 'Initialize Dirichlet'  
  ! Dirichlet
  call mud2(iprmD,fprm,workD,cof,bndc,rhs,zeta(1:nx,:),mgopt,ierror)
  call checkIerror(ierror,'Dirichlet')
 
  write(*,*) 'Initialize Neumman'    
  ! Neumman
  call mud2(iprmN,fprm,workN,cof,bndc,rhs,v(1:nx,:),mgopt,ierror)
  call checkIerror(ierror,'Neumman') 

  ! set initialisation and iguess to 0
  intlD     = 1
  intlN     = 1
  iguessD   = 1
  iguessN   = 1