initial_conditions.f90

subroutine initial_conditions()
use variables
implicit none

!---------------------------------------------------!
!         Initial conditions calculation            !
!---------------------------------------------------!

!$OMP PARALLEL DO PRIVATE(j,i)  
do k=-1,nz+2; do j=-1,ny+2; do i=-1,nx+2

   u(i,j,k) = 0.0
   v(i,j,k) = 0.0
   w(i,j,k) = 0.0
   u1(i,j,k) = 0.0
   v1(i,j,k) = 0.0
   w1(i,j,k) = 0.0
   u2(i,j,k) = 0.0
   v2(i,j,k) = 0.0
   w2(i,j,k) = 0.0
   p(i,j,k) = 0.0
   u_star(i,j,k) = 0.0
   v_star(i,j,k) = 0.0
   w_star(i,j,k) = 0.0
   ETA(i,j,k) = 0.0
   F_tavex(i,j,k) = 0.0
   F_tavey(i,j,k) = 0.0

end do; end do; end do
!$OMP END PARALLEL DO

end subroutine initial_conditions