/FluidMechanics/NavierStokes/RoutineCheck/Static/src/StaticExample.f90

http://github.com/xyan075/examples · Fortran Modern · 785 lines · 455 code · 103 blank · 227 comment · 0 complexity · fbb0fe1843695166385c878676b3cb7b 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. CALL CMISSErrorHandlingModeSet(CMISS_ERRORS_TRAP_ERROR,Err)
  201. !
  202. !================================================================================================================================
  203. !
  204. !CHECK COMPUTATIONAL NODE
  205. !Get the computational nodes information
  206. CALL CMISSComputationalNumberOfNodesGet(NumberOfComputationalNodes,Err)
  207. CALL CMISSComputationalNodeNumberGet(ComputationalNodeNumber,Err)
  208. !
  209. !================================================================================================================================
  210. !
  211. !PROBLEM CONTROL PANEL
  212. !Import cmHeart mesh information
  213. CALL FLUID_MECHANICS_IO_READ_CMHEART(CM,Err)
  214. BASIS_NUMBER_SPACE=CM%ID_M
  215. BASIS_NUMBER_VELOCITY=CM%ID_V
  216. BASIS_NUMBER_PRESSURE=CM%ID_P
  217. NUMBER_OF_DIMENSIONS=CM%D
  218. BASIS_TYPE=CM%IT_T
  219. BASIS_XI_INTERPOLATION_SPACE=CM%IT_M
  220. BASIS_XI_INTERPOLATION_VELOCITY=CM%IT_V
  221. BASIS_XI_INTERPOLATION_PRESSURE=CM%IT_P
  222. NUMBER_OF_NODES_SPACE=CM%N_M
  223. NUMBER_OF_NODES_VELOCITY=CM%N_V
  224. NUMBER_OF_NODES_PRESSURE=CM%N_P
  225. TOTAL_NUMBER_OF_NODES=CM%N_T
  226. TOTAL_NUMBER_OF_ELEMENTS=CM%E_T
  227. NUMBER_OF_ELEMENT_NODES_SPACE=CM%EN_M
  228. NUMBER_OF_ELEMENT_NODES_VELOCITY=CM%EN_V
  229. NUMBER_OF_ELEMENT_NODES_PRESSURE=CM%EN_P
  230. !Set initial values
  231. INITIAL_FIELD_NAVIER_STOKES(1)=0.0_CMISSDP
  232. INITIAL_FIELD_NAVIER_STOKES(2)=0.0_CMISSDP
  233. INITIAL_FIELD_NAVIER_STOKES(3)=0.0_CMISSDP
  234. !Set boundary conditions
  235. FIXED_WALL_NODES_NAVIER_STOKES_FLAG=.TRUE.
  236. INLET_WALL_NODES_NAVIER_STOKES_FLAG=.TRUE.
  237. IF(FIXED_WALL_NODES_NAVIER_STOKES_FLAG) THEN
  238. NUMBER_OF_FIXED_WALL_NODES_NAVIER_STOKES=80
  239. ALLOCATE(FIXED_WALL_NODES_NAVIER_STOKES(NUMBER_OF_FIXED_WALL_NODES_NAVIER_STOKES))
  240. 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, &
  241. & 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, &
  242. & 89,90,91,92,93,94,95,97,99,101,102,103,104,105,106,107,108,111,114,115,116,117,118, &
  243. & 120,122,123,124,125/)
  244. ENDIF
  245. IF(INLET_WALL_NODES_NAVIER_STOKES_FLAG) THEN
  246. NUMBER_OF_INLET_WALL_NODES_NAVIER_STOKES=9
  247. ALLOCATE(INLET_WALL_NODES_NAVIER_STOKES(NUMBER_OF_INLET_WALL_NODES_NAVIER_STOKES))
  248. INLET_WALL_NODES_NAVIER_STOKES=(/6,15,16,23,36,42,81,82,96/)
  249. !Set initial boundary conditions
  250. BOUNDARY_CONDITIONS_NAVIER_STOKES(1)=0.0_CMISSDP
  251. BOUNDARY_CONDITIONS_NAVIER_STOKES(2)=1.0_CMISSDP
  252. BOUNDARY_CONDITIONS_NAVIER_STOKES(3)=0.0_CMISSDP
  253. ENDIF
  254. !Set material parameters
  255. MU_PARAM_NAVIER_STOKES=1.0_CMISSDP
  256. RHO_PARAM_NAVIER_STOKES=1.0_CMISSDP
  257. !Set interpolation parameters
  258. BASIS_XI_GAUSS_SPACE=3
  259. BASIS_XI_GAUSS_VELOCITY=3
  260. BASIS_XI_GAUSS_PRESSURE=3
  261. !Set output parameter
  262. !(NoOutput/ProgressOutput/TimingOutput/SolverOutput/SolverMatrixOutput)
  263. LINEAR_SOLVER_NAVIER_STOKES_OUTPUT_TYPE=CMISS_SOLVER_PROGRESS_OUTPUT
  264. NONLINEAR_SOLVER_NAVIER_STOKES_OUTPUT_TYPE=CMISS_SOLVER_PROGRESS_OUTPUT
  265. !(NoOutput/TimingOutput/MatrixOutput/ElementOutput)
  266. EQUATIONS_NAVIER_STOKES_OUTPUT=CMISS_EQUATIONS_NO_OUTPUT
  267. !Set solver parameters
  268. LINEAR_SOLVER_NAVIER_STOKES_DIRECT_FLAG=.FALSE.
  269. RELATIVE_TOLERANCE=1.0E-10_CMISSDP !default: 1.0E-05_CMISSDP
  270. ABSOLUTE_TOLERANCE=1.0E-10_CMISSDP !default: 1.0E-10_CMISSDP
  271. DIVERGENCE_TOLERANCE=1.0E20 !default: 1.0E5
  272. MAXIMUM_ITERATIONS=100000 !default: 100000
  273. RESTART_VALUE=3000 !default: 30
  274. LINESEARCH_ALPHA=1.0
  275. !
  276. !================================================================================================================================
  277. !
  278. !COORDINATE SYSTEM
  279. !Start the creation of a new RC coordinate system
  280. CALL CMISSCoordinateSystem_Initialise(CoordinateSystem,Err)
  281. CALL CMISSCoordinateSystem_CreateStart(CoordinateSystemUserNumber,CoordinateSystem,Err)
  282. !Set the coordinate system dimension
  283. CALL CMISSCoordinateSystem_DimensionSet(CoordinateSystem,NUMBER_OF_DIMENSIONS,Err)
  284. !Finish the creation of the coordinate system
  285. CALL CMISSCoordinateSystem_CreateFinish(CoordinateSystem,Err)
  286. !
  287. !================================================================================================================================
  288. !
  289. !REGION
  290. !Start the creation of a new region
  291. CALL CMISSRegion_Initialise(Region,Err)
  292. CALL CMISSRegion_CreateStart(RegionUserNumber,WorldRegion,Region,Err)
  293. !Set the regions coordinate system as defined above
  294. CALL CMISSRegion_CoordinateSystemSet(Region,CoordinateSystem,Err)
  295. !Finish the creation of the region
  296. CALL CMISSRegion_CreateFinish(Region,Err)
  297. !
  298. !================================================================================================================================
  299. !
  300. !BASES
  301. !Start the creation of new bases
  302. MESH_NUMBER_OF_COMPONENTS=1
  303. CALL CMISSBasis_Initialise(BasisSpace,Err)
  304. CALL CMISSBasis_CreateStart(BASIS_NUMBER_SPACE,BasisSpace,Err)
  305. !Set the basis type (Lagrange/Simplex)
  306. CALL CMISSBasis_TypeSet(BasisSpace,BASIS_TYPE,Err)
  307. !Set the basis xi number
  308. CALL CMISSBasis_NumberOfXiSet(BasisSpace,NUMBER_OF_DIMENSIONS,Err)
  309. !Set the basis xi interpolation and number of Gauss points
  310. IF(NUMBER_OF_DIMENSIONS==2) THEN
  311. CALL CMISSBasis_InterpolationXiSet(BasisSpace,(/BASIS_XI_INTERPOLATION_SPACE,BASIS_XI_INTERPOLATION_SPACE/),Err)
  312. CALL CMISSBasis_QuadratureNumberOfGaussXiSet(BasisSpace,(/BASIS_XI_GAUSS_SPACE,BASIS_XI_GAUSS_SPACE/),Err)
  313. ELSE IF(NUMBER_OF_DIMENSIONS==3) THEN
  314. CALL CMISSBasis_InterpolationXiSet(BasisSpace,(/BASIS_XI_INTERPOLATION_SPACE,BASIS_XI_INTERPOLATION_SPACE, &
  315. & BASIS_XI_INTERPOLATION_SPACE/),Err)
  316. CALL CMISSBasis_QuadratureNumberOfGaussXiSet(BasisSpace,(/BASIS_XI_GAUSS_SPACE,BASIS_XI_GAUSS_SPACE,BASIS_XI_GAUSS_SPACE/),Err)
  317. ENDIF
  318. !Finish the creation of the basis
  319. CALL CMISSBasis_CreateFinish(BasisSpace,Err)
  320. !Start the creation of another basis
  321. IF(BASIS_XI_INTERPOLATION_VELOCITY==BASIS_XI_INTERPOLATION_SPACE) THEN
  322. BasisVelocity=BasisSpace
  323. ELSE
  324. MESH_NUMBER_OF_COMPONENTS=MESH_NUMBER_OF_COMPONENTS+1
  325. !Initialise a new velocity basis
  326. CALL CMISSBasis_Initialise(BasisVelocity,Err)
  327. !Start the creation of a basis
  328. CALL CMISSBasis_CreateStart(BASIS_NUMBER_VELOCITY,BasisVelocity,Err)
  329. !Set the basis type (Lagrange/Simplex)
  330. CALL CMISSBasis_TypeSet(BasisVelocity,BASIS_TYPE,Err)
  331. !Set the basis xi number
  332. CALL CMISSBasis_NumberOfXiSet(BasisVelocity,NUMBER_OF_DIMENSIONS,Err)
  333. !Set the basis xi interpolation and number of Gauss points
  334. IF(NUMBER_OF_DIMENSIONS==2) THEN
  335. CALL CMISSBasis_InterpolationXiSet(BasisVelocity,(/BASIS_XI_INTERPOLATION_VELOCITY,BASIS_XI_INTERPOLATION_VELOCITY/),Err)
  336. CALL CMISSBasis_QuadratureNumberOfGaussXiSet(BasisVelocity,(/BASIS_XI_GAUSS_VELOCITY,BASIS_XI_GAUSS_VELOCITY/),Err)
  337. ELSE IF(NUMBER_OF_DIMENSIONS==3) THEN
  338. CALL CMISSBasis_InterpolationXiSet(BasisVelocity,(/BASIS_XI_INTERPOLATION_VELOCITY,BASIS_XI_INTERPOLATION_VELOCITY, &
  339. & BASIS_XI_INTERPOLATION_VELOCITY/),Err)
  340. CALL CMISSBasis_QuadratureNumberOfGaussXiSet(BasisVelocity,(/BASIS_XI_GAUSS_VELOCITY,BASIS_XI_GAUSS_VELOCITY, &
  341. & BASIS_XI_GAUSS_VELOCITY/),Err)
  342. ENDIF
  343. !Finish the creation of the basis
  344. CALL CMISSBasis_CreateFinish(BasisVelocity,Err)
  345. ENDIF
  346. !Start the creation of another basis
  347. IF(BASIS_XI_INTERPOLATION_PRESSURE==BASIS_XI_INTERPOLATION_SPACE) THEN
  348. BasisPressure=BasisSpace
  349. ELSE IF(BASIS_XI_INTERPOLATION_PRESSURE==BASIS_XI_INTERPOLATION_VELOCITY) THEN
  350. BasisPressure=BasisVelocity
  351. ELSE
  352. MESH_NUMBER_OF_COMPONENTS=MESH_NUMBER_OF_COMPONENTS+1
  353. !Initialise a new pressure basis
  354. CALL CMISSBasis_Initialise(BasisPressure,Err)
  355. !Start the creation of a basis
  356. CALL CMISSBasis_CreateStart(BASIS_NUMBER_PRESSURE,BasisPressure,Err)
  357. !Set the basis type (Lagrange/Simplex)
  358. CALL CMISSBasis_TypeSet(BasisPressure,BASIS_TYPE,Err)
  359. !Set the basis xi number
  360. CALL CMISSBasis_NumberOfXiSet(BasisPressure,NUMBER_OF_DIMENSIONS,Err)
  361. !Set the basis xi interpolation and number of Gauss points
  362. IF(NUMBER_OF_DIMENSIONS==2) THEN
  363. CALL CMISSBasis_InterpolationXiSet(BasisPressure,(/BASIS_XI_INTERPOLATION_PRESSURE,BASIS_XI_INTERPOLATION_PRESSURE/),Err)
  364. CALL CMISSBasis_QuadratureNumberOfGaussXiSet(BasisPressure,(/BASIS_XI_GAUSS_PRESSURE,BASIS_XI_GAUSS_PRESSURE/),Err)
  365. ELSE IF(NUMBER_OF_DIMENSIONS==3) THEN
  366. CALL CMISSBasis_InterpolationXiSet(BasisPressure,(/BASIS_XI_INTERPOLATION_PRESSURE,BASIS_XI_INTERPOLATION_PRESSURE, &
  367. & BASIS_XI_INTERPOLATION_PRESSURE/),Err)
  368. CALL CMISSBasis_QuadratureNumberOfGaussXiSet(BasisPressure,(/BASIS_XI_GAUSS_PRESSURE,BASIS_XI_GAUSS_PRESSURE, &
  369. & BASIS_XI_GAUSS_PRESSURE/),Err)
  370. ENDIF
  371. !Finish the creation of the basis
  372. CALL CMISSBasis_CreateFinish(BasisPressure,Err)
  373. ENDIF
  374. !
  375. !================================================================================================================================
  376. !
  377. !MESH
  378. !Start the creation of mesh nodes
  379. CALL CMISSNodes_Initialise(Nodes,Err)
  380. CALL CMISSMesh_Initialise(Mesh,Err)
  381. CALL CMISSNodes_CreateStart(Region,TOTAL_NUMBER_OF_NODES,Nodes,Err)
  382. CALL CMISSNodes_CreateFinish(Nodes,Err)
  383. !Start the creation of the mesh
  384. CALL CMISSMesh_CreateStart(MeshUserNumber,Region,NUMBER_OF_DIMENSIONS,Mesh,Err)
  385. !Set number of mesh elements
  386. CALL CMISSMesh_NumberOfElementsSet(Mesh,TOTAL_NUMBER_OF_ELEMENTS,Err)
  387. !Set number of mesh components
  388. CALL CMISSMesh_NumberOfComponentsSet(Mesh,MESH_NUMBER_OF_COMPONENTS,Err)
  389. !Specify spatial mesh component
  390. CALL CMISSMeshElements_Initialise(MeshElementsSpace,Err)
  391. CALL CMISSMeshElements_Initialise(MeshElementsVelocity,Err)
  392. CALL CMISSMeshElements_Initialise(MeshElementsPressure,Err)
  393. MESH_COMPONENT_NUMBER_SPACE=1
  394. MESH_COMPONENT_NUMBER_VELOCITY=1
  395. MESH_COMPONENT_NUMBER_PRESSURE=1
  396. CALL CMISSMeshElements_CreateStart(Mesh,MESH_COMPONENT_NUMBER_SPACE,BasisSpace,MeshElementsSpace,Err)
  397. DO ELEMENT_NUMBER=1,TOTAL_NUMBER_OF_ELEMENTS
  398. CALL CMISSMeshElements_NodesSet(MeshElementsSpace,ELEMENT_NUMBER,CM%M(ELEMENT_NUMBER,1:NUMBER_OF_ELEMENT_NODES_SPACE),Err)
  399. ENDDO
  400. CALL CMISSMeshElements_CreateFinish(MeshElementsSpace,Err)
  401. !Specify velocity mesh component
  402. IF(BASIS_XI_INTERPOLATION_VELOCITY==BASIS_XI_INTERPOLATION_SPACE) THEN
  403. MeshElementsVelocity=MeshElementsSpace
  404. ELSE
  405. MESH_COMPONENT_NUMBER_VELOCITY=MESH_COMPONENT_NUMBER_SPACE+1
  406. CALL CMISSMeshElements_CreateStart(Mesh,MESH_COMPONENT_NUMBER_VELOCITY,BasisVelocity,MeshElementsVelocity,Err)
  407. DO ELEMENT_NUMBER=1,TOTAL_NUMBER_OF_ELEMENTS
  408. CALL CMISSMeshElements_NodesSet(MeshElementsVelocity,ELEMENT_NUMBER,CM%V(ELEMENT_NUMBER, &
  409. & 1:NUMBER_OF_ELEMENT_NODES_VELOCITY),Err)
  410. ENDDO
  411. CALL CMISSMeshElements_CreateFinish(MeshElementsVelocity,Err)
  412. ENDIF
  413. !Specify pressure mesh component
  414. IF(BASIS_XI_INTERPOLATION_PRESSURE==BASIS_XI_INTERPOLATION_SPACE) THEN
  415. MeshElementsPressure=MeshElementsSpace
  416. MESH_COMPONENT_NUMBER_PRESSURE=MESH_COMPONENT_NUMBER_SPACE
  417. ELSE IF(BASIS_XI_INTERPOLATION_PRESSURE==BASIS_XI_INTERPOLATION_VELOCITY) THEN
  418. MeshElementsPressure=MeshElementsVelocity
  419. MESH_COMPONENT_NUMBER_PRESSURE=MESH_COMPONENT_NUMBER_VELOCITY
  420. ELSE
  421. MESH_COMPONENT_NUMBER_PRESSURE=MESH_COMPONENT_NUMBER_VELOCITY+1
  422. CALL CMISSMeshElements_CreateStart(Mesh,MESH_COMPONENT_NUMBER_PRESSURE,BasisPressure,MeshElementsPressure,Err)
  423. DO ELEMENT_NUMBER=1,TOTAL_NUMBER_OF_ELEMENTS
  424. CALL CMISSMeshElements_NodesSet(MeshElementsPressure,ELEMENT_NUMBER,CM%P(ELEMENT_NUMBER, &
  425. & 1:NUMBER_OF_ELEMENT_NODES_PRESSURE),Err)
  426. ENDDO
  427. CALL CMISSMeshElements_CreateFinish(MeshElementsPressure,Err)
  428. ENDIF
  429. !Finish the creation of the mesh
  430. CALL CMISSMesh_CreateFinish(Mesh,Err)
  431. !
  432. !================================================================================================================================
  433. !
  434. !GEOMETRIC FIELD
  435. !Create a decomposition
  436. CALL CMISSDecomposition_Initialise(Decomposition,Err)
  437. CALL CMISSDecomposition_CreateStart(DecompositionUserNumber,Mesh,Decomposition,Err)
  438. !Set the decomposition to be a general decomposition with the specified number of domains
  439. CALL CMISSDecomposition_TypeSet(Decomposition,CMISS_DECOMPOSITION_CALCULATED_TYPE,Err)
  440. CALL CMISSDecomposition_NumberOfDomainsSet(Decomposition,NumberOfComputationalNodes,Err)
  441. !Finish the decomposition
  442. CALL CMISSDecomposition_CreateFinish(Decomposition,Err)
  443. !Start to create a default (geometric) field on the region
  444. CALL CMISSField_Initialise(GeometricField,Err)
  445. CALL CMISSField_CreateStart(GeometricFieldUserNumber,Region,GeometricField,Err)
  446. !Set the field type
  447. CALL CMISSField_TypeSet(GeometricField,CMISS_FIELD_GEOMETRIC_TYPE,Err)
  448. !Set the decomposition to use
  449. CALL CMISSField_MeshDecompositionSet(GeometricField,Decomposition,Err)
  450. !Set the scaling to use
  451. CALL CMISSField_ScalingTypeSet(GeometricField,CMISS_FIELD_NO_SCALING,Err)
  452. !Set the mesh component to be used by the field components.
  453. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS
  454. CALL CMISSField_ComponentMeshComponentSet(GeometricField,CMISS_FIELD_U_VARIABLE_TYPE,COMPONENT_NUMBER, &
  455. & MESH_COMPONENT_NUMBER_SPACE,Err)
  456. ENDDO
  457. !Finish creating the field
  458. CALL CMISSField_CreateFinish(GeometricField,Err)
  459. !Update the geometric field parameters
  460. DO NODE_NUMBER=1,NUMBER_OF_NODES_SPACE
  461. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS
  462. VALUE=CM%N(NODE_NUMBER,COMPONENT_NUMBER)
  463. CALL CMISSDecomposition_NodeDomainGet(Decomposition,NODE_NUMBER,1,BoundaryNodeDomain,Err)
  464. IF(BoundaryNodeDomain==ComputationalNodeNumber) THEN
  465. CALL CMISSField_ParameterSetUpdateNode(GeometricField,CMISS_FIELD_U_VARIABLE_TYPE,CMISS_FIELD_VALUES_SET_TYPE, &
  466. & 1,CMISS_NO_GLOBAL_DERIV,NODE_NUMBER,COMPONENT_NUMBER,VALUE,Err)
  467. ENDIF
  468. ENDDO
  469. ENDDO
  470. CALL CMISSField_ParameterSetUpdateStart(GeometricField,CMISS_FIELD_U_VARIABLE_TYPE,CMISS_FIELD_VALUES_SET_TYPE,Err)
  471. CALL CMISSField_ParameterSetUpdateFinish(GeometricField,CMISS_FIELD_U_VARIABLE_TYPE,CMISS_FIELD_VALUES_SET_TYPE,Err)
  472. !
  473. !================================================================================================================================
  474. !
  475. !EQUATIONS SETS
  476. !Create the equations set for static Navier-Stokes
  477. CALL CMISSEquationsSet_Initialise(EquationsSetNavierStokes,Err)
  478. CALL CMISSField_Initialise(EquationsSetField,Err)
  479. CALL CMISSEquationsSet_CreateStart(EquationsSetUserNumberNavierStokes,Region,GeometricField, &
  480. & CMISS_EQUATIONS_SET_FLUID_MECHANICS_CLASS,CMISS_EQUATIONS_SET_NAVIER_STOKES_EQUATION_TYPE, &
  481. & CMISS_EQUATIONS_SET_STATIC_NAVIER_STOKES_SUBTYPE, &
  482. & EquationsSetFieldUserNumber,EquationsSetField,EquationsSetNavierStokes,Err)
  483. !Set the equations set to be a static Navier-Stokes problem
  484. !Finish creating the equations set
  485. CALL CMISSEquationsSet_CreateFinish(EquationsSetNavierStokes,Err)
  486. !
  487. !================================================================================================================================
  488. !
  489. !DEPENDENT FIELDS
  490. !Create the equations set dependent field variables for static Navier-Stokes
  491. CALL CMISSField_Initialise(DependentFieldNavierStokes,Err)
  492. CALL CMISSEquationsSet_DependentCreateStart(EquationsSetNavierStokes,DependentFieldUserNumberNavierStokes, &
  493. & DependentFieldNavierStokes,Err)
  494. !Set the mesh component to be used by the field components.
  495. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS
  496. CALL CMISSField_ComponentMeshComponentSet(DependentFieldNavierStokes,CMISS_FIELD_U_VARIABLE_TYPE,COMPONENT_NUMBER, &
  497. & MESH_COMPONENT_NUMBER_VELOCITY,Err)
  498. CALL CMISSField_ComponentMeshComponentSet(DependentFieldNavierStokes,CMISS_FIELD_DELUDELN_VARIABLE_TYPE,COMPONENT_NUMBER, &
  499. & MESH_COMPONENT_NUMBER_VELOCITY,Err)
  500. ENDDO
  501. COMPONENT_NUMBER=NUMBER_OF_DIMENSIONS+1
  502. CALL CMISSField_ComponentMeshComponentSet(DependentFieldNavierStokes,CMISS_FIELD_U_VARIABLE_TYPE,COMPONENT_NUMBER, &
  503. & MESH_COMPONENT_NUMBER_PRESSURE,Err)
  504. CALL CMISSField_ComponentMeshComponentSet(DependentFieldNavierStokes,CMISS_FIELD_DELUDELN_VARIABLE_TYPE,COMPONENT_NUMBER, &
  505. & MESH_COMPONENT_NUMBER_PRESSURE,Err)
  506. !Finish the equations set dependent field variables
  507. CALL CMISSEquationsSet_DependentCreateFinish(EquationsSetNavierStokes,Err)
  508. !Initialise dependent field
  509. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS
  510. CALL CMISSField_ComponentValuesInitialise(DependentFieldNavierStokes,CMISS_FIELD_U_VARIABLE_TYPE,CMISS_FIELD_VALUES_SET_TYPE, &
  511. & COMPONENT_NUMBER,INITIAL_FIELD_NAVIER_STOKES(COMPONENT_NUMBER),Err)
  512. ENDDO
  513. !
  514. !================================================================================================================================
  515. !
  516. !MATERIALS FIELDS
  517. !Create the equations set materials field variables for static Navier-Stokes
  518. CALL CMISSField_Initialise(MaterialsFieldNavierStokes,Err)
  519. CALL CMISSEquationsSet_MaterialsCreateStart(EquationsSetNavierStokes,MaterialsFieldUserNumberNavierStokes, &
  520. & MaterialsFieldNavierStokes,Err)
  521. !Finish the equations set materials field variables
  522. CALL CMISSEquationsSet_MaterialsCreateFinish(EquationsSetNavierStokes,Err)
  523. CALL CMISSField_ComponentValuesInitialise(MaterialsFieldNavierStokes,CMISS_FIELD_U_VARIABLE_TYPE,CMISS_FIELD_VALUES_SET_TYPE, &
  524. & MaterialsFieldUserNumberNavierStokesMu,MU_PARAM_NAVIER_STOKES,Err)
  525. CALL CMISSField_ComponentValuesInitialise(MaterialsFieldNavierStokes,CMISS_FIELD_U_VARIABLE_TYPE,CMISS_FIELD_VALUES_SET_TYPE, &
  526. & MaterialsFieldUserNumberNavierStokesRho,RHO_PARAM_NAVIER_STOKES,Err)
  527. !
  528. !================================================================================================================================
  529. !
  530. !EQUATIONS
  531. !Create the equations set equations
  532. CALL CMISSEquations_Initialise(EquationsNavierStokes,Err)
  533. CALL CMISSEquationsSet_EquationsCreateStart(EquationsSetNavierStokes,EquationsNavierStokes,Err)
  534. !Set the equations matrices sparsity type
  535. CALL CMISSEquations_SparsityTypeSet(EquationsNavierStokes,CMISS_EQUATIONS_SPARSE_MATRICES,Err)
  536. !Set the equations set output
  537. CALL CMISSEquations_OutputTypeSet(EquationsNavierStokes,EQUATIONS_NAVIER_STOKES_OUTPUT,Err)
  538. !Finish the equations set equations
  539. CALL CMISSEquationsSet_EquationsCreateFinish(EquationsSetNavierStokes,Err)
  540. !
  541. !================================================================================================================================
  542. !
  543. !PROBLEMS
  544. !Start the creation of a problem.
  545. CALL CMISSProblem_Initialise(Problem,Err)
  546. CALL CMISSControlLoop_Initialise(ControlLoop,Err)
  547. CALL CMISSProblem_CreateStart(ProblemUserNumber,Problem,Err)
  548. !Set the problem to be a static Navier-Stokes problem
  549. CALL CMISSProblem_SpecificationSet(Problem,CMISS_PROBLEM_FLUID_MECHANICS_CLASS,CMISS_PROBLEM_NAVIER_STOKES_EQUATION_TYPE, &
  550. & CMISS_PROBLEM_STATIC_NAVIER_STOKES_SUBTYPE,Err)
  551. !Finish the creation of a problem.
  552. CALL CMISSProblem_CreateFinish(Problem,Err)
  553. !Start the creation of the problem control loop
  554. CALL CMISSProblem_ControlLoopCreateStart(Problem,Err)
  555. !Finish creating the problem control loop
  556. CALL CMISSProblem_ControlLoopCreateFinish(Problem,Err)
  557. !
  558. !================================================================================================================================
  559. !
  560. !SOLVERS
  561. !Start the creation of the problem solvers
  562. CALL CMISSSolver_Initialise(NonlinearSolverNavierStokes,Err)
  563. CALL CMISSSolver_Initialise(LinearSolverNavierStokes,Err)
  564. CALL CMISSProblem_SolversCreateStart(Problem,Err)
  565. !Get the nonlinear static solver
  566. CALL CMISSProblem_SolverGet(Problem,CMISS_CONTROL_LOOP_NODE,SolverNavierStokesUserNumber,NonlinearSolverNavierStokes,Err)
  567. !Set the nonlinear Jacobian type
  568. CALL CMISSSolver_NewtonJacobianCalculationTypeSet(NonlinearSolverNavierStokes,CMISS_SOLVER_NEWTON_JACOBIAN_EQUATIONS_CALCULATED, &
  569. & Err)
  570. !Set the output type
  571. CALL CMISSSolver_OutputTypeSet(NonlinearSolverNavierStokes,NONLINEAR_SOLVER_NAVIER_STOKES_OUTPUT_TYPE,Err)
  572. !Set the solver settings
  573. CALL CMISSSolver_NewtonAbsoluteToleranceSet(NonlinearSolverNavierStokes,ABSOLUTE_TOLERANCE,Err)
  574. CALL CMISSSolver_NewtonRelativeToleranceSet(NonlinearSolverNavierStokes,RELATIVE_TOLERANCE,Err)
  575. !Get the nonlinear linear solver
  576. CALL CMISSSolver_NewtonLinearSolverGet(NonlinearSolverNavierStokes,LinearSolverNavierStokes,Err)
  577. !Set the output type
  578. CALL CMISSSolver_OutputTypeSet(LinearSolverNavierStokes,LINEAR_SOLVER_NAVIER_STOKES_OUTPUT_TYPE,Err)
  579. !Set the solver settings
  580. IF(LINEAR_SOLVER_NAVIER_STOKES_DIRECT_FLAG) THEN
  581. CALL CMISSSolver_LinearTypeSet(LinearSolverNavierStokes,CMISS_SOLVER_LINEAR_DIRECT_SOLVE_TYPE,Err)
  582. CALL CMISSSolver_LibraryTypeSet(LinearSolverNavierStokes,CMISS_SOLVER_MUMPS_LIBRARY,Err)
  583. ELSE
  584. CALL CMISSSolver_LinearTypeSet(LinearSolverNavierStokes,CMISS_SOLVER_LINEAR_ITERATIVE_SOLVE_TYPE,Err)
  585. CALL CMISSSolver_LinearIterativeMaximumIterationsSet(LinearSolverNavierStokes,MAXIMUM_ITERATIONS,Err)
  586. CALL CMISSSolver_LinearIterativeDivergenceToleranceSet(LinearSolverNavierStokes,DIVERGENCE_TOLERANCE,Err)
  587. CALL CMISSSolver_LinearIterativeRelativeToleranceSet(LinearSolverNavierStokes,RELATIVE_TOLERANCE,Err)
  588. CALL CMISSSolver_LinearIterativeAbsoluteToleranceSet(LinearSolverNavierStokes,ABSOLUTE_TOLERANCE,Err)
  589. CALL CMISSSolver_LinearIterativeGMRESRestartSet(LinearSolverNavierStokes,RESTART_VALUE,Err)
  590. ENDIF
  591. !Finish the creation of the problem solver
  592. CALL CMISSProblem_SolversCreateFinish(Problem,Err)
  593. !
  594. !================================================================================================================================
  595. !
  596. !SOLVER EQUATIONS
  597. !Start the creation of the problem solver equations
  598. CALL CMISSSolver_Initialise(LinearSolverNavierStokes,Err)
  599. CALL CMISSSolverEquations_Initialise(SolverEquationsNavierStokes,Err)
  600. CALL CMISSProblem_SolverEquationsCreateStart(Problem,Err)
  601. !Get the linear solver equations
  602. CALL CMISSProblem_SolverGet(Problem,CMISS_CONTROL_LOOP_NODE,SolverNavierStokesUserNumber,LinearSolverNavierStokes,Err)
  603. CALL CMISSSolver_SolverEquationsGet(LinearSolverNavierStokes,SolverEquationsNavierStokes,Err)
  604. !Set the solver equations sparsity
  605. CALL CMISSSolverEquations_SparsityTypeSet(SolverEquationsNavierStokes,CMISS_SOLVER_SPARSE_MATRICES,Err)
  606. !Add in the equations set
  607. CALL CMISSSolverEquations_EquationsSetAdd(SolverEquationsNavierStokes,EquationsSetNavierStokes,EquationsSetIndex,Err)
  608. !Finish the creation of the problem solver equations
  609. CALL CMISSProblem_SolverEquationsCreateFinish(Problem,Err)
  610. !
  611. !================================================================================================================================
  612. !
  613. !BOUNDARY CONDITIONS
  614. !Start the creation of the equations set boundary conditions for Stokes
  615. CALL CMISSBoundaryConditions_Initialise(BoundaryConditionsNavierStokes,Err)
  616. CALL CMISSSolverEquations_BoundaryConditionsCreateStart(SolverEquationsNavierStokes,BoundaryConditionsNavierStokes,Err)
  617. !Set fixed wall nodes
  618. IF(FIXED_WALL_NODES_NAVIER_STOKES_FLAG) THEN
  619. DO NODE_COUNTER=1,NUMBER_OF_FIXED_WALL_NODES_NAVIER_STOKES
  620. NODE_NUMBER=FIXED_WALL_NODES_NAVIER_STOKES(NODE_COUNTER)
  621. CONDITION=CMISS_BOUNDARY_CONDITION_FIXED_WALL
  622. CALL CMISSDecomposition_NodeDomainGet(Decomposition,NODE_NUMBER,1,BoundaryNodeDomain,Err)
  623. IF(BoundaryNodeDomain==ComputationalNodeNumber) THEN
  624. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS
  625. VALUE=0.0_CMISSDP
  626. CALL CMISSBoundaryConditions_SetNode(BoundaryConditionsNavierStokes,DependentFieldNavierStokes, &
  627. & CMISS_FIELD_U_VARIABLE_TYPE,1, &
  628. & CMISS_NO_GLOBAL_DERIV, &
  629. & NODE_NUMBER,COMPONENT_NUMBER,CONDITION,VALUE,Err)
  630. ENDDO
  631. ENDIF
  632. ENDDO
  633. ENDIF
  634. !Set velocity boundary conditions
  635. IF(INLET_WALL_NODES_NAVIER_STOKES_FLAG) THEN
  636. DO NODE_COUNTER=1,NUMBER_OF_INLET_WALL_NODES_NAVIER_STOKES
  637. NODE_NUMBER=INLET_WALL_NODES_NAVIER_STOKES(NODE_COUNTER)
  638. CONDITION=CMISS_BOUNDARY_CONDITION_FIXED_INLET
  639. CALL CMISSDecomposition_NodeDomainGet(Decomposition,NODE_NUMBER,1,BoundaryNodeDomain,Err)
  640. IF(BoundaryNodeDomain==ComputationalNodeNumber) THEN
  641. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS
  642. VALUE=BOUNDARY_CONDITIONS_NAVIER_STOKES(COMPONENT_NUMBER)
  643. CALL CMISSBoundaryConditions_SetNode(BoundaryConditionsNavierStokes,DependentFieldNavierStokes, &
  644. & CMISS_FIELD_U_VARIABLE_TYPE,1, &
  645. & CMISS_NO_GLOBAL_DERIV, &
  646. & NODE_NUMBER,COMPONENT_NUMBER,CONDITION,VALUE,Err)
  647. ENDDO
  648. ENDIF
  649. ENDDO
  650. ENDIF
  651. !Finish the creation of the equations set boundary conditions
  652. CALL CMISSSolverEquations_BoundaryConditionsCreateFinish(SolverEquationsNavierStokes,Err)
  653. !
  654. !================================================================================================================================
  655. !
  656. !RUN SOLVERS
  657. !Turn of PETSc error handling
  658. !CALL PETSC_ERRORHANDLING_SET_ON(ERR,ERROR,*999)
  659. !Solve the problem
  660. WRITE(*,'(A)') "Solving problem..."
  661. CALL CMISSProblem_Solve(Problem,Err)
  662. WRITE(*,'(A)') "Problem solved!"
  663. !
  664. !
  665. !================================================================================================================================
  666. !
  667. !OUTPUT
  668. EXPORT_FIELD_IO=.TRUE.
  669. IF(EXPORT_FIELD_IO) THEN
  670. WRITE(*,'(A)') "Exporting fields..."
  671. CALL CMISSFields_Initialise(Fields,Err)
  672. CALL CMISSFields_Create(Region,Fields,Err)
  673. CALL CMISSFields_NodesExport(Fields,"StaticNavierStokes","FORTRAN",Err)
  674. CALL CMISSFields_ElementsExport(Fields,"StaticNavierStokes","FORTRAN",Err)
  675. CALL CMISSFields_Finalise(Fields,Err)
  676. WRITE(*,'(A)') "Field exported!"
  677. ENDIF
  678. !Finialise CMISS
  679. CALL CMISSFinalise(Err)
  680. WRITE(*,'(A)') "Program successfully completed."
  681. STOP
  682. END PROGRAM NAVIERSTOKESSTATICEXAMPLE