1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv! 2 ! ! 3 ! Subroutine: PHYSICAL_PROP_ROg ! 4 ! Purpose: User hook for calculating the gas phase density. ! 5 ! ! 6 ! ! 7 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^! 8 SUBROUTINE USR_PHYSICAL_PROP_ROg 9 10 ! Fluid grid loop bounds. 11 use compar, only: IJKStart3, IJKEnd3 12 ! Function to identify wall cells 13 use functions, only: WALL_AT 14 15 use error_manager 16 17 implicit none 18 19 ! Local Variables: 20 !---------------------------------------------------------------------// 21 ! Loop indicies 22 INTEGER :: IJK ! Computational cell 23 !......................................................................! 24 25 26 ! The following error message is used to make sure that if a user 27 ! defined gas density is invoked, that this routine has been modified. 28 29 !- REMOVE THE FOLLOWING ---------------------------------------------->> 30 31 CALL INIT_ERR_MSG('USR_PHYSICAL_PROP_ROg') 32 WRITE(ERR_MSG,9999) 33 CALL FLUSH_ERR_MSG(ABORT=.TRUE.) 34 35 9999 FORMAT('ERROR 9999: The user-defined drag routine was invoked ', & 36 'but this',/'generic error message exits. Either choose a ', & 37 'different drag law',/'or correct mfix/model/usr_drag.f') 38 39 !- END REMOVE --------------------------------------------------------<< 40 41 42 IJK_LP: DO IJK = IJKSTART3, IJKEND3 43 IF(WALL_AT(IJK)) cycle IJK_LP 44 45 ! Calculate the fluid density and bulk density 46 ! RO_G(IJK) = 47 ! ROP_G(IJK) = 48 49 ENDDO IJK_LP 50 51 RETURN 52 END SUBROUTINE USR_PHYSICAL_PROP_ROg 53 54 55 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv! 56 ! ! 57 ! Subroutine: PHYSICAL_PROP_ROs ! 58 ! Purpose: User hook for calculating solids phase density. ! 59 ! ! 60 ! Author: J. Musser Date: 28-JUN-13 ! 61 ! Reviewer: Date: ! 62 ! ! 63 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^! 64 SUBROUTINE USR_PHYSICAL_PROP_ROs 65 66 ! Fluid grid loop bounds. 67 use compar, only: IJKStart3, IJKEnd3 68 ! Function to identify wall cells 69 use functions, only: WALL_AT 70 ! Number of solids phases 71 use physprop, only: MMAX 72 73 use error_manager 74 75 76 implicit none 77 78 79 ! Local Variables: 80 !---------------------------------------------------------------------// 81 ! Loop indicies 82 INTEGER :: M, IJK 83 !......................................................................! 84 85 ! The following error message is used to make sure that if a user 86 ! defined gas density is invoked, that this routine has been modified. 87 88 !- REMOVE THE FOLLOWING ---------------------------------------------->> 89 90 CALL INIT_ERR_MSG('USR_PHYSICAL_PROP_ROs') 91 WRITE(ERR_MSG,9999) 92 CALL FLUSH_ERR_MSG(ABORT=.TRUE.) 93 94 9999 FORMAT('ERROR 9999: The user-defined drag routine was invoked ', & 95 'but this',/'generic error message exits. Either choose a ', & 96 'different drag law',/'or correct mfix/model/usr_drag.f') 97 98 !- END REMOVE --------------------------------------------------------<< 99 100 M_LP: DO M=1, MMAX 101 IJK_LP: DO IJK = IJKSTART3, IJKEND3 102 IF(WALL_AT(IJK)) cycle IJK_LP 103 104 ! Calculate the solids density. 105 ! RO_S(IJK,M) = 106 107 ENDDO IJK_LP 108 ENDDO M_LP 109 110 RETURN 111 END SUBROUTINE USR_PHYSICAL_PROP_ROs 112 113 114 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv! 115 ! ! 116 ! Subroutine: PHYSICAL_PROP_CPg ! 117 ! Purpose: User hook for calculating the gas phase constant pressure ! 118 ! specific heat. ! 119 ! ! 120 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^! 121 SUBROUTINE USR_PHYSICAL_PROP_CPg 122 123 ! Fluid grid loop bounds. 124 use compar, only: IJKStart3, IJKEnd3 125 ! Function to identify wall cells 126 use functions, only: WALL_AT 127 128 use error_manager 129 130 implicit none 131 132 ! Local Variables: 133 !---------------------------------------------------------------------// 134 ! Loop indicies 135 INTEGER :: IJK 136 !......................................................................! 137 138 ! The following error message is used to make sure that if a user 139 ! defined gas density is invoked, that this routine has been modified. 140 141 !- REMOVE THE FOLLOWING ---------------------------------------------->> 142 143 CALL INIT_ERR_MSG('USR_PHYSICAL_PROP_CPg') 144 WRITE(ERR_MSG,9999) 145 CALL FLUSH_ERR_MSG(ABORT=.TRUE.) 146 147 9999 FORMAT('ERROR 9999: The user-defined drag routine was invoked ', & 148 'but this',/'generic error message exits. Either choose a ', & 149 'different drag law',/'or correct mfix/model/usr_drag.f') 150 151 !- END REMOVE --------------------------------------------------------<< 152 153 IJK_LP: DO IJK = IJKSTART3, IJKEND3 154 IF(WALL_AT(IJK)) cycle IJK_LP 155 156 ! Calculate the fluid density and bulk density 157 ! C_PG(IJK) = 158 159 ENDDO IJK_LP 160 161 RETURN 162 END SUBROUTINE USR_PHYSICAL_PROP_CPg 163 164 165 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv! 166 ! ! 167 ! Subroutine: PHYSICAL_PROP_CPs ! 168 ! Purpose: User hook for calculating solids phase constant pressure ! 169 ! specific heat. ! 170 ! ! 171 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^! 172 SUBROUTINE USR_PHYSICAL_PROP_CPs 173 174 ! Fluid grid loop bounds. 175 use compar, only: IJKStart3, IJKEnd3 176 ! Function to identify wall cells 177 use functions, only: WALL_AT 178 ! Number of solids phases 179 use physprop, only: MMAX 180 181 use error_manager 182 183 implicit none 184 185 ! Local Variables: 186 !---------------------------------------------------------------------// 187 ! Loop indicies 188 INTEGER :: M, IJK 189 !......................................................................! 190 191 ! The following error message is used to make sure that if a user 192 ! defined gas density is invoked, that this routine has been modified. 193 194 !- REMOVE THE FOLLOWING ---------------------------------------------->> 195 196 CALL INIT_ERR_MSG('USR_PHYSICAL_PROP_CPs') 197 WRITE(ERR_MSG,9999) 198 CALL FLUSH_ERR_MSG(ABORT=.TRUE.) 199 200 9999 FORMAT('ERROR 9999: The user-defined drag routine was invoked ', & 201 'but this',/'generic error message exits. Either choose a ', & 202 'different drag law',/'or correct mfix/model/usr_drag.f') 203 204 !- END REMOVE --------------------------------------------------------<< 205 206 M_LP: DO M=1, MMAX 207 IJK_LP: DO IJK = IJKSTART3, IJKEND3 208 IF(WALL_AT(IJK)) cycle IJK_LP 209 210 ! Calculate the solids phase specific heat. 211 ! C_PS(IJK,M) = 212 213 ENDDO IJK_LP 214 ENDDO M_LP 215 216 RETURN 217 END SUBROUTINE USR_PHYSICAL_PROP_CPs 218