mahmood19227
/
SimulationCore2
1
0
Fork 0
Code Issues 0 Pull Requests 0 Releases 0 Wiki Activity
Simulation Core
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
54 Commits
2 Branches
134 MiB
 
 
 
 
 
 
Tree: c8c7b4c604
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from 'c8c7b4c604'
${ noResults }
SimulationCore2/Equipments/BopStack/LOSS_INPUTS.i90

1270 lines
67 KiB
Raw Blame History 

  1. # 1 "/mnt/c/Projects/VSIM/SimulationCore2/Equipments/BopStack/LOSS_INPUTS.f90"

  2.  SUBROUTINE DEALLOCATE_ARRAYS()

  3. use SimulationVariables

  4. use UnityModule

  5. implicit none

  6. !===========================================================================

  7. !                   	RAMLINE MINOR LOSSES INPUT

  8. !===========================================================================

  9.  DEALLOCATE (data%State%BopStackInput%MINORS1,data%State%RamLine%MINORDIAMETER_RAMLINE,data%State%RamLine%AREAMINOR_RAMLINE  &

  10.             ,data%State%RamLine%LF_RAMLINE,data%State%RamLine%CV_RAMLINE,data%State%RamLine%NOTE_RAMLINE,data%State%RAMS%minlosspa,data%State%RAMS%minloss)

  11. !===========================================================================

  12. !                   	RAMLINE PIPNING LOSSES INPUT

  13. !===========================================================================

  14.  DEALLOCATE (data%State%BopStackInput%PIPINGS_RAMLINE,data%State%RamLine%DIAM_RAMLINE_INCH,  &

  15.             data%State%RamLine%AREA_RAMLINE,data%State%RamLine%LENGT_RAMLINE,data%State%RamLine%ROUGHNESS_RAMLINE,data%State%RamLine%RELROUGH_RAMLINE   &

  16.             ,data%State%RAMS%Re_ramline,data%State%RAMS%fric,data%State%RAMS%fricloss)

  17. !===========================================================================

  18. !                   	ANNULAR MINOR LOSSES INPUT

  19. !===========================================================================

  20.  DEALLOCATE (data%State%BopStackInput%MINORS_ANNULAR,data%State%AnnularComputational%MINORDIAMETER_ANNULARLINE,data%State%AnnularComputational%AREAMINOR_ANNULARLINE  &

  21.             ,data%State%AnnularComputational%LF_ANNULARLINE,data%State%AnnularComputational%CV_ANNULARLINE,data%State%AnnularComputational%NOTE_ANNULARLINE,data%State%AnnularComputational%minlosspa_ANNULAR,data%State%AnnularComputational%minloss_ANNULAR)

  22. !===========================================================================

  23. !                   	ANNULAR PIPNING LOSSES INPUT

  24. !===========================================================================

  25.  DEALLOCATE (data%State%BopStackInput%PIPINGS_ANNULAR,data%State%AnnularComputational%DIAM_ANNULARLINE_INCH,data%State%AnnularComputational%AREA_ANNULARLINE, &

  26.             data%State%AnnularComputational%LENGT_ANNULARLINE,data%State%AnnularComputational%ROUGHNESS_ANNULARLINE,data%State%AnnularComputational%RELROUGH_ANNULARLINE   &

  27.             ,data%State%AnnularComputational%Re_ANNULARline,data%State%AnnularComputational%fricANNULAR,data%State%AnnularComputational%friclossANNULAR)

  28. !===========================================================================

  29. !                   	AIR PUMP LOSSES INPUT

  30. !===========================================================================

  31. DEALLOCATE (data%State%BopStackInput%PIPINGS_AIRPUMP,data%State%RamLine%DIAM_AIR_INCH,  &

  32.             data%State%RamLine%Re_air,data%State%RamLine%AREA_AIR,data%State%RamLine%LENGT_AIR,data%State%RamLine%ROUGHNESS_AIRPLINE,data%State%RamLine%REL_ROUGHAIR,  &

  33.            data%State%RamLine%fric_air,data%State%RamLine%fricloss_air)

  34. !================================================================

  35. DEALLOCATE (data%State%BopStackInput%MINORS_AIRPUMP,data%State%RamLine%MINORDIAM_AIR_INCH,  &

  36.             data%State%RamLine%MINORAREA_AIR,data%State%RamLine%LF_AIR,data%State%RamLine%CV_AIR,data%State%RamLine%NOTE_AIR &

  37.             ,data%State%RamLine%minlosspa_air,data%State%RamLine%minloss_air)  

  38.  

  39. !===========================================================================

  40. !                   	DELAY ARRAYS

  41. !===========================================================================

  42. call data%State%AnnularComputational%Pannular_regDelay%Empty()

  43. 

  44. 

  45. 

  46.     END

  47. 

  48. 

  49. 

  50. 

  51. 

  52. 

  53. 

  54.     

  55.     

  56.     

  57.     

  58. SUBROUTINE LOSS_INPUTS()

  59. use SimulationVariables

  60. implicit none

  61. INTEGER I

  62. 

  63. !===========================================================================

  64. !                   	RAMLINE MINOR LOSSES INPUT

  65. !===========================================================================

  66. data%State%RamLine%NO_MINORSRAMLINE=34

  67. 

  68. ALLOCATE (data%State%BopStackInput%MINORS1(data%State%RamLine%NO_MINORSRAMLINE,4))

  69. 

  70. !               ID(INCH)     LF	      CV	    NOTE(BAR)   DESCRIPTION

  71. data%State%BopStackInput%MINORS1(1,1:4)=     (/2.,        2.,      0.,	    0./)      !Acc.tee

  72. data%State%BopStackInput%MINORS1(2,1:4)=     (/2.,        0.9,     0.,	    0./)      !Avg.acc.tee

  73. data%State%BopStackInput%MINORS1(3,1:4)=     (/2.,        0.9,     0.,	    0./)      !Avg.acc.tee

  74. data%State%BopStackInput%MINORS1(4,1:4)=     (/2.,        0.9,     0.,	    0./)      !Avg.acc.tee

  75. data%State%BopStackInput%MINORS1(5,1:4)=     (/2.,        0.9,     0.,	    0./)      !tee

  76. data%State%BopStackInput%MINORS1(6,1:4)=     (/2.,        2.,      0.,	    0./)      !tee

  77. data%State%BopStackInput%MINORS1(7,1:4)=     (/2.,        0.,      105.,	    0./)      !valve

  78. data%State%BopStackInput%MINORS1(8,1:4)=     (/2.,        0.9,     0.,	    0./)      !tee

  79. data%State%BopStackInput%MINORS1(9,1:4)=     (/2.,        0.,      105.,	    0./)      !valve

  80. data%State%BopStackInput%MINORS1(10,1:4)=    (/2.,        0.42,    0.,	    0./)      !elbow

  81. data%State%BopStackInput%MINORS1(11,1:4)=    (/2.,        0.42,    0.,	    0./)      !elbow

  82. data%State%BopStackInput%MINORS1(12,1:4)=    (/2.,        0.8,     0.,	    0./)      !unionA

  83. data%State%BopStackInput%MINORS1(13,1:4)=    (/2.,        0.8,     0.,	    0./)      !unionA

  84. data%State%BopStackInput%MINORS1(14,1:4)=    (/2.,        1.5,     0.,	    0./)      !elbow

  85. data%State%BopStackInput%MINORS1(15,1:4)=    (/2.,        0.,      425.,	    0./)      !valve

  86. data%State%BopStackInput%MINORS1(16,1:4)=    (/2.,        2.,      0.,	    0./)      !tee

  87. data%State%BopStackInput%MINORS1(17,1:4)=    (/0.75,      0.,      1.5,	    0./)      !REGULATOR

  88. data%State%BopStackInput%MINORS1(18,1:4)=    (/1.,        2.,      0.,	    0./)      !tee

  89. data%State%BopStackInput%MINORS1(19,1:4)=    (/1.,        1.5,     0.,	    0./)      !elbow

  90. data%State%BopStackInput%MINORS1(20,1:4)=    (/1.,        0.42,    0.,	    0./)      !elbow

  91. data%State%BopStackInput%MINORS1(21,1:4)=    (/1.,        0.42,    0.,	    0./)      !elbow

  92. data%State%BopStackInput%MINORS1(22,1:4)=    (/1.,        1.5,     0.,	    0./)      !elbow

  93. data%State%BopStackInput%MINORS1(23,1:4)=    (/1.,        0.,      105.,     0./)      !valve

  94. data%State%BopStackInput%MINORS1(24,1:4)=    (/1.,        0.9,     0.,	    0./)      !tee

  95. data%State%BopStackInput%MINORS1(25,1:4)=    (/1.,        0.,      0.,	    0.5/)     !FT

  96. data%State%BopStackInput%MINORS1(26,1:4)=    (/1.,        0.,      0.,       3.4/)     !filter

  97. data%State%BopStackInput%MINORS1(27,1:4)=    (/1.,        0.,      105.,	    0./)      !valve

  98. data%State%BopStackInput%MINORS1(28,1:4)=    (/1.,        0.9,     0.,	    0./)      !tee

  99. data%State%BopStackInput%MINORS1(29,1:4)=    (/1.,        1.5,     0.,	    0./)      !elbow

  100. data%State%BopStackInput%MINORS1(30,1:4)=    (/1.,        1.5,     0.,	    0./)      !elbow

  101. data%State%BopStackInput%MINORS1(31,1:4)=    (/1.,        0.,      9.2,	    0./)      !valve

  102. data%State%BopStackInput%MINORS1(32,1:4)=    (/1.,        0.8,     0.,	    0./)      !unionA

  103. data%State%BopStackInput%MINORS1(33,1:4)=    (/1.,        0.8,     0.,	    0./)      !unionA

  104. data%State%BopStackInput%MINORS1(34,1:4)=    (/0.75,      0.35,     0.,	    0./)      !contraction

  105. 

  106. 

  107. ALLOCATE (data%State%RamLine%MINORDIAMETER_RAMLINE(data%State%RamLine%NO_MINORSRAMLINE),data%State%RamLine%AREAMINOR_RAMLINE(data%State%RamLine%NO_MINORSRAMLINE),data%State%RamLine%LF_RAMLINE(data%State%RamLine%NO_MINORSRAMLINE),data%State%RamLine%CV_RAMLINE(data%State%RamLine%NO_MINORSRAMLINE)  &

  108.             ,data%State%RamLine%NOTE_RAMLINE(data%State%RamLine%NO_MINORSRAMLINE),data%State%RAMS%minlosspa(6,data%State%RamLine%NO_MINORSRAMLINE),data%State%RAMS%minloss(6,data%State%RamLine%NO_MINORSRAMLINE))

  109. 

  110. 

  111. 

  112. DO I=1,data%State%RamLine%NO_MINORSRAMLINE

  113.     data%State%RamLine%MINORDIAMETER_RAMLINE(I)=data%State%BopStackInput%MINORS1(I,1)

  114.     data%State%RamLine%LF_RAMLINE(I)=data%State%BopStackInput%MINORS1(I,2)

  115.     data%State%RamLine%CV_RAMLINE(I)=data%State%BopStackInput%MINORS1(I,3)

  116.     data%State%RamLine%NOTE_RAMLINE(I)=data%State%BopStackInput%MINORS1(I,4)

  117.     

  118.     

  119.     data%State%RamLine%AREAMINOR_RAMLINE(I)=PI*(data%State%RamLine%MINORDIAMETER_RAMLINE(I)*0.0254)**2/4.   !D(in), AREA(m)

  120. ENDDO

  121.     

  122. 

  123. !===========================================================================

  124. !                   	RAMLINE PIPNING LOSSES INPUT

  125. !===========================================================================

  126. data%State%RamLine%NO_PIPINGSRAMLINE=15

  127. 

  128. ALLOCATE (data%State%BopStackInput%PIPINGS_RAMLINE(data%State%RamLine%NO_PIPINGSRAMLINE,3))

  129. 

  130. !               ID(INCH)	L(MM)	ROUGHNESS(MM)=e  DESCRIPTION

  131. data%State%BopStackInput%PIPINGS_RAMLINE(1,1:3)=		 (/2.,		1035.,		0.03/)      !Avg.acc.distance

  132. data%State%BopStackInput%PIPINGS_RAMLINE(2,1:3)=		 (/2.,		730.,       0.03/)      !Acc.end.horizontal

  133. data%State%BopStackInput%PIPINGS_RAMLINE(3,1:3)=		 (/2.,		2000.,		0.03/)      !Acc.end.vertical

  134. data%State%BopStackInput%PIPINGS_RAMLINE(4,1:3)=      (/2.,		6000.,		0.05/)      !Hyd.hose

  135. data%State%BopStackInput%PIPINGS_RAMLINE(5,1:3)=      (/2.,		2370.,		0.03/)      !Corner.vertical1

  136. data%State%BopStackInput%PIPINGS_RAMLINE(6,1:3)=      (/2.,		210.,       0.03/)      !Add.from.bend

  137. data%State%BopStackInput%PIPINGS_RAMLINE(7,1:3)=      (/1.,		780.,       0.03/)      !Corner.horizontal1

  138. data%State%BopStackInput%PIPINGS_RAMLINE(8,1:3)=      (/1.,		780.,       0.03/)      !Corner.horizontal2

  139. data%State%BopStackInput%PIPINGS_RAMLINE(9,1:3)=      (/1.,		750.,       0.03/)      !Extra.length.back

  140. data%State%BopStackInput%PIPINGS_RAMLINE(10,1:3)=     (/1.,		800.,       0.03/)      !Corner.horizontal3

  141. data%State%BopStackInput%PIPINGS_RAMLINE(11,1:3)=     (/1.,		1650.,		0.03/)      !Corner.vertical2

  142. data%State%BopStackInput%PIPINGS_RAMLINE(12,1:3)=     (/1.,		340.,       0.03/)      !12.Valves.horizontal

  143. data%State%BopStackInput%PIPINGS_RAMLINE(13,1:3)=     (/1.,		1650.,		0.03/)      !Valves.vertical

  144. data%State%BopStackInput%PIPINGS_RAMLINE(14,1:3)=     (/1.,		31000.,		10./)       !Hyd.hose

  145. data%State%BopStackInput%PIPINGS_RAMLINE(15,1:3)=     (/1.,		526.,       0.03/)      !Add.from.bend

  146. 

  147. 

  148. 

  149.  ALLOCATE (data%State%RamLine%DIAM_RAMLINE_INCH(data%State%RamLine%NO_PIPINGSRAMLINE),  &

  150.             data%State%RamLine%AREA_RAMLINE(data%State%RamLine%NO_PIPINGSRAMLINE),data%State%RamLine%LENGT_RAMLINE(data%State%RamLine%NO_PIPINGSRAMLINE),data%State%RamLine%ROUGHNESS_RAMLINE(data%State%RamLine%NO_PIPINGSRAMLINE),data%State%RamLine%RELROUGH_RAMLINE(data%State%RamLine%NO_PIPINGSRAMLINE)   &

  151.             ,data%State%RAMS%Re_ramline(6,data%State%RamLine%NO_PIPINGSRAMLINE),data%State%RAMS%fric(6,data%State%RamLine%NO_PIPINGSRAMLINE),data%State%RAMS%fricloss(6,data%State%RamLine%NO_PIPINGSRAMLINE))

  152. 

  153. DO I=1,data%State%RamLine%NO_PIPINGSRAMLINE

  154.     data%State%RamLine%DIAM_RAMLINE_INCH(I)=data%State%BopStackInput%PIPINGS_RAMLINE(I,1)

  155.     data%State%RamLine%LENGT_RAMLINE(I)=data%State%BopStackInput%PIPINGS_RAMLINE(I,2)

  156.     data%State%RamLine%ROUGHNESS_RAMLINE(I)=data%State%BopStackInput%PIPINGS_RAMLINE(I,3)

  157.     

  158.     

  159.     

  160.     data%State%RamLine%AREA_RAMLINE(I)=PI*(data%State%RamLine%DIAM_RAMLINE_INCH(I)*0.0254)**2/4   !D(in), AREA(m)

  161.     data%State%RamLine%RELROUGH_RAMLINE(I)=data%State%RamLine%ROUGHNESS_RAMLINE(I)/(data%State%RamLine%DIAM_RAMLINE_INCH(I)*25.4)    !e/D

  162. !DIAM_RAMLINE_MM(I)=DIAM_RAMLINE_MM(I)*.001  !  (m)

  163.     data%State%RamLine%LENGT_RAMLINE(I)=data%State%RamLine%LENGT_RAMLINE(I)*.001      !  (m)

  164. ENDDO

  165.   

  166. 

  167. 

  168. !===========================================================================

  169. !                   	ANNULAR PREVENTER MINOR LOSSES INPUT

  170. !===========================================================================

  171. data%State%AnnularComputational%NO_MinorsAnnularLine=29

  172. 

  173. ALLOCATE (data%State%BopStackInput%MINORS_ANNULAR(data%State%AnnularComputational%NO_MinorsAnnularLine,4))

  174. 

  175. !			            ID(INCH)	LF	    CV	NOTE(BAR)     DESCRIPTION

  176. data%State%BopStackInput%MINORS_ANNULAR(1,1:4)=      (/2.,       2.,	    0.,	    0./)      !Acc.tee

  177. data%State%BopStackInput%MINORS_ANNULAR(2,1:4)=      (/2.,       0.9,	0.,	    0./)      !Avg.acc.tee

  178. data%State%BopStackInput%MINORS_ANNULAR(3,1:4)=      (/2.,       0.9,	0.,	    0./)      !Avg.acc.tee

  179. data%State%BopStackInput%MINORS_ANNULAR(4,1:4)=      (/2.,       0.9,	0.,	    0./)      !Avg.acc.tee

  180. data%State%BopStackInput%MINORS_ANNULAR(5,1:4)=      (/2.,       0.9,	0.,	    0./)      !tee

  181. data%State%BopStackInput%MINORS_ANNULAR(6,1:4)=      (/2.,       2.,	    0.,	    0./)      !tee

  182. data%State%BopStackInput%MINORS_ANNULAR(7,1:4)=      (/2.,       0.,	    105.,	0./)      !valve

  183. data%State%BopStackInput%MINORS_ANNULAR(8,1:4)=      (/2.,       0.9,	0.,	    0./)      !tee

  184. data%State%BopStackInput%MINORS_ANNULAR(9,1:4)=      (/2.,       0.,	    105.,	0./)      !valve

  185. data%State%BopStackInput%MINORS_ANNULAR(10,1:4)=     (/2.,       0.42,	0.,	    0./)      !elbow

  186. data%State%BopStackInput%MINORS_ANNULAR(11,1:4)=     (/2.,       0.42,	0.,	    0./)      !elbow

  187. data%State%BopStackInput%MINORS_ANNULAR(12,1:4)=     (/2.,       0.8,	0.,	    0./)      !unionA

  188. data%State%BopStackInput%MINORS_ANNULAR(13,1:4)=     (/2.,       0.8,	0.,	    0./)      !unionA

  189. data%State%BopStackInput%MINORS_ANNULAR(14,1:4)=     (/2.,       1.5,	0.,	    0./)      !elbow

  190. data%State%BopStackInput%MINORS_ANNULAR(15,1:4)=     (/2.,       0.,	    425.,	0./)      !valve

  191. data%State%BopStackInput%MINORS_ANNULAR(16,1:4)=     (/2.,       2.,	    0.,	    0./)      !tee

  192. data%State%BopStackInput%MINORS_ANNULAR(17,1:4)=     (/0.75,     0.,	    1.5,	0./)      !REGULATOR

  193. data%State%BopStackInput%MINORS_ANNULAR(18,1:4)=     (/1.,       2.,	    0.,	    0./)      !tee

  194. data%State%BopStackInput%MINORS_ANNULAR(19,1:4)=     (/1.,       1.5,	0.,	    0./)      !elbow

  195. data%State%BopStackInput%MINORS_ANNULAR(20,1:4)=     (/1.,       0.42,	0.,	    0./)      !elbow

  196. data%State%BopStackInput%MINORS_ANNULAR(21,1:4)=     (/1.,       0.42,	0.,	    0./)      !elbow

  197. data%State%BopStackInput%MINORS_ANNULAR(22,1:4)=     (/1.,       1.5,	0.,	    0./)      !elbow

  198. data%State%BopStackInput%MINORS_ANNULAR(23,1:4)=     (/1.,       0.,	    3.2,	0./)      !valve

  199. data%State%BopStackInput%MINORS_ANNULAR(24,1:4)=     (/1.,       2.,	    0.,	    0./)      !tee

  200. data%State%BopStackInput%MINORS_ANNULAR(25,1:4)=     (/1.,       1.5,	0.,	    0./)      !elbow

  201. data%State%BopStackInput%MINORS_ANNULAR(26,1:4)=     (/1.,       0.42,	0.,	    0./)      !elbow

  202. data%State%BopStackInput%MINORS_ANNULAR(27,1:4)=     (/1.,       0.42,	0.,	    0./)      !elbow

  203. data%State%BopStackInput%MINORS_ANNULAR(28,1:4)=     (/1.,       1.5,	0.,	    0./)      !elbow

  204. data%State%BopStackInput%MINORS_ANNULAR(29,1:4)=     (/1.,       0.,	    3.2,	0./)      !valve

  205. 

  206. 

  207. 

  208.  ALLOCATE (data%State%AnnularComputational%MINORDIAMETER_ANNULARLINE(data%State%AnnularComputational%NO_MinorsAnnularLine),data%State%AnnularComputational%AREAMINOR_ANNULARLINE(data%State%AnnularComputational%NO_MinorsAnnularLine),data%State%AnnularComputational%LF_ANNULARLINE(data%State%AnnularComputational%NO_MinorsAnnularLine)  &

  209.             ,data%State%AnnularComputational%CV_ANNULARLINE(data%State%AnnularComputational%NO_MinorsAnnularLine),data%State%AnnularComputational%NOTE_ANNULARLINE(data%State%AnnularComputational%NO_MinorsAnnularLine),data%State%AnnularComputational%minlosspa_ANNULAR(data%State%AnnularComputational%NO_MinorsAnnularLine),data%State%AnnularComputational%minloss_ANNULAR(data%State%AnnularComputational%NO_MinorsAnnularLine))

  210. 

  211. 

  212. 

  213. DO I=1,data%State%AnnularComputational%NO_MinorsAnnularLine

  214.     data%State%AnnularComputational%MINORDIAMETER_ANNULARLINE(I)=data%State%BopStackInput%MINORS_ANNULAR(I,1)

  215.     data%State%AnnularComputational%LF_ANNULARLINE(I)=data%State%BopStackInput%MINORS_ANNULAR(I,2)

  216.     data%State%AnnularComputational%CV_ANNULARLINE(I)=data%State%BopStackInput%MINORS_ANNULAR(I,3)

  217.     data%State%AnnularComputational%NOTE_ANNULARLINE(I)=data%State%BopStackInput%MINORS_ANNULAR(I,4)

  218.     

  219.     data%State%AnnularComputational%AREAMINOR_ANNULARLINE(I)=PI*(data%State%AnnularComputational%MINORDIAMETER_ANNULARLINE(I)*0.0254)**2/4.   !D(in), AREA(m)

  220. ENDDO

  221.    

  222. 

  223. !===========================================================================

  224. !                   ANNULAR PREVENTER PIPNING LOSSES INPUT

  225. !===========================================================================

  226. data%State%AnnularComputational%NO_PipingsAnnularLine=10

  227. 

  228. ALLOCATE (data%State%BopStackInput%PIPINGS_ANNULAR(data%State%AnnularComputational%NO_PipingsAnnularLine,3))

  229. 

  230. !            ID(INCH)   L(MM)	ROUGHNESS(MM)=e		 DESCRIPTION

  231. data%State%BopStackInput%PIPINGS_ANNULAR(1,1:3)=		 (/2.,		1035.,		0.03/)		     !Avg.acc.distance

  232. data%State%BopStackInput%PIPINGS_ANNULAR(2,1:3)=		 (/2.,		730.,		0.03/)		     !Acc.endhorizontal

  233. data%State%BopStackInput%PIPINGS_ANNULAR(3,1:3)=		 (/2.,		2000.,		0.03/)		     !Acc.endvertical

  234. data%State%BopStackInput%PIPINGS_ANNULAR(4,1:3)=		 (/2.,		6000.,		0.03/)		     !Hyd.hose

  235. data%State%BopStackInput%PIPINGS_ANNULAR(5,1:3)=		 (/2.,		2370.,		0.03/)		     !Corner.vertical1

  236. data%State%BopStackInput%PIPINGS_ANNULAR(6,1:3)=		 (/2.,		210.,		0.03/)		     !Add.frombend

  237. data%State%BopStackInput%PIPINGS_ANNULAR(7,1:3)=		 (/2.,		1000.,		0.03/)		     !manifold

  238. data%State%BopStackInput%PIPINGS_ANNULAR(8,1:3)=		 (/1.,		46000.,		0.03/)		     !pipe

  239. data%State%BopStackInput%PIPINGS_ANNULAR(9,1:3)=		 (/2.,		1000.,		0.03/)		     !manifold

  240. data%State%BopStackInput%PIPINGS_ANNULAR(10,1:3)=     (/1.,		46000.,		0.03/)           !pipe

  241. 

  242. 

  243. 

  244. 

  245.  ALLOCATE (data%State%AnnularComputational%DIAM_ANNULARLINE_INCH(data%State%AnnularComputational%NO_PipingsAnnularLine),data%State%AnnularComputational%AREA_ANNULARLINE(data%State%AnnularComputational%NO_PipingsAnnularLine),data%State%AnnularComputational%LENGT_ANNULARLINE(data%State%AnnularComputational%NO_PipingsAnnularLine)  &

  246.             ,data%State%AnnularComputational%ROUGHNESS_ANNULARLINE(data%State%AnnularComputational%NO_PipingsAnnularLine),data%State%AnnularComputational%RELROUGH_ANNULARLINE(data%State%AnnularComputational%NO_PipingsAnnularLine)   &

  247.             ,data%State%AnnularComputational%Re_ANNULARline(data%State%AnnularComputational%NO_PipingsAnnularLine),data%State%AnnularComputational%fricANNULAR(data%State%AnnularComputational%NO_PipingsAnnularLine),data%State%AnnularComputational%friclossANNULAR(data%State%AnnularComputational%NO_PipingsAnnularLine))

  248. 

  249. 

  250. DO I=1,data%State%AnnularComputational%NO_PipingsAnnularLine

  251.     data%State%AnnularComputational%DIAM_ANNULARLINE_INCH(I)=data%State%BopStackInput%PIPINGS_ANNULAR(I,1)

  252.     data%State%AnnularComputational%LENGT_ANNULARLINE(I)=data%State%BopStackInput%PIPINGS_ANNULAR(I,2)

  253.     data%State%AnnularComputational%ROUGHNESS_ANNULARLINE(I)=data%State%BopStackInput%PIPINGS_ANNULAR(I,3)

  254.     

  255.     

  256.     data%State%AnnularComputational%AREA_ANNULARLINE(I)=PI*(data%State%AnnularComputational%DIAM_ANNULARLINE_INCH(I)*0.0254)**2/4.   !D(in), AREA(m)

  257.     data%State%AnnularComputational%RELROUGH_ANNULARLINE(I)=data%State%AnnularComputational%ROUGHNESS_ANNULARLINE(I)/(data%State%AnnularComputational%DIAM_ANNULARLINE_INCH(I)*25.4)

  258. !DIAM_ANNULARLINE_MM(I)=DIAM_ANNULARLINE_MM(I)*.001  !  (m)

  259.     data%State%AnnularComputational%LENGT_ANNULARLINE(I)=data%State%AnnularComputational%LENGT_ANNULARLINE(I)*.001      !  (m)

  260.     ENDDO

  261. 

  262. 

  263. !===========================================================================

  264. !                   	AIR PUMP LOSSES INPUT

  265. !===========================================================================

  266. data%State%RamLine%NO_PIPINGS_AIRPLINE=1

  267. 

  268. ALLOCATE (data%State%BopStackInput%PIPINGS_AIRPUMP(data%State%RamLine%NO_PIPINGS_AIRPLINE,3))

  269. 

  270. !			ID(INCH)    L(MM)	  ROUGHNESS(MM)=e   DESCRIPTION

  271. data%State%BopStackInput%PIPINGS_AIRPUMP(1,1:3)=		 (/2.,      10000.,   0.03/)            !Avg.acc.distance

  272.  

  273. 

  274. 

  275.   

  276.   ALLOCATE (data%State%RamLine%DIAM_AIR_INCH(data%State%RamLine%NO_PIPINGS_AIRPLINE),data%State%RamLine%Re_air(data%State%RamLine%NO_PIPINGS_AIRPLINE),data%State%RamLine%AREA_AIR(data%State%RamLine%NO_PIPINGS_AIRPLINE),  &

  277.             data%State%RamLine%LENGT_AIR(data%State%RamLine%NO_PIPINGS_AIRPLINE),data%State%RamLine%ROUGHNESS_AIRPLINE(data%State%RamLine%NO_PIPINGS_AIRPLINE),data%State%RamLine%REL_ROUGHAIR(data%State%RamLine%NO_PIPINGS_AIRPLINE),  &

  278.            data%State%RamLine%fric_air(data%State%RamLine%NO_PIPINGS_AIRPLINE),data%State%RamLine%fricloss_air(data%State%RamLine%NO_PIPINGS_AIRPLINE))

  279. 

  280. 

  281. 

  282. DO I=1,data%State%RamLine%NO_PIPINGS_AIRPLINE

  283.     data%State%RamLine%DIAM_AIR_INCH(I)=data%State%BopStackInput%PIPINGS_AIRPUMP(I,1)

  284.     data%State%RamLine%LENGT_AIR(I)=data%State%BopStackInput%PIPINGS_AIRPUMP(I,2)

  285.     data%State%RamLine%ROUGHNESS_AIRPLINE(I)=data%State%BopStackInput%PIPINGS_AIRPUMP(I,3)

  286.     

  287.     

  288.     data%State%RamLine%AREA_AIR(I)=PI*(data%State%RamLine%DIAM_AIR_INCH(I)*0.0254)**2/4   !D(in), AREA(m)

  289.     data%State%RamLine%REL_ROUGHAIR(I)=data%State%RamLine%ROUGHNESS_AIRPLINE(I)/(data%State%RamLine%DIAM_AIR_INCH(I)*25.4)

  290. !DIAM_RAMLINE_MM(I)=DIAM_RAMLINE_MM(I)*.001  !  (m)

  291.     data%State%RamLine%LENGT_AIR(I)=data%State%RamLine%LENGT_AIR(I)*.001      !  (m)

  292.     ENDDO

  293. 

  294. 

  295.     

  296. !================================================================

  297. data%State%RamLine%NO_MINORS_AIRPLINE=6

  298. 

  299. ALLOCATE (data%State%BopStackInput%MINORS_AIRPUMP(data%State%RamLine%NO_MINORS_AIRPLINE,4))

  300. 

  301. !		ID(INCH)		LF	    CV	NOTE(BAR)      DESCRIPTION

  302. data%State%BopStackInput%MINORS_AIRPUMP(1,1:4)=      (/2.,			10.,	0.,	    0./)       !Acc.tee

  303. data%State%BopStackInput%MINORS_AIRPUMP(2,1:4)=      (/2.,			11.,	0.,	    0./)       !elbow

  304. data%State%BopStackInput%MINORS_AIRPUMP(3,1:4)=      (/1.,			0.,	    0.,	    3.4/)      !filter

  305. data%State%BopStackInput%MINORS_AIRPUMP(4,1:4)=      (/2.,			0.,	    105.,	0./)       !valve

  306. data%State%BopStackInput%MINORS_AIRPUMP(5,1:4)=      (/1.,			0.,	    9.2,	0./)       !valve

  307. data%State%BopStackInput%MINORS_AIRPUMP(6,1:4)=      (/2.,			6.4,	0.,	    0./)       !unionA

  308.   

  309. 

  310.     

  311.   ALLOCATE (data%State%RamLine%MINORDIAM_AIR_INCH(data%State%RamLine%NO_MINORS_AIRPLINE),data%State%RamLine%MINORAREA_AIR(data%State%RamLine%NO_MINORS_AIRPLINE),  &

  312.             data%State%RamLine%LF_AIR(data%State%RamLine%NO_MINORS_AIRPLINE),data%State%RamLine%CV_AIR(data%State%RamLine%NO_MINORS_AIRPLINE),data%State%RamLine%NOTE_AIR(data%State%RamLine%NO_MINORS_AIRPLINE) &

  313.             ,data%State%RamLine%minlosspa_air(data%State%RamLine%NO_MINORS_AIRPLINE),data%State%RamLine%minloss_air(data%State%RamLine%NO_MINORS_AIRPLINE))   

  314.  

  315. 

  316.  

  317. DO I=1,data%State%RamLine%NO_MINORS_AIRPLINE

  318.     data%State%RamLine%MINORDIAM_AIR_INCH(I)=data%State%BopStackInput%MINORS_AIRPUMP(I,1)

  319.     data%State%RamLine%LF_AIR(I)=data%State%BopStackInput%MINORS_AIRPUMP(I,2)

  320.     data%State%RamLine%CV_AIR(I)=data%State%BopStackInput%MINORS_AIRPUMP(I,3)

  321.     data%State%RamLine%NOTE_AIR(I)=data%State%BopStackInput%MINORS_AIRPUMP(I,4)

  322.     

  323.     

  324.     data%State%RamLine%MINORAREA_AIR(I)=PI*(data%State%RamLine%MINORDIAM_AIR_INCH(I)*0.0254)**2/4.   !D(in), AREA(m)

  325. ENDDO    

  326.     

  327. 

  328.  

  329.     END

  330.     

  331. SUBROUTINE pumps_charge_bottle()

  332. use SimulationVariables

  333. USE CAccumulatorVariables

  334. USE CBopStackVariables

  335. use CBopControlPanel

  336.     use SimulationVariables

  337. USE CEquipmentsConstants

  338. ! use CSimulationVariables

  339. implicit none

  340. 

  341. !Pannular_regset=min(AnnularRegulatorSetControl,1700.)  ! for changing its set conditions instantaneously

  342. !write(*,*) 'pumps_charge_bottle'

  343.     

  344.     if(data%Equipments%BopControlPanel%ByePassValve == -1.0) then

  345.         data%State%BopStackAcc%ByPassOld= -1.0

  346.     elseif(data%Equipments%BopControlPanel%ByePassValve == 1.0) then

  347.         data%State%BopStackAcc%ByPassOld= 1.0

  348.     endif

  349. 

  350.     

  351. !=====================================================================

  352. !                       ACCUMULATOR

  353. !=====for a 10 gal bottle,precharge=1000psig curve BOSCH-isotherm=====

  354. !for charging bottles by the pump

  355. !((((((((IN OUTER LOOP))))))

  356. ! ba1=1003; ba2=.03375; ba3=4.014; ba4=.2458;

  357. if (data%State%RamLine%AIRP_SWITCH==0) data%State%Pumps%DELTAV_AIR=0

  358. data%State%RamLine%FVR=data%State%RamLine%FVR+data%State%Pumps%DELTAV_AIR+data%State%Pumps%DELTAV_ELECP

  359. 

  360. data%State%RamLine%pacc_before=data%State%RamLine%P_ACC

  361. data%State%RamLine%P_ACC=data%State%RamLine%B1*exp(data%State%RamLine%B2*data%State%RamLine%FVR/data%State%BopStackAcc%NOBOTTLES)+data%State%RamLine%B3*exp(data%State%RamLine%B4*data%State%RamLine%FVR/data%State%BopStackAcc%NOBOTTLES) ! adiabatic(psig)<<<< 8=no. of bottles

  362. 

  363. !=====================================================================

  364. 

  365. 

  366.     if(data%State%BopStackAcc%ByPassOld == 1.0) then

  367. 

  368.         if (data%State%BopStackAcc%pram_reg<data%State%BopStackAcc%PRAMS_REGSET) then

  369. !pram_reg=pram_reg+(p_acc-pacc_before)

  370.             data%State%BopStackAcc%PressureDifference= (data%State%RamLine%P_ACC-data%State%RamLine%pacc_before)

  371.         else

  372. !pram_reg=PRAMS_REGSET

  373.             data%State%BopStackAcc%PressureDifference= data%State%BopStackAcc%PRAMS_REGSET - data%State%BopStackAcc%pram_reg

  374.         endif

  375.         

  376.         

  377. !old=pram_reg

  378. !difference= new-old

  379.          if (ABS(data%State%BopStackAcc%PressureDifference) > BaseDifferenceP) then

  380.              data%State%BopStackAcc%pram_reg= data%State%BopStackAcc%pram_reg + (data%State%BopStackAcc%PressureDifference/PressureDifferenceSteps)      ! PressureDifferenceSteps = 20.

  381.          else

  382.              

  383.                 if (data%State%BopStackAcc%pram_reg<data%State%BopStackAcc%PRAMS_REGSET) then

  384.                     data%State%BopStackAcc%pram_reg=data%State%BopStackAcc%pram_reg+(data%State%RamLine%P_ACC-data%State%RamLine%pacc_before)

  385.                 else

  386.                     data%State%BopStackAcc%pram_reg=data%State%BopStackAcc%PRAMS_REGSET

  387.                 endif

  388.         

  389.          endif        

  390.         

  391.         

  392.         

  393.         

  394. 

  395.     else        !if(ByPassOld == -1.0) then

  396. 

  397. 

  398. !old=pram_reg

  399. !new=p_acc- MAXVAL(data%State%RAM%loss_before)

  400. !difference= new-old

  401.          data%State%BopStackAcc%PressureDifference= data%State%RamLine%P_ACC- MAXVAL(data%State%RAM%loss_before) - data%State%BopStackAcc%pram_reg

  402.          if (ABS(data%State%BopStackAcc%PressureDifference) > BaseDifferenceP) then

  403.              data%State%BopStackAcc%pram_reg= data%State%BopStackAcc%pram_reg + (data%State%BopStackAcc%PressureDifference/PressureDifferenceSteps)      ! PressureDifferenceSteps = 20.

  404.          else

  405.             data%State%BopStackAcc%pram_reg= data%State%RamLine%P_ACC- MAXVAL(data%State%RAM%loss_before)

  406.          endif         

  407. 

  408.     endif

  409. 

  410. 

  411. !Pannular_reg= min(p_acc,Pannular_regset)

  412.     

  413. 

  414.     end

  415.     

  416.     

  417. SUBROUTINE airpump_code()    

  418. use SimulationVariables

  419. use CSounds

  420. implicit none

  421. INTEGER I

  422. 

  423. 

  424. 

  425. data%State%Pumps%QAIR=data%State%Pumps%Qiter+.1   !(gpm)  maximum flow for the start

  426. data%State%RamLine%diffp_air=-10

  427. data%State%RamLine%losses_air=10

  428. 

  429. !===================================================================

  430. !                       AIR OPERATED PUMP

  431. !=================for air consumption at 8 bar======================

  432. do while (data%State%RamLine%diffp_air<0)

  433.  data%State%Pumps%QAIR=data%State%Pumps%QAIR-.1

  434. ! Qup=QAIR_PUMP;

  435. ! bba1 =31.8; bba2 =-725.7 ; bba3 =4154;

  436.  data%State%Pumps%P_AIRP=data%State%RamLine%BBA1*data%State%Pumps%QAIR**2+data%State%RamLine%BBA2*data%State%Pumps%QAIR+data%State%RamLine%BBA3 !(psig)

  437.  data%State%RamLine%kinetic_air=sg*wdens*(data%State%Pumps%QAIR*6.30902e-5/((1/4.)*pi*(2*0.254e-1)**2))**2/(2*6895) !(psi)

  438. 

  439. 

  440.  data%State%RamLine%diffp_air=data%State%Pumps%P_AIRP+data%State%RamLine%kinetic_air-data%State%RamLine%P_ACC

  441. 

  442. end do    !returns Qup

  443. 

  444. 

  445.  do while (abs((data%State%RamLine%diffp_air-data%State%RamLine%losses_air)/data%State%RamLine%diffp_air)>data%State%Pumps%TOL_AIR) !finding correct QAIR_pump for 1 timecounter_ram

  446. 

  447. if (data%State%RamLine%diffp_air-data%State%RamLine%losses_air>0) then

  448.     data%State%Pumps%QAIR=data%State%Pumps%QAIR+.01

  449. else

  450.     data%State%Pumps%QAIR=data%State%Pumps%QAIR-.01

  451. endif

  452. 

  453. !===================================================================

  454. !                       AIR OPERATED PUMP

  455. !                   Maximator - Model: GX (35)

  456. !=================for air consumption at 8 bar======================

  457. 

  458.  data%State%Pumps%P_AIRP=data%State%RamLine%BBA1*data%State%Pumps%QAIR**2+data%State%RamLine%BBA2*data%State%Pumps%QAIR+data%State%RamLine%BBA3 !(psig)

  459.  data%State%RamLine%kinetic_air=sg*wdens*(data%State%Pumps%QAIR*6.30902e-005/((1/4.)*pi*(2*0.254e-1)**2))**2/(2*6895) !(psi)

  460. 

  461.  data%State%RamLine%diffp_air=data%State%Pumps%P_AIRP+data%State%RamLine%kinetic_air-data%State%RamLine%P_ACC

  462. 

  463. 

  464. !===========================LOSSES====================================

  465.  do i=1,data%State%RamLine%NO_PIPINGS_AIRPLINE

  466. data%State%RamLine%Re_air(i)=data%State%Pumps%QAIR*6.30902e-005*data%State%RamLine%DIAM_AIR_INCH(I)*0.0254/(data%State%RamLine%AREA_AIR(i)*nu)

  467. enddo

  468. 

  469. do i=1,data%State%RamLine%NO_PIPINGS_AIRPLINE

  470.      if (data%State%RamLine%Re_air(i)<Re_cr) then

  471.          data%State%RamLine%fric_air(i)=64/data%State%RamLine%Re_air(i)

  472.      else

  473.         data%State%RamLine%fric_air(i)=1/(-1.8*log10((data%State%RamLine%REL_ROUGHAIR(i)/3.7)**1.11+6.9/data%State%RamLine%Re_air(i)))**2

  474.      endif

  475.      

  476. data%State%RamLine%fricloss_air(i)=((data%State%RamLine%fric_air(i)*(wdens*sg*data%State%RamLine%LENGT_AIR(i)*(data%State%Pumps%QAIR*6.30902e-005/data%State%RamLine%AREA_AIR(i))**2))/(2*data%State%RamLine%DIAM_AIR_INCH(I)*0.0254))/6895

  477. 

  478.  enddo

  479. 

  480. 

  481. 

  482. 

  483. data%State%RamLine%pipe_loss1air=sum(data%State%RamLine%fricloss_air)  !loss before regulator(psi)

  484. 

  485. do i=1,data%State%RamLine%NO_MINORS_AIRPLINE

  486.     if (data%State%RamLine%LF_AIR(i)/=0) then

  487.         data%State%RamLine%minlosspa_air(i)=data%State%RamLine%LF_AIR(i)*wdens*sg*(data%State%Pumps%QAIR*6.30902e-005/data%State%RamLine%MINORAREA_AIR(i))**2/2  !(Pa)

  488.         data%State%RamLine%minloss_air(i)=data%State%RamLine%minlosspa_air(i)/6895 !(psi)

  489.     elseif (data%State%RamLine%CV_AIR(i)/=0) then

  490.         data%State%RamLine%minlosspa_air(i)=1000*sg*((11.7*data%State%Pumps%QAIR*6.30902e-005*3600)/(data%State%RamLine%CV_AIR(i)))**2 !(pa)

  491.         data%State%RamLine%minloss_air(i)=data%State%RamLine%minlosspa_air(i)/6895 !(psi)

  492.     else

  493.         data%State%RamLine%minlosspa_air(i)=data%State%RamLine%NOTE_AIR(i)*1e5 !(pa)

  494.         data%State%RamLine%minloss_air(i)=data%State%RamLine%minlosspa_air(i)/6895 !(psi)

  495.     endif

  496. enddo

  497. 

  498.         

  499. data%State%RamLine%minor_loss1air=sum(data%State%RamLine%minloss_air)  !loss before regulator(psi)

  500. 

  501. data%State%RamLine%static_loss1air=0

  502. 

  503. data%State%RamLine%losses_air=data%State%RamLine%pipe_loss1air+data%State%RamLine%minor_loss1air+data%State%RamLine%static_loss1air !(psi)

  504. !=========================================================================

  505. 

  506. 

  507. 

  508.  end do  !returns correct QAIR_PUMP for 1 counter

  509. 

  510. data%State%RamLine%counter_airp=data%State%RamLine%counter_airp+1 

  511.    data%State%Pumps%Qiter=data%State%Pumps%QAIR

  512.    

  513.    data%State%BopStackAcc%Cumulative_AirVolume= data%State%BopStackAcc%Cumulative_AirVolume + ((data%State%Pumps%QAIR * data%State%RamLine%DeltaT_BOP / 60.0 ) / (6.71*0.004329004) )        ! =strokes

  514. !write(*,*) 'Cumulative_AirVolume=' , Cumulative_AirVolume

  515.    

  516.    if ( data%State%BopStackAcc%Cumulative_AirVolume > 1.0 ) then

  517.        data%State%BopStackAcc%SoundKoomeyAirPump= 120.0

  518.        data%State%BopStackAcc%Cumulative_AirVolume= data%State%BopStackAcc%Cumulative_AirVolume - 1.0

  519.    else

  520.        data%State%BopStackAcc%SoundKoomeyAirPump= 0.0

  521.    endif

  522.    

  523.    

  524. !SoundKoomeyAirPump= int ( (QAIR_PUMP/ (6.71*0.004329004))  )   !  6.71: Volume per stroke (in^3)  , 0.004329004: in^3 to gal

  525. !  2.0 : two air pumps   60:seconds

  526. !write(*,*) 'sedaaaa=' ,(QAIR_PUMP/ (6.71*0.004329004)) *DeltaT_BOP/(60.0*2.0)

  527.    call SetSoundKoomeyAirPump(data%State%BopStackAcc%SoundKoomeyAirPump)

  528.    

  529. 

  530. data%State%Pumps%DELTAV_AIR=data%State%Pumps%QAIR*data%State%RamLine%DeltaT_BOP/60.  !(galon) delta_t=1sec , Q(gpm)

  531. 

  532. !write(*,*) 'QAIR_PUMP=' , QAIR_PUMP

  533. 

  534. 

  535. !((((((((IN OUTER LOOP))))))

  536. !fvr_air=fvr_air+deltav_air  !=fvr_after (galon) for next timecounter_ram

  537. 

  538. 

  539. !========================================================

  540. 

  541. !WRITE(50,50) QAIR_PUMP,p_airp, &

  542. !p_acc,diffp_air,losses_air,deltav_air

  543.     50     FORMAT(6(f15.5))

  544. 

  545.     

  546.     end

  547.     

  548.     

  549.     SUBROUTINE bop_code(ramtype,H_RAMBOP,RNUMBER) 

  550.     

  551.     use SimulationVariables

  552.     USE CAccumulatorVariables

  553.     use CBopControlPanel

  554.     use SimulationVariables

  555.     USE CEquipmentsConstants

  556. ! use CSimulationVariables

  557.     use SimulationVariables !@

  558.     use UnityModule

  559.    implicit none 

  560.     integer ramtype, I

  561.     real H_RAMBOP

  562.     INTEGER RNUMBER

  563.    

  564. !Pannular_regset=min(AnnularRegulatorSetControl,1700.)  ! for changing its set conditions instantaneously

  565.     

  566. !====================================================

  567. !       BOP back pressure IN PRESENCE OF DP

  568. !====================================================

  569. if (data%State%RAM(RNUMBER)%bop_type==2) then

  570.      data%State%RAM(RNUMBER)%checkp= data%State%RAM(RNUMBER)%p_bop

  571. if ( data%State%RAM(RNUMBER)%vdis_tot<data%State%ShearRAM%VA) then

  572.      data%State%RAM(RNUMBER)%p_bop=data%State%ShearRAM%PA

  573. elseif (data%State%ShearRAM%VA<= data%State%RAM(RNUMBER)%vdis_tot .and.  data%State%RAM(RNUMBER)%vdis_tot<=data%State%ShearRAM%VB) then

  574.      data%State%RAM(RNUMBER)%p_bop=data%State%ShearRAM%PA+data%State%ShearRAM%PB*(( data%State%RAM(RNUMBER)%vdis_tot-data%State%ShearRAM%VA)/(data%State%ShearRAM%VB-data%State%ShearRAM%VA))

  575. else

  576.      data%State%RAM(RNUMBER)%p_bop=data%State%ShearRAM%PA

  577. endif

  578.  

  579. if ( data%State%RAM(RNUMBER)%checkp/=data%State%ShearRAM%PA .and.  data%State%RAM(RNUMBER)%p_bop==data%State%ShearRAM%PA) then

  580.     data%State%RAM(RNUMBER)%Q=data%State%RAM(RNUMBER)%Qzero

  581. else

  582.  data%State%RAM(RNUMBER)%Q=data%State%RAM(RNUMBER)%flow

  583. endif

  584. 

  585.     endif    

  586.     

  587. !====================================================

  588. !       BOP back pressure without DP

  589. !====================================================

  590. if (data%State%RAM(RNUMBER)%bop_type==3) then 

  591.      data%State%RAM(RNUMBER)%p_bop=data%State%ShearRAM%PA

  592. ! Q=flow

  593.     endif 

  594.     

  595.       

  596.     

  597.  data%State%RAM(RNUMBER)%clock=0

  598. !======================Losses============================

  599.  data%State%RAM(RNUMBER)%loss_after=0 !initial value

  600.  data%State%RAM(RNUMBER)%diffp_ram=1000  !initial value

  601.  data%State%RAM(RNUMBER)%loss_before=0

  602. 

  603. !Q=0.0055;  %initial flow rate (m^3/s)

  604. 

  605.  data%State%RAM(RNUMBER)%Q=data%State%RAM(RNUMBER)%flow

  606. 

  607. 

  608. do while (abs( data%State%RAM(RNUMBER)%diffp_ram- data%State%RAM(RNUMBER)%loss_after)/ data%State%RAM(RNUMBER)%diffp_ram>data%State%RAM(RNUMBER)%tol)

  609.     if (data%State%RAM(RNUMBER)%Bottles_Charged_MalfActive==.true.) exit

  610. ! while abs( data%State%RAM(RNUMBER)%diffp_ram- data%State%RAM(RNUMBER)%loss_after)>10

  611.  data%State%RAM(RNUMBER)%clock= data%State%RAM(RNUMBER)%clock+1

  612. 

  613. if ( data%State%RAM(RNUMBER)%clock>20) then

  614. !    tclock=clock

  615.     data%State%RAM(RNUMBER)%tol= data%State%RAM(RNUMBER)%TOLZERO+(floor( data%State%RAM(RNUMBER)%clock/10)-1)*.001

  616. endif

  617. !      if (clock==1) continue

  618.        

  619.     if ( data%State%RAM(RNUMBER)%clock/=1 .and.  data%State%RAM(RNUMBER)%loss_after> data%State%RAM(RNUMBER)%diffp_ram) then

  620.         data%State%RAM(RNUMBER)%Q=data%State%RAM(RNUMBER)%Q-.1

  621.     elseif ( data%State%RAM(RNUMBER)%clock/=1 .and.  data%State%RAM(RNUMBER)%loss_after<= data%State%RAM(RNUMBER)%diffp_ram) then

  622.         data%State%RAM(RNUMBER)%Q=data%State%RAM(RNUMBER)%Q+.01

  623. 

  624.     endif

  625. 

  626. !====================Before Regulator=========================

  627. !pipe_loss1=3.2;  % should be replaced by formula

  628. do i=1,data%State%RamLine%NO_PIPINGSRAMLINE

  629.  data%State%RAMS%Re_ramline(RNUMBER,i)=data%State%RAM(RNUMBER)%Q*6.30902e-5*data%State%RamLine%DIAM_RAMLINE_INCH(i)*0.0254/(data%State%RamLine%AREA_RAMLINE(i)*nu)

  630. 

  631. enddo

  632. 

  633. 

  634. do i=1,data%State%RamLine%NO_PIPINGSRAMLINE

  635.      if ( data%State%RAMS%Re_ramline(RNUMBER,i)<Re_cr) then

  636.           data%State%RAMS%fric(RNUMBER,i)=64/ data%State%RAMS%Re_ramline(RNUMBER,i)

  637.      else

  638. 

  639.          data%State%RAMS%fric(RNUMBER,i)=1/(-1.8*log10((data%State%RamLine%RELROUGH_RAMLINE(i)/3.7)**1.11+6.9/ data%State%RAMS%Re_ramline(RNUMBER,i)))**2

  640.      endif

  641. 

  642.  data%State%RAMS%fricloss(RNUMBER,i)=(( data%State%RAMS%fric(RNUMBER,i)*(wdens*sg*data%State%RamLine%LENGT_RAMLINE(i)*(data%State%RAM(RNUMBER)%Q*6.30902e-5/data%State%RamLine%AREA_RAMLINE(i))**2))/(2*data%State%RamLine%DIAM_RAMLINE_INCH(i)*0.0254))/6895

  643. 

  644. enddo

  645. data%State%RAM(RNUMBER)%pipe_loss1=sum( data%State%RAMS%fricloss(RNUMBER,1:6))  !loss before regulator(psi)

  646. 

  647. !   Q*6.30902e-5   for (gpm) to (m^3/sec)

  648. do i=1,data%State%RamLine%NO_MINORSRAMLINE

  649.     if (data%State%RamLine%LF_RAMLINE(i)/=0) then

  650.          data%State%RAMS%minlosspa(RNUMBER,i)=data%State%RamLine%LF_RAMLINE(i)*wdens*sg*(data%State%RAM(RNUMBER)%Q*6.30902e-5/data%State%RamLine%AREAMINOR_RAMLINE(i))**2/2  !(Pa)

  651.          data%State%RAMS%minloss(RNUMBER,i)= data%State%RAMS%minlosspa(RNUMBER,i)/6895 !(psi)

  652.     elseif (data%State%RamLine%CV_RAMLINE(i)/=0) then

  653.          data%State%RAMS%minlosspa(RNUMBER,i)=1000*sg*((11.7*data%State%RAM(RNUMBER)%Q*6.30902e-5*3600)/(data%State%RamLine%CV_RAMLINE(i)))**2 !(pa)

  654.          data%State%RAMS%minloss(RNUMBER,i)= data%State%RAMS%minlosspa(RNUMBER,i)/6895 !(psi)

  655.     else

  656.          data%State%RAMS%minlosspa(RNUMBER,i)=data%State%RamLine%NOTE_RAMLINE(i)*1e5 !(pa)

  657.          data%State%RAMS%minloss(RNUMBER,i)= data%State%RAMS%minlosspa(RNUMBER,i)/6895 !(psi)

  658.     endif

  659. enddo

  660. 

  661.  data%State%RAMS%minloss(RNUMBER,17)= 5.*data%State%RAM(RNUMBER)%Q       ! <==== regulator loss (psig)

  662.         

  663.  data%State%RAM(RNUMBER)%minor_loss1=sum( data%State%RAMS%minloss(RNUMBER,1:16))  !loss before regulator(psi)

  664. 

  665.  data%State%RAM(RNUMBER)%static_loss1=0

  666. 

  667.  data%State%RAM(RNUMBER)%kinetic_loss1=sg*wdens*(data%State%RAM(RNUMBER)%Q*6.30902e-5/((1/4.)*pi*(.72*0.254e-1)**2))**2/(2*6895) !(psi)

  668. !.72 for regulator diameter)

  669. 

  670.  data%State%RAM(RNUMBER)%loss_before= data%State%RAM(RNUMBER)%pipe_loss1+ data%State%RAM(RNUMBER)%minor_loss1+ data%State%RAM(RNUMBER)%static_loss1+ data%State%RAM(RNUMBER)%kinetic_loss1 !(psi)

  671. !===================after Regulator======================

  672. ! all in (psi)

  673.  data%State%RAM(RNUMBER)%pipe_loss2=sum( data%State%RAMS%fricloss(RNUMBER,7:data%State%RamLine%NO_PIPINGSRAMLINE))  !7=after regulator<<<<<<<

  674.  data%State%RAM(RNUMBER)%minor_loss2=sum( data%State%RAMS%minloss(RNUMBER,17:data%State%RamLine%NO_MINORSRAMLINE))  !17=after regulator<<<<<<<

  675. 

  676.  data%State%RAM(RNUMBER)%deltah=H_RAMBOP-data%State%ShearRAM%H_REGRAM

  677. ! data%State%RAM(RNUMBER)%static_loss2pa= data%State%RAM(RNUMBER)%deltah*sg*wdens*gravity; %(pa)

  678.  data%State%RAM(RNUMBER)%static_loss2= data%State%RAM(RNUMBER)%deltah*sg*wdens*gravity/6895 !(psi)

  679.  data%State%RAM(RNUMBER)%kinetic_loss2=0

  680. 

  681.  data%State%RAM(RNUMBER)%loss_after= data%State%RAM(RNUMBER)%pipe_loss2+ data%State%RAM(RNUMBER)%minor_loss2+ data%State%RAM(RNUMBER)%static_loss2 !(psi)

  682. !=======================================================

  683. 

  684.     if(data%Equipments%BopControlPanel%ByePassValve == -1.0) then

  685.         data%State%BopStackAcc%ByPassOld= -1.0

  686.     elseif(data%Equipments%BopControlPanel%ByePassValve == 1.0) then

  687.         data%State%BopStackAcc%ByPassOld= 1.0

  688.     endif

  689. !   write(*,*) 'ByPassOld=', ByPassOld

  690.  

  691. 

  692. !==========================p_regulator======================

  693. if (data%State%BopStackAcc%PRAMS_REGSET<data%State%RamLine%P_ACC- data%State%RAM(RNUMBER)%loss_before) then

  694.     

  695.     if(data%State%BopStackAcc%ByPassOld == 1.0)  then  

  696. !old=pram_reg

  697. !new=PRAMS_REGSET- sum(data%State%RAMS%minloss(2:6,17))     ! 2:6 = RNUMBER

  698. !difference= new-old

  699.          data%State%BopStackAcc%PressureDifference= data%State%BopStackAcc%PRAMS_REGSET- sum(data%State%RAMS%minloss(2:6,17)) - data%State%BopStackAcc%pram_reg

  700.          if (ABS(data%State%BopStackAcc%PressureDifference) > BaseDifferenceP) then

  701.              data%State%BopStackAcc%pram_reg= data%State%BopStackAcc%pram_reg + (data%State%BopStackAcc%PressureDifference/PressureDifferenceSteps)      ! PressureDifferenceSteps = 20.

  702.          else

  703.             data%State%BopStackAcc%pram_reg= data%State%BopStackAcc%PRAMS_REGSET- sum(data%State%RAMS%minloss(2:6,17))     ! 2:6 = RNUMBER

  704.          endif

  705. !pram_reg=PRAMS_REGSET- sum(data%State%RAMS%minloss(2:6,17))     ! 2:6 = RNUMBER

  706.     else    !(ByPassOld == -1.0)

  707. !old=pram_reg

  708. !new=p_acc- data%State%RAM(RNUMBER)%loss_before

  709. !difference= new-old

  710.          data%State%BopStackAcc%PressureDifference= data%State%RamLine%P_ACC- data%State%RAM(RNUMBER)%loss_before - data%State%BopStackAcc%pram_reg

  711.          if (ABS(data%State%BopStackAcc%PressureDifference) > BaseDifferenceP) then

  712.              data%State%BopStackAcc%pram_reg= data%State%BopStackAcc%pram_reg + (data%State%BopStackAcc%PressureDifference/PressureDifferenceSteps)      ! PressureDifferenceSteps = 20.

  713.          else

  714.             data%State%BopStackAcc%pram_reg= data%State%RamLine%P_ACC- data%State%RAM(RNUMBER)%loss_before

  715.          endif    

  716.     

  717.     endif

  718. !pram_reg=p_acc- data%State%RAM(RNUMBER)%loss_before

  719. 

  720. else

  721.     data%State%BopStackAcc%pram_reg=data%State%RamLine%P_ACC- data%State%RAM(RNUMBER)%loss_before- sum(data%State%RAMS%minloss(2:6,17))     ! 2:6 = RNUMBER

  722. endif

  723. !===========================================================

  724. 

  725. 

  726. !====================Flow Rate calculation==================

  727.  data%State%RAM(RNUMBER)%diffp_ram=data%State%BopStackAcc%pram_reg+data%State%RAMS%minloss(RNUMBER,17)- data%State%RAM(RNUMBER)%p_bop

  728. !write(*,*) 'Q calculation'

  729. 

  730. 

  731. end do  ! returns correct Q for one timecounter_ram

  732. !write(*,*)'data%State%RAM(RNUMBER)%Q= ' , data%State%RAM(RNUMBER)%Q

  733. 

  734. !write(*,*) 'UpperRamsLeakMalf , Bottles_Charged_MalfActive=' , UpperRamsLeakMalf , Bottles_Charged_MalfActive

  735. if (data%State%RAM(RNUMBER)%Bottles_Charged_MalfActive==.true.) then

  736.     data%State%RAM(RNUMBER)%Q= 5.0 ! gpm

  737. !write(*,*) 'malffffffffff'

  738. endif

  739. 

  740. data%State%RAM(RNUMBER)%flow=data%State%RAM(RNUMBER)%Q

  741. 

  742. 

  743. 

  744. !====================p_accumulator for next iterate==================

  745. ! deltav=Q*1*1000;   %delta_t=1sec (lit)

  746. ! fvr=fvr-deltav;  %=fvr_after (lit)

  747. ! p_acc=0.001178*fvr^3+0.02019*fvr^2+16.1*fvr+2379; %(psia)

  748. ! vdis_bottles=fvr_tot-fvr;   %from the start***needed for  p_bop calculating (lit)

  749. 

  750. !========for a 10 gal bottle,precharge=1000psig curve BOSCH======

  751. data%State%RAM(RNUMBER)%deltav_bottles=data%State%RAM(RNUMBER)%Q*data%State%RamLine%DeltaT_BOP/60.   !(galon) delta_t=1sec , Q(gpm)

  752. data%State%RamLine%FVR=data%State%RamLine%FVR- data%State%RAM(RNUMBER)%deltav_bottles  !=fvr_after (galon)

  753. ! b1=993.7; b2=0.164; b3=5.492; b4=0.9796;

  754.  data%State%RAM(RNUMBER)%p_acccheck=data%State%RamLine%P_ACC

  755. data%State%RamLine%P_ACC=data%State%RamLine%B1*exp(data%State%RamLine%B2*data%State%RamLine%FVR/data%State%BopStackAcc%NOBOTTLES)+data%State%RamLine%B3*exp(data%State%RamLine%B4*data%State%RamLine%FVR/data%State%BopStackAcc%NOBOTTLES) !(psig)<<<< 8=no. of bottles

  756. !formula is derived for fvr (galon)

  757. !vdis_bottles=fvr_tot-fvr

  758.  data%State%RAM(RNUMBER)%vdis_bottles= data%State%RAM(RNUMBER)%vdis_bottles+ data%State%RAM(RNUMBER)%deltav_bottles

  759. 

  760.  data%State%RamLine%NoActiveRmas=6.-real(data%State%Annular%Annular_closed+data%State%RamLine%ShearBop_closed+data%State%PipeRam1%closed+data%State%PipeRam2%closed+data%State%ChokeLine%closed+data%State%KillLine%closed)

  761.  

  762.  data%State%RAM(RNUMBER)%vdis_elecp= data%State%RAM(RNUMBER)%vdis_elecp+(data%State%Pumps%DELTAV_ELECP/data%State%RamLine%NoActiveRmas)  !overal discharged

  763.   data%State%RAM(RNUMBER)%fvr_air= data%State%RAM(RNUMBER)%fvr_air+(data%State%Pumps%DELTAV_AIR/data%State%RamLine%NoActiveRmas)  !=fvr_after (galon) for next timecounter_ram

  764. !write(*,*) 'deltav_air , NoActiveRmas' , deltav_air , NoActiveRmas

  765. 

  766.  data%State%RAM(RNUMBER)%vdis_tot= data%State%RAM(RNUMBER)%vdis_bottles+ data%State%RAM(RNUMBER)%fvr_air+ data%State%RAM(RNUMBER)%vdis_elecp      !overal

  767. !  write(*,*)  data%State%RAM(RNUMBER)%vdis_tot, vdis_bottles, data%State%RAM(RNUMBER)%fvr_air, vdis_elecp      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

  768. 

  769. if (ramtype==1) then    !for shear ram

  770.      if (data%State%RamLine%ShearRamIsClosing) then

  771.         data%State%ShearRAM%IDshearBop=(2.*(data%State%ShearRAM%NeededVolumeShearRams- data%State%RAM(RNUMBER)%vdis_tot)*231./data%State%ShearRAM%AbopShearRam) +data%State%ShearRAM%ShearIsNotAllowed*data%State%ShearRAM%ODDrillpipe_inShearRam

  772. !write(*,*) 'stop now' !, IDshearBop, NeededVolumeShearRams, data%State%RAM(RNUMBER)%vdis_tot, AbopShearRam

  773.      endif

  774.      

  775.      if (data%State%RamLine%ShearRamIsOpening) then

  776.         data%State%ShearRAM%IDshearBop=data%State%ShearRAM%IDshearBopBase-2.*(data%State%ShearRAM%NeededVolumeShearRams- data%State%RAM(RNUMBER)%vdis_tot)*231./data%State%ShearRAM%AbopShearRam    

  777.      endif

  778.      

  779. 

  780.         

  781.       if ( data%State%RAM(RNUMBER)%vdis_tot>data%State%ShearRAM%NeededVolumeShearRams) then

  782. 

  783.              data%State%RamLine%ShearBop_closed=1

  784. !ShearBop_closed_withPossibility= ShearBop_closed * TD_BOPConnectionPossibility(3)

  785.              

  786.              data%State%BopStackAcc%pram_reg=data%State%BopStackAcc%pram_reg+data%State%RAMS%minloss(RNUMBER,17)

  787.              data%State%RamLine%P_ACC= data%State%RAM(RNUMBER)%p_acccheck

  788. 

  789.            if (data%State%RamLine%ShearRamIsClosing) then

  790.               data%State%ShearRAM%IDshearBop=0. + data%State%ShearRAM%ShearIsNotAllowed*data%State%ShearRAM%ODDrillpipe_inShearRam      

  791.               data%Equipments%BopControlPanel%MiddleRamsCloseLED = LedOn

  792.               data%State%BopStackInput%MiddleRamsCloseLEDMine = LedOn

  793.               data%Equipments%BopControlPanel%MiddleRamsOpenLED = LedOff

  794.               data%State%BopStackInput%MiddleRamsOpenLEDMine = LedOff

  795.               if (data%State%TD_BOP%BOPConnectionPossibility(3) == 1 .and. data%State%ShearRAM%ShearIsNotAllowed==0) then

  796.                   CALL CloseMiddleRams

  797.                   data%State%RamLine%ShearBop_Situation_forTD= 1       ! closed - for TD code

  798.               endif

  799.           endif

  800.      

  801.           if (data%State%RamLine%ShearRamIsOpening) then

  802.              data%State%ShearRAM%IDshearBop=data%State%ShearRAM%IDshearBopBase

  803.              data%Equipments%BopControlPanel%MiddleRamsOpenLED = LedOn

  804.              data%State%BopStackInput%MiddleRamsOpenLEDMine = LedOn

  805.              data%Equipments%BopControlPanel%MiddleRamsCloseLED = LedOff

  806.              data%State%BopStackInput%MiddleRamsCloseLEDMine = LedOff

  807.              CALL OpenMiddleRams

  808.              data%State%RamLine%ShearBop_Situation_forTD= 0       ! open - for TD code

  809.          endif

  810.     

  811.       endif  

  812.       

  813.       data%State%ShearRAM%IDshearBopFinal= data%State%ShearRAM%IDshearBop       ! for output data

  814.      

  815. endif

  816. 

  817. if (ramtype==2) then    !for pipe ram1

  818.      if (data%State%PipeRam1%IsClosing) then

  819.         data%State%PipeRam1%ID=(2.*(data%State%PipeRam1%NeededVolume- data%State%RAM(RNUMBER)%vdis_tot)*231./data%State%PipeRam1%A)+max(data%State%PipeRam1%ODDrillpipe_in,data%State%PipeRam1%ODDrillpipe_inBase)   

  820.      endif

  821.      

  822.      if (data%State%PipeRam1%IsOpening) then

  823.         data%State%PipeRam1%ID=data%State%PipeRam1%IDBase-2.*(data%State%PipeRam1%NeededVolume- data%State%RAM(RNUMBER)%vdis_tot)*231./data%State%PipeRam1%A    

  824.      endif

  825. 

  826.      if ( data%State%RAM(RNUMBER)%vdis_tot>data%State%PipeRam1%NeededVolume) then

  827. 

  828.               data%State%PipeRam1%closed=1

  829. !PipeRam1_Situation_forTD= PipeRam1_closed * TD_BOPConnectionPossibility(2)

  830.               data%State%BopStackAcc%pram_reg=data%State%BopStackAcc%pram_reg+data%State%RAMS%minloss(RNUMBER,17)

  831.               data%State%RamLine%P_ACC= data%State%RAM(RNUMBER)%p_acccheck

  832. 

  833.            if (data%State%PipeRam1%IsClosing) then

  834.               data%State%PipeRam1%ID=max(data%State%PipeRam1%ODDrillpipe_in,data%State%PipeRam1%ODDrillpipe_inBase)       

  835.               data%Equipments%BopControlPanel%UpperRamsCloseLED = LedOn

  836.               data%State%BopStackInput%UpperRamsCloseLEDMine = LedOn

  837.               data%Equipments%BopControlPanel%UpperRamsOpenLED = LedOff

  838.               data%State%BopStackInput%UpperRamsOpenLEDMine = LedOff

  839.               if (data%State%TD_BOP%BOPConnectionPossibility(2) == 1) then

  840.                   CALL CloseUpperRams       ! for C code

  841.                   call Set_BlowoutFromAnnular(.true.) 

  842. 

  843.  

  844.                   data%State%PipeRam1%Situation_forTD= 1       ! closed - for TD code

  845.               endif

  846.           endif

  847.      

  848.           if (data%State%PipeRam1%IsOpening) then

  849.               data%State%PipeRam1%ID=data%State%PipeRam1%IDBase

  850.               data%Equipments%BopControlPanel%UpperRamsOpenLED = LedOn

  851.               data%State%BopStackInput%UpperRamsOpenLEDMine = LedOn

  852.               data%Equipments%BopControlPanel%UpperRamsCloseLED = LedOff

  853.               data%State%BopStackInput%UpperRamsCloseLEDMine = LedOff

  854.               Call OpenUpperRams       ! for C code

  855.               data%State%PipeRam1%Situation_forTD= 0       ! open - for TD code

  856.          endif

  857.     

  858.      endif 

  859.      

  860.      data%State%ShearRAM%IDPipeRam1Final= data%State%PipeRam1%ID       ! for output data

  861. 

  862. endif

  863. 

  864. 

  865. if (ramtype==3) then    !for pipe ram2

  866.      if (data%State%PipeRam2%IsClosing) then

  867.         data%State%PipeRam2%ID=(2.*(data%State%PipeRam2%NeededVolume- data%State%RAM(RNUMBER)%vdis_tot)*231./data%State%PipeRam1%A)+max(data%State%PipeRam2%ODDrillpipe_in,data%State%PipeRam1%ODDrillpipe_inBase)   

  868.      endif

  869.      

  870.      if (data%State%PipeRam2%IsOpening) then

  871.         data%State%PipeRam2%ID=data%State%PipeRam1%IDBase-2.*(data%State%PipeRam2%NeededVolume- data%State%RAM(RNUMBER)%vdis_tot)*231./data%State%PipeRam1%A    

  872.      endif

  873. 

  874.      if ( data%State%RAM(RNUMBER)%vdis_tot>data%State%PipeRam2%NeededVolume) then

  875. 

  876.               data%State%PipeRam2%closed=1

  877. !PipeRam2_closed_withPossibility= PipeRam2_closed * TD_BOPConnectionPossibility(4)

  878.               data%State%BopStackAcc%pram_reg=data%State%BopStackAcc%pram_reg+data%State%RAMS%minloss(RNUMBER,17)

  879.               data%State%RamLine%P_ACC= data%State%RAM(RNUMBER)%p_acccheck

  880. 

  881.            if (data%State%PipeRam2%IsClosing) then

  882.               data%State%PipeRam2%ID=max(data%State%PipeRam2%ODDrillpipe_in,data%State%PipeRam1%ODDrillpipe_inBase)       

  883.               data%Equipments%BopControlPanel%LowerRamsCloseLED = LedOn

  884.               data%State%BopStackInput%LowerRamsCloseLEDMine = LedOn

  885.               data%Equipments%BopControlPanel%LowerRamsOpenLED = LedOff

  886.               data%State%BopStackInput%LowerRamsOpenLEDMine = LedOff

  887.               if (data%State%TD_BOP%BOPConnectionPossibility(4) == 1) then

  888.                   CALL CloseLowerRams

  889.                   data%State%PipeRam2%Situation_forTD= 1       ! closed - for TD code

  890.               endif

  891.           endif

  892.      

  893.           if (data%State%PipeRam2%IsOpening) then

  894.               data%State%PipeRam2%ID=data%State%PipeRam1%IDBase

  895.               data%Equipments%BopControlPanel%LowerRamsOpenLED = LedOn

  896.               data%State%BopStackInput%LowerRamsOpenLEDMine = LedOn

  897.               data%Equipments%BopControlPanel%LowerRamsCloseLED = LedOff

  898.               data%State%BopStackInput%LowerRamsCloseLEDMine = LedOff

  899.               CALL OpenLowerRams

  900.               data%State%PipeRam2%Situation_forTD= 0       ! open - for TD code

  901.          endif

  902.     

  903.      endif  

  904.      

  905.      data%State%ShearRAM%IDPipeRam2Final= data%State%PipeRam2%ID       ! for output data

  906. 

  907. endif

  908. 

  909. 

  910. if (ramtype==4) then    !for Choke Line

  911.      if (data%State%ChokeLine%IsClosing) then

  912.         data%State%ChokeLine%ID=(2.*(data%State%ChokeLine%NeededVolume- data%State%RAM(RNUMBER)%vdis_tot)*231./data%State%ChokeLine%Abop)+max(data%State%ChokeLine%ODDrillpipe_in,data%State%ChokeLine%ODDrillpipe_inBase)  

  913.      endif

  914.      

  915.      if (data%State%ChokeLine%IsOpening) then

  916.         data%State%ChokeLine%ID=data%State%ChokeLine%IDBase-2.*(data%State%ChokeLine%NeededVolume- data%State%RAM(RNUMBER)%vdis_tot)*231./data%State%ChokeLine%Abop    

  917.      endif

  918. 

  919.      if ( data%State%RAM(RNUMBER)%vdis_tot>data%State%ChokeLine%NeededVolume) then

  920. 

  921.               data%State%ChokeLine%closed=1

  922.               data%State%BopStackAcc%pram_reg=data%State%BopStackAcc%pram_reg+data%State%RAMS%minloss(RNUMBER,17)

  923.               data%State%RamLine%P_ACC= data%State%RAM(RNUMBER)%p_acccheck

  924. 

  925.            if (data%State%ChokeLine%IsClosing) then

  926.               data%State%ChokeLine%ID=max(data%State%ChokeLine%ODDrillpipe_in,data%State%ChokeLine%ODDrillpipe_inBase)      

  927.               data%Equipments%BopControlPanel%ChokeLineCloseLED = LedOn

  928.               data%State%BopStackInput%ChokeLineCloseLEDMine = LedOn

  929.               data%Equipments%BopControlPanel%ChokeLineOpenLED = LedOff

  930.               data%State%BopStackInput%ChokeLineOpenLEDMine = LedOff

  931.               CALL CloseChokeLine

  932.           endif

  933.      

  934.           if (data%State%ChokeLine%IsOpening) then

  935.               data%State%ChokeLine%ID=data%State%ChokeLine%IDBase

  936.               data%Equipments%BopControlPanel%ChokeLineOpenLED = LedOn

  937.               data%State%BopStackInput%ChokeLineOpenLEDMine = LedOn

  938.               data%Equipments%BopControlPanel%ChokeLineCloseLED = LedOff

  939.               data%State%BopStackInput%ChokeLineCloseLEDMine = LedOff

  940.               CALL OpenChokeLine

  941.          endif

  942.     

  943.      endif  

  944. 

  945. endif

  946. 

  947. 

  948. 

  949. if (ramtype==5) then    !for Kill Line

  950.      if (data%State%KillLine%IsClosing) then

  951.         data%State%KillLine%ID=(2.*(data%State%KillLine%NeededVolume- data%State%RAM(RNUMBER)%vdis_tot)*231./data%State%KillLine%Abop)+max(data%State%KillLine%ODDrillpipe_in,data%State%KillLine%ODDrillpipe_inBase)  

  952.      endif

  953.      

  954.      if (data%State%KillLine%IsOpening) then

  955.         data%State%KillLine%ID=data%State%KillLine%IDBase-2.*(data%State%KillLine%NeededVolume- data%State%RAM(RNUMBER)%vdis_tot)*231./data%State%KillLine%Abop    

  956.      endif

  957. 

  958.      if ( data%State%RAM(RNUMBER)%vdis_tot>data%State%KillLine%NeededVolume) then

  959. 

  960.               data%State%KillLine%closed=1

  961.               data%State%BopStackAcc%pram_reg=data%State%BopStackAcc%pram_reg+data%State%RAMS%minloss(RNUMBER,17)

  962.               data%State%RamLine%P_ACC= data%State%RAM(RNUMBER)%p_acccheck

  963. 

  964.            if (data%State%KillLine%IsClosing) then

  965.               data%State%KillLine%ID=max(data%State%KillLine%ODDrillpipe_in,data%State%KillLine%ODDrillpipe_inBase)      

  966.               data%Equipments%BopControlPanel%KillLineCloseLED = LedOn

  967.               data%State%BopStackInput%KillLineCloseLedMine = LedOn

  968.               data%Equipments%BopControlPanel%KillLineOpenLED = LedOff

  969.               data%State%BopStackInput%KillLineOpenLedMine = LedOff

  970.               CALL CloseKillLine

  971.           endif

  972.      

  973.           if (data%State%KillLine%IsOpening) then

  974.               data%State%KillLine%ID=data%State%KillLine%IDBase

  975.               data%Equipments%BopControlPanel%KillLineOpenLED = LedOn

  976.               data%State%BopStackInput%KillLineOpenLedMine = LedOn

  977.               data%Equipments%BopControlPanel%KillLineCloseLED = LedOff

  978.               data%State%BopStackInput%KillLineCloseLedMine = LedOff

  979.               CALL OpenKillLine

  980.          endif

  981.     

  982.      endif  

  983. 

  984. endif

  985. 

  986. 

  987. 

  988. 

  989. !if (ramtype==4) then    !for annular

  990. !     if (AnnularIsClosing) then

  991. !        IDAnnular=((NeededVolumeAnnular-vdis_tot)*231./AbopAnnular)+ODDrillpipe_inAnnular

  992. !     endif

  993. !

  994. !     if (AnnularIsOpening) then

  995. !        IDAnnular=IDAnnularBase-(NeededVolumeAnnular-vdis_tot)*231./AbopAnnular

  996. !     endif

  997. !

  998. !     if (vdis_tot>NeededVolumeAnnular) then

  999. !

  1000. !              Annular_closed=1

  1001. !              p_acc= data%State%RAM(RNUMBER)%p_acccheck

  1002. !

  1003. !           if (AnnularIsClosing) then

  1004. !              IDAnnular=ODDrillpipe_inAnnular

  1005. !              AnnularCloseLed = LedOn

  1006. !              AnnularOpenLed = LedOff

  1007. !          endif

  1008. !

  1009. !          if (AnnularIsOpening) then

  1010. !              IDAnnular=IDAnnularBase

  1011. !              AnnularOpenLed = LedOn

  1012. !              AnnularCloseLed = LedOff

  1013. !         endif

  1014. !

  1015. !     endif

  1016. !

  1017. !endif

  1018. 

  1019. 

  1020.     end

  1021.     

  1022.     

  1023.    

  1024.     SUBROUTINE bop_codeAnnular(RNUMBER) 

  1025.     

  1026.     use SimulationVariables 

  1027.     use CBopControlPanel

  1028.     use SimulationVariables

  1029.     USE CEquipmentsConstants

  1030.     use SimulationVariables !@

  1031. implicit none

  1032. INTEGER RNUMBER, I

  1033. 

  1034. 

  1035. !Pannular_regset=min(AnnularRegulatorSetControl,1700.)  ! for changing its set conditions instantaneously

  1036.     

  1037. !====================================================

  1038. !       BOP back pressure without DP

  1039. !====================================================

  1040. !if (bop_type==3) then

  1041. !p_annular=510.725-(30.145*IDAnnular)

  1042.     data%State%AnnularComputational%p_annular=448-(19.7*data%State%Annular%IDAnnular)

  1043. !    Q=flow

  1044. !    endif

  1045.     

  1046.       

  1047.     

  1048.  data%State%RAM(RNUMBER)%clock=0

  1049. !======================Losses============================

  1050.  data%State%RAM(RNUMBER)%loss_after=0 !initial value

  1051.  data%State%RAM(RNUMBER)%diffp_ram=1000  !initial value

  1052.  data%State%RAM(RNUMBER)%loss_before=0

  1053. 

  1054. !Q=0.0055;  %initial flow rate (m^3/s)

  1055. 

  1056.  data%State%RAM(RNUMBER)%Q=data%State%RAM(RNUMBER)%flow

  1057. 

  1058. !write(*,*) 'Q1=' , Q

  1059. !write(*,*) 'tol=' , tol

  1060. do while (abs( data%State%RAM(RNUMBER)%diffp_ram- data%State%RAM(RNUMBER)%loss_after)/ data%State%RAM(RNUMBER)%diffp_ram>data%State%Annular%tolAnnular)

  1061.     if (data%State%RAM(RNUMBER)%Bottles_Charged_MalfActive==.true.) exit

  1062. ! while abs( data%State%RAM(RNUMBER)%diffp_ram- data%State%RAM(RNUMBER)%loss_after)>10

  1063.  data%State%RAM(RNUMBER)%clock= data%State%RAM(RNUMBER)%clock+1

  1064. 

  1065. if ( data%State%RAM(RNUMBER)%clock>20) then

  1066. !    tclock=clock

  1067.     data%State%Annular%tolAnnular=data%State%Annular%tolzeroAnnular+(floor( data%State%RAM(RNUMBER)%clock/10)-1)*.001

  1068. endif

  1069. !      if (clock==1) continue

  1070.        

  1071.     if ( data%State%RAM(RNUMBER)%clock/=1 .and.  data%State%RAM(RNUMBER)%loss_after> data%State%RAM(RNUMBER)%diffp_ram) then

  1072.         data%State%RAM(RNUMBER)%Q=data%State%RAM(RNUMBER)%Q-.1

  1073.     elseif ( data%State%RAM(RNUMBER)%clock/=1 .and.  data%State%RAM(RNUMBER)%loss_after<= data%State%RAM(RNUMBER)%diffp_ram) then

  1074.         data%State%RAM(RNUMBER)%Q=data%State%RAM(RNUMBER)%Q+.01

  1075. 

  1076.     endif

  1077. 

  1078. !====================Before Regulator=========================

  1079. do i=1,data%State%AnnularComputational%NO_PipingsAnnularLine

  1080. data%State%AnnularComputational%Re_ANNULARline(i)=data%State%RAM(RNUMBER)%Q*6.30902e-5*data%State%AnnularComputational%DIAM_ANNULARLINE_INCH(i)*0.0254/(data%State%AnnularComputational%AREA_ANNULARLINE(i)*nu)

  1081. 

  1082. enddo

  1083. 

  1084. 

  1085. 

  1086. do i=1,data%State%AnnularComputational%NO_PipingsAnnularLine

  1087.      if (data%State%AnnularComputational%Re_ANNULARline(i)<Re_cr) then

  1088.          data%State%AnnularComputational%fricANNULAR(i)=64/data%State%AnnularComputational%Re_ANNULARline(i)

  1089.      else

  1090. 

  1091.         data%State%AnnularComputational%fricANNULAR(i)=1/(-1.8*log10((data%State%AnnularComputational%RELROUGH_ANNULARLINE(i)/3.7)**1.11+6.9/data%State%AnnularComputational%Re_ANNULARline(i)))**2

  1092.      endif

  1093. 

  1094. data%State%AnnularComputational%friclossANNULAR(i)=((data%State%AnnularComputational%fricANNULAR(i)*(wdens*sg*data%State%AnnularComputational%LENGT_ANNULARLINE(i)*(data%State%RAM(RNUMBER)%Q*6.30902e-5/data%State%AnnularComputational%AREA_ANNULARLINE(i))**2))/(2*data%State%AnnularComputational%DIAM_ANNULARLINE_INCH(i)*0.0254))/6895

  1095. 

  1096. enddo

  1097. 

  1098. 

  1099. 

  1100. 

  1101.  data%State%RAM(RNUMBER)%pipe_loss1=sum(data%State%AnnularComputational%friclossANNULAR(1:6))  !loss before regulator(psi)

  1102. 

  1103. 

  1104. do i=1,data%State%AnnularComputational%NO_MinorsAnnularLine

  1105.     if (data%State%AnnularComputational%LF_ANNULARLINE(i)/=0) then

  1106.         data%State%AnnularComputational%minlosspa_ANNULAR(i)=data%State%AnnularComputational%LF_ANNULARLINE(i)*wdens*sg*(data%State%RAM(RNUMBER)%Q*6.30902e-5/data%State%AnnularComputational%AREAMINOR_ANNULARLINE(i))**2/2  !(Pa)

  1107.         data%State%AnnularComputational%minloss_ANNULAR(i)=data%State%AnnularComputational%minlosspa_ANNULAR(i)/6895 !(psi)

  1108.     elseif (data%State%AnnularComputational%CV_ANNULARLINE(i)/=0) then

  1109.         data%State%AnnularComputational%minlosspa_ANNULAR(i)=1000*sg*((11.7*data%State%RAM(RNUMBER)%Q*6.30902e-5*3600)/(data%State%AnnularComputational%CV_ANNULARLINE(i)))**2 !(pa)

  1110.         data%State%AnnularComputational%minloss_ANNULAR(i)=data%State%AnnularComputational%minlosspa_ANNULAR(i)/6895 !(psi)

  1111.     else

  1112.         data%State%AnnularComputational%minlosspa_ANNULAR(i)=data%State%AnnularComputational%NOTE_ANNULARLINE(i)*1e5 !(pa)

  1113.         data%State%AnnularComputational%minloss_ANNULAR(i)=data%State%AnnularComputational%minlosspa_ANNULAR(i)/6895 !(psi)

  1114.     endif

  1115. enddo

  1116. 

  1117. data%State%AnnularComputational%minloss_ANNULAR(17)= 3.4*data%State%RAM(RNUMBER)%Q       ! <==== regulator loss (psig)

  1118.         

  1119.  data%State%RAM(RNUMBER)%minor_loss1=sum(data%State%AnnularComputational%minloss_ANNULAR(1:16))  !loss before regulator(psi)

  1120. 

  1121.  data%State%RAM(RNUMBER)%static_loss1=0

  1122. 

  1123.  data%State%RAM(RNUMBER)%kinetic_loss1=sg*wdens*(data%State%RAM(RNUMBER)%Q*6.30902e-5/((1/4.)*pi*(.72*0.254e-1)**2))**2/(2*6895) !(psi)

  1124. !.72 for regulator diameter)

  1125. 

  1126.  data%State%RAM(RNUMBER)%loss_before= data%State%RAM(RNUMBER)%pipe_loss1+ data%State%RAM(RNUMBER)%minor_loss1+ data%State%RAM(RNUMBER)%static_loss1+ data%State%RAM(RNUMBER)%kinetic_loss1 !(psi)

  1127. !===================after Regulator======================

  1128. ! all in (psi)

  1129.  data%State%RAM(RNUMBER)%pipe_loss2=sum(data%State%AnnularComputational%friclossANNULAR(7:data%State%AnnularComputational%NO_PipingsAnnularLine))  !7=after regulator<<<<<<<

  1130.  data%State%RAM(RNUMBER)%minor_loss2=sum(data%State%AnnularComputational%minloss_ANNULAR(17:data%State%AnnularComputational%NO_MinorsAnnularLine))  !17=after regulator<<<<<<<

  1131. 

  1132.  data%State%RAM(RNUMBER)%deltah=data%State%Annular%H_AnnularBop-data%State%ShearRAM%H_REGRAM

  1133. ! data%State%RAM(RNUMBER)%static_loss2pa= data%State%RAM(RNUMBER)%deltah*sg*wdens*gravity; %(pa)

  1134.  data%State%RAM(RNUMBER)%static_loss2= data%State%RAM(RNUMBER)%deltah*sg*wdens*gravity/6895 !(psi)

  1135.  data%State%RAM(RNUMBER)%kinetic_loss2=0

  1136. 

  1137.  data%State%RAM(RNUMBER)%loss_after= data%State%RAM(RNUMBER)%pipe_loss2+ data%State%RAM(RNUMBER)%minor_loss2+ data%State%RAM(RNUMBER)%static_loss2 !(psi)

  1138. !=======================================================

  1139. 

  1140. 

  1141. 

  1142. !==========================p_regulator======================

  1143. if (data%State%Annular%Pannular_regset<data%State%RamLine%P_ACC- data%State%RAM(RNUMBER)%loss_before) then

  1144.     

  1145. !!new= Pannular_regset-minloss_ANNULAR(17)

  1146. !

  1147. !CALL Pannular_regDelay%AddToFirst( Pannular_regset-minloss_ANNULAR(17) )

  1148. !CALL Pannular_regDelay%Remove(Pannular_regDelay%Length())

  1149. !Pannular_reg= Pannular_regDelay%Last()

  1150.             

  1151.             data%State%Annular%Pannular_reg= data%State%Annular%Pannular_regset-data%State%AnnularComputational%minloss_ANNULAR(17)     

  1152. 

  1153. ! data%State%RAM(RNUMBER)%loss_after= data%State%RAM(RNUMBER)%loss_after-minloss(17);   %minloss(17)=regulator loss

  1154. else

  1155.     

  1156. !!new= p_acc- data%State%RAM(RNUMBER)%loss_before-minloss_ANNULAR(17)

  1157. !

  1158. !CALL Pannular_regDelay%AddToFirst( p_acc- data%State%RAM(RNUMBER)%loss_before-minloss_ANNULAR(17) )

  1159. !CALL Pannular_regDelay%Remove(Pannular_regDelay%Length())

  1160. !Pannular_reg= Pannular_regDelay%Last()

  1161.             

  1162.             data%State%Annular%Pannular_reg= data%State%RamLine%P_ACC- data%State%RAM(RNUMBER)%loss_before-data%State%AnnularComputational%minloss_ANNULAR(17)    

  1163. 

  1164. endif

  1165. !===========================================================

  1166. 

  1167. 

  1168. !====================Flow Rate calculation==================

  1169.  data%State%RAM(RNUMBER)%diffp_ram=data%State%Annular%Pannular_reg+data%State%AnnularComputational%minloss_ANNULAR(17)-data%State%AnnularComputational%p_annular

  1170. 

  1171. 

  1172. 

  1173. 

  1174. 

  1175. end do  ! returns correct Q for one timecounter_ram

  1176. 

  1177. if (data%State%RAM(RNUMBER)%Bottles_Charged_MalfActive==.true.) then

  1178.     data%State%RAM(RNUMBER)%Q= 5.0 ! gpm

  1179. !write(*,*) 'malffffffffff'

  1180. endif

  1181. 

  1182. data%State%RAM(RNUMBER)%flow=data%State%RAM(RNUMBER)%Q

  1183. !write(*,*) 'minloss_ANNULAR(17)=' , minloss_ANNULAR(17)

  1184. !write(*,*) '%loss_after=' ,data%State%RAM(RNUMBER)%loss_after

  1185. !====================p_accumulator for next iterate==================

  1186. ! deltav=Q*1*1000;   %delta_t=1sec (lit)

  1187. ! fvr=fvr-deltav;  %=fvr_after (lit)

  1188. ! p_acc=0.001178*fvr^3+0.02019*fvr^2+16.1*fvr+2379; %(psia)

  1189. ! vdis_bottles=fvr_tot-fvr;   %from the start***needed for p_BOP calculating (lit)

  1190. 

  1191. !========for a 10 gal bottle,precharge=1000psig curve BOSCH======

  1192.  data%State%RAM(RNUMBER)%deltav_bottles=data%State%RAM(RNUMBER)%Q*data%State%RamLine%DeltaT_BOP/60.   !(galon) delta_t=1sec , Q(gpm)

  1193. data%State%RamLine%FVR=data%State%RamLine%FVR- data%State%RAM(RNUMBER)%deltav_bottles  !=fvr_after (galon)

  1194. ! b1=993.7; b2=0.164; b3=5.492; b4=0.9796;

  1195.  data%State%RAM(RNUMBER)%p_acccheck=data%State%RamLine%P_ACC

  1196. data%State%RamLine%P_ACC=data%State%RamLine%B1*exp(data%State%RamLine%B2*data%State%RamLine%FVR/data%State%BopStackAcc%NOBOTTLES)+data%State%RamLine%B3*exp(data%State%RamLine%B4*data%State%RamLine%FVR/data%State%BopStackAcc%NOBOTTLES) !(psig)<<<< 8=no. of bottles

  1197. !formula is derived for fvr (galon)

  1198. !vdis_bottles=fvr_tot-fvr

  1199.  data%State%RAM(RNUMBER)%vdis_bottles= data%State%RAM(RNUMBER)%vdis_bottles+ data%State%RAM(RNUMBER)%deltav_bottles

  1200. 

  1201.  data%State%RamLine%NoActiveRmas=6.-real(data%State%Annular%Annular_closed+data%State%RamLine%ShearBop_closed+data%State%PipeRam1%closed+data%State%PipeRam2%closed+data%State%ChokeLine%closed+data%State%KillLine%closed)

  1202.  data%State%RAM(RNUMBER)%vdis_elecp= data%State%RAM(RNUMBER)%vdis_elecp+(data%State%Pumps%DELTAV_ELECP/data%State%RamLine%NoActiveRmas)  !overal discharged

  1203.   data%State%RAM(RNUMBER)%fvr_air= data%State%RAM(RNUMBER)%fvr_air+(data%State%Pumps%DELTAV_AIR/data%State%RamLine%NoActiveRmas)  !=fvr_after (galon) for next timecounter_ram

  1204. 

  1205. 

  1206.  data%State%RAM(RNUMBER)%vdis_tot= data%State%RAM(RNUMBER)%vdis_bottles+ data%State%RAM(RNUMBER)%fvr_air+ data%State%RAM(RNUMBER)%vdis_elecp      !overal

  1207. !write(*,*)  'data%State%RAM(RNUMBER)%vdis_tot=' , data%State%RAM(RNUMBER)%vdis_tot!, vdis_bottles,fvr_air, vdis_elecp      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<

  1208. 

  1209. !write(*,*) 'a)IDAnnular=' , IDAnnular

  1210. !    WRITE(*,*) 'a)NeededVolumeAnnular=' , NeededVolumeAnnular

  1211. !

  1212. !    write(*,*) 'AbopAnnular= ' , AbopAnnular

  1213. !    write(*,*) 'max=' , max(ODDrillpipe_inAnnular,ODDrillpipe_inAnnularBase)

  1214. 

  1215. !if (ramtype==4) then    !for annular

  1216.      if (data%State%Annular%AnnularIsClosing) then

  1217.         data%State%Annular%IDAnnular=(2.*(data%State%Annular%NeededVolumeAnnular- data%State%RAM(RNUMBER)%vdis_tot)*231./data%State%Annular%AbopAnnular)+max(data%State%Annular%ODDrillpipe_inAnnular,data%State%Annular%ODDrillpipe_inAnnularBase)  

  1218. !write(*,*) 'b)IDAnnular=' , IDAnnular

  1219.      endif

  1220.      

  1221.      if (data%State%Annular%AnnularIsOpening) then

  1222.         data%State%Annular%IDAnnular=data%State%Annular%IDAnnularBase-2.*(data%State%Annular%NeededVolumeAnnular- data%State%RAM(RNUMBER)%vdis_tot)*231./data%State%Annular%AbopAnnular

  1223.      endif

  1224. 

  1225.      if ( data%State%RAM(RNUMBER)%vdis_tot>data%State%Annular%NeededVolumeAnnular) then

  1226. 

  1227.               data%State%Annular%Annular_closed=1

  1228. !Annular_closed_withPossibility= Annular_closed * TD_BOPConnectionPossibility(1)

  1229. !pram_reg=pram_reg+minloss_ANNULAR(17)

  1230.             data%State%Annular%Pannular_reg= data%State%Annular%Pannular_reg+data%State%AnnularComputational%minloss_ANNULAR(17)               

  1231. 

  1232.               data%State%RamLine%P_ACC= data%State%RAM(RNUMBER)%p_acccheck

  1233.            if (data%State%Annular%AnnularIsClosing) then

  1234.                

  1235.               data%State%Annular%IDAnnular=max(data%State%Annular%ODDrillpipe_inAnnular,data%State%Annular%ODDrillpipe_inAnnularBase)      

  1236.               data%Equipments%BopControlPanel%AnnularCloseLED = LedOn

  1237.               data%State%BopStackInput%AnnularCloseLedMine = LedOn

  1238.               

  1239.               data%Equipments%BopControlPanel%AnnularOpenLED = LedOff

  1240.               data%State%BopStackInput%AnnularOpenLedMine = LedOff

  1241.               

  1242.               if (data%State%TD_BOP%BOPConnectionPossibility(1) == 1) then

  1243.                   CALL CloseAnnular

  1244.                   data%State%Annular%Annular_Situation_forTD= 1       ! closed - for TD code

  1245.               endif

  1246.            endif

  1247.            

  1248.            

  1249.      

  1250.           if (data%State%Annular%AnnularIsOpening) then

  1251.               

  1252.               data%State%Annular%IDAnnular=data%State%Annular%IDAnnularBase

  1253.               data%Equipments%BopControlPanel%AnnularOpenLED = LedOn

  1254.               data%State%BopStackInput%AnnularOpenLedMine = LedOn

  1255.               

  1256.               data%Equipments%BopControlPanel%AnnularCloseLED = LedOff

  1257.               data%State%BopStackInput%AnnularCloseLedMine = LedOff

  1258.               

  1259.               CALL OpenAnnular

  1260.               data%State%Annular%Annular_Situation_forTD= 0       ! open - for TD code

  1261.          endif

  1262.     

  1263.      endif 

  1264.      

  1265.         data%State%ShearRAM%IDAnnularFinal= data%State%Annular%IDAnnular       ! for output data

  1266. !endif

  1267. 

  1268. 

  1269.     end