C *********************************************************** C Fortran file translated from WHIRL Wed Sep 29 15:51:40 2004 C *********************************************************** C *********************************************************** SUBROUTINE box_forward(ILEV1) use w2f__types use active_module use all_globals_mod IMPLICIT NONE C C **** Parameters and Result **** C INTEGER(w2f__i4) ILEV1 C C **** Local Variables and functions **** C EXTERNAL box_cycle_fields EXTERNAL box_density EXTERNAL box_robert_filter EXTERNAL box_timestep EXTERNAL box_transport INTEGER(w2f__i4) IKEY INTEGER(w2f__i4) ISBYTE PARAMETER ( ISBYTE = 8) INTEGER(w2f__i4) L INTEGER(w2f__i4) NLEV1 PARAMETER ( NLEV1 = 73) INTEGER(w2f__i4) NLEV2 PARAMETER ( NLEV2 = 50) C C **** statements **** C IKEY = ILEV1 CALL box_density(TNOW,SNOW,RHO) CALL box_transport(RHO,UVEL) CALL box_timestep(GAMMA_T,TSTAR,NULLFORCE,UVEL,TNOW,TOLD,TNEW) CALL box_timestep(GAMMA_S,SSTAR,FW,UVEL,SNOW,SOLD,SNEW) CALL box_robert_filter(TNOW,TOLD,TNEW) CALL box_robert_filter(SNOW,SOLD,SNEW) CALL box_cycle_fields() DO L = 1, 3, 1 IF (TNOW(L)%v .LT. (-2.0D00)) THEN TNOW(INT(L))%v = 2.0D00 CALL zero_deriv(TNOW(INT(L))) ENDIF END DO RETURN END SUBROUTINE SUBROUTINE box_final_state() use w2f__types use active_module use all_globals_mod IMPLICIT NONE C C **** Local Variables and functions **** C INTEGER(w2f__i4) ISBYTE PARAMETER ( ISBYTE = 8) INTEGER(w2f__i4) L INTEGER(w2f__i4) NLEV1 PARAMETER ( NLEV1 = 73) INTEGER(w2f__i4) NLEV2 PARAMETER ( NLEV2 = 50) C C **** statements **** C DO L = 1, 6, 1 TSVEC(INT(L))%v = 0.0D00 CALL zero_deriv(TSVEC(INT(L))) END DO DO L = 1, 3, 1 TSVEC(INT(L))%v = TNOW(L)%v CALL setderiv(TSVEC(L),TNOW(L)) TSVEC(INT(L+3))%v = SNOW(L)%v CALL setderiv(TSVEC(L+3),SNOW(L)) END DO RETURN END SUBROUTINE SUBROUTINE box_ini_fields() use w2f__types use active_module use all_globals_mod IMPLICIT NONE C C **** Global Variables **** C REAL(w2f__8) OpenAD_Symbol_140 REAL(w2f__8) OpenAD_Symbol_141 REAL(w2f__8) OpenAD_Symbol_142 REAL(w2f__8) OpenAD_Symbol_138 REAL(w2f__8) OpenAD_Symbol_144 REAL(w2f__8) OpenAD_Symbol_139 REAL(w2f__8) OpenAD_Symbol_150 REAL(w2f__8) OpenAD_Symbol_145 REAL(w2f__8) OpenAD_Symbol_151 REAL(w2f__8) OpenAD_Symbol_146 REAL(w2f__8) OpenAD_Symbol_152 REAL(w2f__8) OpenAD_Symbol_147 REAL(w2f__8) OpenAD_Symbol_153 REAL(w2f__8) OpenAD_Symbol_148 REAL(w2f__8) OpenAD_Symbol_154 type(active) :: OpenAD_Symbol_440 REAL(w2f__8) OpenAD_Symbol_441 REAL(w2f__8) OpenAD_Symbol_436 type(active) :: OpenAD_Symbol_442 REAL(w2f__8) OpenAD_Symbol_437 REAL(w2f__8) OpenAD_Symbol_443 type(active) :: OpenAD_Symbol_438 type(active) :: OpenAD_Symbol_444 REAL(w2f__8) OpenAD_Symbol_439 type(active) :: OpenAD_Symbol_450 REAL(w2f__8) OpenAD_Symbol_445 REAL(w2f__8) OpenAD_Symbol_451 type(active) :: OpenAD_Symbol_446 type(active) :: OpenAD_Symbol_452 REAL(w2f__8) OpenAD_Symbol_447 REAL(w2f__8) OpenAD_Symbol_453 type(active) :: OpenAD_Symbol_448 type(active) :: OpenAD_Symbol_454 REAL(w2f__8) OpenAD_Symbol_449 C C **** Local Variables and functions **** C INTEGER(w2f__i4) ISBYTE PARAMETER ( ISBYTE = 8) INTEGER(w2f__i4) L INTEGER(w2f__i4) NLEV1 PARAMETER ( NLEV1 = 73) INTEGER(w2f__i4) NLEV2 PARAMETER ( NLEV2 = 50) C C **** statements **** C THC_TOT%v = 0.0D00 THC_T%v = 0.0D00 THC_S%v = 0.0D00 METRIC1%v = 0.0D00 METRIC2%v = 0.0D00 METRIC%v = 0.0D00 NULLFORCE(1)%v = 0.0D00 NULLFORCE(2)%v = 0.0D00 CALL zero_deriv(THC_TOT) CALL zero_deriv(THC_T) CALL zero_deriv(THC_S) CALL zero_deriv(METRIC1) CALL zero_deriv(METRIC2) CALL zero_deriv(METRIC) CALL zero_deriv(NULLFORCE(1)) CALL zero_deriv(NULLFORCE(2)) OpenAD_Symbol_141 = (HUNDRED%v/YEAR%v) OpenAD_Symbol_138 = (OpenAD_Symbol_141 * 3.5D+01) OpenAD_Symbol_146 = (AREA(1)%v*OpenAD_Symbol_138) OpenAD_Symbol_139 = OpenAD_Symbol_138 OpenAD_Symbol_144 = (INT(1 _w2f__i8)/YEAR%v) OpenAD_Symbol_145 = (-(HUNDRED%v/(YEAR%v*YEAR%v))) OpenAD_Symbol_142 = 3.5D+01 OpenAD_Symbol_140 = AREA(1)%v FW(1)%v = OpenAD_Symbol_146 OpenAD_Symbol_436 = (OpenAD_Symbol_142 * OpenAD_Symbol_140) OpenAD_Symbol_437 = (OpenAD_Symbol_144 * OpenAD_Symbol_436) OpenAD_Symbol_439 = (OpenAD_Symbol_145 * OpenAD_Symbol_436) OpenAD_Symbol_441 = OpenAD_Symbol_139 CALL setderiv(OpenAD_Symbol_442,AREA(1)) CALL setderiv(OpenAD_Symbol_440,YEAR) CALL setderiv(OpenAD_Symbol_438,HUNDRED) CALL sax(OpenAD_Symbol_437,OpenAD_Symbol_438,FW(1)) CALL saxpy(OpenAD_Symbol_439,OpenAD_Symbol_440,FW(1)) CALL saxpy(OpenAD_Symbol_441,OpenAD_Symbol_442,FW(1)) FW(2)%v = (-FW(1)%v) OpenAD_Symbol_147 = (-1_w2f__i8) TSTAR(1)%v = 2.2D+01 TSTAR(2)%v = 0.0D00 SSTAR(1)%v = 3.6D+01 SSTAR(2)%v = 3.4D+01 OpenAD_Symbol_443 = OpenAD_Symbol_147 CALL setderiv(OpenAD_Symbol_444,FW(1)) CALL zero_deriv(TSTAR(1)) CALL zero_deriv(TSTAR(2)) CALL zero_deriv(SSTAR(1)) CALL zero_deriv(SSTAR(2)) CALL sax(OpenAD_Symbol_443,OpenAD_Symbol_444,FW(2)) OpenAD_Symbol_150 = (SV%v*2.0D+01) OpenAD_Symbol_148 = 2.0D+01 UBAR%v = OpenAD_Symbol_150 T(1)%v = 2.0D+01 T(2)%v = 1.0D00 T(3)%v = 1.0D00 S(1)%v = 3.55D+01 S(2)%v = 3.45D+01 S(3)%v = 3.45D+01 OpenAD_Symbol_445 = OpenAD_Symbol_148 CALL setderiv(OpenAD_Symbol_446,SV) CALL zero_deriv(T(1)) CALL zero_deriv(T(2)) CALL zero_deriv(T(3)) CALL zero_deriv(S(1)) CALL zero_deriv(S(2)) CALL zero_deriv(S(3)) CALL sax(OpenAD_Symbol_445,OpenAD_Symbol_446,UBAR) DO L = 1, 6, 1 TSVEC(INT(L))%v = 0.0D00 CALL zero_deriv(TSVEC(INT(L))) END DO DO L = 1, 3, 1 T(INT(L))%v = (T(L)%v+XX(L)%v) OpenAD_Symbol_151 = 1_w2f__i8 OpenAD_Symbol_152 = 1_w2f__i8 OpenAD_Symbol_447 = OpenAD_Symbol_151 OpenAD_Symbol_449 = OpenAD_Symbol_152 CALL setderiv(OpenAD_Symbol_450,XX(L)) CALL setderiv(OpenAD_Symbol_448,T(L)) CALL sax(OpenAD_Symbol_447,OpenAD_Symbol_448,T(L)) CALL saxpy(OpenAD_Symbol_449,OpenAD_Symbol_450,T(L)) S(INT(L))%v = (S(L)%v+XX(L+3)%v) OpenAD_Symbol_153 = 1_w2f__i8 OpenAD_Symbol_154 = 1_w2f__i8 OpenAD_Symbol_451 = OpenAD_Symbol_153 OpenAD_Symbol_453 = OpenAD_Symbol_154 CALL setderiv(OpenAD_Symbol_454,XX(L+3)) CALL setderiv(OpenAD_Symbol_452,S(L)) CALL sax(OpenAD_Symbol_451,OpenAD_Symbol_452,S(L)) CALL saxpy(OpenAD_Symbol_453,OpenAD_Symbol_454,S(L)) END DO DO L = 1, 3, 1 TNEW(INT(L))%v = T(L)%v CALL setderiv(TNEW(L),T(L)) SNEW(INT(L))%v = S(L)%v CALL setderiv(SNEW(L),S(L)) TOLD(INT(L))%v = T(L)%v CALL setderiv(TOLD(L),T(L)) SOLD(INT(L))%v = S(L)%v CALL setderiv(SOLD(L),S(L)) TNOW(INT(L))%v = T(L)%v CALL setderiv(TNOW(L),T(L)) SNOW(INT(L))%v = S(L)%v CALL setderiv(SNOW(L),S(L)) END DO UVEL%v = UBAR%v CALL setderiv(UVEL,UBAR) RETURN END SUBROUTINE SUBROUTINE box_ini_params() use w2f__types use active_module use all_globals_mod IMPLICIT NONE C C **** Global Variables **** C REAL(w2f__8) OpenAD_Symbol_201 REAL(w2f__8) OpenAD_Symbol_202 REAL(w2f__8) OpenAD_Symbol_203 REAL(w2f__8) OpenAD_Symbol_204 REAL(w2f__8) OpenAD_Symbol_210 REAL(w2f__8) OpenAD_Symbol_205 REAL(w2f__8) OpenAD_Symbol_160 REAL(w2f__8) OpenAD_Symbol_155 REAL(w2f__8) OpenAD_Symbol_211 REAL(w2f__8) OpenAD_Symbol_206 REAL(w2f__8) OpenAD_Symbol_161 REAL(w2f__8) OpenAD_Symbol_156 REAL(w2f__8) OpenAD_Symbol_162 REAL(w2f__8) OpenAD_Symbol_157 REAL(w2f__8) OpenAD_Symbol_213 REAL(w2f__8) OpenAD_Symbol_208 REAL(w2f__8) OpenAD_Symbol_163 REAL(w2f__8) OpenAD_Symbol_158 REAL(w2f__8) OpenAD_Symbol_214 REAL(w2f__8) OpenAD_Symbol_209 REAL(w2f__8) OpenAD_Symbol_164 REAL(w2f__8) OpenAD_Symbol_159 REAL(w2f__8) OpenAD_Symbol_220 REAL(w2f__8) OpenAD_Symbol_170 REAL(w2f__8) OpenAD_Symbol_165 REAL(w2f__8) OpenAD_Symbol_216 REAL(w2f__8) OpenAD_Symbol_171 REAL(w2f__8) OpenAD_Symbol_166 REAL(w2f__8) OpenAD_Symbol_222 REAL(w2f__8) OpenAD_Symbol_217 REAL(w2f__8) OpenAD_Symbol_172 REAL(w2f__8) OpenAD_Symbol_167 REAL(w2f__8) OpenAD_Symbol_223 REAL(w2f__8) OpenAD_Symbol_173 REAL(w2f__8) OpenAD_Symbol_168 REAL(w2f__8) OpenAD_Symbol_219 REAL(w2f__8) OpenAD_Symbol_174 REAL(w2f__8) OpenAD_Symbol_169 REAL(w2f__8) OpenAD_Symbol_225 REAL(w2f__8) OpenAD_Symbol_180 REAL(w2f__8) OpenAD_Symbol_175 REAL(w2f__8) OpenAD_Symbol_231 REAL(w2f__8) OpenAD_Symbol_226 REAL(w2f__8) OpenAD_Symbol_181 REAL(w2f__8) OpenAD_Symbol_176 REAL(w2f__8) OpenAD_Symbol_232 REAL(w2f__8) OpenAD_Symbol_182 REAL(w2f__8) OpenAD_Symbol_177 REAL(w2f__8) OpenAD_Symbol_228 REAL(w2f__8) OpenAD_Symbol_183 REAL(w2f__8) OpenAD_Symbol_178 REAL(w2f__8) OpenAD_Symbol_234 REAL(w2f__8) OpenAD_Symbol_229 REAL(w2f__8) OpenAD_Symbol_184 REAL(w2f__8) OpenAD_Symbol_179 REAL(w2f__8) OpenAD_Symbol_240 REAL(w2f__8) OpenAD_Symbol_235 REAL(w2f__8) OpenAD_Symbol_190 REAL(w2f__8) OpenAD_Symbol_241 REAL(w2f__8) OpenAD_Symbol_186 REAL(w2f__8) OpenAD_Symbol_242 REAL(w2f__8) OpenAD_Symbol_237 REAL(w2f__8) OpenAD_Symbol_192 REAL(w2f__8) OpenAD_Symbol_187 REAL(w2f__8) OpenAD_Symbol_243 REAL(w2f__8) OpenAD_Symbol_238 REAL(w2f__8) OpenAD_Symbol_193 REAL(w2f__8) OpenAD_Symbol_188 REAL(w2f__8) OpenAD_Symbol_244 REAL(w2f__8) OpenAD_Symbol_239 REAL(w2f__8) OpenAD_Symbol_250 REAL(w2f__8) OpenAD_Symbol_195 REAL(w2f__8) OpenAD_Symbol_251 REAL(w2f__8) OpenAD_Symbol_246 REAL(w2f__8) OpenAD_Symbol_196 REAL(w2f__8) OpenAD_Symbol_252 REAL(w2f__8) OpenAD_Symbol_247 REAL(w2f__8) OpenAD_Symbol_253 REAL(w2f__8) OpenAD_Symbol_248 REAL(w2f__8) OpenAD_Symbol_198 REAL(w2f__8) OpenAD_Symbol_254 REAL(w2f__8) OpenAD_Symbol_249 REAL(w2f__8) OpenAD_Symbol_199 REAL(w2f__8) OpenAD_Symbol_260 REAL(w2f__8) OpenAD_Symbol_255 REAL(w2f__8) OpenAD_Symbol_261 REAL(w2f__8) OpenAD_Symbol_262 REAL(w2f__8) OpenAD_Symbol_257 REAL(w2f__8) OpenAD_Symbol_263 REAL(w2f__8) OpenAD_Symbol_258 REAL(w2f__8) OpenAD_Symbol_264 REAL(w2f__8) OpenAD_Symbol_259 REAL(w2f__8) OpenAD_Symbol_270 REAL(w2f__8) OpenAD_Symbol_265 REAL(w2f__8) OpenAD_Symbol_271 REAL(w2f__8) OpenAD_Symbol_266 REAL(w2f__8) OpenAD_Symbol_272 REAL(w2f__8) OpenAD_Symbol_267 REAL(w2f__8) OpenAD_Symbol_273 REAL(w2f__8) OpenAD_Symbol_268 REAL(w2f__8) OpenAD_Symbol_274 REAL(w2f__8) OpenAD_Symbol_269 REAL(w2f__8) OpenAD_Symbol_275 type(active) :: OpenAD_Symbol_500 REAL(w2f__8) OpenAD_Symbol_501 type(active) :: OpenAD_Symbol_502 REAL(w2f__8) OpenAD_Symbol_503 type(active) :: OpenAD_Symbol_504 REAL(w2f__8) OpenAD_Symbol_510 REAL(w2f__8) OpenAD_Symbol_505 type(active) :: OpenAD_Symbol_460 REAL(w2f__8) OpenAD_Symbol_455 type(active) :: OpenAD_Symbol_511 type(active) :: OpenAD_Symbol_506 REAL(w2f__8) OpenAD_Symbol_461 type(active) :: OpenAD_Symbol_456 REAL(w2f__8) OpenAD_Symbol_512 REAL(w2f__8) OpenAD_Symbol_507 type(active) :: OpenAD_Symbol_462 REAL(w2f__8) OpenAD_Symbol_457 type(active) :: OpenAD_Symbol_513 REAL(w2f__8) OpenAD_Symbol_508 REAL(w2f__8) OpenAD_Symbol_463 type(active) :: OpenAD_Symbol_458 REAL(w2f__8) OpenAD_Symbol_514 type(active) :: OpenAD_Symbol_509 type(active) :: OpenAD_Symbol_464 REAL(w2f__8) OpenAD_Symbol_459 REAL(w2f__8) OpenAD_Symbol_520 type(active) :: OpenAD_Symbol_515 type(active) :: OpenAD_Symbol_470 REAL(w2f__8) OpenAD_Symbol_465 type(active) :: OpenAD_Symbol_521 REAL(w2f__8) OpenAD_Symbol_516 REAL(w2f__8) OpenAD_Symbol_471 type(active) :: OpenAD_Symbol_466 REAL(w2f__8) OpenAD_Symbol_522 type(active) :: OpenAD_Symbol_517 type(active) :: OpenAD_Symbol_472 REAL(w2f__8) OpenAD_Symbol_467 type(active) :: OpenAD_Symbol_523 REAL(w2f__8) OpenAD_Symbol_518 REAL(w2f__8) OpenAD_Symbol_473 type(active) :: OpenAD_Symbol_468 REAL(w2f__8) OpenAD_Symbol_524 type(active) :: OpenAD_Symbol_519 type(active) :: OpenAD_Symbol_474 REAL(w2f__8) OpenAD_Symbol_469 REAL(w2f__8) OpenAD_Symbol_530 type(active) :: OpenAD_Symbol_525 type(active) :: OpenAD_Symbol_480 REAL(w2f__8) OpenAD_Symbol_475 type(active) :: OpenAD_Symbol_531 REAL(w2f__8) OpenAD_Symbol_526 REAL(w2f__8) OpenAD_Symbol_481 type(active) :: OpenAD_Symbol_476 REAL(w2f__8) OpenAD_Symbol_532 type(active) :: OpenAD_Symbol_527 type(active) :: OpenAD_Symbol_482 REAL(w2f__8) OpenAD_Symbol_477 type(active) :: OpenAD_Symbol_533 REAL(w2f__8) OpenAD_Symbol_528 REAL(w2f__8) OpenAD_Symbol_483 type(active) :: OpenAD_Symbol_478 REAL(w2f__8) OpenAD_Symbol_534 type(active) :: OpenAD_Symbol_529 type(active) :: OpenAD_Symbol_484 REAL(w2f__8) OpenAD_Symbol_479 type(active) :: OpenAD_Symbol_540 type(active) :: OpenAD_Symbol_535 type(active) :: OpenAD_Symbol_490 REAL(w2f__8) OpenAD_Symbol_485 REAL(w2f__8) OpenAD_Symbol_541 REAL(w2f__8) OpenAD_Symbol_536 REAL(w2f__8) OpenAD_Symbol_491 type(active) :: OpenAD_Symbol_486 type(active) :: OpenAD_Symbol_542 REAL(w2f__8) OpenAD_Symbol_537 type(active) :: OpenAD_Symbol_492 REAL(w2f__8) OpenAD_Symbol_487 REAL(w2f__8) OpenAD_Symbol_543 type(active) :: OpenAD_Symbol_538 REAL(w2f__8) OpenAD_Symbol_493 type(active) :: OpenAD_Symbol_488 type(active) :: OpenAD_Symbol_544 REAL(w2f__8) OpenAD_Symbol_539 type(active) :: OpenAD_Symbol_494 REAL(w2f__8) OpenAD_Symbol_489 type(active) :: OpenAD_Symbol_550 REAL(w2f__8) OpenAD_Symbol_545 REAL(w2f__8) OpenAD_Symbol_495 REAL(w2f__8) OpenAD_Symbol_551 type(active) :: OpenAD_Symbol_546 type(active) :: OpenAD_Symbol_496 type(active) :: OpenAD_Symbol_552 REAL(w2f__8) OpenAD_Symbol_547 REAL(w2f__8) OpenAD_Symbol_497 REAL(w2f__8) OpenAD_Symbol_553 type(active) :: OpenAD_Symbol_548 REAL(w2f__8) OpenAD_Symbol_498 type(active) :: OpenAD_Symbol_554 REAL(w2f__8) OpenAD_Symbol_549 REAL(w2f__8) OpenAD_Symbol_499 type(active) :: OpenAD_Symbol_560 REAL(w2f__8) OpenAD_Symbol_555 REAL(w2f__8) OpenAD_Symbol_561 type(active) :: OpenAD_Symbol_556 type(active) :: OpenAD_Symbol_562 REAL(w2f__8) OpenAD_Symbol_557 type(active) :: OpenAD_Symbol_558 REAL(w2f__8) OpenAD_Symbol_559 C C **** Local Variables and functions **** C INTEGER(w2f__i4) I INTEGER(w2f__i4) ISBYTE PARAMETER ( ISBYTE = 8) INTEGER(w2f__i4) J INTEGER(w2f__i4) L INTEGER(w2f__i4) NLEV1 PARAMETER ( NLEV1 = 73) INTEGER(w2f__i4) NLEV2 PARAMETER ( NLEV2 = 50) C C **** statements **** C VERBMODE = .FALSE. BLENGTH(1)%v = 5.0D+08 BLENGTH(2)%v = 1.0D+08 CALL zero_deriv(BLENGTH(1)) CALL zero_deriv(BLENGTH(2)) BLENGTH(3)%v = BLENGTH(1)%v BHEIGHT(1)%v = 1.0D+05 BHEIGHT(3)%v = 4.0D+05 CALL setderiv(BLENGTH(3),BLENGTH(1)) CALL zero_deriv(BHEIGHT(1)) CALL zero_deriv(BHEIGHT(3)) BHEIGHT(2)%v = (BHEIGHT(1)%v+BHEIGHT(3)%v) OpenAD_Symbol_155 = 1_w2f__i8 OpenAD_Symbol_156 = 1_w2f__i8 OpenAD_Symbol_455 = OpenAD_Symbol_155 OpenAD_Symbol_457 = OpenAD_Symbol_156 CALL setderiv(OpenAD_Symbol_458,BHEIGHT(3)) CALL setderiv(OpenAD_Symbol_456,BHEIGHT(1)) CALL sax(OpenAD_Symbol_455,OpenAD_Symbol_456,BHEIGHT(2)) CALL saxpy(OpenAD_Symbol_457,OpenAD_Symbol_458,BHEIGHT(2)) OpenAD_Symbol_157 = (BHEIGHT(1)%v+BHEIGHT(3)%v) OpenAD_Symbol_162 = (BHEIGHT(1)%v/OpenAD_Symbol_157) OpenAD_Symbol_158 = (INT(1_w2f__i8) / OpenAD_Symbol_157) OpenAD_Symbol_160 = 1_w2f__i8 OpenAD_Symbol_161 = 1_w2f__i8 OpenAD_Symbol_159 = (-(BHEIGHT(1)%v/(OpenAD_Symbol_157*OpenAD_Symb +ol_157))) DELTA%v = OpenAD_Symbol_162 BWIDTH%v = 4.0D+08 OpenAD_Symbol_459 = (OpenAD_Symbol_160 * OpenAD_Symbol_159) OpenAD_Symbol_461 = (OpenAD_Symbol_161 * OpenAD_Symbol_159) OpenAD_Symbol_463 = OpenAD_Symbol_158 CALL setderiv(OpenAD_Symbol_464,BHEIGHT(1)) CALL setderiv(OpenAD_Symbol_462,BHEIGHT(3)) CALL setderiv(OpenAD_Symbol_460,BHEIGHT(1)) CALL zero_deriv(BWIDTH) CALL sax(OpenAD_Symbol_459,OpenAD_Symbol_460,DELTA) CALL saxpy(OpenAD_Symbol_461,OpenAD_Symbol_462,DELTA) CALL saxpy(OpenAD_Symbol_463,OpenAD_Symbol_464,DELTA) OpenAD_Symbol_165 = (BLENGTH(1)%v*BWIDTH%v) OpenAD_Symbol_163 = BWIDTH%v OpenAD_Symbol_164 = BLENGTH(1)%v AREA(1)%v = OpenAD_Symbol_165 OpenAD_Symbol_465 = OpenAD_Symbol_163 OpenAD_Symbol_467 = OpenAD_Symbol_164 CALL setderiv(OpenAD_Symbol_468,BWIDTH) CALL setderiv(OpenAD_Symbol_466,BLENGTH(1)) CALL sax(OpenAD_Symbol_465,OpenAD_Symbol_466,AREA(1)) CALL saxpy(OpenAD_Symbol_467,OpenAD_Symbol_468,AREA(1)) OpenAD_Symbol_168 = (BLENGTH(2)%v*BWIDTH%v) OpenAD_Symbol_166 = BWIDTH%v OpenAD_Symbol_167 = BLENGTH(2)%v AREA(2)%v = OpenAD_Symbol_168 OpenAD_Symbol_469 = OpenAD_Symbol_166 OpenAD_Symbol_471 = OpenAD_Symbol_167 CALL setderiv(OpenAD_Symbol_472,BWIDTH) CALL setderiv(OpenAD_Symbol_470,BLENGTH(2)) CALL sax(OpenAD_Symbol_469,OpenAD_Symbol_470,AREA(2)) CALL saxpy(OpenAD_Symbol_471,OpenAD_Symbol_472,AREA(2)) OpenAD_Symbol_171 = (BLENGTH(3)%v*BWIDTH%v) OpenAD_Symbol_169 = BWIDTH%v OpenAD_Symbol_170 = BLENGTH(3)%v AREA(3)%v = OpenAD_Symbol_171 OpenAD_Symbol_473 = OpenAD_Symbol_169 OpenAD_Symbol_475 = OpenAD_Symbol_170 CALL setderiv(OpenAD_Symbol_476,BWIDTH) CALL setderiv(OpenAD_Symbol_474,BLENGTH(3)) CALL sax(OpenAD_Symbol_473,OpenAD_Symbol_474,AREA(3)) CALL saxpy(OpenAD_Symbol_475,OpenAD_Symbol_476,AREA(3)) OpenAD_Symbol_174 = (AREA(1)%v*BHEIGHT(1)%v) OpenAD_Symbol_172 = BHEIGHT(1)%v OpenAD_Symbol_173 = AREA(1)%v VOL(1)%v = OpenAD_Symbol_174 OpenAD_Symbol_477 = OpenAD_Symbol_172 OpenAD_Symbol_479 = OpenAD_Symbol_173 CALL setderiv(OpenAD_Symbol_480,BHEIGHT(1)) CALL setderiv(OpenAD_Symbol_478,AREA(1)) CALL sax(OpenAD_Symbol_477,OpenAD_Symbol_478,VOL(1)) CALL saxpy(OpenAD_Symbol_479,OpenAD_Symbol_480,VOL(1)) OpenAD_Symbol_175 = (BHEIGHT(1)%v+BHEIGHT(3)%v) OpenAD_Symbol_180 = (AREA(2)%v*OpenAD_Symbol_175) OpenAD_Symbol_176 = OpenAD_Symbol_175 OpenAD_Symbol_178 = 1_w2f__i8 OpenAD_Symbol_179 = 1_w2f__i8 OpenAD_Symbol_177 = AREA(2)%v VOL(2)%v = OpenAD_Symbol_180 OpenAD_Symbol_481 = (OpenAD_Symbol_178 * OpenAD_Symbol_177) OpenAD_Symbol_483 = (OpenAD_Symbol_179 * OpenAD_Symbol_177) OpenAD_Symbol_485 = OpenAD_Symbol_176 CALL setderiv(OpenAD_Symbol_486,AREA(2)) CALL setderiv(OpenAD_Symbol_484,BHEIGHT(3)) CALL setderiv(OpenAD_Symbol_482,BHEIGHT(1)) CALL sax(OpenAD_Symbol_481,OpenAD_Symbol_482,VOL(2)) CALL saxpy(OpenAD_Symbol_483,OpenAD_Symbol_484,VOL(2)) CALL saxpy(OpenAD_Symbol_485,OpenAD_Symbol_486,VOL(2)) OpenAD_Symbol_183 = (AREA(3)%v*BHEIGHT(3)%v) OpenAD_Symbol_181 = BHEIGHT(3)%v OpenAD_Symbol_182 = AREA(3)%v VOL(3)%v = OpenAD_Symbol_183 ROBERT_FILTER_COEFF%v = 2.5D-01 HUNDRED%v = 1.0D+02 THOUSAND%v = 1.0D+03 DAY%v = 8.64D+04 OpenAD_Symbol_487 = OpenAD_Symbol_181 OpenAD_Symbol_489 = OpenAD_Symbol_182 CALL setderiv(OpenAD_Symbol_490,BHEIGHT(3)) CALL setderiv(OpenAD_Symbol_488,AREA(3)) CALL zero_deriv(ROBERT_FILTER_COEFF) CALL zero_deriv(HUNDRED) CALL zero_deriv(THOUSAND) CALL zero_deriv(DAY) CALL sax(OpenAD_Symbol_487,OpenAD_Symbol_488,VOL(3)) CALL saxpy(OpenAD_Symbol_489,OpenAD_Symbol_490,VOL(3)) OpenAD_Symbol_186 = (DAY%v*3.65D+02) OpenAD_Symbol_184 = 3.65D+02 YEAR%v = OpenAD_Symbol_186 SV%v = 1.0D+12 OpenAD_Symbol_491 = OpenAD_Symbol_184 CALL setderiv(OpenAD_Symbol_492,DAY) CALL zero_deriv(SV) CALL sax(OpenAD_Symbol_491,OpenAD_Symbol_492,YEAR) OpenAD_Symbol_187 = (SV%v*1.6D+01) OpenAD_Symbol_192 = (OpenAD_Symbol_187 / > 4.00000000000000019169D-04) OpenAD_Symbol_190 = 1.6D+01 OpenAD_Symbol_188 = (INT(1_w2f__i8) / 4.00000000000000019169D-04) U0%v = OpenAD_Symbol_192 ALPHA%v = 1.6679999999999998864D-04 BETA%v = 7.81000000000000010186D-04 DAYS_PER_50M_MIXED_LAYER%v = 5.0D+01 OpenAD_Symbol_493 = (OpenAD_Symbol_190 * OpenAD_Symbol_188) CALL setderiv(OpenAD_Symbol_494,SV) CALL zero_deriv(ALPHA) CALL zero_deriv(BETA) CALL zero_deriv(DAYS_PER_50M_MIXED_LAYER) CALL sax(OpenAD_Symbol_493,OpenAD_Symbol_494,U0) OpenAD_Symbol_193 = (DAY%v*3.0D+02) OpenAD_Symbol_198 = (1.0D00 / OpenAD_Symbol_193) OpenAD_Symbol_196 = 3.0D+02 OpenAD_Symbol_195 = (-(1.0D00 /(OpenAD_Symbol_193 * > OpenAD_Symbol_193))) GAMMA_T%v = OpenAD_Symbol_198 OpenAD_Symbol_495 = (OpenAD_Symbol_196 * OpenAD_Symbol_195) CALL setderiv(OpenAD_Symbol_496,DAY) CALL sax(OpenAD_Symbol_495,OpenAD_Symbol_496,GAMMA_T) OpenAD_Symbol_202 = (BHEIGHT(1)%v/5.0D+03) OpenAD_Symbol_203 = (DAY%v*DAYS_PER_50M_MIXED_LAYER%v) OpenAD_Symbol_199 = (OpenAD_Symbol_202 * OpenAD_Symbol_203) OpenAD_Symbol_210 = (0.0D00 / OpenAD_Symbol_199) OpenAD_Symbol_206 = (INT(1_w2f__i8) / 5.0D+03) OpenAD_Symbol_204 = OpenAD_Symbol_203 OpenAD_Symbol_208 = DAYS_PER_50M_MIXED_LAYER%v OpenAD_Symbol_209 = DAY%v OpenAD_Symbol_205 = OpenAD_Symbol_202 OpenAD_Symbol_201 = (-(0.0D00 /(OpenAD_Symbol_199 * > OpenAD_Symbol_199))) GAMMA_S%v = OpenAD_Symbol_210 EPSILON_IC%v = (-1.00000000000000004792D-04) OpenAD_Symbol_497 = (OpenAD_Symbol_204 * OpenAD_Symbol_201) OpenAD_Symbol_498 = (OpenAD_Symbol_205 * OpenAD_Symbol_201) OpenAD_Symbol_499 = (OpenAD_Symbol_208 * OpenAD_Symbol_498) OpenAD_Symbol_501 = (OpenAD_Symbol_209 * OpenAD_Symbol_498) OpenAD_Symbol_503 = (OpenAD_Symbol_206 * OpenAD_Symbol_497) CALL setderiv(OpenAD_Symbol_504,BHEIGHT(1)) CALL setderiv(OpenAD_Symbol_502,DAYS_PER_50M_MIXED_LAYER) CALL setderiv(OpenAD_Symbol_500,DAY) CALL zero_deriv(EPSILON_IC) CALL sax(OpenAD_Symbol_499,OpenAD_Symbol_500,GAMMA_S) CALL saxpy(OpenAD_Symbol_501,OpenAD_Symbol_502,GAMMA_S) CALL saxpy(OpenAD_Symbol_503,OpenAD_Symbol_504,GAMMA_S) OpenAD_Symbol_213 = (DAY%v*1.5D+01) OpenAD_Symbol_211 = 1.5D+01 NOISE_CORRELATION_TIME%v = OpenAD_Symbol_213 OpenAD_Symbol_216 = (DAY%v*5.0D00) OpenAD_Symbol_214 = 5.0D00 DELTA_T%v = OpenAD_Symbol_216 OpenAD_Symbol_219 = (YEAR%v*5.0D+01) OpenAD_Symbol_217 = 5.0D+01 INTEGRATION_TIME%v = OpenAD_Symbol_219 OpenAD_Symbol_505 = OpenAD_Symbol_211 OpenAD_Symbol_507 = OpenAD_Symbol_214 OpenAD_Symbol_508 = OpenAD_Symbol_217 CALL setderiv(OpenAD_Symbol_509,YEAR) CALL setderiv(OpenAD_Symbol_506,DAY) CALL sax(OpenAD_Symbol_505,OpenAD_Symbol_506,NOISE_CORRELATION_TIM +E) CALL sax(OpenAD_Symbol_507,OpenAD_Symbol_506,DELTA_T) CALL sax(OpenAD_Symbol_508,OpenAD_Symbol_509,INTEGRATION_TIME) N_MAX = INT((INTEGRATION_TIME%v/DELTA_T%v)) FDEPS%v = 9.99999999999999954748D-07 CALL zero_deriv(FDEPS) OpenAD_Symbol_222 = (THOUSAND%v/3.2680000000000000715D-02) OpenAD_Symbol_220 = (INT(1_w2f__i8) / 3.2680000000000000715D-02) Y(1)%v = OpenAD_Symbol_222 OpenAD_Symbol_510 = OpenAD_Symbol_220 CALL setderiv(OpenAD_Symbol_511,THOUSAND) CALL sax(OpenAD_Symbol_510,OpenAD_Symbol_511,Y(1)) OpenAD_Symbol_225 = (THOUSAND%v/7.9399999999999991418D-03) OpenAD_Symbol_223 = (INT(1_w2f__i8) / 7.9399999999999991418D-03) Y(2)%v = OpenAD_Symbol_225 OpenAD_Symbol_512 = OpenAD_Symbol_223 CALL setderiv(OpenAD_Symbol_513,THOUSAND) CALL sax(OpenAD_Symbol_512,OpenAD_Symbol_513,Y(2)) OpenAD_Symbol_228 = (THOUSAND%v/1.39999999999999998578D-03) OpenAD_Symbol_226 = (INT(1_w2f__i8) / 1.39999999999999998578D-03) Y(3)%v = OpenAD_Symbol_228 OpenAD_Symbol_514 = OpenAD_Symbol_226 CALL setderiv(OpenAD_Symbol_515,THOUSAND) CALL sax(OpenAD_Symbol_514,OpenAD_Symbol_515,Y(3)) OpenAD_Symbol_231 = (THOUSAND%v/1.41699999999999992628D-01) OpenAD_Symbol_229 = (INT(1_w2f__i8) / 1.41699999999999992628D-01) Y(4)%v = OpenAD_Symbol_231 OpenAD_Symbol_516 = OpenAD_Symbol_229 CALL setderiv(OpenAD_Symbol_517,THOUSAND) CALL sax(OpenAD_Symbol_516,OpenAD_Symbol_517,Y(4)) OpenAD_Symbol_234 = (THOUSAND%v/1.28599999999999992095D-01) OpenAD_Symbol_232 = (INT(1_w2f__i8) / 1.28599999999999992095D-01) Y(5)%v = OpenAD_Symbol_234 OpenAD_Symbol_518 = OpenAD_Symbol_232 CALL setderiv(OpenAD_Symbol_519,THOUSAND) CALL sax(OpenAD_Symbol_518,OpenAD_Symbol_519,Y(5)) OpenAD_Symbol_237 = (THOUSAND%v/8.7800000000000003042D-02) OpenAD_Symbol_235 = (INT(1_w2f__i8) / 8.7800000000000003042D-02) Y(6)%v = OpenAD_Symbol_237 OpenAD_Symbol_520 = OpenAD_Symbol_235 CALL setderiv(OpenAD_Symbol_521,THOUSAND) CALL sax(OpenAD_Symbol_520,OpenAD_Symbol_521,Y(6)) OpenAD_Symbol_240 = (ALPHA%v*DELTA%v) OpenAD_Symbol_238 = DELTA%v OpenAD_Symbol_239 = ALPHA%v R(1)%v = OpenAD_Symbol_240 OpenAD_Symbol_522 = OpenAD_Symbol_238 OpenAD_Symbol_524 = OpenAD_Symbol_239 CALL setderiv(OpenAD_Symbol_525,DELTA) CALL setderiv(OpenAD_Symbol_523,ALPHA) CALL sax(OpenAD_Symbol_522,OpenAD_Symbol_523,R(1)) CALL saxpy(OpenAD_Symbol_524,OpenAD_Symbol_525,R(1)) R(2)%v = (-ALPHA%v) OpenAD_Symbol_241 = (-1_w2f__i8) OpenAD_Symbol_526 = OpenAD_Symbol_241 CALL setderiv(OpenAD_Symbol_527,ALPHA) CALL sax(OpenAD_Symbol_526,OpenAD_Symbol_527,R(2)) OpenAD_Symbol_242 = (1.0D00-DELTA%v) OpenAD_Symbol_247 = (ALPHA%v*OpenAD_Symbol_242) OpenAD_Symbol_243 = OpenAD_Symbol_242 OpenAD_Symbol_246 = (-1_w2f__i8) OpenAD_Symbol_244 = ALPHA%v R(3)%v = OpenAD_Symbol_247 OpenAD_Symbol_528 = (OpenAD_Symbol_246 * OpenAD_Symbol_244) OpenAD_Symbol_530 = OpenAD_Symbol_243 CALL setderiv(OpenAD_Symbol_531,ALPHA) CALL setderiv(OpenAD_Symbol_529,DELTA) CALL sax(OpenAD_Symbol_528,OpenAD_Symbol_529,R(3)) CALL saxpy(OpenAD_Symbol_530,OpenAD_Symbol_531,R(3)) OpenAD_Symbol_251 = (-(BETA%v*DELTA%v)) OpenAD_Symbol_249 = DELTA%v OpenAD_Symbol_250 = BETA%v OpenAD_Symbol_248 = (-1_w2f__i8) R(4)%v = OpenAD_Symbol_251 OpenAD_Symbol_532 = (OpenAD_Symbol_249 * OpenAD_Symbol_248) OpenAD_Symbol_534 = (OpenAD_Symbol_250 * OpenAD_Symbol_248) CALL setderiv(OpenAD_Symbol_535,DELTA) CALL setderiv(OpenAD_Symbol_533,BETA) CALL sax(OpenAD_Symbol_532,OpenAD_Symbol_533,R(4)) CALL saxpy(OpenAD_Symbol_534,OpenAD_Symbol_535,R(4)) R(5)%v = BETA%v CALL setderiv(R(5),BETA) OpenAD_Symbol_253 = (1.0D00-DELTA%v) OpenAD_Symbol_258 = (-(BETA%v*OpenAD_Symbol_253)) OpenAD_Symbol_254 = OpenAD_Symbol_253 OpenAD_Symbol_257 = (-1_w2f__i8) OpenAD_Symbol_255 = BETA%v OpenAD_Symbol_252 = (-1_w2f__i8) R(6)%v = OpenAD_Symbol_258 OpenAD_Symbol_536 = (OpenAD_Symbol_257 * OpenAD_Symbol_255) OpenAD_Symbol_537 = (OpenAD_Symbol_254 * OpenAD_Symbol_252) OpenAD_Symbol_539 = (OpenAD_Symbol_536 * OpenAD_Symbol_252) CALL setderiv(OpenAD_Symbol_540,DELTA) CALL setderiv(OpenAD_Symbol_538,BETA) CALL sax(OpenAD_Symbol_537,OpenAD_Symbol_538,R(6)) CALL saxpy(OpenAD_Symbol_539,OpenAD_Symbol_540,R(6)) DO L = 1, 6, 1 OpenAD_Symbol_259 = (R(L)%v*U0%v) OpenAD_Symbol_264 = (OpenAD_Symbol_259/SV%v) OpenAD_Symbol_262 = U0%v OpenAD_Symbol_263 = R(L)%v OpenAD_Symbol_260 = (INT(1 _w2f__i8)/SV%v) OpenAD_Symbol_261 = (-(OpenAD_Symbol_259/(SV%v*SV%v))) R(INT(L))%v = OpenAD_Symbol_264 OpenAD_Symbol_541 = (OpenAD_Symbol_262 * OpenAD_Symbol_260) OpenAD_Symbol_543 = (OpenAD_Symbol_263 * OpenAD_Symbol_260) OpenAD_Symbol_545 = OpenAD_Symbol_261 CALL setderiv(OpenAD_Symbol_546,SV) CALL setderiv(OpenAD_Symbol_544,U0) CALL setderiv(OpenAD_Symbol_542,R(L)) CALL sax(OpenAD_Symbol_541,OpenAD_Symbol_542,R(L)) CALL saxpy(OpenAD_Symbol_543,OpenAD_Symbol_544,R(L)) CALL saxpy(OpenAD_Symbol_545,OpenAD_Symbol_546,R(L)) END DO DO L = 1, 6, 1 IF(L .LE. INT(3_w2f__i8)) THEN PROJ_T(INT(L))%v = 1.0D00 PROJ_S(INT(L))%v = 0.0D00 CALL zero_deriv(PROJ_T(INT(L))) CALL zero_deriv(PROJ_S(INT(L))) ELSE PROJ_T(INT(L))%v = 0.0D00 PROJ_S(INT(L))%v = 1.0D00 CALL zero_deriv(PROJ_T(INT(L))) CALL zero_deriv(PROJ_S(INT(L))) ENDIF END DO DO L = 1, 6, 1 OpenAD_Symbol_267 = (PROJ_T(L)%v*R(L)%v) OpenAD_Symbol_265 = R(L)%v OpenAD_Symbol_266 = PROJ_T(L)%v R_T(INT(L))%v = OpenAD_Symbol_267 OpenAD_Symbol_547 = OpenAD_Symbol_265 OpenAD_Symbol_549 = OpenAD_Symbol_266 CALL setderiv(OpenAD_Symbol_550,R(L)) CALL setderiv(OpenAD_Symbol_548,PROJ_T(L)) CALL sax(OpenAD_Symbol_547,OpenAD_Symbol_548,R_T(L)) CALL saxpy(OpenAD_Symbol_549,OpenAD_Symbol_550,R_T(L)) OpenAD_Symbol_270 = (PROJ_S(L)%v*R(L)%v) OpenAD_Symbol_268 = R(L)%v OpenAD_Symbol_269 = PROJ_S(L)%v R_S(INT(L))%v = OpenAD_Symbol_270 OpenAD_Symbol_551 = OpenAD_Symbol_268 OpenAD_Symbol_553 = OpenAD_Symbol_269 CALL setderiv(OpenAD_Symbol_554,R(L)) CALL setderiv(OpenAD_Symbol_552,PROJ_S(L)) CALL sax(OpenAD_Symbol_551,OpenAD_Symbol_552,R_S(L)) CALL saxpy(OpenAD_Symbol_553,OpenAD_Symbol_554,R_S(L)) END DO R1(1)%v = 1.0D00 R1(2)%v = (-1.0D00) R1(3)%v = 0.0D00 R1(4)%v = 0.0D00 R1(5)%v = 0.0D00 R1(6)%v = 0.0D00 CALL zero_deriv(R1(1)) CALL zero_deriv(R1(2)) CALL zero_deriv(R1(3)) CALL zero_deriv(R1(4)) CALL zero_deriv(R1(5)) CALL zero_deriv(R1(6)) DO J = 1, 6, 1 DO I = 1, 6, 1 OpenAD_Symbol_273 = (R(I)%v*R(J)%v) OpenAD_Symbol_271 = R(J)%v OpenAD_Symbol_272 = R(I)%v X(INT(I),INT(J))%v = OpenAD_Symbol_273 OpenAD_Symbol_555 = OpenAD_Symbol_271 OpenAD_Symbol_557 = OpenAD_Symbol_272 CALL setderiv(OpenAD_Symbol_558,R(J)) CALL setderiv(OpenAD_Symbol_556,R(I)) CALL sax(OpenAD_Symbol_555,OpenAD_Symbol_556,X(I,J)) CALL saxpy(OpenAD_Symbol_557,OpenAD_Symbol_558,X(I,J)) END DO END DO EPSILON_REGULARIZE%v = 1.00000000000000002092D-08 CALL zero_deriv(EPSILON_REGULARIZE) DO J = 1, 6, 1 DO I = 1, 6, 1 X(INT(I),INT(J))%v = (X(I,J)%v+EPSILON_REGULARIZE%v) OpenAD_Symbol_274 = 1_w2f__i8 OpenAD_Symbol_275 = 1_w2f__i8 OpenAD_Symbol_559 = OpenAD_Symbol_274 OpenAD_Symbol_561 = OpenAD_Symbol_275 CALL setderiv(OpenAD_Symbol_562,EPSILON_REGULARIZE) CALL setderiv(OpenAD_Symbol_560,X(I,J)) CALL sax(OpenAD_Symbol_559,OpenAD_Symbol_560,X(I,J)) CALL saxpy(OpenAD_Symbol_561,OpenAD_Symbol_562,X(I,J)) END DO END DO RETURN END SUBROUTINE SUBROUTINE box_model_body() use w2f__types use active_module use all_globals_mod IMPLICIT NONE C C **** Local Variables and functions **** C EXTERNAL box_final_state EXTERNAL box_forward EXTERNAL box_ini_fields INTEGER(w2f__i4) ILEV1 INTEGER(w2f__i4) ILEV2 INTEGER(w2f__i4) ILOOP INTEGER(w2f__i4) ISBYTE PARAMETER ( ISBYTE = 8) INTEGER(w2f__i4) MAXLEV2 INTEGER(w2f__i4) NLEV1 PARAMETER ( NLEV1 = 73) INTEGER(w2f__i4) NLEV2 PARAMETER ( NLEV2 = 50) C C **** statements **** C CALL box_ini_fields() MAXLEV2 = (N_MAX / 73) + 1 IF(N_MAX .GT. INT(3650_w2f__i8)) THEN WRITE(*, *) 'NEED TO SET nlev1*nlev2 >= n_max ' ELSE DO ILEV2 = 1, 50, 1 IF(ILEV2 .LE. MAXLEV2) THEN DO ILEV1 = 1, 73, 1 ILOOP = ILEV1 + ILEV2 * 73 +(-73) IF(ILOOP .LE. N_MAX) THEN CALL box_forward(ILEV1) ENDIF END DO ENDIF END DO CALL box_final_state() ENDIF RETURN END SUBROUTINE SUBROUTINE box_density(TLOC, SLOC, RHOLOC) use w2f__types use active_module use all_globals_mod IMPLICIT NONE C C **** Global Variables **** C REAL(w2f__8) OpenAD_Symbol_284 REAL(w2f__8) OpenAD_Symbol_290 type(active) :: OpenAD_Symbol_285 type(active) :: OpenAD_Symbol_291 REAL(w2f__8) OpenAD_Symbol_286 type(active) :: OpenAD_Symbol_287 REAL(w2f__8) OpenAD_Symbol_288 type(active) :: OpenAD_Symbol_289 REAL(w2f__8) OpenAD_Symbol_0 REAL(w2f__8) OpenAD_Symbol_1 REAL(w2f__8) OpenAD_Symbol_2 REAL(w2f__8) OpenAD_Symbol_3 REAL(w2f__8) OpenAD_Symbol_4 REAL(w2f__8) OpenAD_Symbol_5 REAL(w2f__8) OpenAD_Symbol_6 C C **** Parameters and Result **** C type(active) :: TLOC(1 : 3) type(active) :: SLOC(1 : 3) type(active) :: RHOLOC(1 : 3) C C **** Local Variables and functions **** C INTEGER(w2f__i4) ISBYTE PARAMETER ( ISBYTE = 8) INTEGER(w2f__i4) L INTEGER(w2f__i4) NLEV1 PARAMETER ( NLEV1 = 73) INTEGER(w2f__i4) NLEV2 PARAMETER ( NLEV2 = 50) C C **** statements **** C DO L = 1, 3, 1 OpenAD_Symbol_6 = (SLOC(L)%v*BETA%v-TLOC(L)%v*ALPHA%v) OpenAD_Symbol_2 = BETA%v OpenAD_Symbol_3 = SLOC(L)%v OpenAD_Symbol_0 = 1_w2f__i8 OpenAD_Symbol_4 = ALPHA%v OpenAD_Symbol_5 = TLOC(L)%v OpenAD_Symbol_1 = (-1_w2f__i8) RHOLOC(INT(L))%v = OpenAD_Symbol_6 OpenAD_Symbol_284 = (OpenAD_Symbol_2 * OpenAD_Symbol_0) OpenAD_Symbol_286 = (OpenAD_Symbol_3 * OpenAD_Symbol_0) OpenAD_Symbol_288 = (OpenAD_Symbol_4 * OpenAD_Symbol_1) OpenAD_Symbol_290 = (OpenAD_Symbol_5 * OpenAD_Symbol_1) CALL setderiv(OpenAD_Symbol_291,ALPHA) CALL setderiv(OpenAD_Symbol_289,TLOC(L)) CALL setderiv(OpenAD_Symbol_287,BETA) CALL setderiv(OpenAD_Symbol_285,SLOC(L)) CALL sax(OpenAD_Symbol_284,OpenAD_Symbol_285,RHOLOC(L)) CALL saxpy(OpenAD_Symbol_286,OpenAD_Symbol_287,RHOLOC(L)) CALL saxpy(OpenAD_Symbol_288,OpenAD_Symbol_289,RHOLOC(L)) CALL saxpy(OpenAD_Symbol_290,OpenAD_Symbol_291,RHOLOC(L)) END DO RETURN END SUBROUTINE SUBROUTINE box_transport(RHOLOC, UVELLOC) use w2f__types use active_module use all_globals_mod IMPLICIT NONE C C **** Global Variables **** C REAL(w2f__8) OpenAD_Symbol_300 type(active) :: OpenAD_Symbol_301 REAL(w2f__8) OpenAD_Symbol_302 REAL(w2f__8) OpenAD_Symbol_303 type(active) :: OpenAD_Symbol_304 REAL(w2f__8) OpenAD_Symbol_305 REAL(w2f__8) OpenAD_Symbol_306 type(active) :: OpenAD_Symbol_307 REAL(w2f__8) OpenAD_Symbol_292 type(active) :: OpenAD_Symbol_293 REAL(w2f__8) OpenAD_Symbol_294 REAL(w2f__8) OpenAD_Symbol_295 REAL(w2f__8) OpenAD_Symbol_296 type(active) :: OpenAD_Symbol_297 REAL(w2f__8) OpenAD_Symbol_298 REAL(w2f__8) OpenAD_Symbol_299 REAL(w2f__8) OpenAD_Symbol_10 REAL(w2f__8) OpenAD_Symbol_11 REAL(w2f__8) OpenAD_Symbol_12 REAL(w2f__8) OpenAD_Symbol_13 REAL(w2f__8) OpenAD_Symbol_14 REAL(w2f__8) OpenAD_Symbol_15 REAL(w2f__8) OpenAD_Symbol_16 REAL(w2f__8) OpenAD_Symbol_21 REAL(w2f__8) OpenAD_Symbol_17 REAL(w2f__8) OpenAD_Symbol_22 REAL(w2f__8) OpenAD_Symbol_18 REAL(w2f__8) OpenAD_Symbol_19 REAL(w2f__8) OpenAD_Symbol_7 REAL(w2f__8) OpenAD_Symbol_8 REAL(w2f__8) OpenAD_Symbol_9 C C **** Parameters and Result **** C type(active) :: RHOLOC(1 : 3) type(active) :: UVELLOC C C **** Local Variables and functions **** C INTEGER(w2f__i4) ISBYTE PARAMETER ( ISBYTE = 8) INTEGER(w2f__i4) NLEV1 PARAMETER ( NLEV1 = 73) INTEGER(w2f__i4) NLEV2 PARAMETER ( NLEV2 = 50) C C **** statements **** C OpenAD_Symbol_17 = (1.0D00-DELTA%v) OpenAD_Symbol_8 = (RHOLOC(1)%v*DELTA%v+RHOLOC(3)%v*OpenAD_Symbol_1 +7-RHOLOC(2)%v) OpenAD_Symbol_22 = (-(U0%v*OpenAD_Symbol_8)) OpenAD_Symbol_9 = OpenAD_Symbol_8 OpenAD_Symbol_15 = DELTA%v OpenAD_Symbol_16 = RHOLOC(1)%v OpenAD_Symbol_13 = 1_w2f__i8 OpenAD_Symbol_18 = OpenAD_Symbol_17 OpenAD_Symbol_21 = (-1_w2f__i8) OpenAD_Symbol_19 = RHOLOC(3)%v OpenAD_Symbol_14 = 1_w2f__i8 OpenAD_Symbol_11 = 1_w2f__i8 OpenAD_Symbol_12 = (-1_w2f__i8) OpenAD_Symbol_10 = U0%v OpenAD_Symbol_7 = (-1_w2f__i8) UVELLOC%v = OpenAD_Symbol_22 OpenAD_Symbol_292 = (OpenAD_Symbol_9 * OpenAD_Symbol_7) OpenAD_Symbol_294 = (OpenAD_Symbol_10 * OpenAD_Symbol_7) OpenAD_Symbol_295 = (OpenAD_Symbol_11 * OpenAD_Symbol_294) OpenAD_Symbol_296 = (OpenAD_Symbol_12 * OpenAD_Symbol_294) OpenAD_Symbol_298 = (OpenAD_Symbol_13 * OpenAD_Symbol_295) OpenAD_Symbol_299 = (OpenAD_Symbol_14 * OpenAD_Symbol_295) OpenAD_Symbol_300 = (OpenAD_Symbol_18 * OpenAD_Symbol_299) OpenAD_Symbol_302 = (OpenAD_Symbol_19 * OpenAD_Symbol_299) OpenAD_Symbol_303 = (OpenAD_Symbol_15 * OpenAD_Symbol_298) OpenAD_Symbol_305 = (OpenAD_Symbol_16 * OpenAD_Symbol_298) OpenAD_Symbol_306 = (OpenAD_Symbol_305 + OpenAD_Symbol_21 * > OpenAD_Symbol_302) CALL setderiv(OpenAD_Symbol_307,DELTA) CALL setderiv(OpenAD_Symbol_304,RHOLOC(1)) CALL setderiv(OpenAD_Symbol_301,RHOLOC(3)) CALL setderiv(OpenAD_Symbol_297,RHOLOC(2)) CALL setderiv(OpenAD_Symbol_293,U0) CALL sax(OpenAD_Symbol_292,OpenAD_Symbol_293,UVELLOC) CALL saxpy(OpenAD_Symbol_296,OpenAD_Symbol_297,UVELLOC) CALL saxpy(OpenAD_Symbol_300,OpenAD_Symbol_301,UVELLOC) CALL saxpy(OpenAD_Symbol_303,OpenAD_Symbol_304,UVELLOC) CALL saxpy(OpenAD_Symbol_306,OpenAD_Symbol_307,UVELLOC) RETURN END SUBROUTINE SUBROUTINE box_robert_filter(FLDNOW, FLDOLD, FLDNEW) use w2f__types use active_module use all_globals_mod IMPLICIT NONE C C **** Global Variables **** C REAL(w2f__8) OpenAD_Symbol_130 REAL(w2f__8) OpenAD_Symbol_131 REAL(w2f__8) OpenAD_Symbol_126 REAL(w2f__8) OpenAD_Symbol_132 REAL(w2f__8) OpenAD_Symbol_127 REAL(w2f__8) OpenAD_Symbol_133 REAL(w2f__8) OpenAD_Symbol_128 REAL(w2f__8) OpenAD_Symbol_134 REAL(w2f__8) OpenAD_Symbol_129 REAL(w2f__8) OpenAD_Symbol_135 REAL(w2f__8) OpenAD_Symbol_137 REAL(w2f__8) OpenAD_Symbol_423 type(active) :: OpenAD_Symbol_424 type(active) :: OpenAD_Symbol_430 REAL(w2f__8) OpenAD_Symbol_425 REAL(w2f__8) OpenAD_Symbol_431 REAL(w2f__8) OpenAD_Symbol_426 REAL(w2f__8) OpenAD_Symbol_432 type(active) :: OpenAD_Symbol_427 type(active) :: OpenAD_Symbol_433 REAL(w2f__8) OpenAD_Symbol_428 REAL(w2f__8) OpenAD_Symbol_434 REAL(w2f__8) OpenAD_Symbol_429 type(active) :: OpenAD_Symbol_435 C C **** Parameters and Result **** C type(active) :: FLDNOW(1 : 3) type(active) :: FLDOLD(1 : 3) type(active) :: FLDNEW(1 : 3) C C **** Local Variables and functions **** C INTEGER(w2f__i4) ISBYTE PARAMETER ( ISBYTE = 8) INTEGER(w2f__i4) L INTEGER(w2f__i4) NLEV1 PARAMETER ( NLEV1 = 73) INTEGER(w2f__i4) NLEV2 PARAMETER ( NLEV2 = 50) C C **** statements **** C DO L = 1, 3, 1 OpenAD_Symbol_128 = (FLDOLD(L)%v+FLDNEW(L)%v-FLDNOW(L)%v*2.0D00) OpenAD_Symbol_137 = (FLDNOW(L)%v+ROBERT_FILTER_COEFF%v*OpenAD_Sy +mbol_128) OpenAD_Symbol_126 = 1_w2f__i8 OpenAD_Symbol_129 = OpenAD_Symbol_128 OpenAD_Symbol_131 = 1_w2f__i8 OpenAD_Symbol_133 = 1_w2f__i8 OpenAD_Symbol_135 = 2.0D00 OpenAD_Symbol_134 = (-1_w2f__i8) OpenAD_Symbol_132 = 1_w2f__i8 OpenAD_Symbol_130 = ROBERT_FILTER_COEFF%v OpenAD_Symbol_127 = 1_w2f__i8 FLDNOW(INT(L))%v = OpenAD_Symbol_137 OpenAD_Symbol_423 = (OpenAD_Symbol_129 * OpenAD_Symbol_127) OpenAD_Symbol_425 = (OpenAD_Symbol_130 * OpenAD_Symbol_127) OpenAD_Symbol_426 = (OpenAD_Symbol_131 * OpenAD_Symbol_425) OpenAD_Symbol_428 = (OpenAD_Symbol_132 * OpenAD_Symbol_425) OpenAD_Symbol_429 = (OpenAD_Symbol_133 * OpenAD_Symbol_428) OpenAD_Symbol_431 = (OpenAD_Symbol_134 * OpenAD_Symbol_428) OpenAD_Symbol_432 = (OpenAD_Symbol_135 * OpenAD_Symbol_431) OpenAD_Symbol_434 = OpenAD_Symbol_126 CALL setderiv(OpenAD_Symbol_435,FLDNOW(L)) CALL setderiv(OpenAD_Symbol_433,FLDNOW(L)) CALL setderiv(OpenAD_Symbol_430,FLDNEW(L)) CALL setderiv(OpenAD_Symbol_427,FLDOLD(L)) CALL setderiv(OpenAD_Symbol_424,ROBERT_FILTER_COEFF) CALL sax(OpenAD_Symbol_423,OpenAD_Symbol_424,FLDNOW(L)) CALL saxpy(OpenAD_Symbol_426,OpenAD_Symbol_427,FLDNOW(L)) CALL saxpy(OpenAD_Symbol_429,OpenAD_Symbol_430,FLDNOW(L)) CALL saxpy(OpenAD_Symbol_432,OpenAD_Symbol_433,FLDNOW(L)) CALL saxpy(OpenAD_Symbol_434,OpenAD_Symbol_435,FLDNOW(L)) END DO RETURN END SUBROUTINE SUBROUTINE box_cycle_fields() use w2f__types use active_module use all_globals_mod IMPLICIT NONE C C **** Local Variables and functions **** C INTEGER(w2f__i4) ISBYTE PARAMETER ( ISBYTE = 8) INTEGER(w2f__i4) L INTEGER(w2f__i4) NLEV1 PARAMETER ( NLEV1 = 73) INTEGER(w2f__i4) NLEV2 PARAMETER ( NLEV2 = 50) C C **** statements **** C DO L = 1, 3, 1 TOLD(INT(L))%v = TNOW(L)%v CALL setderiv(TOLD(L),TNOW(L)) TNOW(INT(L))%v = TNEW(L)%v CALL setderiv(TNOW(L),TNEW(L)) SOLD(INT(L))%v = SNOW(L)%v CALL setderiv(SOLD(L),SNOW(L)) SNOW(INT(L))%v = SNEW(L)%v CALL setderiv(SNOW(L),SNEW(L)) END DO RETURN END SUBROUTINE SUBROUTINE box_update(FLDNEW, FLDOLD, DFLDDT) use w2f__types use active_module use all_globals_mod IMPLICIT NONE C C **** Global Variables **** C REAL(w2f__8) OpenAD_Symbol_280 REAL(w2f__8) OpenAD_Symbol_281 REAL(w2f__8) OpenAD_Symbol_276 REAL(w2f__8) OpenAD_Symbol_277 REAL(w2f__8) OpenAD_Symbol_283 REAL(w2f__8) OpenAD_Symbol_278 REAL(w2f__8) OpenAD_Symbol_279 REAL(w2f__8) OpenAD_Symbol_563 REAL(w2f__8) OpenAD_Symbol_564 type(active) :: OpenAD_Symbol_565 REAL(w2f__8) OpenAD_Symbol_566 type(active) :: OpenAD_Symbol_567 REAL(w2f__8) OpenAD_Symbol_568 type(active) :: OpenAD_Symbol_569 C C **** Parameters and Result **** C type(active) :: FLDNEW(1 : 3) type(active) :: FLDOLD(1 : 3) type(active) :: DFLDDT(1 : 3) C C **** Local Variables and functions **** C INTEGER(w2f__i4) ISBYTE PARAMETER ( ISBYTE = 8) INTEGER(w2f__i4) L INTEGER(w2f__i4) NLEV1 PARAMETER ( NLEV1 = 73) INTEGER(w2f__i4) NLEV2 PARAMETER ( NLEV2 = 50) C C **** statements **** C DO L = 1, 3, 1 OpenAD_Symbol_278 = (DELTA_T%v*2.0D00) OpenAD_Symbol_283 = (FLDOLD(L)%v+DFLDDT(L)%v*OpenAD_Symbol_278) OpenAD_Symbol_276 = 1_w2f__i8 OpenAD_Symbol_279 = OpenAD_Symbol_278 OpenAD_Symbol_281 = 2.0D00 OpenAD_Symbol_280 = DFLDDT(L)%v OpenAD_Symbol_277 = 1_w2f__i8 FLDNEW(INT(L))%v = OpenAD_Symbol_283 OpenAD_Symbol_563 = (OpenAD_Symbol_281 * OpenAD_Symbol_280) OpenAD_Symbol_564 = (OpenAD_Symbol_279 * OpenAD_Symbol_277) OpenAD_Symbol_566 = (OpenAD_Symbol_563 * OpenAD_Symbol_277) OpenAD_Symbol_568 = OpenAD_Symbol_276 CALL setderiv(OpenAD_Symbol_569,FLDOLD(L)) CALL setderiv(OpenAD_Symbol_567,DELTA_T) CALL setderiv(OpenAD_Symbol_565,DFLDDT(L)) CALL sax(OpenAD_Symbol_564,OpenAD_Symbol_565,FLDNEW(L)) CALL saxpy(OpenAD_Symbol_566,OpenAD_Symbol_567,FLDNEW(L)) CALL saxpy(OpenAD_Symbol_568,OpenAD_Symbol_569,FLDNEW(L)) END DO RETURN END SUBROUTINE SUBROUTINE box_timestep(GAMMALOC, FLDSTAR, EXTFORLOC, UVELLOC, > FLDNOW, FLDOLD, FLDNEW) use w2f__types use active_module use all_globals_mod IMPLICIT NONE C C **** Global Variables **** C REAL(w2f__8) OpenAD_Symbol_100 REAL(w2f__8) OpenAD_Symbol_101 REAL(w2f__8) OpenAD_Symbol_102 REAL(w2f__8) OpenAD_Symbol_103 REAL(w2f__8) OpenAD_Symbol_104 REAL(w2f__8) OpenAD_Symbol_105 REAL(w2f__8) OpenAD_Symbol_110 REAL(w2f__8) OpenAD_Symbol_106 REAL(w2f__8) OpenAD_Symbol_111 REAL(w2f__8) OpenAD_Symbol_107 REAL(w2f__8) OpenAD_Symbol_112 REAL(w2f__8) OpenAD_Symbol_113 REAL(w2f__8) OpenAD_Symbol_108 REAL(w2f__8) OpenAD_Symbol_114 REAL(w2f__8) OpenAD_Symbol_109 REAL(w2f__8) OpenAD_Symbol_120 REAL(w2f__8) OpenAD_Symbol_115 REAL(w2f__8) OpenAD_Symbol_121 REAL(w2f__8) OpenAD_Symbol_116 REAL(w2f__8) OpenAD_Symbol_122 REAL(w2f__8) OpenAD_Symbol_117 REAL(w2f__8) OpenAD_Symbol_123 REAL(w2f__8) OpenAD_Symbol_118 REAL(w2f__8) OpenAD_Symbol_124 REAL(w2f__8) OpenAD_Symbol_119 REAL(w2f__8) OpenAD_Symbol_125 REAL(w2f__8) OpenAD_Symbol_310 type(active) :: OpenAD_Symbol_311 REAL(w2f__8) OpenAD_Symbol_312 REAL(w2f__8) OpenAD_Symbol_313 REAL(w2f__8) OpenAD_Symbol_308 type(active) :: OpenAD_Symbol_314 REAL(w2f__8) OpenAD_Symbol_309 type(active) :: OpenAD_Symbol_320 REAL(w2f__8) OpenAD_Symbol_315 REAL(w2f__8) OpenAD_Symbol_321 REAL(w2f__8) OpenAD_Symbol_316 REAL(w2f__8) OpenAD_Symbol_322 type(active) :: OpenAD_Symbol_317 type(active) :: OpenAD_Symbol_323 REAL(w2f__8) OpenAD_Symbol_318 REAL(w2f__8) OpenAD_Symbol_324 REAL(w2f__8) OpenAD_Symbol_319 REAL(w2f__8) OpenAD_Symbol_330 type(active) :: OpenAD_Symbol_325 type(active) :: OpenAD_Symbol_331 REAL(w2f__8) OpenAD_Symbol_326 REAL(w2f__8) OpenAD_Symbol_332 type(active) :: OpenAD_Symbol_327 REAL(w2f__8) OpenAD_Symbol_333 REAL(w2f__8) OpenAD_Symbol_328 REAL(w2f__8) OpenAD_Symbol_334 type(active) :: OpenAD_Symbol_329 REAL(w2f__8) OpenAD_Symbol_340 type(active) :: OpenAD_Symbol_335 type(active) :: OpenAD_Symbol_341 REAL(w2f__8) OpenAD_Symbol_336 REAL(w2f__8) OpenAD_Symbol_342 REAL(w2f__8) OpenAD_Symbol_337 REAL(w2f__8) OpenAD_Symbol_343 type(active) :: OpenAD_Symbol_338 type(active) :: OpenAD_Symbol_344 REAL(w2f__8) OpenAD_Symbol_339 REAL(w2f__8) OpenAD_Symbol_400 REAL(w2f__8) OpenAD_Symbol_350 REAL(w2f__8) OpenAD_Symbol_345 type(active) :: OpenAD_Symbol_401 type(active) :: OpenAD_Symbol_351 REAL(w2f__8) OpenAD_Symbol_346 REAL(w2f__8) OpenAD_Symbol_402 REAL(w2f__8) OpenAD_Symbol_352 type(active) :: OpenAD_Symbol_347 REAL(w2f__8) OpenAD_Symbol_403 type(active) :: OpenAD_Symbol_353 REAL(w2f__8) OpenAD_Symbol_348 type(active) :: OpenAD_Symbol_404 REAL(w2f__8) OpenAD_Symbol_354 type(active) :: OpenAD_Symbol_349 type(active) :: OpenAD_Symbol_410 REAL(w2f__8) OpenAD_Symbol_405 type(active) :: OpenAD_Symbol_360 type(active) :: OpenAD_Symbol_355 REAL(w2f__8) OpenAD_Symbol_411 type(active) :: OpenAD_Symbol_406 REAL(w2f__8) OpenAD_Symbol_361 REAL(w2f__8) OpenAD_Symbol_356 type(active) :: OpenAD_Symbol_412 REAL(w2f__8) OpenAD_Symbol_407 type(active) :: OpenAD_Symbol_362 type(active) :: OpenAD_Symbol_357 REAL(w2f__8) OpenAD_Symbol_413 type(active) :: OpenAD_Symbol_408 REAL(w2f__8) OpenAD_Symbol_363 REAL(w2f__8) OpenAD_Symbol_358 REAL(w2f__8) OpenAD_Symbol_414 REAL(w2f__8) OpenAD_Symbol_409 type(active) :: OpenAD_Symbol_364 REAL(w2f__8) OpenAD_Symbol_359 type(active) :: OpenAD_Symbol_420 type(active) :: OpenAD_Symbol_415 REAL(w2f__8) OpenAD_Symbol_370 REAL(w2f__8) OpenAD_Symbol_365 REAL(w2f__8) OpenAD_Symbol_421 REAL(w2f__8) OpenAD_Symbol_416 type(active) :: OpenAD_Symbol_371 REAL(w2f__8) OpenAD_Symbol_366 type(active) :: OpenAD_Symbol_422 type(active) :: OpenAD_Symbol_417 REAL(w2f__8) OpenAD_Symbol_372 REAL(w2f__8) OpenAD_Symbol_367 REAL(w2f__8) OpenAD_Symbol_418 REAL(w2f__8) OpenAD_Symbol_373 type(active) :: OpenAD_Symbol_368 REAL(w2f__8) OpenAD_Symbol_419 type(active) :: OpenAD_Symbol_374 REAL(w2f__8) OpenAD_Symbol_369 type(active) :: OpenAD_Symbol_380 REAL(w2f__8) OpenAD_Symbol_375 REAL(w2f__8) OpenAD_Symbol_381 REAL(w2f__8) OpenAD_Symbol_376 type(active) :: OpenAD_Symbol_382 type(active) :: OpenAD_Symbol_377 REAL(w2f__8) OpenAD_Symbol_383 REAL(w2f__8) OpenAD_Symbol_378 type(active) :: OpenAD_Symbol_384 REAL(w2f__8) OpenAD_Symbol_379 REAL(w2f__8) OpenAD_Symbol_390 REAL(w2f__8) OpenAD_Symbol_385 REAL(w2f__8) OpenAD_Symbol_391 type(active) :: OpenAD_Symbol_386 type(active) :: OpenAD_Symbol_392 REAL(w2f__8) OpenAD_Symbol_387 REAL(w2f__8) OpenAD_Symbol_393 type(active) :: OpenAD_Symbol_388 REAL(w2f__8) OpenAD_Symbol_394 REAL(w2f__8) OpenAD_Symbol_389 type(active) :: OpenAD_Symbol_395 REAL(w2f__8) OpenAD_Symbol_396 REAL(w2f__8) OpenAD_Symbol_397 type(active) :: OpenAD_Symbol_398 REAL(w2f__8) OpenAD_Symbol_399 REAL(w2f__8) OpenAD_Symbol_23 REAL(w2f__8) OpenAD_Symbol_24 REAL(w2f__8) OpenAD_Symbol_25 REAL(w2f__8) OpenAD_Symbol_30 REAL(w2f__8) OpenAD_Symbol_26 REAL(w2f__8) OpenAD_Symbol_31 REAL(w2f__8) OpenAD_Symbol_27 REAL(w2f__8) OpenAD_Symbol_32 REAL(w2f__8) OpenAD_Symbol_28 REAL(w2f__8) OpenAD_Symbol_33 REAL(w2f__8) OpenAD_Symbol_29 REAL(w2f__8) OpenAD_Symbol_34 REAL(w2f__8) OpenAD_Symbol_35 REAL(w2f__8) OpenAD_Symbol_40 REAL(w2f__8) OpenAD_Symbol_36 REAL(w2f__8) OpenAD_Symbol_41 REAL(w2f__8) OpenAD_Symbol_37 REAL(w2f__8) OpenAD_Symbol_42 REAL(w2f__8) OpenAD_Symbol_38 REAL(w2f__8) OpenAD_Symbol_43 REAL(w2f__8) OpenAD_Symbol_39 REAL(w2f__8) OpenAD_Symbol_44 REAL(w2f__8) OpenAD_Symbol_45 REAL(w2f__8) OpenAD_Symbol_50 REAL(w2f__8) OpenAD_Symbol_46 REAL(w2f__8) OpenAD_Symbol_51 REAL(w2f__8) OpenAD_Symbol_47 REAL(w2f__8) OpenAD_Symbol_52 REAL(w2f__8) OpenAD_Symbol_48 REAL(w2f__8) OpenAD_Symbol_53 REAL(w2f__8) OpenAD_Symbol_49 REAL(w2f__8) OpenAD_Symbol_54 REAL(w2f__8) OpenAD_Symbol_55 REAL(w2f__8) OpenAD_Symbol_60 REAL(w2f__8) OpenAD_Symbol_56 REAL(w2f__8) OpenAD_Symbol_61 REAL(w2f__8) OpenAD_Symbol_57 REAL(w2f__8) OpenAD_Symbol_62 REAL(w2f__8) OpenAD_Symbol_58 REAL(w2f__8) OpenAD_Symbol_63 REAL(w2f__8) OpenAD_Symbol_59 REAL(w2f__8) OpenAD_Symbol_64 REAL(w2f__8) OpenAD_Symbol_65 REAL(w2f__8) OpenAD_Symbol_70 REAL(w2f__8) OpenAD_Symbol_66 REAL(w2f__8) OpenAD_Symbol_71 REAL(w2f__8) OpenAD_Symbol_67 REAL(w2f__8) OpenAD_Symbol_72 REAL(w2f__8) OpenAD_Symbol_68 REAL(w2f__8) OpenAD_Symbol_73 REAL(w2f__8) OpenAD_Symbol_69 REAL(w2f__8) OpenAD_Symbol_74 REAL(w2f__8) OpenAD_Symbol_75 REAL(w2f__8) OpenAD_Symbol_80 REAL(w2f__8) OpenAD_Symbol_76 REAL(w2f__8) OpenAD_Symbol_81 REAL(w2f__8) OpenAD_Symbol_77 REAL(w2f__8) OpenAD_Symbol_82 REAL(w2f__8) OpenAD_Symbol_78 REAL(w2f__8) OpenAD_Symbol_83 REAL(w2f__8) OpenAD_Symbol_79 REAL(w2f__8) OpenAD_Symbol_84 REAL(w2f__8) OpenAD_Symbol_85 REAL(w2f__8) OpenAD_Symbol_90 REAL(w2f__8) OpenAD_Symbol_86 REAL(w2f__8) OpenAD_Symbol_91 REAL(w2f__8) OpenAD_Symbol_87 REAL(w2f__8) OpenAD_Symbol_92 REAL(w2f__8) OpenAD_Symbol_88 REAL(w2f__8) OpenAD_Symbol_93 REAL(w2f__8) OpenAD_Symbol_89 REAL(w2f__8) OpenAD_Symbol_94 REAL(w2f__8) OpenAD_Symbol_95 REAL(w2f__8) OpenAD_Symbol_96 REAL(w2f__8) OpenAD_Symbol_97 REAL(w2f__8) OpenAD_Symbol_98 REAL(w2f__8) OpenAD_Symbol_99 C C **** Parameters and Result **** C type(active) :: GAMMALOC type(active) :: FLDSTAR(1 : 2) type(active) :: EXTFORLOC(1 : 2) type(active) :: UVELLOC type(active) :: FLDNOW(1 : 3) type(active) :: FLDOLD(1 : 3) type(active) :: FLDNEW(1 : 3) C C **** Local Variables and functions **** C EXTERNAL box_update type(active) :: DFLDDT(1 : 3) INTEGER(w2f__i4) ISBYTE PARAMETER ( ISBYTE = 8) INTEGER(w2f__i4) NLEV1 PARAMETER ( NLEV1 = 73) INTEGER(w2f__i4) NLEV2 PARAMETER ( NLEV2 = 50) C C **** statements **** C IF (UVELLOC%v .GE. 0.0D00) THEN OpenAD_Symbol_33 = (FLDSTAR(1)%v-FLDNOW(1)%v) OpenAD_Symbol_30 = (GAMMALOC%v*OpenAD_Symbol_33) OpenAD_Symbol_38 = (FLDNOW(3)%v-FLDNOW(1)%v) OpenAD_Symbol_23 = (EXTFORLOC(1)%v+VOL(1)%v*OpenAD_Symbol_30+UVE +LLOC%v*OpenAD_Symbol_38) OpenAD_Symbol_43 = (OpenAD_Symbol_23/VOL(1)%v) OpenAD_Symbol_28 = 1_w2f__i8 OpenAD_Symbol_31 = OpenAD_Symbol_30 OpenAD_Symbol_34 = OpenAD_Symbol_33 OpenAD_Symbol_36 = 1_w2f__i8 OpenAD_Symbol_37 = (-1_w2f__i8) OpenAD_Symbol_35 = GAMMALOC%v OpenAD_Symbol_32 = VOL(1)%v OpenAD_Symbol_29 = 1_w2f__i8 OpenAD_Symbol_26 = 1_w2f__i8 OpenAD_Symbol_39 = OpenAD_Symbol_38 OpenAD_Symbol_41 = 1_w2f__i8 OpenAD_Symbol_42 = (-1_w2f__i8) OpenAD_Symbol_40 = UVELLOC%v OpenAD_Symbol_27 = 1_w2f__i8 OpenAD_Symbol_24 = (INT(1 _w2f__i8)/VOL(1)%v) OpenAD_Symbol_25 = (-(OpenAD_Symbol_23/(VOL(1)%v*VOL(1)%v))) DFLDDT(1)%v = OpenAD_Symbol_43 OpenAD_Symbol_308 = (OpenAD_Symbol_26 * OpenAD_Symbol_24) OpenAD_Symbol_309 = (OpenAD_Symbol_27 * OpenAD_Symbol_24) OpenAD_Symbol_310 = (OpenAD_Symbol_28 * OpenAD_Symbol_308) OpenAD_Symbol_312 = (OpenAD_Symbol_29 * OpenAD_Symbol_308) OpenAD_Symbol_313 = (OpenAD_Symbol_31 * OpenAD_Symbol_312) OpenAD_Symbol_315 = (OpenAD_Symbol_32 * OpenAD_Symbol_312) OpenAD_Symbol_316 = (OpenAD_Symbol_39 * OpenAD_Symbol_309) OpenAD_Symbol_318 = (OpenAD_Symbol_40 * OpenAD_Symbol_309) OpenAD_Symbol_319 = (OpenAD_Symbol_34 * OpenAD_Symbol_315) OpenAD_Symbol_321 = (OpenAD_Symbol_35 * OpenAD_Symbol_315) OpenAD_Symbol_322 = (OpenAD_Symbol_41 * OpenAD_Symbol_318) OpenAD_Symbol_324 = (OpenAD_Symbol_42 * OpenAD_Symbol_318) OpenAD_Symbol_326 = (OpenAD_Symbol_36 * OpenAD_Symbol_321) OpenAD_Symbol_328 = (OpenAD_Symbol_37 * OpenAD_Symbol_321) OpenAD_Symbol_330 = OpenAD_Symbol_25 CALL setderiv(OpenAD_Symbol_331,VOL(1)) CALL setderiv(OpenAD_Symbol_329,FLDNOW(1)) CALL setderiv(OpenAD_Symbol_327,FLDSTAR(1)) CALL setderiv(OpenAD_Symbol_325,FLDNOW(1)) CALL setderiv(OpenAD_Symbol_323,FLDNOW(3)) CALL setderiv(OpenAD_Symbol_320,GAMMALOC) CALL setderiv(OpenAD_Symbol_317,UVELLOC) CALL setderiv(OpenAD_Symbol_314,VOL(1)) CALL setderiv(OpenAD_Symbol_311,EXTFORLOC(1)) CALL sax(OpenAD_Symbol_310,OpenAD_Symbol_311,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_313,OpenAD_Symbol_314,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_316,OpenAD_Symbol_317,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_319,OpenAD_Symbol_320,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_322,OpenAD_Symbol_323,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_324,OpenAD_Symbol_325,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_326,OpenAD_Symbol_327,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_328,OpenAD_Symbol_329,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_330,OpenAD_Symbol_331,DFLDDT(1)) OpenAD_Symbol_54 = (FLDSTAR(2)%v-FLDNOW(2)%v) OpenAD_Symbol_51 = (GAMMALOC%v*OpenAD_Symbol_54) OpenAD_Symbol_59 = (FLDNOW(1)%v-FLDNOW(2)%v) OpenAD_Symbol_44 = (EXTFORLOC(2)%v+VOL(2)%v*OpenAD_Symbol_51+UVE +LLOC%v*OpenAD_Symbol_59) OpenAD_Symbol_64 = (OpenAD_Symbol_44/VOL(2)%v) OpenAD_Symbol_49 = 1_w2f__i8 OpenAD_Symbol_52 = OpenAD_Symbol_51 OpenAD_Symbol_55 = OpenAD_Symbol_54 OpenAD_Symbol_57 = 1_w2f__i8 OpenAD_Symbol_58 = (-1_w2f__i8) OpenAD_Symbol_56 = GAMMALOC%v OpenAD_Symbol_53 = VOL(2)%v OpenAD_Symbol_50 = 1_w2f__i8 OpenAD_Symbol_47 = 1_w2f__i8 OpenAD_Symbol_60 = OpenAD_Symbol_59 OpenAD_Symbol_62 = 1_w2f__i8 OpenAD_Symbol_63 = (-1_w2f__i8) OpenAD_Symbol_61 = UVELLOC%v OpenAD_Symbol_48 = 1_w2f__i8 OpenAD_Symbol_45 = (INT(1 _w2f__i8)/VOL(2)%v) OpenAD_Symbol_46 = (-(OpenAD_Symbol_44/(VOL(2)%v*VOL(2)%v))) DFLDDT(2)%v = OpenAD_Symbol_64 OpenAD_Symbol_332 = (OpenAD_Symbol_47 * OpenAD_Symbol_45) OpenAD_Symbol_333 = (OpenAD_Symbol_48 * OpenAD_Symbol_45) OpenAD_Symbol_334 = (OpenAD_Symbol_49 * OpenAD_Symbol_332) OpenAD_Symbol_336 = (OpenAD_Symbol_50 * OpenAD_Symbol_332) OpenAD_Symbol_337 = (OpenAD_Symbol_52 * OpenAD_Symbol_336) OpenAD_Symbol_339 = (OpenAD_Symbol_53 * OpenAD_Symbol_336) OpenAD_Symbol_340 = (OpenAD_Symbol_60 * OpenAD_Symbol_333) OpenAD_Symbol_342 = (OpenAD_Symbol_61 * OpenAD_Symbol_333) OpenAD_Symbol_343 = (OpenAD_Symbol_55 * OpenAD_Symbol_339) OpenAD_Symbol_345 = (OpenAD_Symbol_56 * OpenAD_Symbol_339) OpenAD_Symbol_346 = (OpenAD_Symbol_62 * OpenAD_Symbol_342) OpenAD_Symbol_348 = (OpenAD_Symbol_63 * OpenAD_Symbol_342) OpenAD_Symbol_350 = (OpenAD_Symbol_57 * OpenAD_Symbol_345) OpenAD_Symbol_352 = (OpenAD_Symbol_58 * OpenAD_Symbol_345) OpenAD_Symbol_354 = OpenAD_Symbol_46 CALL setderiv(OpenAD_Symbol_355,VOL(2)) CALL setderiv(OpenAD_Symbol_353,FLDNOW(2)) CALL setderiv(OpenAD_Symbol_351,FLDSTAR(2)) CALL setderiv(OpenAD_Symbol_349,FLDNOW(2)) CALL setderiv(OpenAD_Symbol_347,FLDNOW(1)) CALL setderiv(OpenAD_Symbol_344,GAMMALOC) CALL setderiv(OpenAD_Symbol_341,UVELLOC) CALL setderiv(OpenAD_Symbol_338,VOL(2)) CALL setderiv(OpenAD_Symbol_335,EXTFORLOC(2)) CALL sax(OpenAD_Symbol_334,OpenAD_Symbol_335,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_337,OpenAD_Symbol_338,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_340,OpenAD_Symbol_341,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_343,OpenAD_Symbol_344,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_346,OpenAD_Symbol_347,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_348,OpenAD_Symbol_349,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_350,OpenAD_Symbol_351,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_352,OpenAD_Symbol_353,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_354,OpenAD_Symbol_355,DFLDDT(2)) OpenAD_Symbol_68 = (FLDNOW(2)%v-FLDNOW(3)%v) OpenAD_Symbol_65 = (UVELLOC%v*OpenAD_Symbol_68) OpenAD_Symbol_73 = (OpenAD_Symbol_65/VOL(3)%v) OpenAD_Symbol_69 = OpenAD_Symbol_68 OpenAD_Symbol_71 = 1_w2f__i8 OpenAD_Symbol_72 = (-1_w2f__i8) OpenAD_Symbol_70 = UVELLOC%v OpenAD_Symbol_66 = (INT(1 _w2f__i8)/VOL(3)%v) OpenAD_Symbol_67 = (-(OpenAD_Symbol_65/(VOL(3)%v*VOL(3)%v))) DFLDDT(3)%v = OpenAD_Symbol_73 OpenAD_Symbol_356 = (OpenAD_Symbol_69 * OpenAD_Symbol_66) OpenAD_Symbol_358 = (OpenAD_Symbol_70 * OpenAD_Symbol_66) OpenAD_Symbol_359 = (OpenAD_Symbol_71 * OpenAD_Symbol_358) OpenAD_Symbol_361 = (OpenAD_Symbol_72 * OpenAD_Symbol_358) OpenAD_Symbol_363 = OpenAD_Symbol_67 CALL setderiv(OpenAD_Symbol_364,VOL(3)) CALL setderiv(OpenAD_Symbol_362,FLDNOW(3)) CALL setderiv(OpenAD_Symbol_360,FLDNOW(2)) CALL setderiv(OpenAD_Symbol_357,UVELLOC) CALL sax(OpenAD_Symbol_356,OpenAD_Symbol_357,DFLDDT(3)) CALL saxpy(OpenAD_Symbol_359,OpenAD_Symbol_360,DFLDDT(3)) CALL saxpy(OpenAD_Symbol_361,OpenAD_Symbol_362,DFLDDT(3)) CALL saxpy(OpenAD_Symbol_363,OpenAD_Symbol_364,DFLDDT(3)) ELSE OpenAD_Symbol_84 = (FLDSTAR(1)%v-FLDNOW(1)%v) OpenAD_Symbol_81 = (GAMMALOC%v*OpenAD_Symbol_84) OpenAD_Symbol_89 = (FLDNOW(2)%v-FLDNOW(1)%v) OpenAD_Symbol_74 = (EXTFORLOC(1)%v+VOL(1)%v*OpenAD_Symbol_81-UVE +LLOC%v*OpenAD_Symbol_89) OpenAD_Symbol_94 = (OpenAD_Symbol_74/VOL(1)%v) OpenAD_Symbol_79 = 1_w2f__i8 OpenAD_Symbol_82 = OpenAD_Symbol_81 OpenAD_Symbol_85 = OpenAD_Symbol_84 OpenAD_Symbol_87 = 1_w2f__i8 OpenAD_Symbol_88 = (-1_w2f__i8) OpenAD_Symbol_86 = GAMMALOC%v OpenAD_Symbol_83 = VOL(1)%v OpenAD_Symbol_80 = 1_w2f__i8 OpenAD_Symbol_77 = 1_w2f__i8 OpenAD_Symbol_90 = OpenAD_Symbol_89 OpenAD_Symbol_92 = 1_w2f__i8 OpenAD_Symbol_93 = (-1_w2f__i8) OpenAD_Symbol_91 = UVELLOC%v OpenAD_Symbol_78 = (-1_w2f__i8) OpenAD_Symbol_75 = (INT(1 _w2f__i8)/VOL(1)%v) OpenAD_Symbol_76 = (-(OpenAD_Symbol_74/(VOL(1)%v*VOL(1)%v))) DFLDDT(1)%v = OpenAD_Symbol_94 OpenAD_Symbol_365 = (OpenAD_Symbol_77 * OpenAD_Symbol_75) OpenAD_Symbol_366 = (OpenAD_Symbol_78 * OpenAD_Symbol_75) OpenAD_Symbol_367 = (OpenAD_Symbol_79 * OpenAD_Symbol_365) OpenAD_Symbol_369 = (OpenAD_Symbol_80 * OpenAD_Symbol_365) OpenAD_Symbol_370 = (OpenAD_Symbol_82 * OpenAD_Symbol_369) OpenAD_Symbol_372 = (OpenAD_Symbol_83 * OpenAD_Symbol_369) OpenAD_Symbol_373 = (OpenAD_Symbol_90 * OpenAD_Symbol_366) OpenAD_Symbol_375 = (OpenAD_Symbol_91 * OpenAD_Symbol_366) OpenAD_Symbol_376 = (OpenAD_Symbol_85 * OpenAD_Symbol_372) OpenAD_Symbol_378 = (OpenAD_Symbol_86 * OpenAD_Symbol_372) OpenAD_Symbol_379 = (OpenAD_Symbol_92 * OpenAD_Symbol_375) OpenAD_Symbol_381 = (OpenAD_Symbol_93 * OpenAD_Symbol_375) OpenAD_Symbol_383 = (OpenAD_Symbol_87 * OpenAD_Symbol_378) OpenAD_Symbol_385 = (OpenAD_Symbol_88 * OpenAD_Symbol_378) OpenAD_Symbol_387 = OpenAD_Symbol_76 CALL setderiv(OpenAD_Symbol_388,VOL(1)) CALL setderiv(OpenAD_Symbol_386,FLDNOW(1)) CALL setderiv(OpenAD_Symbol_384,FLDSTAR(1)) CALL setderiv(OpenAD_Symbol_382,FLDNOW(1)) CALL setderiv(OpenAD_Symbol_380,FLDNOW(2)) CALL setderiv(OpenAD_Symbol_377,GAMMALOC) CALL setderiv(OpenAD_Symbol_374,UVELLOC) CALL setderiv(OpenAD_Symbol_371,VOL(1)) CALL setderiv(OpenAD_Symbol_368,EXTFORLOC(1)) CALL sax(OpenAD_Symbol_367,OpenAD_Symbol_368,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_370,OpenAD_Symbol_371,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_373,OpenAD_Symbol_374,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_376,OpenAD_Symbol_377,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_379,OpenAD_Symbol_380,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_381,OpenAD_Symbol_382,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_383,OpenAD_Symbol_384,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_385,OpenAD_Symbol_386,DFLDDT(1)) CALL saxpy(OpenAD_Symbol_387,OpenAD_Symbol_388,DFLDDT(1)) OpenAD_Symbol_105 = (FLDSTAR(2)%v-FLDNOW(2)%v) OpenAD_Symbol_102 = (GAMMALOC%v*OpenAD_Symbol_105) OpenAD_Symbol_110 = (FLDNOW(3)%v-FLDNOW(2)%v) OpenAD_Symbol_95 = (EXTFORLOC(2)%v+VOL(2)%v*OpenAD_Symbol_102-UV +ELLOC%v*OpenAD_Symbol_110) OpenAD_Symbol_115 = (OpenAD_Symbol_95/VOL(2)%v) OpenAD_Symbol_100 = 1_w2f__i8 OpenAD_Symbol_103 = OpenAD_Symbol_102 OpenAD_Symbol_106 = OpenAD_Symbol_105 OpenAD_Symbol_108 = 1_w2f__i8 OpenAD_Symbol_109 = (-1_w2f__i8) OpenAD_Symbol_107 = GAMMALOC%v OpenAD_Symbol_104 = VOL(2)%v OpenAD_Symbol_101 = 1_w2f__i8 OpenAD_Symbol_98 = 1_w2f__i8 OpenAD_Symbol_111 = OpenAD_Symbol_110 OpenAD_Symbol_113 = 1_w2f__i8 OpenAD_Symbol_114 = (-1_w2f__i8) OpenAD_Symbol_112 = UVELLOC%v OpenAD_Symbol_99 = (-1_w2f__i8) OpenAD_Symbol_96 = (INT(1 _w2f__i8)/VOL(2)%v) OpenAD_Symbol_97 = (-(OpenAD_Symbol_95/(VOL(2)%v*VOL(2)%v))) DFLDDT(2)%v = OpenAD_Symbol_115 OpenAD_Symbol_389 = (OpenAD_Symbol_98 * OpenAD_Symbol_96) OpenAD_Symbol_390 = (OpenAD_Symbol_99 * OpenAD_Symbol_96) OpenAD_Symbol_391 = (OpenAD_Symbol_100 * OpenAD_Symbol_389) OpenAD_Symbol_393 = (OpenAD_Symbol_101 * OpenAD_Symbol_389) OpenAD_Symbol_394 = (OpenAD_Symbol_103 * OpenAD_Symbol_393) OpenAD_Symbol_396 = (OpenAD_Symbol_104 * OpenAD_Symbol_393) OpenAD_Symbol_397 = (OpenAD_Symbol_111 * OpenAD_Symbol_390) OpenAD_Symbol_399 = (OpenAD_Symbol_112 * OpenAD_Symbol_390) OpenAD_Symbol_400 = (OpenAD_Symbol_106 * OpenAD_Symbol_396) OpenAD_Symbol_402 = (OpenAD_Symbol_107 * OpenAD_Symbol_396) OpenAD_Symbol_403 = (OpenAD_Symbol_113 * OpenAD_Symbol_399) OpenAD_Symbol_405 = (OpenAD_Symbol_114 * OpenAD_Symbol_399) OpenAD_Symbol_407 = (OpenAD_Symbol_108 * OpenAD_Symbol_402) OpenAD_Symbol_409 = (OpenAD_Symbol_109 * OpenAD_Symbol_402) OpenAD_Symbol_411 = OpenAD_Symbol_97 CALL setderiv(OpenAD_Symbol_412,VOL(2)) CALL setderiv(OpenAD_Symbol_410,FLDNOW(2)) CALL setderiv(OpenAD_Symbol_408,FLDSTAR(2)) CALL setderiv(OpenAD_Symbol_406,FLDNOW(2)) CALL setderiv(OpenAD_Symbol_404,FLDNOW(3)) CALL setderiv(OpenAD_Symbol_401,GAMMALOC) CALL setderiv(OpenAD_Symbol_398,UVELLOC) CALL setderiv(OpenAD_Symbol_395,VOL(2)) CALL setderiv(OpenAD_Symbol_392,EXTFORLOC(2)) CALL sax(OpenAD_Symbol_391,OpenAD_Symbol_392,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_394,OpenAD_Symbol_395,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_397,OpenAD_Symbol_398,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_400,OpenAD_Symbol_401,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_403,OpenAD_Symbol_404,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_405,OpenAD_Symbol_406,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_407,OpenAD_Symbol_408,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_409,OpenAD_Symbol_410,DFLDDT(2)) CALL saxpy(OpenAD_Symbol_411,OpenAD_Symbol_412,DFLDDT(2)) OpenAD_Symbol_120 = (FLDNOW(1)%v-FLDNOW(3)%v) OpenAD_Symbol_117 = (UVELLOC%v*OpenAD_Symbol_120) OpenAD_Symbol_125 = (-(OpenAD_Symbol_117/VOL(3)%v)) OpenAD_Symbol_121 = OpenAD_Symbol_120 OpenAD_Symbol_123 = 1_w2f__i8 OpenAD_Symbol_124 = (-1_w2f__i8) OpenAD_Symbol_122 = UVELLOC%v OpenAD_Symbol_118 = (INT(1 _w2f__i8)/VOL(3)%v) OpenAD_Symbol_119 = (-(OpenAD_Symbol_117/(VOL(3)%v*VOL(3)%v))) OpenAD_Symbol_116 = (-1_w2f__i8) DFLDDT(3)%v = OpenAD_Symbol_125 OpenAD_Symbol_413 = (OpenAD_Symbol_118 * OpenAD_Symbol_116) OpenAD_Symbol_414 = (OpenAD_Symbol_119 * OpenAD_Symbol_116) OpenAD_Symbol_416 = (OpenAD_Symbol_121 * OpenAD_Symbol_413) OpenAD_Symbol_418 = (OpenAD_Symbol_122 * OpenAD_Symbol_413) OpenAD_Symbol_419 = (OpenAD_Symbol_123 * OpenAD_Symbol_418) OpenAD_Symbol_421 = (OpenAD_Symbol_124 * OpenAD_Symbol_418) CALL setderiv(OpenAD_Symbol_422,FLDNOW(3)) CALL setderiv(OpenAD_Symbol_420,FLDNOW(1)) CALL setderiv(OpenAD_Symbol_417,UVELLOC) CALL setderiv(OpenAD_Symbol_415,VOL(3)) CALL sax(OpenAD_Symbol_414,OpenAD_Symbol_415,DFLDDT(3)) CALL saxpy(OpenAD_Symbol_416,OpenAD_Symbol_417,DFLDDT(3)) CALL saxpy(OpenAD_Symbol_419,OpenAD_Symbol_420,DFLDDT(3)) CALL saxpy(OpenAD_Symbol_421,OpenAD_Symbol_422,DFLDDT(3)) ENDIF CALL box_update(FLDNEW,FLDOLD,DFLDDT) RETURN END SUBROUTINE