/FluidMechanics/NavierStokes/RoutineCheck/ParallelStatic/src/ParallelStaticExample.f90

http://github.com/xyan075/examples · Fortran Modern · 793 lines · 460 code · 105 blank · 228 comment · 0 complexity · 8ef41b9505041bfde45ee3a3f2f7e6e1 MD5 · raw file

  1. !> \file
  2. !> \author Sebastian Krittian
  3. !> \brief This is an example program to solve a static Navier-Stokes equation using OpenCMISS calls.
  4. !>
  5. !> \section LICENSE
  6. !>
  7. !> Version: MPL 1.1/GPL 2.0/LGPL 2.1
  8. !>
  9. !> The contents of this file are subject to the Mozilla Public License
  10. !> Version 1.1 (the "License"); you may not use this file except in
  11. !> compliance with the License. You may obtain a copy of the License at
  12. !> http://www.mozilla.org/MPL/
  13. !>
  14. !> Software distributed under the License is distributed on an "AS IS"
  15. !> basis, WITHOUT WARRANTY OF ANY KIND, either express or implied. See the
  16. !> License for the specific language governing rights and limitations
  17. !> under the License.
  18. !>
  19. !> The Original Code is OpenCMISS
  20. !>
  21. !> The Initial Developer of the Original Code is University of Auckland,
  22. !> Auckland, New Zealand and University of Oxford, Oxford, United
  23. !> Kingdom. Portions created by the University of Auckland and University
  24. !> of Oxford are Copyright (C) 2007 by the University of Auckland and
  25. !> the University of Oxford. All Rights Reserved.
  26. !>
  27. !> Contributor(s):
  28. !>
  29. !> Alternatively, the contents of this file may be used under the terms of
  30. !> either the GNU General Public License Version 2 or later (the "GPL"), or
  31. !> the GNU Lesser General Public License Version 2.1 or later (the "LGPL"),
  32. !> in which case the provisions of the GPL or the LGPL are applicable instead
  33. !> of those above. If you wish to allow use of your version of this file only
  34. !> under the terms of either the GPL or the LGPL, and not to allow others to
  35. !> use your version of this file under the terms of the MPL, indicate your
  36. !> decision by deleting the provisions above and replace them with the notice
  37. !> and other provisions required by the GPL or the LGPL. If you do not delete
  38. !> the provisions above, a recipient may use your version of this file under
  39. !> the terms of any one of the MPL, the GPL or the LGPL.
  40. !>
  41. !> \example FluidMechanics/NavierStokes/RoutineCheck/Static/src/StaticExample.f90
  42. !! Example program to solve a static Navier-Stokes equation using OpenCMISS calls.
  43. !! \htmlinclude FluidMechanics/NavierStokes/RoutineCheck/Static/history.html
  44. !!
  45. !<
  46. !> Main program
  47. PROGRAM NAVIERSTOKESSTATICEXAMPLE
  48. !
  49. !================================================================================================================================
  50. !
  51. !PROGRAM LIBRARIES
  52. USE OPENCMISS
  53. USE FLUID_MECHANICS_IO_ROUTINES
  54. USE MPI
  55. #ifdef WIN32
  56. USE IFQWINCMISS
  57. #endif
  58. !
  59. !================================================================================================================================
  60. !
  61. !PROGRAM VARIABLES AND TYPES
  62. IMPLICIT NONE
  63. INTEGER(CMISSIntg), PARAMETER :: EquationsSetFieldUserNumber=1337
  64. TYPE(CMISSFieldType) :: EquationsSetField
  65. !Test program parameters
  66. INTEGER(CMISSIntg), PARAMETER :: CoordinateSystemUserNumber=1
  67. INTEGER(CMISSIntg), PARAMETER :: RegionUserNumber=2
  68. INTEGER(CMISSIntg), PARAMETER :: MeshUserNumber=3
  69. INTEGER(CMISSIntg), PARAMETER :: DecompositionUserNumber=4
  70. INTEGER(CMISSIntg), PARAMETER :: GeometricFieldUserNumber=5
  71. INTEGER(CMISSIntg), PARAMETER :: DependentFieldUserNumberNavierStokes=6
  72. INTEGER(CMISSIntg), PARAMETER :: MaterialsFieldUserNumberNavierStokes=7
  73. INTEGER(CMISSIntg), PARAMETER :: IndependentFieldUserNumberNavierStokes=8
  74. INTEGER(CMISSIntg), PARAMETER :: EquationsSetUserNumberNavierStokes=9
  75. INTEGER(CMISSIntg), PARAMETER :: ProblemUserNumber=10
  76. INTEGER(CMISSIntg), PARAMETER :: DomainUserNumber=2
  77. INTEGER(CMISSIntg), PARAMETER :: SolverNavierStokesUserNumber=1
  78. INTEGER(CMISSIntg), PARAMETER :: MaterialsFieldUserNumberNavierStokesMu=1
  79. INTEGER(CMISSIntg), PARAMETER :: MaterialsFieldUserNumberNavierStokesRho=2
  80. !Program types
  81. TYPE(EXPORT_CONTAINER):: CM
  82. !Program variables
  83. INTEGER(CMISSIntg) :: NUMBER_OF_DIMENSIONS
  84. INTEGER(CMISSIntg) :: BASIS_TYPE
  85. INTEGER(CMISSIntg) :: BASIS_NUMBER_SPACE
  86. INTEGER(CMISSIntg) :: BASIS_NUMBER_VELOCITY
  87. INTEGER(CMISSIntg) :: BASIS_NUMBER_PRESSURE
  88. INTEGER(CMISSIntg) :: BASIS_XI_GAUSS_SPACE
  89. INTEGER(CMISSIntg) :: BASIS_XI_GAUSS_VELOCITY
  90. INTEGER(CMISSIntg) :: BASIS_XI_GAUSS_PRESSURE
  91. INTEGER(CMISSIntg) :: BASIS_XI_INTERPOLATION_SPACE
  92. INTEGER(CMISSIntg) :: BASIS_XI_INTERPOLATION_VELOCITY
  93. INTEGER(CMISSIntg) :: BASIS_XI_INTERPOLATION_PRESSURE
  94. INTEGER(CMISSIntg) :: MESH_NUMBER_OF_COMPONENTS
  95. INTEGER(CMISSIntg) :: MESH_COMPONENT_NUMBER_SPACE
  96. INTEGER(CMISSIntg) :: MESH_COMPONENT_NUMBER_VELOCITY
  97. INTEGER(CMISSIntg) :: MESH_COMPONENT_NUMBER_PRESSURE
  98. INTEGER(CMISSIntg) :: NUMBER_OF_NODES_SPACE
  99. INTEGER(CMISSIntg) :: NUMBER_OF_NODES_VELOCITY
  100. INTEGER(CMISSIntg) :: NUMBER_OF_NODES_PRESSURE
  101. INTEGER(CMISSIntg) :: NUMBER_OF_ELEMENT_NODES_SPACE
  102. INTEGER(CMISSIntg) :: NUMBER_OF_ELEMENT_NODES_VELOCITY
  103. INTEGER(CMISSIntg) :: NUMBER_OF_ELEMENT_NODES_PRESSURE
  104. INTEGER(CMISSIntg) :: TOTAL_NUMBER_OF_NODES
  105. INTEGER(CMISSIntg) :: TOTAL_NUMBER_OF_ELEMENTS
  106. INTEGER(CMISSIntg) :: MAXIMUM_ITERATIONS
  107. INTEGER(CMISSIntg) :: RESTART_VALUE
  108. ! INTEGER(CMISSIntg) :: MPI_IERROR
  109. INTEGER(CMISSIntg) :: NUMBER_OF_FIXED_WALL_NODES_NAVIER_STOKES
  110. INTEGER(CMISSIntg) :: NUMBER_OF_INLET_WALL_NODES_NAVIER_STOKES
  111. INTEGER(CMISSIntg) :: EQUATIONS_NAVIER_STOKES_OUTPUT
  112. INTEGER(CMISSIntg) :: COMPONENT_NUMBER
  113. INTEGER(CMISSIntg) :: NODE_NUMBER
  114. INTEGER(CMISSIntg) :: ELEMENT_NUMBER
  115. INTEGER(CMISSIntg) :: NODE_COUNTER
  116. INTEGER(CMISSIntg) :: CONDITION
  117. INTEGER(CMISSIntg) :: LINEAR_SOLVER_NAVIER_STOKES_OUTPUT_TYPE
  118. INTEGER(CMISSIntg) :: NONLINEAR_SOLVER_NAVIER_STOKES_OUTPUT_TYPE
  119. INTEGER, ALLOCATABLE, DIMENSION(:):: FIXED_WALL_NODES_NAVIER_STOKES
  120. INTEGER, ALLOCATABLE, DIMENSION(:):: INLET_WALL_NODES_NAVIER_STOKES
  121. REAL(CMISSDP) :: INITIAL_FIELD_NAVIER_STOKES(3)
  122. REAL(CMISSDP) :: BOUNDARY_CONDITIONS_NAVIER_STOKES(3)
  123. REAL(CMISSDP) :: DIVERGENCE_TOLERANCE
  124. REAL(CMISSDP) :: RELATIVE_TOLERANCE
  125. REAL(CMISSDP) :: ABSOLUTE_TOLERANCE
  126. REAL(CMISSDP) :: LINESEARCH_ALPHA
  127. REAL(CMISSDP) :: VALUE
  128. REAL(CMISSDP) :: MU_PARAM_NAVIER_STOKES
  129. REAL(CMISSDP) :: RHO_PARAM_NAVIER_STOKES
  130. LOGICAL :: EXPORT_FIELD_IO
  131. LOGICAL :: LINEAR_SOLVER_NAVIER_STOKES_DIRECT_FLAG
  132. LOGICAL :: FIXED_WALL_NODES_NAVIER_STOKES_FLAG
  133. LOGICAL :: INLET_WALL_NODES_NAVIER_STOKES_FLAG
  134. !CMISS variables
  135. !Regions
  136. TYPE(CMISSRegionType) :: Region
  137. TYPE(CMISSRegionType) :: WorldRegion
  138. !Coordinate systems
  139. TYPE(CMISSCoordinateSystemType) :: CoordinateSystem
  140. TYPE(CMISSCoordinateSystemType) :: WorldCoordinateSystem
  141. !Basis
  142. TYPE(CMISSBasisType) :: BasisSpace
  143. TYPE(CMISSBasisType) :: BasisVelocity
  144. TYPE(CMISSBasisType) :: BasisPressure
  145. !Nodes
  146. TYPE(CMISSNodesType) :: Nodes
  147. !Elements
  148. TYPE(CMISSMeshElementsType) :: MeshElementsSpace
  149. TYPE(CMISSMeshElementsType) :: MeshElementsVelocity
  150. TYPE(CMISSMeshElementsType) :: MeshElementsPressure
  151. !Meshes
  152. TYPE(CMISSMeshType) :: Mesh
  153. !Decompositions
  154. TYPE(CMISSDecompositionType) :: Decomposition
  155. !Fields
  156. TYPE(CMISSFieldsType) :: Fields
  157. !Field types
  158. TYPE(CMISSFieldType) :: GeometricField
  159. TYPE(CMISSFieldType) :: DependentFieldNavierStokes
  160. TYPE(CMISSFieldType) :: MaterialsFieldNavierStokes
  161. !Boundary conditions
  162. TYPE(CMISSBoundaryConditionsType) :: BoundaryConditionsNavierStokes
  163. !Equations sets
  164. TYPE(CMISSEquationsSetType) :: EquationsSetNavierStokes
  165. !Equations
  166. TYPE(CMISSEquationsType) :: EquationsNavierStokes
  167. !Problems
  168. TYPE(CMISSProblemType) :: Problem
  169. !Control loops
  170. TYPE(CMISSControlLoopType) :: ControlLoop
  171. !Solvers
  172. TYPE(CMISSSolverType) :: NonlinearSolverNavierStokes
  173. TYPE(CMISSSolverType) :: LinearSolverNavierStokes
  174. !Solver equations
  175. TYPE(CMISSSolverEquationsType) :: SolverEquationsNavierStokes
  176. #ifdef WIN32
  177. !Quickwin type
  178. LOGICAL :: QUICKWIN_STATUS=.FALSE.
  179. TYPE(WINDOWCONFIG) :: QUICKWIN_WINDOW_CONFIG
  180. #endif
  181. !Generic CMISS variables
  182. INTEGER(CMISSIntg) :: NumberOfComputationalNodes,ComputationalNodeNumber,BoundaryNodeDomain
  183. INTEGER(CMISSIntg) :: EquationsSetIndex
  184. INTEGER(CMISSIntg) :: Err
  185. #ifdef WIN32
  186. !Initialise QuickWin
  187. QUICKWIN_WINDOW_CONFIG%TITLE="General Output" !Window title
  188. QUICKWIN_WINDOW_CONFIG%NUMTEXTROWS=-1 !Max possible number of rows
  189. QUICKWIN_WINDOW_CONFIG%MODE=QWIN$SCROLLDOWN
  190. !Set the window parameters
  191. QUICKWIN_STATUS=SETWINDOWCONFIG(QUICKWIN_WINDOW_CONFIG)
  192. !If attempt fails set with system estimated values
  193. IF(.NOT.QUICKWIN_STATUS) QUICKWIN_STATUS=SETWINDOWCONFIG(QUICKWIN_WINDOW_CONFIG)
  194. #endif
  195. !
  196. !================================================================================================================================
  197. !
  198. !INITIALISE OPENCMISS
  199. CALL CMISSInitialise(WorldCoordinateSystem,WorldRegion,Err)
  200. !
  201. !================================================================================================================================
  202. !
  203. !CHECK COMPUTATIONAL NODE
  204. !Get the computational nodes information
  205. CALL CMISSComputationalNumberOfNodesGet(NumberOfComputationalNodes,Err)
  206. CALL CMISSComputationalNodeNumberGet(ComputationalNodeNumber,Err)
  207. !
  208. !================================================================================================================================
  209. !
  210. !PROBLEM CONTROL PANEL
  211. !Import cmHeart mesh information
  212. CALL FLUID_MECHANICS_IO_READ_CMHEART(CM,Err)
  213. BASIS_NUMBER_SPACE=CM%ID_M
  214. BASIS_NUMBER_VELOCITY=CM%ID_V
  215. BASIS_NUMBER_PRESSURE=CM%ID_P
  216. NUMBER_OF_DIMENSIONS=CM%D
  217. BASIS_TYPE=CM%IT_T
  218. BASIS_XI_INTERPOLATION_SPACE=CM%IT_M
  219. BASIS_XI_INTERPOLATION_VELOCITY=CM%IT_V
  220. BASIS_XI_INTERPOLATION_PRESSURE=CM%IT_P
  221. NUMBER_OF_NODES_SPACE=CM%N_M
  222. NUMBER_OF_NODES_VELOCITY=CM%N_V
  223. NUMBER_OF_NODES_PRESSURE=CM%N_P
  224. TOTAL_NUMBER_OF_NODES=CM%N_T
  225. TOTAL_NUMBER_OF_ELEMENTS=CM%E_T
  226. NUMBER_OF_ELEMENT_NODES_SPACE=CM%EN_M
  227. NUMBER_OF_ELEMENT_NODES_VELOCITY=CM%EN_V
  228. NUMBER_OF_ELEMENT_NODES_PRESSURE=CM%EN_P
  229. !Set initial values
  230. INITIAL_FIELD_NAVIER_STOKES(1)=0.0_CMISSDP
  231. INITIAL_FIELD_NAVIER_STOKES(2)=0.0_CMISSDP
  232. INITIAL_FIELD_NAVIER_STOKES(3)=0.0_CMISSDP
  233. !Set boundary conditions
  234. FIXED_WALL_NODES_NAVIER_STOKES_FLAG=.TRUE.
  235. INLET_WALL_NODES_NAVIER_STOKES_FLAG=.TRUE.
  236. IF(FIXED_WALL_NODES_NAVIER_STOKES_FLAG) THEN
  237. NUMBER_OF_FIXED_WALL_NODES_NAVIER_STOKES=80
  238. ALLOCATE(FIXED_WALL_NODES_NAVIER_STOKES(NUMBER_OF_FIXED_WALL_NODES_NAVIER_STOKES))
  239. FIXED_WALL_NODES_NAVIER_STOKES=(/1,2,3,4,5,7,9,10,11,12,13,14,17,20,24,28,29,30,31,32,33,34,35,37,39, &
  240. & 41,44,46,47,48,50,51,52,53,54,57,60,64,65,66,67,68,70,72,74,76,77,78,79,80,83,86, &
  241. & 89,90,91,92,93,94,95,97,99,101,102,103,104,105,106,107,108,111,114,115,116,117,118, &
  242. & 120,122,123,124,125/)
  243. ENDIF
  244. IF(INLET_WALL_NODES_NAVIER_STOKES_FLAG) THEN
  245. NUMBER_OF_INLET_WALL_NODES_NAVIER_STOKES=9
  246. ALLOCATE(INLET_WALL_NODES_NAVIER_STOKES(NUMBER_OF_INLET_WALL_NODES_NAVIER_STOKES))
  247. INLET_WALL_NODES_NAVIER_STOKES=(/6,15,16,23,36,42,81,82,96/)
  248. !Set initial boundary conditions
  249. BOUNDARY_CONDITIONS_NAVIER_STOKES(1)=0.0_CMISSDP
  250. BOUNDARY_CONDITIONS_NAVIER_STOKES(2)=1.0_CMISSDP
  251. BOUNDARY_CONDITIONS_NAVIER_STOKES(3)=0.0_CMISSDP
  252. ENDIF
  253. !Set material parameters
  254. MU_PARAM_NAVIER_STOKES=1.0_CMISSDP
  255. RHO_PARAM_NAVIER_STOKES=0.0_CMISSDP
  256. !Set interpolation parameters
  257. BASIS_XI_GAUSS_SPACE=3
  258. BASIS_XI_GAUSS_VELOCITY=3
  259. BASIS_XI_GAUSS_PRESSURE=3
  260. !Set output parameter
  261. !(NoOutput/ProgressOutput/TimingOutput/SolverOutput/SolverMatrixOutput)
  262. LINEAR_SOLVER_NAVIER_STOKES_OUTPUT_TYPE=CMISS_SOLVER_PROGRESS_OUTPUT
  263. NONLINEAR_SOLVER_NAVIER_STOKES_OUTPUT_TYPE=CMISS_SOLVER_MATRIX_OUTPUT
  264. !(NoOutput/TimingOutput/MatrixOutput/ElementOutput)
  265. EQUATIONS_NAVIER_STOKES_OUTPUT=CMISS_EQUATIONS_NO_OUTPUT
  266. !Set solver parameters
  267. LINEAR_SOLVER_NAVIER_STOKES_DIRECT_FLAG=.TRUE.
  268. RELATIVE_TOLERANCE=1.0E-5_CMISSDP !default: 1.0E-05_CMISSDP
  269. ABSOLUTE_TOLERANCE=1.0E-10_CMISSDP !default: 1.0E-10_CMISSDP
  270. DIVERGENCE_TOLERANCE=1.0E20 !default: 1.0E5
  271. MAXIMUM_ITERATIONS=100000 !default: 100000
  272. RESTART_VALUE=3000 !default: 30
  273. LINESEARCH_ALPHA=1.0
  274. !
  275. !================================================================================================================================
  276. !
  277. !COORDINATE SYSTEM
  278. !Start the creation of a new RC coordinate system
  279. CALL CMISSCoordinateSystem_Initialise(CoordinateSystem,Err)
  280. CALL CMISSCoordinateSystem_CreateStart(CoordinateSystemUserNumber,CoordinateSystem,Err)
  281. !Set the coordinate system dimension
  282. CALL CMISSCoordinateSystem_DimensionSet(CoordinateSystem,NUMBER_OF_DIMENSIONS,Err)
  283. !Finish the creation of the coordinate system
  284. CALL CMISSCoordinateSystem_CreateFinish(CoordinateSystem,Err)
  285. !
  286. !================================================================================================================================
  287. !
  288. !REGION
  289. !Start the creation of a new region
  290. CALL CMISSRegion_Initialise(Region,Err)
  291. CALL CMISSRegion_CreateStart(RegionUserNumber,WorldRegion,Region,Err)
  292. !Set the regions coordinate system as defined above
  293. CALL CMISSRegion_CoordinateSystemSet(Region,CoordinateSystem,Err)
  294. !Finish the creation of the region
  295. CALL CMISSRegion_CreateFinish(Region,Err)
  296. !
  297. !================================================================================================================================
  298. !
  299. !BASES
  300. !Start the creation of new bases
  301. MESH_NUMBER_OF_COMPONENTS=1
  302. CALL CMISSBasis_Initialise(BasisSpace,Err)
  303. CALL CMISSBasis_CreateStart(BASIS_NUMBER_SPACE,BasisSpace,Err)
  304. !Set the basis type (Lagrange/Simplex)
  305. CALL CMISSBasis_TypeSet(BasisSpace,BASIS_TYPE,Err)
  306. !Set the basis xi number
  307. CALL CMISSBasis_NumberOfXiSet(BasisSpace,NUMBER_OF_DIMENSIONS,Err)
  308. !Set the basis xi interpolation and number of Gauss points
  309. IF(NUMBER_OF_DIMENSIONS==2) THEN
  310. CALL CMISSBasis_InterpolationXiSet(BasisSpace,(/BASIS_XI_INTERPOLATION_SPACE,BASIS_XI_INTERPOLATION_SPACE/),Err)
  311. CALL CMISSBasis_QuadratureNumberOfGaussXiSet(BasisSpace,(/BASIS_XI_GAUSS_SPACE,BASIS_XI_GAUSS_SPACE/),Err)
  312. ELSE IF(NUMBER_OF_DIMENSIONS==3) THEN
  313. CALL CMISSBasis_InterpolationXiSet(BasisSpace,(/BASIS_XI_INTERPOLATION_SPACE,BASIS_XI_INTERPOLATION_SPACE, &
  314. & BASIS_XI_INTERPOLATION_SPACE/),Err)
  315. CALL CMISSBasis_QuadratureNumberOfGaussXiSet(BasisSpace,(/BASIS_XI_GAUSS_SPACE,BASIS_XI_GAUSS_SPACE,BASIS_XI_GAUSS_SPACE/),Err)
  316. ENDIF
  317. !Finish the creation of the basis
  318. CALL CMISSBasis_CreateFinish(BasisSpace,Err)
  319. !Start the creation of another basis
  320. IF(BASIS_XI_INTERPOLATION_VELOCITY==BASIS_XI_INTERPOLATION_SPACE) THEN
  321. BasisVelocity=BasisSpace
  322. ELSE
  323. MESH_NUMBER_OF_COMPONENTS=MESH_NUMBER_OF_COMPONENTS+1
  324. !Initialise a new velocity basis
  325. CALL CMISSBasis_Initialise(BasisVelocity,Err)
  326. !Start the creation of a basis
  327. CALL CMISSBasis_CreateStart(BASIS_NUMBER_VELOCITY,BasisVelocity,Err)
  328. !Set the basis type (Lagrange/Simplex)
  329. CALL CMISSBasis_TypeSet(BasisVelocity,BASIS_TYPE,Err)
  330. !Set the basis xi number
  331. CALL CMISSBasis_NumberOfXiSet(BasisVelocity,NUMBER_OF_DIMENSIONS,Err)
  332. !Set the basis xi interpolation and number of Gauss points
  333. IF(NUMBER_OF_DIMENSIONS==2) THEN
  334. CALL CMISSBasis_InterpolationXiSet(BasisVelocity,(/BASIS_XI_INTERPOLATION_VELOCITY,BASIS_XI_INTERPOLATION_VELOCITY/),Err)
  335. CALL CMISSBasis_QuadratureNumberOfGaussXiSet(BasisVelocity,(/BASIS_XI_GAUSS_VELOCITY,BASIS_XI_GAUSS_VELOCITY/),Err)
  336. ELSE IF(NUMBER_OF_DIMENSIONS==3) THEN
  337. CALL CMISSBasis_InterpolationXiSet(BasisVelocity,(/BASIS_XI_INTERPOLATION_VELOCITY,BASIS_XI_INTERPOLATION_VELOCITY, &
  338. & BASIS_XI_INTERPOLATION_VELOCITY/),Err)
  339. CALL CMISSBasis_QuadratureNumberOfGaussXiSet(BasisVelocity,(/BASIS_XI_GAUSS_VELOCITY,BASIS_XI_GAUSS_VELOCITY, &
  340. & BASIS_XI_GAUSS_VELOCITY/),Err)
  341. ENDIF
  342. !Finish the creation of the basis
  343. CALL CMISSBasis_CreateFinish(BasisVelocity,Err)
  344. ENDIF
  345. !Start the creation of another basis
  346. IF(BASIS_XI_INTERPOLATION_PRESSURE==BASIS_XI_INTERPOLATION_SPACE) THEN
  347. BasisPressure=BasisSpace
  348. ELSE IF(BASIS_XI_INTERPOLATION_PRESSURE==BASIS_XI_INTERPOLATION_VELOCITY) THEN
  349. BasisPressure=BasisVelocity
  350. ELSE
  351. MESH_NUMBER_OF_COMPONENTS=MESH_NUMBER_OF_COMPONENTS+1
  352. !Initialise a new pressure basis
  353. CALL CMISSBasis_Initialise(BasisPressure,Err)
  354. !Start the creation of a basis
  355. CALL CMISSBasis_CreateStart(BASIS_NUMBER_PRESSURE,BasisPressure,Err)
  356. !Set the basis type (Lagrange/Simplex)
  357. CALL CMISSBasis_TypeSet(BasisPressure,BASIS_TYPE,Err)
  358. !Set the basis xi number
  359. CALL CMISSBasis_NumberOfXiSet(BasisPressure,NUMBER_OF_DIMENSIONS,Err)
  360. !Set the basis xi interpolation and number of Gauss points
  361. IF(NUMBER_OF_DIMENSIONS==2) THEN
  362. CALL CMISSBasis_InterpolationXiSet(BasisPressure,(/BASIS_XI_INTERPOLATION_PRESSURE,BASIS_XI_INTERPOLATION_PRESSURE/),Err)
  363. CALL CMISSBasis_QuadratureNumberOfGaussXiSet(BasisPressure,(/BASIS_XI_GAUSS_PRESSURE,BASIS_XI_GAUSS_PRESSURE/),Err)
  364. ELSE IF(NUMBER_OF_DIMENSIONS==3) THEN
  365. CALL CMISSBasis_InterpolationXiSet(BasisPressure,(/BASIS_XI_INTERPOLATION_PRESSURE,BASIS_XI_INTERPOLATION_PRESSURE, &
  366. & BASIS_XI_INTERPOLATION_PRESSURE/),Err)
  367. CALL CMISSBasis_QuadratureNumberOfGaussXiSet(BasisPressure,(/BASIS_XI_GAUSS_PRESSURE,BASIS_XI_GAUSS_PRESSURE, &
  368. & BASIS_XI_GAUSS_PRESSURE/),Err)
  369. ENDIF
  370. !Finish the creation of the basis
  371. CALL CMISSBasis_CreateFinish(BasisPressure,Err)
  372. ENDIF
  373. !
  374. !================================================================================================================================
  375. !
  376. !MESH
  377. !Start the creation of mesh nodes
  378. CALL CMISSNodes_Initialise(Nodes,Err)
  379. CALL CMISSMesh_Initialise(Mesh,Err)
  380. CALL CMISSNodes_CreateStart(Region,TOTAL_NUMBER_OF_NODES,Nodes,Err)
  381. CALL CMISSNodes_CreateFinish(Nodes,Err)
  382. !Start the creation of the mesh
  383. CALL CMISSMesh_CreateStart(MeshUserNumber,Region,NUMBER_OF_DIMENSIONS,Mesh,Err)
  384. !Set number of mesh elements
  385. CALL CMISSMesh_NumberOfElementsSet(Mesh,TOTAL_NUMBER_OF_ELEMENTS,Err)
  386. !Set number of mesh components
  387. CALL CMISSMesh_NumberOfComponentsSet(Mesh,MESH_NUMBER_OF_COMPONENTS,Err)
  388. !Specify spatial mesh component
  389. CALL CMISSMeshElements_Initialise(MeshElementsSpace,Err)
  390. CALL CMISSMeshElements_Initialise(MeshElementsVelocity,Err)
  391. CALL CMISSMeshElements_Initialise(MeshElementsPressure,Err)
  392. MESH_COMPONENT_NUMBER_SPACE=1
  393. MESH_COMPONENT_NUMBER_VELOCITY=1
  394. MESH_COMPONENT_NUMBER_PRESSURE=1
  395. CALL CMISSMeshElements_CreateStart(Mesh,MESH_COMPONENT_NUMBER_SPACE,BasisSpace,MeshElementsSpace,Err)
  396. DO ELEMENT_NUMBER=1,TOTAL_NUMBER_OF_ELEMENTS
  397. CALL CMISSMeshElements_NodesSet(MeshElementsSpace,ELEMENT_NUMBER,CM%M(ELEMENT_NUMBER,1:NUMBER_OF_ELEMENT_NODES_SPACE),Err)
  398. ENDDO
  399. CALL CMISSMeshElements_CreateFinish(MeshElementsSpace,Err)
  400. !Specify velocity mesh component
  401. IF(BASIS_XI_INTERPOLATION_VELOCITY==BASIS_XI_INTERPOLATION_SPACE) THEN
  402. MeshElementsVelocity=MeshElementsSpace
  403. ELSE
  404. MESH_COMPONENT_NUMBER_VELOCITY=MESH_COMPONENT_NUMBER_SPACE+1
  405. CALL CMISSMeshElements_CreateStart(Mesh,MESH_COMPONENT_NUMBER_VELOCITY,BasisVelocity,MeshElementsVelocity,Err)
  406. DO ELEMENT_NUMBER=1,TOTAL_NUMBER_OF_ELEMENTS
  407. CALL CMISSMeshElements_NodesSet(MeshElementsVelocity,ELEMENT_NUMBER,CM%V(ELEMENT_NUMBER, &
  408. & 1:NUMBER_OF_ELEMENT_NODES_VELOCITY),Err)
  409. ENDDO
  410. CALL CMISSMeshElements_CreateFinish(MeshElementsVelocity,Err)
  411. ENDIF
  412. !Specify pressure mesh component
  413. IF(BASIS_XI_INTERPOLATION_PRESSURE==BASIS_XI_INTERPOLATION_SPACE) THEN
  414. MeshElementsPressure=MeshElementsSpace
  415. MESH_COMPONENT_NUMBER_PRESSURE=MESH_COMPONENT_NUMBER_SPACE
  416. ELSE IF(BASIS_XI_INTERPOLATION_PRESSURE==BASIS_XI_INTERPOLATION_VELOCITY) THEN
  417. MeshElementsPressure=MeshElementsVelocity
  418. MESH_COMPONENT_NUMBER_PRESSURE=MESH_COMPONENT_NUMBER_VELOCITY
  419. ELSE
  420. MESH_COMPONENT_NUMBER_PRESSURE=MESH_COMPONENT_NUMBER_VELOCITY+1
  421. CALL CMISSMeshElements_CreateStart(Mesh,MESH_COMPONENT_NUMBER_PRESSURE,BasisPressure,MeshElementsPressure,Err)
  422. DO ELEMENT_NUMBER=1,TOTAL_NUMBER_OF_ELEMENTS
  423. CALL CMISSMeshElements_NodesSet(MeshElementsPressure,ELEMENT_NUMBER,CM%P(ELEMENT_NUMBER, &
  424. & 1:NUMBER_OF_ELEMENT_NODES_PRESSURE),Err)
  425. ENDDO
  426. CALL CMISSMeshElements_CreateFinish(MeshElementsPressure,Err)
  427. ENDIF
  428. !Finish the creation of the mesh
  429. CALL CMISSMesh_CreateFinish(Mesh,Err)
  430. !
  431. !================================================================================================================================
  432. !
  433. !GEOMETRIC FIELD
  434. !Create a decomposition
  435. CALL CMISSDecomposition_Initialise(Decomposition,Err)
  436. CALL CMISSDecomposition_CreateStart(DecompositionUserNumber,Mesh,Decomposition,Err)
  437. !Set the decomposition to be a general decomposition with the specified number of domains
  438. CALL CMISSDecomposition_TypeSet(Decomposition,CMISS_DECOMPOSITION_USER_DEFINED_TYPE,Err)
  439. CALL CMISSDecomposition_NumberOfDomainsSet(Decomposition,NumberOfComputationalNodes,Err)
  440. CALL CMISSDecomposition_ElementDomainSet(Decomposition,1,0,Err)
  441. CALL CMISSDecomposition_ElementDomainSet(Decomposition,2,0,Err)
  442. CALL CMISSDecomposition_ElementDomainSet(Decomposition,3,0,Err)
  443. CALL CMISSDecomposition_ElementDomainSet(Decomposition,4,1,Err)
  444. CALL CMISSDecomposition_ElementDomainSet(Decomposition,5,1,Err)
  445. CALL CMISSDecomposition_ElementDomainSet(Decomposition,6,1,Err)
  446. CALL CMISSDecomposition_ElementDomainSet(Decomposition,7,1,Err)
  447. CALL CMISSDecomposition_ElementDomainSet(Decomposition,8,1,Err)
  448. !Finish the decomposition
  449. CALL CMISSDecomposition_CreateFinish(Decomposition,Err)
  450. !Start to create a default (geometric) field on the region
  451. CALL CMISSField_Initialise(GeometricField,Err)
  452. CALL CMISSField_CreateStart(GeometricFieldUserNumber,Region,GeometricField,Err)
  453. !Set the field type
  454. CALL CMISSField_TypeSet(GeometricField,CMISS_FIELD_GEOMETRIC_TYPE,Err)
  455. !Set the decomposition to use
  456. CALL CMISSField_MeshDecompositionSet(GeometricField,Decomposition,Err)
  457. !Set the scaling to use
  458. CALL CMISSField_ScalingTypeSet(GeometricField,CMISS_FIELD_NO_SCALING,Err)
  459. !Set the mesh component to be used by the field components.
  460. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS
  461. CALL CMISSField_ComponentMeshComponentSet(GeometricField,CMISS_FIELD_U_VARIABLE_TYPE,COMPONENT_NUMBER, &
  462. & MESH_COMPONENT_NUMBER_SPACE,Err)
  463. ENDDO
  464. !Finish creating the field
  465. CALL CMISSField_CreateFinish(GeometricField,Err)
  466. !Update the geometric field parameters
  467. DO NODE_NUMBER=1,NUMBER_OF_NODES_SPACE
  468. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS
  469. VALUE=CM%N(NODE_NUMBER,COMPONENT_NUMBER)
  470. CALL CMISSDecomposition_NodeDomainGet(Decomposition,NODE_NUMBER,1,BoundaryNodeDomain,Err)
  471. IF(BoundaryNodeDomain==ComputationalNodeNumber) THEN
  472. CALL CMISSField_ParameterSetUpdateNode(GeometricField,CMISS_FIELD_U_VARIABLE_TYPE,CMISS_FIELD_VALUES_SET_TYPE, &
  473. & 1,CMISS_NO_GLOBAL_DERIV,NODE_NUMBER,COMPONENT_NUMBER,VALUE,Err)
  474. ENDIF
  475. ENDDO
  476. ENDDO
  477. CALL CMISSField_ParameterSetUpdateStart(GeometricField,CMISS_FIELD_U_VARIABLE_TYPE,CMISS_FIELD_VALUES_SET_TYPE,Err)
  478. CALL CMISSField_ParameterSetUpdateFinish(GeometricField,CMISS_FIELD_U_VARIABLE_TYPE,CMISS_FIELD_VALUES_SET_TYPE,Err)
  479. !
  480. !================================================================================================================================
  481. !
  482. !EQUATIONS SETS
  483. !Create the equations set for static Navier-Stokes
  484. CALL CMISSEquationsSet_Initialise(EquationsSetNavierStokes,Err)
  485. CALL CMISSField_Initialise(EquationsSetField,Err)
  486. CALL CMISSEquationsSet_CreateStart(EquationsSetUserNumberNavierStokes,Region,GeometricField, &
  487. & CMISS_EQUATIONS_SET_FLUID_MECHANICS_CLASS,CMISS_EQUATIONS_SET_NAVIER_STOKES_EQUATION_TYPE, &
  488. & CMISS_EQUATIONS_SET_STATIC_NAVIER_STOKES_SUBTYPE, &
  489. & EquationsSetFieldUserNumber,EquationsSetField,EquationsSetNavierStokes,Err)
  490. !Set the equations set to be a static Navier-Stokes problem
  491. !Finish creating the equations set
  492. CALL CMISSEquationsSet_CreateFinish(EquationsSetNavierStokes,Err)
  493. !
  494. !================================================================================================================================
  495. !
  496. !DEPENDENT FIELDS
  497. !Create the equations set dependent field variables for static Navier-Stokes
  498. CALL CMISSField_Initialise(DependentFieldNavierStokes,Err)
  499. CALL CMISSEquationsSet_DependentCreateStart(EquationsSetNavierStokes,DependentFieldUserNumberNavierStokes, &
  500. & DependentFieldNavierStokes,Err)
  501. !Set the mesh component to be used by the field components.
  502. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS
  503. CALL CMISSField_ComponentMeshComponentSet(DependentFieldNavierStokes,CMISS_FIELD_U_VARIABLE_TYPE,COMPONENT_NUMBER, &
  504. & MESH_COMPONENT_NUMBER_VELOCITY,Err)
  505. CALL CMISSField_ComponentMeshComponentSet(DependentFieldNavierStokes,CMISS_FIELD_DELUDELN_VARIABLE_TYPE,COMPONENT_NUMBER, &
  506. & MESH_COMPONENT_NUMBER_VELOCITY,Err)
  507. ENDDO
  508. COMPONENT_NUMBER=NUMBER_OF_DIMENSIONS+1
  509. CALL CMISSField_ComponentMeshComponentSet(DependentFieldNavierStokes,CMISS_FIELD_U_VARIABLE_TYPE,COMPONENT_NUMBER, &
  510. & MESH_COMPONENT_NUMBER_PRESSURE,Err)
  511. CALL CMISSField_ComponentMeshComponentSet(DependentFieldNavierStokes,CMISS_FIELD_DELUDELN_VARIABLE_TYPE,COMPONENT_NUMBER, &
  512. & MESH_COMPONENT_NUMBER_PRESSURE,Err)
  513. !Finish the equations set dependent field variables
  514. CALL CMISSEquationsSet_DependentCreateFinish(EquationsSetNavierStokes,Err)
  515. !Initialise dependent field
  516. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS
  517. CALL CMISSField_ComponentValuesInitialise(DependentFieldNavierStokes,CMISS_FIELD_U_VARIABLE_TYPE,CMISS_FIELD_VALUES_SET_TYPE, &
  518. & COMPONENT_NUMBER,INITIAL_FIELD_NAVIER_STOKES(COMPONENT_NUMBER),Err)
  519. ENDDO
  520. !
  521. !================================================================================================================================
  522. !
  523. !MATERIALS FIELDS
  524. !Create the equations set materials field variables for static Navier-Stokes
  525. CALL CMISSField_Initialise(MaterialsFieldNavierStokes,Err)
  526. CALL CMISSEquationsSet_MaterialsCreateStart(EquationsSetNavierStokes,MaterialsFieldUserNumberNavierStokes, &
  527. & MaterialsFieldNavierStokes,Err)
  528. !Finish the equations set materials field variables
  529. CALL CMISSEquationsSet_MaterialsCreateFinish(EquationsSetNavierStokes,Err)
  530. CALL CMISSField_ComponentValuesInitialise(MaterialsFieldNavierStokes,CMISS_FIELD_U_VARIABLE_TYPE,CMISS_FIELD_VALUES_SET_TYPE, &
  531. & MaterialsFieldUserNumberNavierStokesMu,MU_PARAM_NAVIER_STOKES,Err)
  532. CALL CMISSField_ComponentValuesInitialise(MaterialsFieldNavierStokes,CMISS_FIELD_U_VARIABLE_TYPE,CMISS_FIELD_VALUES_SET_TYPE, &
  533. & MaterialsFieldUserNumberNavierStokesRho,RHO_PARAM_NAVIER_STOKES,Err)
  534. !
  535. !================================================================================================================================
  536. !
  537. !EQUATIONS
  538. !Create the equations set equations
  539. CALL CMISSEquations_Initialise(EquationsNavierStokes,Err)
  540. CALL CMISSEquationsSet_EquationsCreateStart(EquationsSetNavierStokes,EquationsNavierStokes,Err)
  541. !Set the equations matrices sparsity type
  542. CALL CMISSEquations_SparsityTypeSet(EquationsNavierStokes,CMISS_EQUATIONS_SPARSE_MATRICES,Err)
  543. !Set the equations set output
  544. CALL CMISSEquations_OutputTypeSet(EquationsNavierStokes,EQUATIONS_NAVIER_STOKES_OUTPUT,Err)
  545. !Finish the equations set equations
  546. CALL CMISSEquationsSet_EquationsCreateFinish(EquationsSetNavierStokes,Err)
  547. !
  548. !================================================================================================================================
  549. !
  550. !PROBLEMS
  551. !Start the creation of a problem.
  552. CALL CMISSProblem_Initialise(Problem,Err)
  553. CALL CMISSControlLoop_Initialise(ControlLoop,Err)
  554. CALL CMISSProblem_CreateStart(ProblemUserNumber,Problem,Err)
  555. !Set the problem to be a static Navier-Stokes problem
  556. CALL CMISSProblem_SpecificationSet(Problem,CMISS_PROBLEM_FLUID_MECHANICS_CLASS,CMISS_PROBLEM_NAVIER_STOKES_EQUATION_TYPE, &
  557. & CMISS_PROBLEM_STATIC_NAVIER_STOKES_SUBTYPE,Err)
  558. !Finish the creation of a problem.
  559. CALL CMISSProblem_CreateFinish(Problem,Err)
  560. !Start the creation of the problem control loop
  561. CALL CMISSProblem_ControlLoopCreateStart(Problem,Err)
  562. !Finish creating the problem control loop
  563. CALL CMISSProblem_ControlLoopCreateFinish(Problem,Err)
  564. !
  565. !================================================================================================================================
  566. !
  567. !SOLVERS
  568. !Start the creation of the problem solvers
  569. CALL CMISSSolver_Initialise(NonlinearSolverNavierStokes,Err)
  570. CALL CMISSSolver_Initialise(LinearSolverNavierStokes,Err)
  571. CALL CMISSProblem_SolversCreateStart(Problem,Err)
  572. !Get the nonlinear static solver
  573. CALL CMISSProblem_SolverGet(Problem,CMISS_CONTROL_LOOP_NODE,SolverNavierStokesUserNumber,NonlinearSolverNavierStokes,Err)
  574. !Set the nonlinear Jacobian type
  575. ! CALL CMISSSolver_NewtonJacobianCalculationTypeSet(NonlinearSolverNavierStokes,CMISS_SOLVER_NEWTON_JACOBIAN_EQUATIONS_CALCULATED,Err)
  576. !Set the output type
  577. CALL CMISSSolver_OutputTypeSet(NonlinearSolverNavierStokes,NONLINEAR_SOLVER_NAVIER_STOKES_OUTPUT_TYPE,Err)
  578. !Set the solver settings
  579. CALL CMISSSolver_NewtonAbsoluteToleranceSet(NonlinearSolverNavierStokes,ABSOLUTE_TOLERANCE,Err)
  580. CALL CMISSSolver_NewtonRelativeToleranceSet(NonlinearSolverNavierStokes,RELATIVE_TOLERANCE,Err)
  581. !Get the nonlinear linear solver
  582. CALL CMISSSolver_NewtonLinearSolverGet(NonlinearSolverNavierStokes,LinearSolverNavierStokes,Err)
  583. !Set the output type
  584. CALL CMISSSolver_OutputTypeSet(LinearSolverNavierStokes,LINEAR_SOLVER_NAVIER_STOKES_OUTPUT_TYPE,Err)
  585. !Set the solver settings
  586. IF(LINEAR_SOLVER_NAVIER_STOKES_DIRECT_FLAG) THEN
  587. CALL CMISSSolver_LinearTypeSet(LinearSolverNavierStokes,CMISS_SOLVER_LINEAR_DIRECT_SOLVE_TYPE,Err)
  588. CALL CMISSSolver_LibraryTypeSet(LinearSolverNavierStokes,CMISS_SOLVER_MUMPS_LIBRARY,Err)
  589. ELSE
  590. CALL CMISSSolver_LinearTypeSet(LinearSolverNavierStokes,CMISS_SOLVER_LINEAR_ITERATIVE_SOLVE_TYPE,Err)
  591. CALL CMISSSolver_LinearIterativeMaximumIterationsSet(LinearSolverNavierStokes,MAXIMUM_ITERATIONS,Err)
  592. CALL CMISSSolver_LinearIterativeDivergenceToleranceSet(LinearSolverNavierStokes,DIVERGENCE_TOLERANCE,Err)
  593. CALL CMISSSolver_LinearIterativeRelativeToleranceSet(LinearSolverNavierStokes,RELATIVE_TOLERANCE,Err)
  594. CALL CMISSSolver_LinearIterativeAbsoluteToleranceSet(LinearSolverNavierStokes,ABSOLUTE_TOLERANCE,Err)
  595. CALL CMISSSolver_LinearIterativeGMRESRestartSet(LinearSolverNavierStokes,RESTART_VALUE,Err)
  596. ENDIF
  597. !Finish the creation of the problem solver
  598. CALL CMISSProblem_SolversCreateFinish(Problem,Err)
  599. !
  600. !================================================================================================================================
  601. !
  602. !SOLVER EQUATIONS
  603. !Start the creation of the problem solver equations
  604. CALL CMISSSolver_Initialise(LinearSolverNavierStokes,Err)
  605. CALL CMISSSolverEquations_Initialise(SolverEquationsNavierStokes,Err)
  606. CALL CMISSProblem_SolverEquationsCreateStart(Problem,Err)
  607. !Get the linear solver equations
  608. CALL CMISSProblem_SolverGet(Problem,CMISS_CONTROL_LOOP_NODE,SolverNavierStokesUserNumber,LinearSolverNavierStokes,Err)
  609. CALL CMISSSolver_SolverEquationsGet(LinearSolverNavierStokes,SolverEquationsNavierStokes,Err)
  610. !Set the solver equations sparsity
  611. CALL CMISSSolverEquations_SparsityTypeSet(SolverEquationsNavierStokes,CMISS_SOLVER_SPARSE_MATRICES,Err)
  612. !Add in the equations set
  613. CALL CMISSSolverEquations_EquationsSetAdd(SolverEquationsNavierStokes,EquationsSetNavierStokes,EquationsSetIndex,Err)
  614. !Finish the creation of the problem solver equations
  615. CALL CMISSProblem_SolverEquationsCreateFinish(Problem,Err)
  616. !
  617. !================================================================================================================================
  618. !
  619. !BOUNDARY CONDITIONS
  620. !Start the creation of the equations set boundary conditions for Stokes
  621. CALL CMISSBoundaryConditions_Initialise(BoundaryConditionsNavierStokes,Err)
  622. CALL CMISSSolverEquations_BoundaryConditionsCreateStart(SolverEquationsNavierStokes,BoundaryConditionsNavierStokes,Err)
  623. !Set fixed wall nodes
  624. IF(FIXED_WALL_NODES_NAVIER_STOKES_FLAG) THEN
  625. DO NODE_COUNTER=1,NUMBER_OF_FIXED_WALL_NODES_NAVIER_STOKES
  626. NODE_NUMBER=FIXED_WALL_NODES_NAVIER_STOKES(NODE_COUNTER)
  627. CONDITION=CMISS_BOUNDARY_CONDITION_FIXED_WALL
  628. CALL CMISSDecomposition_NodeDomainGet(Decomposition,NODE_NUMBER,1,BoundaryNodeDomain,Err)
  629. IF(BoundaryNodeDomain==ComputationalNodeNumber) THEN
  630. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS
  631. VALUE=0.0_CMISSDP
  632. CALL CMISSBoundaryConditions_SetNode(BoundaryConditionsNavierStokes,DependentFieldNavierStokes, &
  633. & CMISS_FIELD_U_VARIABLE_TYPE,1, &
  634. & CMISS_NO_GLOBAL_DERIV, &
  635. & NODE_NUMBER,COMPONENT_NUMBER,CONDITION,VALUE,Err)
  636. ENDDO
  637. ENDIF
  638. ENDDO
  639. ENDIF
  640. !Set velocity boundary conditions
  641. IF(INLET_WALL_NODES_NAVIER_STOKES_FLAG) THEN
  642. DO NODE_COUNTER=1,NUMBER_OF_INLET_WALL_NODES_NAVIER_STOKES
  643. NODE_NUMBER=INLET_WALL_NODES_NAVIER_STOKES(NODE_COUNTER)
  644. CONDITION=CMISS_BOUNDARY_CONDITION_FIXED_INLET
  645. CALL CMISSDecomposition_NodeDomainGet(Decomposition,NODE_NUMBER,1,BoundaryNodeDomain,Err)
  646. IF(BoundaryNodeDomain==ComputationalNodeNumber) THEN
  647. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS
  648. VALUE=BOUNDARY_CONDITIONS_NAVIER_STOKES(COMPONENT_NUMBER)
  649. CALL CMISSBoundaryConditions_SetNode(BoundaryConditionsNavierStokes,DependentFieldNavierStokes, &
  650. & CMISS_FIELD_U_VARIABLE_TYPE,1, &
  651. & CMISS_NO_GLOBAL_DERIV, &
  652. & NODE_NUMBER,COMPONENT_NUMBER,CONDITION,VALUE,Err)
  653. ENDDO
  654. ENDIF
  655. ENDDO
  656. ENDIF
  657. !Finish the creation of the equations set boundary conditions
  658. CALL CMISSSolverEquations_BoundaryConditionsCreateFinish(SolverEquationsNavierStokes,Err)
  659. !
  660. !================================================================================================================================
  661. !
  662. !RUN SOLVERS
  663. !Turn of PETSc error handling
  664. !CALL PETSC_ERRORHANDLING_SET_ON(ERR,ERROR,*999)
  665. !Solve the problem
  666. WRITE(*,'(A)') "Solving problem..."
  667. CALL CMISSProblem_Solve(Problem,Err)
  668. WRITE(*,'(A)') "Problem solved!"
  669. !
  670. !
  671. !================================================================================================================================
  672. !
  673. !OUTPUT
  674. EXPORT_FIELD_IO=.TRUE.
  675. IF(EXPORT_FIELD_IO) THEN
  676. WRITE(*,'(A)') "Exporting fields..."
  677. CALL CMISSFields_Initialise(Fields,Err)
  678. CALL CMISSFields_Create(Region,Fields,Err)
  679. CALL CMISSFields_NodesExport(Fields,"StaticNavierStokes","FORTRAN",Err)
  680. CALL CMISSFields_ElementsExport(Fields,"StaticNavierStokes","FORTRAN",Err)
  681. CALL CMISSFields_Finalise(Fields,Err)
  682. WRITE(*,'(A)') "Field exported!"
  683. ENDIF
  684. !Finialise CMISS
  685. CALL CMISSFinalise(Err)
  686. WRITE(*,'(A)') "Program successfully completed."
  687. STOP
  688. END PROGRAM NAVIERSTOKESSTATICEXAMPLE