Blame view

highEkmanNLin/initialize.f90 803 Bytes
af19620a   Kévin Duquette   Add all files
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
  ! 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