C *********************************************************** C Fortran file translated from WHIRL Wed Sep 29 14:46:31 2004 C *********************************************************** C *********************************************************** SUBROUTINE head(NX, NY, X, FVEC, R) use w2f__types use active_module IMPLICIT NONE C C **** Global Variables **** C REAL(w2f__8) OpenAD_Symbol_4 REAL(w2f__8) OpenAD_Symbol_5 REAL(w2f__8) OpenAD_Symbol_6 REAL(w2f__8) OpenAD_Symbol_7 REAL(w2f__8) OpenAD_Symbol_8 REAL(w2f__8) OpenAD_Symbol_9 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_17 REAL(w2f__8) OpenAD_Symbol_200 REAL(w2f__8) OpenAD_Symbol_201 REAL(w2f__8) OpenAD_Symbol_202 type(active) :: OpenAD_Symbol_203 REAL(w2f__8) OpenAD_Symbol_204 type(active) :: OpenAD_Symbol_205 REAL(w2f__8) OpenAD_Symbol_206 type(active) :: OpenAD_Symbol_207 REAL(w2f__8) OpenAD_Symbol_196 type(active) :: OpenAD_Symbol_197 REAL(w2f__8) OpenAD_Symbol_198 type(active) :: OpenAD_Symbol_199 REAL(w2f__8) OpenAD_Symbol_1 REAL(w2f__8) OpenAD_Symbol_2 C C **** Parameters and Result **** C INTEGER(w2f__i4) NX INTEGER(w2f__i4) NY type(active) :: X(1 : INT((NX * NY))) type(active) :: FVEC(1 : INT((NX * NY))) type(active) :: R C C **** Local Variables and functions **** C INTEGER(w2f__i8) t__1 INTEGER(w2f__i8) t__2 EXTERNAL foo REAL(w2f__8) FOUR PARAMETER ( FOUR = 4.000000000000000) type(active) :: HX type(active) :: HX2 type(active) :: HY type(active) :: HY2 INTEGER(w2f__i4) K INTEGER(w2f__i4) N type(active) :: NXP1 type(active) :: NYP1 REAL(w2f__8) ONE PARAMETER ( ONE = 1.000000000000000) REAL(w2f__8) THREE PARAMETER ( THREE = 3.000000000000000) REAL(w2f__8) TWO PARAMETER ( TWO = 2.000000000000000) REAL(w2f__8) ZERO PARAMETER ( ZERO = 0.000000000000000) C C **** Temporary variables **** C INTEGER(w2f__i4) doloop_ub C C **** statements **** C t__1 = INT(NX * NY) t__2 = MAX(t__1, 0_w2f__i8) N = NX * NY NXP1%v = (NX+INT(1 _w2f__i8)) NYP1%v = (NY+INT(1 _w2f__i8)) CALL zero_deriv(NXP1) CALL zero_deriv(NYP1) OpenAD_Symbol_2 = (1.0D00/NXP1%v) OpenAD_Symbol_1 = (-(1.0D00/(NXP1%v*NXP1%v))) HX%v = OpenAD_Symbol_2 OpenAD_Symbol_5 = (1.0D00/NYP1%v) OpenAD_Symbol_4 = (-(1.0D00/(NYP1%v*NYP1%v))) HY%v = OpenAD_Symbol_5 OpenAD_Symbol_8 = (HY%v*HY%v) OpenAD_Symbol_6 = HY%v OpenAD_Symbol_7 = HY%v HY2%v = OpenAD_Symbol_8 OpenAD_Symbol_11 = (HX%v*HX%v) OpenAD_Symbol_9 = HX%v OpenAD_Symbol_10 = HX%v HX2%v = OpenAD_Symbol_11 OpenAD_Symbol_196 = (OpenAD_Symbol_1 *(OpenAD_Symbol_10 + > OpenAD_Symbol_9)) OpenAD_Symbol_198 = (OpenAD_Symbol_4 *(OpenAD_Symbol_7 + > OpenAD_Symbol_6)) OpenAD_Symbol_200 = OpenAD_Symbol_1 OpenAD_Symbol_201 = OpenAD_Symbol_4 CALL setderiv(OpenAD_Symbol_199,NYP1) CALL setderiv(OpenAD_Symbol_197,NXP1) CALL sax(OpenAD_Symbol_196,OpenAD_Symbol_197,HX2) CALL sax(OpenAD_Symbol_198,OpenAD_Symbol_199,HY2) CALL sax(OpenAD_Symbol_200,OpenAD_Symbol_197,HX) CALL sax(OpenAD_Symbol_201,OpenAD_Symbol_199,HY) CALL foo(NX,NY,X,FVEC,R,N,HX,HY,HY2,HX2) doloop_ub = N DO K = 1, doloop_ub, 1 OpenAD_Symbol_12 = (HX2%v*HY2%v) OpenAD_Symbol_17 = (FVEC(K)%v*OpenAD_Symbol_12) OpenAD_Symbol_13 = OpenAD_Symbol_12 OpenAD_Symbol_15 = HY2%v OpenAD_Symbol_16 = HX2%v OpenAD_Symbol_14 = FVEC(K)%v FVEC(INT(K))%v = OpenAD_Symbol_17 OpenAD_Symbol_202 = (OpenAD_Symbol_15 * OpenAD_Symbol_14) OpenAD_Symbol_204 = (OpenAD_Symbol_16 * OpenAD_Symbol_14) OpenAD_Symbol_206 = OpenAD_Symbol_13 CALL setderiv(OpenAD_Symbol_207,FVEC(K)) CALL setderiv(OpenAD_Symbol_205,HY2) CALL setderiv(OpenAD_Symbol_203,HX2) CALL sax(OpenAD_Symbol_202,OpenAD_Symbol_203,FVEC(K)) CALL saxpy(OpenAD_Symbol_204,OpenAD_Symbol_205,FVEC(K)) CALL saxpy(OpenAD_Symbol_206,OpenAD_Symbol_207,FVEC(K)) END DO RETURN END SUBROUTINE SUBROUTINE foo(NX, NY, X, FVEC, R, N, HX, HY, HY2, HX2) use w2f__types use active_module IMPLICIT NONE C C **** Global Variables **** C REAL(w2f__8) OpenAD_Symbol_20 REAL(w2f__8) OpenAD_Symbol_21 REAL(w2f__8) OpenAD_Symbol_22 REAL(w2f__8) OpenAD_Symbol_18 REAL(w2f__8) OpenAD_Symbol_23 REAL(w2f__8) OpenAD_Symbol_19 REAL(w2f__8) OpenAD_Symbol_24 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_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_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_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_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_85 REAL(w2f__8) OpenAD_Symbol_90 REAL(w2f__8) OpenAD_Symbol_86 REAL(w2f__8) OpenAD_Symbol_87 REAL(w2f__8) OpenAD_Symbol_88 REAL(w2f__8) OpenAD_Symbol_210 type(active) :: OpenAD_Symbol_211 REAL(w2f__8) OpenAD_Symbol_212 type(active) :: OpenAD_Symbol_213 REAL(w2f__8) OpenAD_Symbol_208 REAL(w2f__8) OpenAD_Symbol_214 type(active) :: OpenAD_Symbol_209 REAL(w2f__8) OpenAD_Symbol_220 type(active) :: OpenAD_Symbol_215 type(active) :: OpenAD_Symbol_221 REAL(w2f__8) OpenAD_Symbol_216 REAL(w2f__8) OpenAD_Symbol_222 type(active) :: OpenAD_Symbol_217 REAL(w2f__8) OpenAD_Symbol_223 REAL(w2f__8) OpenAD_Symbol_218 REAL(w2f__8) OpenAD_Symbol_224 type(active) :: OpenAD_Symbol_219 REAL(w2f__8) OpenAD_Symbol_230 REAL(w2f__8) OpenAD_Symbol_225 REAL(w2f__8) OpenAD_Symbol_231 REAL(w2f__8) OpenAD_Symbol_226 type(active) :: OpenAD_Symbol_232 REAL(w2f__8) OpenAD_Symbol_227 REAL(w2f__8) OpenAD_Symbol_233 REAL(w2f__8) OpenAD_Symbol_228 REAL(w2f__8) OpenAD_Symbol_234 type(active) :: OpenAD_Symbol_229 REAL(w2f__8) OpenAD_Symbol_240 REAL(w2f__8) OpenAD_Symbol_235 REAL(w2f__8) OpenAD_Symbol_241 REAL(w2f__8) OpenAD_Symbol_236 type(active) :: OpenAD_Symbol_242 REAL(w2f__8) OpenAD_Symbol_237 REAL(w2f__8) OpenAD_Symbol_243 REAL(w2f__8) OpenAD_Symbol_238 REAL(w2f__8) OpenAD_Symbol_244 REAL(w2f__8) OpenAD_Symbol_239 REAL(w2f__8) OpenAD_Symbol_250 type(active) :: OpenAD_Symbol_245 REAL(w2f__8) OpenAD_Symbol_251 REAL(w2f__8) OpenAD_Symbol_246 type(active) :: OpenAD_Symbol_252 REAL(w2f__8) OpenAD_Symbol_247 REAL(w2f__8) OpenAD_Symbol_253 REAL(w2f__8) OpenAD_Symbol_248 type(active) :: OpenAD_Symbol_254 REAL(w2f__8) OpenAD_Symbol_249 type(active) :: OpenAD_Symbol_260 REAL(w2f__8) OpenAD_Symbol_255 REAL(w2f__8) OpenAD_Symbol_261 type(active) :: OpenAD_Symbol_256 type(active) :: OpenAD_Symbol_262 REAL(w2f__8) OpenAD_Symbol_257 REAL(w2f__8) OpenAD_Symbol_263 type(active) :: OpenAD_Symbol_258 REAL(w2f__8) OpenAD_Symbol_264 REAL(w2f__8) OpenAD_Symbol_259 type(active) :: OpenAD_Symbol_265 REAL(w2f__8) OpenAD_Symbol_266 type(active) :: OpenAD_Symbol_267 REAL(w2f__8) OpenAD_Symbol_268 type(active) :: OpenAD_Symbol_269 C C **** Parameters and Result **** C INTEGER(w2f__i4) NX INTEGER(w2f__i4) NY type(active) :: X(1 : INT((NX * NY))) type(active) :: FVEC(1 : INT((NX * NY))) type(active) :: R INTEGER(w2f__i4) N type(active) :: HX type(active) :: HY type(active) :: HY2 type(active) :: HX2 C C **** Local Variables and functions **** C INTEGER(w2f__i8) t__3 INTEGER(w2f__i8) t__4 EXTERNAL bar type(active) :: DPDX type(active) :: DPDY REAL(w2f__8) FOUR PARAMETER ( FOUR = 4.000000000000000) INTEGER(w2f__i4) I INTEGER(w2f__i4) J INTEGER(w2f__i4) K REAL(w2f__8) ONE PARAMETER ( ONE = 1.000000000000000) type(active) :: P type(active) :: PB type(active) :: PBB type(active) :: PBL type(active) :: PBLAP type(active) :: PBR type(active) :: PL type(active) :: PLAP type(active) :: PLL type(active) :: PLLAP type(active) :: PR type(active) :: PRLAP type(active) :: PRR type(active) :: PT type(active) :: PTL type(active) :: PTLAP type(active) :: PTR type(active) :: PTT REAL(w2f__8) THREE PARAMETER ( THREE = 3.000000000000000) REAL(w2f__8) TWO PARAMETER ( TWO = 2.000000000000000) REAL(w2f__8) ZERO PARAMETER ( ZERO = 0.000000000000000) C C **** Temporary variables **** C INTEGER(w2f__i4) doloop_ub INTEGER(w2f__i4) doloop_ub0 C C **** statements **** C t__3 = INT(NX * NY) t__4 = MAX(t__3, 0_w2f__i8) doloop_ub = NY DO I = 1, doloop_ub, 1 doloop_ub0 = NX DO J = 1, doloop_ub0, 1 K = J + NX *(I +(-1)) IF((I .eq. INT(1_w2f__i8)) .OR.(J .eq. INT(1_w2f__i8))) THEN PBL%v = 0.0D00 CALL zero_deriv(PBL) ELSE PBL%v = X(K-NX+(-1))%v CALL setderiv(PBL,X(K-NX+(-1))) ENDIF IF(I .eq. INT(1_w2f__i8)) THEN PB%v = 0.0D00 CALL zero_deriv(PB) PBB%v = X(K)%v CALL setderiv(PBB,X(K)) ELSE IF(I .eq. INT(2_w2f__i8)) THEN PB%v = X(K-NX)%v PBB%v = 0.0D00 CALL setderiv(PB,X(K-NX)) CALL zero_deriv(PBB) ELSE PB%v = X(K-NX)%v CALL setderiv(PB,X(K-NX)) PBB%v = X(K-NX*2)%v CALL setderiv(PBB,X(K-NX*2)) ENDIF ENDIF IF((NX .eq. J) .OR.(I .eq. INT(1_w2f__i8))) THEN PBR%v = 0.0D00 CALL zero_deriv(PBR) ELSE PBR%v = X(K-NX+1)%v CALL setderiv(PBR,X(K-NX+1)) ENDIF IF(J .eq. INT(1_w2f__i8)) THEN PL%v = 0.0D00 CALL zero_deriv(PL) PLL%v = X(K)%v CALL setderiv(PLL,X(K)) ELSE IF(J .eq. INT(2_w2f__i8)) THEN PL%v = X(K+(-1))%v PLL%v = 0.0D00 CALL setderiv(PL,X(K+(-1))) CALL zero_deriv(PLL) ELSE PL%v = X(K+(-1))%v CALL setderiv(PL,X(K+(-1))) PLL%v = X(K+(-2))%v CALL setderiv(PLL,X(K+(-2))) ENDIF ENDIF P%v = X(K)%v CALL setderiv(P,X(K)) IF(J .eq.(NX + INT((-1_w2f__i8)))) THEN PR%v = X(K+1)%v PRR%v = 0.0D00 CALL setderiv(PR,X(K+1)) CALL zero_deriv(PRR) ELSE IF(NX .eq. J) THEN PR%v = 0.0D00 CALL zero_deriv(PR) PRR%v = X(K)%v CALL setderiv(PRR,X(K)) ELSE PR%v = X(K+1)%v CALL setderiv(PR,X(K+1)) PRR%v = X(K+2)%v CALL setderiv(PRR,X(K+2)) ENDIF ENDIF IF((NY .eq. I) .OR.(J .eq. INT(1_w2f__i8))) THEN PTL%v = 0.0D00 CALL zero_deriv(PTL) ELSE PTL%v = X(NX+K+(-1))%v CALL setderiv(PTL,X(NX+K+(-1))) ENDIF IF(I .eq.(NY + INT((-1_w2f__i8)))) THEN PT%v = X(NX+K)%v PTT%v = 0.0D00 CALL setderiv(PT,X(NX+K)) CALL zero_deriv(PTT) ELSE IF(NY .eq. I) THEN PT%v = 0.0D00 CALL zero_deriv(PT) OpenAD_Symbol_90 = (X(K)%v+HY%v*2.0D00) OpenAD_Symbol_86 = 1_w2f__i8 OpenAD_Symbol_88 = 2.0D00 OpenAD_Symbol_87 = 1_w2f__i8 PTT%v = OpenAD_Symbol_90 OpenAD_Symbol_266 = (OpenAD_Symbol_88 * OpenAD_Symbol_87) OpenAD_Symbol_268 = OpenAD_Symbol_86 CALL setderiv(OpenAD_Symbol_269,X(K)) CALL setderiv(OpenAD_Symbol_267,HY) CALL sax(OpenAD_Symbol_266,OpenAD_Symbol_267,PTT) CALL saxpy(OpenAD_Symbol_268,OpenAD_Symbol_269,PTT) ELSE PT%v = X(NX+K)%v CALL setderiv(PT,X(NX+K)) PTT%v = X(K+NX*2)%v CALL setderiv(PTT,X(K+NX*2)) ENDIF ENDIF IF((NX .eq. J) .OR.(NY .eq. I)) THEN PTR%v = 0.0D00 CALL zero_deriv(PTR) ELSE PTR%v = X(NX+K+1)%v CALL setderiv(PTR,X(NX+K+1)) ENDIF OpenAD_Symbol_18 = (PT%v-PB%v) OpenAD_Symbol_19 = (HY%v*2.0D00) OpenAD_Symbol_26 = (OpenAD_Symbol_18 / OpenAD_Symbol_19) OpenAD_Symbol_22 = 1_w2f__i8 OpenAD_Symbol_23 = (-1_w2f__i8) OpenAD_Symbol_20 = (INT(1_w2f__i8) / OpenAD_Symbol_19) OpenAD_Symbol_24 = 2.0D00 OpenAD_Symbol_21 = (-(OpenAD_Symbol_18 /(OpenAD_Symbol_19 * > OpenAD_Symbol_19))) DPDY%v = OpenAD_Symbol_26 OpenAD_Symbol_27 = (PR%v-PL%v) OpenAD_Symbol_28 = (HX%v*2.0D00) OpenAD_Symbol_35 = (OpenAD_Symbol_27 / OpenAD_Symbol_28) OpenAD_Symbol_31 = 1_w2f__i8 OpenAD_Symbol_32 = (-1_w2f__i8) OpenAD_Symbol_29 = (INT(1_w2f__i8) / OpenAD_Symbol_28) OpenAD_Symbol_33 = 2.0D00 OpenAD_Symbol_30 = (-(OpenAD_Symbol_27 /(OpenAD_Symbol_28 * > OpenAD_Symbol_28))) DPDX%v = OpenAD_Symbol_35 OpenAD_Symbol_208 = (OpenAD_Symbol_22 * OpenAD_Symbol_20) OpenAD_Symbol_210 = (OpenAD_Symbol_23 * OpenAD_Symbol_20) OpenAD_Symbol_212 = (OpenAD_Symbol_24 * OpenAD_Symbol_21) OpenAD_Symbol_214 = (OpenAD_Symbol_31 * OpenAD_Symbol_29) OpenAD_Symbol_216 = (OpenAD_Symbol_32 * OpenAD_Symbol_29) OpenAD_Symbol_218 = (OpenAD_Symbol_33 * OpenAD_Symbol_30) CALL setderiv(OpenAD_Symbol_219,HX) CALL setderiv(OpenAD_Symbol_217,PL) CALL setderiv(OpenAD_Symbol_215,PR) CALL setderiv(OpenAD_Symbol_213,HY) CALL setderiv(OpenAD_Symbol_211,PB) CALL setderiv(OpenAD_Symbol_209,PT) CALL sax(OpenAD_Symbol_208,OpenAD_Symbol_209,DPDY) CALL saxpy(OpenAD_Symbol_210,OpenAD_Symbol_211,DPDY) CALL saxpy(OpenAD_Symbol_212,OpenAD_Symbol_213,DPDY) CALL sax(OpenAD_Symbol_214,OpenAD_Symbol_215,DPDX) CALL saxpy(OpenAD_Symbol_216,OpenAD_Symbol_217,DPDX) CALL saxpy(OpenAD_Symbol_218,OpenAD_Symbol_219,DPDX) CALL bar(PBB,PBL,PB,PBR,PLL,PL,P,PR,PRR,PTL,PT,PTR,PTT,HY2,HX2 +,PLAP,PBLAP,PLLAP,PRLAP,PTLAP) OpenAD_Symbol_40 = (PBLAP%v+PTLAP%v-PLAP%v*2.0D00) OpenAD_Symbol_49 = (PLLAP%v+PRLAP%v-PLAP%v*2.0D00) OpenAD_Symbol_67 = (PRLAP%v-PLLAP%v) OpenAD_Symbol_63 = (DPDY%v*OpenAD_Symbol_67) OpenAD_Symbol_64 = (HX%v*2.0D00) OpenAD_Symbol_78 = (PTLAP%v-PBLAP%v) OpenAD_Symbol_74 = (DPDX%v*OpenAD_Symbol_78) OpenAD_Symbol_75 = (HY%v*2.0D00) OpenAD_Symbol_58 = ((OpenAD_Symbol_63 / OpenAD_Symbol_64) -( > OpenAD_Symbol_74 / OpenAD_Symbol_75)) OpenAD_Symbol_85 = ((OpenAD_Symbol_40/HY2%v)+(OpenAD_Symbol_49 +/HX2%v)-R%v*OpenAD_Symbol_58) OpenAD_Symbol_43 = 1_w2f__i8 OpenAD_Symbol_45 = 1_w2f__i8 OpenAD_Symbol_47 = 2.0D00 OpenAD_Symbol_46 = (-1_w2f__i8) OpenAD_Symbol_44 = 1_w2f__i8 OpenAD_Symbol_41 = (INT(1 _w2f__i8)/HY2%v) OpenAD_Symbol_42 = (-(OpenAD_Symbol_40/(HY2%v*HY2%v))) OpenAD_Symbol_38 = 1_w2f__i8 OpenAD_Symbol_52 = 1_w2f__i8 OpenAD_Symbol_54 = 1_w2f__i8 OpenAD_Symbol_56 = 2.0D00 OpenAD_Symbol_55 = (-1_w2f__i8) OpenAD_Symbol_53 = 1_w2f__i8 OpenAD_Symbol_50 = (INT(1 _w2f__i8)/HX2%v) OpenAD_Symbol_51 = (-(OpenAD_Symbol_49/(HX2%v*HX2%v))) OpenAD_Symbol_39 = 1_w2f__i8 OpenAD_Symbol_36 = 1_w2f__i8 OpenAD_Symbol_59 = OpenAD_Symbol_58 OpenAD_Symbol_68 = OpenAD_Symbol_67 OpenAD_Symbol_70 = 1_w2f__i8 OpenAD_Symbol_71 = (-1_w2f__i8) OpenAD_Symbol_69 = DPDY%v OpenAD_Symbol_65 = (INT(1_w2f__i8) / OpenAD_Symbol_64) OpenAD_Symbol_72 = 2.0D00 OpenAD_Symbol_66 = (-(OpenAD_Symbol_63 /(OpenAD_Symbol_64 * > OpenAD_Symbol_64))) OpenAD_Symbol_61 = 1_w2f__i8 OpenAD_Symbol_79 = OpenAD_Symbol_78 OpenAD_Symbol_81 = 1_w2f__i8 OpenAD_Symbol_82 = (-1_w2f__i8) OpenAD_Symbol_80 = DPDX%v OpenAD_Symbol_76 = (INT(1_w2f__i8) / OpenAD_Symbol_75) OpenAD_Symbol_83 = 2.0D00 OpenAD_Symbol_77 = (-(OpenAD_Symbol_74 /(OpenAD_Symbol_75 * > OpenAD_Symbol_75))) OpenAD_Symbol_62 = (-1_w2f__i8) OpenAD_Symbol_60 = R%v OpenAD_Symbol_37 = (-1_w2f__i8) FVEC(INT(K))%v = OpenAD_Symbol_85 OpenAD_Symbol_220 = (OpenAD_Symbol_59 * OpenAD_Symbol_37) OpenAD_Symbol_222 = (OpenAD_Symbol_60 * OpenAD_Symbol_37) OpenAD_Symbol_223 = (OpenAD_Symbol_38 * OpenAD_Symbol_36) OpenAD_Symbol_224 = (OpenAD_Symbol_39 * OpenAD_Symbol_36) OpenAD_Symbol_225 = (OpenAD_Symbol_61 * OpenAD_Symbol_222) OpenAD_Symbol_226 = (OpenAD_Symbol_62 * OpenAD_Symbol_222) OpenAD_Symbol_227 = (OpenAD_Symbol_50 * OpenAD_Symbol_224) OpenAD_Symbol_228 = (OpenAD_Symbol_51 * OpenAD_Symbol_224) OpenAD_Symbol_230 = (OpenAD_Symbol_41 * OpenAD_Symbol_223) OpenAD_Symbol_231 = (OpenAD_Symbol_42 * OpenAD_Symbol_223) OpenAD_Symbol_233 = (OpenAD_Symbol_76 * OpenAD_Symbol_226) OpenAD_Symbol_234 = (OpenAD_Symbol_77 * OpenAD_Symbol_226) OpenAD_Symbol_235 = (OpenAD_Symbol_65 * OpenAD_Symbol_225) OpenAD_Symbol_236 = (OpenAD_Symbol_66 * OpenAD_Symbol_225) OpenAD_Symbol_237 = (OpenAD_Symbol_52 * OpenAD_Symbol_227) OpenAD_Symbol_238 = (OpenAD_Symbol_53 * OpenAD_Symbol_227) OpenAD_Symbol_239 = (OpenAD_Symbol_43 * OpenAD_Symbol_230) OpenAD_Symbol_240 = (OpenAD_Symbol_44 * OpenAD_Symbol_230) OpenAD_Symbol_241 = (OpenAD_Symbol_79 * OpenAD_Symbol_233) OpenAD_Symbol_243 = (OpenAD_Symbol_80 * OpenAD_Symbol_233) OpenAD_Symbol_244 = (OpenAD_Symbol_68 * OpenAD_Symbol_235) OpenAD_Symbol_246 = (OpenAD_Symbol_69 * OpenAD_Symbol_235) OpenAD_Symbol_247 = (OpenAD_Symbol_54 * OpenAD_Symbol_238) OpenAD_Symbol_248 = (OpenAD_Symbol_55 * OpenAD_Symbol_238) OpenAD_Symbol_249 = (OpenAD_Symbol_45 * OpenAD_Symbol_240) OpenAD_Symbol_250 = (OpenAD_Symbol_46 * OpenAD_Symbol_240) OpenAD_Symbol_251 = (OpenAD_Symbol_249 + OpenAD_Symbol_81 * > OpenAD_Symbol_243) OpenAD_Symbol_253 = (OpenAD_Symbol_239 + OpenAD_Symbol_82 * > OpenAD_Symbol_243) OpenAD_Symbol_255 = (OpenAD_Symbol_247 + OpenAD_Symbol_70 * > OpenAD_Symbol_246) OpenAD_Symbol_257 = (OpenAD_Symbol_237 + OpenAD_Symbol_71 * > OpenAD_Symbol_246) OpenAD_Symbol_259 = (OpenAD_Symbol_83 * OpenAD_Symbol_234) OpenAD_Symbol_261 = (OpenAD_Symbol_72 * OpenAD_Symbol_236) OpenAD_Symbol_263 = (OpenAD_Symbol_56 * OpenAD_Symbol_248) OpenAD_Symbol_264 = (OpenAD_Symbol_263 + OpenAD_Symbol_47 * > OpenAD_Symbol_250) CALL setderiv(OpenAD_Symbol_265,PLAP) CALL setderiv(OpenAD_Symbol_262,HX) CALL setderiv(OpenAD_Symbol_260,HY) CALL setderiv(OpenAD_Symbol_258,PLLAP) CALL setderiv(OpenAD_Symbol_256,PRLAP) CALL setderiv(OpenAD_Symbol_254,PBLAP) CALL setderiv(OpenAD_Symbol_252,PTLAP) CALL setderiv(OpenAD_Symbol_245,DPDY) CALL setderiv(OpenAD_Symbol_242,DPDX) CALL setderiv(OpenAD_Symbol_232,HY2) CALL setderiv(OpenAD_Symbol_229,HX2) CALL setderiv(OpenAD_Symbol_221,R) CALL sax(OpenAD_Symbol_220,OpenAD_Symbol_221,FVEC(K)) CALL saxpy(OpenAD_Symbol_228,OpenAD_Symbol_229,FVEC(K)) CALL saxpy(OpenAD_Symbol_231,OpenAD_Symbol_232,FVEC(K)) CALL saxpy(OpenAD_Symbol_241,OpenAD_Symbol_242,FVEC(K)) CALL saxpy(OpenAD_Symbol_244,OpenAD_Symbol_245,FVEC(K)) CALL saxpy(OpenAD_Symbol_251,OpenAD_Symbol_252,FVEC(K)) CALL saxpy(OpenAD_Symbol_253,OpenAD_Symbol_254,FVEC(K)) CALL saxpy(OpenAD_Symbol_255,OpenAD_Symbol_256,FVEC(K)) CALL saxpy(OpenAD_Symbol_257,OpenAD_Symbol_258,FVEC(K)) CALL saxpy(OpenAD_Symbol_259,OpenAD_Symbol_260,FVEC(K)) CALL saxpy(OpenAD_Symbol_261,OpenAD_Symbol_262,FVEC(K)) CALL saxpy(OpenAD_Symbol_264,OpenAD_Symbol_265,FVEC(K)) END DO END DO RETURN END SUBROUTINE SUBROUTINE bar(PBB, PBL, PB, PBR, PLL, PL, P, PR, PRR, PTL, PT, > PTR, PTT, HY2, HX2, PLAP, PBLAP, PLLAP, PRLAP, PTLAP) use w2f__types use active_module IMPLICIT NONE C C **** Global Variables **** C REAL(w2f__8) OpenAD_Symbol_91 REAL(w2f__8) OpenAD_Symbol_92 REAL(w2f__8) OpenAD_Symbol_93 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 REAL(w2f__8) OpenAD_Symbol_100 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_106 REAL(w2f__8) OpenAD_Symbol_111 REAL(w2f__8) OpenAD_Symbol_107 REAL(w2f__8) OpenAD_Symbol_112 REAL(w2f__8) OpenAD_Symbol_108 REAL(w2f__8) OpenAD_Symbol_113 REAL(w2f__8) OpenAD_Symbol_109 REAL(w2f__8) OpenAD_Symbol_114 REAL(w2f__8) OpenAD_Symbol_115 REAL(w2f__8) OpenAD_Symbol_120 REAL(w2f__8) OpenAD_Symbol_116 REAL(w2f__8) OpenAD_Symbol_121 REAL(w2f__8) OpenAD_Symbol_117 REAL(w2f__8) OpenAD_Symbol_118 REAL(w2f__8) OpenAD_Symbol_123 REAL(w2f__8) OpenAD_Symbol_119 REAL(w2f__8) OpenAD_Symbol_124 REAL(w2f__8) OpenAD_Symbol_125 REAL(w2f__8) OpenAD_Symbol_130 REAL(w2f__8) OpenAD_Symbol_126 REAL(w2f__8) OpenAD_Symbol_127 REAL(w2f__8) OpenAD_Symbol_132 REAL(w2f__8) OpenAD_Symbol_128 REAL(w2f__8) OpenAD_Symbol_133 REAL(w2f__8) OpenAD_Symbol_129 REAL(w2f__8) OpenAD_Symbol_134 REAL(w2f__8) OpenAD_Symbol_135 REAL(w2f__8) OpenAD_Symbol_140 REAL(w2f__8) OpenAD_Symbol_136 REAL(w2f__8) OpenAD_Symbol_141 REAL(w2f__8) OpenAD_Symbol_137 REAL(w2f__8) OpenAD_Symbol_142 REAL(w2f__8) OpenAD_Symbol_138 REAL(w2f__8) OpenAD_Symbol_139 REAL(w2f__8) OpenAD_Symbol_144 REAL(w2f__8) OpenAD_Symbol_145 REAL(w2f__8) OpenAD_Symbol_150 REAL(w2f__8) OpenAD_Symbol_146 REAL(w2f__8) OpenAD_Symbol_151 REAL(w2f__8) OpenAD_Symbol_147 REAL(w2f__8) OpenAD_Symbol_148 REAL(w2f__8) OpenAD_Symbol_153 REAL(w2f__8) OpenAD_Symbol_149 REAL(w2f__8) OpenAD_Symbol_154 REAL(w2f__8) OpenAD_Symbol_155 REAL(w2f__8) OpenAD_Symbol_160 REAL(w2f__8) OpenAD_Symbol_156 REAL(w2f__8) OpenAD_Symbol_161 REAL(w2f__8) OpenAD_Symbol_157 REAL(w2f__8) OpenAD_Symbol_162 REAL(w2f__8) OpenAD_Symbol_158 REAL(w2f__8) OpenAD_Symbol_163 REAL(w2f__8) OpenAD_Symbol_159 REAL(w2f__8) OpenAD_Symbol_165 REAL(w2f__8) OpenAD_Symbol_170 REAL(w2f__8) OpenAD_Symbol_166 REAL(w2f__8) OpenAD_Symbol_171 REAL(w2f__8) OpenAD_Symbol_167 REAL(w2f__8) OpenAD_Symbol_172 REAL(w2f__8) OpenAD_Symbol_168 REAL(w2f__8) OpenAD_Symbol_174 REAL(w2f__8) OpenAD_Symbol_169 REAL(w2f__8) OpenAD_Symbol_180 REAL(w2f__8) OpenAD_Symbol_175 REAL(w2f__8) OpenAD_Symbol_181 REAL(w2f__8) OpenAD_Symbol_176 REAL(w2f__8) OpenAD_Symbol_182 REAL(w2f__8) OpenAD_Symbol_177 REAL(w2f__8) OpenAD_Symbol_183 REAL(w2f__8) OpenAD_Symbol_178 REAL(w2f__8) OpenAD_Symbol_184 REAL(w2f__8) OpenAD_Symbol_179 REAL(w2f__8) OpenAD_Symbol_190 REAL(w2f__8) OpenAD_Symbol_191 REAL(w2f__8) OpenAD_Symbol_186 REAL(w2f__8) OpenAD_Symbol_192 REAL(w2f__8) OpenAD_Symbol_187 REAL(w2f__8) OpenAD_Symbol_193 REAL(w2f__8) OpenAD_Symbol_188 REAL(w2f__8) OpenAD_Symbol_189 REAL(w2f__8) OpenAD_Symbol_300 REAL(w2f__8) OpenAD_Symbol_195 type(active) :: OpenAD_Symbol_301 REAL(w2f__8) OpenAD_Symbol_302 REAL(w2f__8) OpenAD_Symbol_303 type(active) :: OpenAD_Symbol_304 type(active) :: OpenAD_Symbol_310 REAL(w2f__8) OpenAD_Symbol_305 REAL(w2f__8) OpenAD_Symbol_311 REAL(w2f__8) OpenAD_Symbol_306 REAL(w2f__8) OpenAD_Symbol_312 type(active) :: OpenAD_Symbol_307 type(active) :: OpenAD_Symbol_313 REAL(w2f__8) OpenAD_Symbol_308 REAL(w2f__8) OpenAD_Symbol_314 REAL(w2f__8) OpenAD_Symbol_309 REAL(w2f__8) OpenAD_Symbol_320 REAL(w2f__8) OpenAD_Symbol_315 REAL(w2f__8) OpenAD_Symbol_270 type(active) :: OpenAD_Symbol_321 REAL(w2f__8) OpenAD_Symbol_316 REAL(w2f__8) OpenAD_Symbol_271 REAL(w2f__8) OpenAD_Symbol_322 REAL(w2f__8) OpenAD_Symbol_317 type(active) :: OpenAD_Symbol_272 REAL(w2f__8) OpenAD_Symbol_323 type(active) :: OpenAD_Symbol_318 REAL(w2f__8) OpenAD_Symbol_273 type(active) :: OpenAD_Symbol_324 REAL(w2f__8) OpenAD_Symbol_319 REAL(w2f__8) OpenAD_Symbol_274 REAL(w2f__8) OpenAD_Symbol_330 REAL(w2f__8) OpenAD_Symbol_325 REAL(w2f__8) OpenAD_Symbol_280 type(active) :: OpenAD_Symbol_275 REAL(w2f__8) OpenAD_Symbol_331 REAL(w2f__8) OpenAD_Symbol_326 REAL(w2f__8) OpenAD_Symbol_281 REAL(w2f__8) OpenAD_Symbol_276 type(active) :: OpenAD_Symbol_332 type(active) :: OpenAD_Symbol_327 REAL(w2f__8) OpenAD_Symbol_282 REAL(w2f__8) OpenAD_Symbol_277 REAL(w2f__8) OpenAD_Symbol_333 REAL(w2f__8) OpenAD_Symbol_328 REAL(w2f__8) OpenAD_Symbol_283 REAL(w2f__8) OpenAD_Symbol_278 REAL(w2f__8) OpenAD_Symbol_334 REAL(w2f__8) OpenAD_Symbol_329 REAL(w2f__8) OpenAD_Symbol_284 REAL(w2f__8) OpenAD_Symbol_279 REAL(w2f__8) OpenAD_Symbol_340 type(active) :: OpenAD_Symbol_335 REAL(w2f__8) OpenAD_Symbol_290 REAL(w2f__8) OpenAD_Symbol_285 REAL(w2f__8) OpenAD_Symbol_341 REAL(w2f__8) OpenAD_Symbol_336 REAL(w2f__8) OpenAD_Symbol_291 REAL(w2f__8) OpenAD_Symbol_286 REAL(w2f__8) OpenAD_Symbol_342 REAL(w2f__8) OpenAD_Symbol_337 REAL(w2f__8) OpenAD_Symbol_292 REAL(w2f__8) OpenAD_Symbol_287 REAL(w2f__8) OpenAD_Symbol_343 REAL(w2f__8) OpenAD_Symbol_338 type(active) :: OpenAD_Symbol_293 REAL(w2f__8) OpenAD_Symbol_288 REAL(w2f__8) OpenAD_Symbol_344 REAL(w2f__8) OpenAD_Symbol_339 REAL(w2f__8) OpenAD_Symbol_294 REAL(w2f__8) OpenAD_Symbol_289 REAL(w2f__8) OpenAD_Symbol_350 REAL(w2f__8) OpenAD_Symbol_345 REAL(w2f__8) OpenAD_Symbol_295 REAL(w2f__8) OpenAD_Symbol_351 REAL(w2f__8) OpenAD_Symbol_346 type(active) :: OpenAD_Symbol_296 REAL(w2f__8) OpenAD_Symbol_352 REAL(w2f__8) OpenAD_Symbol_347 REAL(w2f__8) OpenAD_Symbol_297 REAL(w2f__8) OpenAD_Symbol_353 REAL(w2f__8) OpenAD_Symbol_348 REAL(w2f__8) OpenAD_Symbol_298 REAL(w2f__8) OpenAD_Symbol_354 REAL(w2f__8) OpenAD_Symbol_349 REAL(w2f__8) OpenAD_Symbol_299 C C **** Parameters and Result **** C type(active) :: PBB type(active) :: PBL type(active) :: PB type(active) :: PBR type(active) :: PLL type(active) :: PL type(active) :: P type(active) :: PR type(active) :: PRR type(active) :: PTL type(active) :: PT type(active) :: PTR type(active) :: PTT type(active) :: HY2 type(active) :: HX2 type(active) :: PLAP type(active) :: PBLAP type(active) :: PLLAP type(active) :: PRLAP type(active) :: PTLAP C C **** Local Variables and functions **** C REAL(w2f__8) FOUR PARAMETER ( FOUR = 4.000000000000000) REAL(w2f__8) ONE PARAMETER ( ONE = 1.000000000000000) REAL(w2f__8) THREE PARAMETER ( THREE = 3.000000000000000) REAL(w2f__8) TWO PARAMETER ( TWO = 2.000000000000000) REAL(w2f__8) ZERO PARAMETER ( ZERO = 0.000000000000000) C C **** statements **** C OpenAD_Symbol_93 = (PBB%v+P%v-PB%v*2.0D00) OpenAD_Symbol_102 = (PBL%v+PBR%v-PB%v*2.0D00) OpenAD_Symbol_111 = ((OpenAD_Symbol_93/HY2%v)+(OpenAD_Symbol_102/H +X2%v)) OpenAD_Symbol_96 = 1_w2f__i8 OpenAD_Symbol_98 = 1_w2f__i8 OpenAD_Symbol_100 = 2.0D00 OpenAD_Symbol_99 = (-1_w2f__i8) OpenAD_Symbol_97 = 1_w2f__i8 OpenAD_Symbol_94 = (INT(1 _w2f__i8)/HY2%v) OpenAD_Symbol_95 = (-(OpenAD_Symbol_93/(HY2%v*HY2%v))) OpenAD_Symbol_91 = 1_w2f__i8 OpenAD_Symbol_105 = 1_w2f__i8 OpenAD_Symbol_107 = 1_w2f__i8 OpenAD_Symbol_109 = 2.0D00 OpenAD_Symbol_108 = (-1_w2f__i8) OpenAD_Symbol_106 = 1_w2f__i8 OpenAD_Symbol_103 = (INT(1 _w2f__i8)/HX2%v) OpenAD_Symbol_104 = (-(OpenAD_Symbol_102/(HX2%v*HX2%v))) OpenAD_Symbol_92 = 1_w2f__i8 PBLAP%v = OpenAD_Symbol_111 OpenAD_Symbol_114 = (PBL%v+PTL%v-PL%v*2.0D00) OpenAD_Symbol_123 = (PLL%v+P%v-PL%v*2.0D00) OpenAD_Symbol_132 = ((OpenAD_Symbol_114/HY2%v)+(OpenAD_Symbol_123/ +HX2%v)) OpenAD_Symbol_117 = 1_w2f__i8 OpenAD_Symbol_119 = 1_w2f__i8 OpenAD_Symbol_121 = 2.0D00 OpenAD_Symbol_120 = (-1_w2f__i8) OpenAD_Symbol_118 = 1_w2f__i8 OpenAD_Symbol_115 = (INT(1 _w2f__i8)/HY2%v) OpenAD_Symbol_116 = (-(OpenAD_Symbol_114/(HY2%v*HY2%v))) OpenAD_Symbol_112 = 1_w2f__i8 OpenAD_Symbol_126 = 1_w2f__i8 OpenAD_Symbol_128 = 1_w2f__i8 OpenAD_Symbol_130 = 2.0D00 OpenAD_Symbol_129 = (-1_w2f__i8) OpenAD_Symbol_127 = 1_w2f__i8 OpenAD_Symbol_124 = (INT(1 _w2f__i8)/HX2%v) OpenAD_Symbol_125 = (-(OpenAD_Symbol_123/(HX2%v*HX2%v))) OpenAD_Symbol_113 = 1_w2f__i8 PLLAP%v = OpenAD_Symbol_132 OpenAD_Symbol_135 = (PB%v+PT%v-P%v*2.0D00) OpenAD_Symbol_144 = (PL%v+PR%v-P%v*2.0D00) OpenAD_Symbol_153 = ((OpenAD_Symbol_135/HY2%v)+(OpenAD_Symbol_144/ +HX2%v)) OpenAD_Symbol_138 = 1_w2f__i8 OpenAD_Symbol_140 = 1_w2f__i8 OpenAD_Symbol_142 = 2.0D00 OpenAD_Symbol_141 = (-1_w2f__i8) OpenAD_Symbol_139 = 1_w2f__i8 OpenAD_Symbol_136 = (INT(1 _w2f__i8)/HY2%v) OpenAD_Symbol_137 = (-(OpenAD_Symbol_135/(HY2%v*HY2%v))) OpenAD_Symbol_133 = 1_w2f__i8 OpenAD_Symbol_147 = 1_w2f__i8 OpenAD_Symbol_149 = 1_w2f__i8 OpenAD_Symbol_151 = 2.0D00 OpenAD_Symbol_150 = (-1_w2f__i8) OpenAD_Symbol_148 = 1_w2f__i8 OpenAD_Symbol_145 = (INT(1 _w2f__i8)/HX2%v) OpenAD_Symbol_146 = (-(OpenAD_Symbol_144/(HX2%v*HX2%v))) OpenAD_Symbol_134 = 1_w2f__i8 PLAP%v = OpenAD_Symbol_153 OpenAD_Symbol_156 = (PBR%v+PTR%v-PR%v*2.0D00) OpenAD_Symbol_165 = (P%v+PRR%v-PR%v*2.0D00) OpenAD_Symbol_174 = ((OpenAD_Symbol_156/HY2%v)+(OpenAD_Symbol_165/ +HX2%v)) OpenAD_Symbol_159 = 1_w2f__i8 OpenAD_Symbol_161 = 1_w2f__i8 OpenAD_Symbol_163 = 2.0D00 OpenAD_Symbol_162 = (-1_w2f__i8) OpenAD_Symbol_160 = 1_w2f__i8 OpenAD_Symbol_157 = (INT(1 _w2f__i8)/HY2%v) OpenAD_Symbol_158 = (-(OpenAD_Symbol_156/(HY2%v*HY2%v))) OpenAD_Symbol_154 = 1_w2f__i8 OpenAD_Symbol_168 = 1_w2f__i8 OpenAD_Symbol_170 = 1_w2f__i8 OpenAD_Symbol_172 = 2.0D00 OpenAD_Symbol_171 = (-1_w2f__i8) OpenAD_Symbol_169 = 1_w2f__i8 OpenAD_Symbol_166 = (INT(1 _w2f__i8)/HX2%v) OpenAD_Symbol_167 = (-(OpenAD_Symbol_165/(HX2%v*HX2%v))) OpenAD_Symbol_155 = 1_w2f__i8 PRLAP%v = OpenAD_Symbol_174 OpenAD_Symbol_177 = (P%v+PTT%v-PT%v*2.0D00) OpenAD_Symbol_186 = (PTL%v+PTR%v-PT%v*2.0D00) OpenAD_Symbol_195 = ((OpenAD_Symbol_177/HY2%v)+(OpenAD_Symbol_186/ +HX2%v)) OpenAD_Symbol_180 = 1_w2f__i8 OpenAD_Symbol_182 = 1_w2f__i8 OpenAD_Symbol_184 = 2.0D00 OpenAD_Symbol_183 = (-1_w2f__i8) OpenAD_Symbol_181 = 1_w2f__i8 OpenAD_Symbol_178 = (INT(1 _w2f__i8)/HY2%v) OpenAD_Symbol_179 = (-(OpenAD_Symbol_177/(HY2%v*HY2%v))) OpenAD_Symbol_175 = 1_w2f__i8 OpenAD_Symbol_189 = 1_w2f__i8 OpenAD_Symbol_191 = 1_w2f__i8 OpenAD_Symbol_193 = 2.0D00 OpenAD_Symbol_192 = (-1_w2f__i8) OpenAD_Symbol_190 = 1_w2f__i8 OpenAD_Symbol_187 = (INT(1 _w2f__i8)/HX2%v) OpenAD_Symbol_188 = (-(OpenAD_Symbol_186/(HX2%v*HX2%v))) OpenAD_Symbol_176 = 1_w2f__i8 PTLAP%v = OpenAD_Symbol_195 OpenAD_Symbol_270 = (OpenAD_Symbol_187 * OpenAD_Symbol_176) OpenAD_Symbol_271 = (OpenAD_Symbol_188 * OpenAD_Symbol_176) OpenAD_Symbol_273 = (OpenAD_Symbol_178 * OpenAD_Symbol_175) OpenAD_Symbol_274 = (OpenAD_Symbol_179 * OpenAD_Symbol_175) OpenAD_Symbol_276 = (OpenAD_Symbol_166 * OpenAD_Symbol_155) OpenAD_Symbol_277 = (OpenAD_Symbol_167 * OpenAD_Symbol_155) OpenAD_Symbol_278 = (OpenAD_Symbol_157 * OpenAD_Symbol_154) OpenAD_Symbol_279 = (OpenAD_Symbol_158 * OpenAD_Symbol_154) OpenAD_Symbol_280 = (OpenAD_Symbol_145 * OpenAD_Symbol_134) OpenAD_Symbol_281 = (OpenAD_Symbol_146 * OpenAD_Symbol_134) OpenAD_Symbol_282 = (OpenAD_Symbol_136 * OpenAD_Symbol_133) OpenAD_Symbol_283 = (OpenAD_Symbol_137 * OpenAD_Symbol_133) OpenAD_Symbol_284 = (OpenAD_Symbol_124 * OpenAD_Symbol_113) OpenAD_Symbol_285 = (OpenAD_Symbol_125 * OpenAD_Symbol_113) OpenAD_Symbol_286 = (OpenAD_Symbol_115 * OpenAD_Symbol_112) OpenAD_Symbol_287 = (OpenAD_Symbol_116 * OpenAD_Symbol_112) OpenAD_Symbol_288 = (OpenAD_Symbol_103 * OpenAD_Symbol_92) OpenAD_Symbol_289 = (OpenAD_Symbol_104 * OpenAD_Symbol_92) OpenAD_Symbol_290 = (OpenAD_Symbol_94 * OpenAD_Symbol_91) OpenAD_Symbol_291 = (OpenAD_Symbol_95 * OpenAD_Symbol_91) OpenAD_Symbol_292 = (OpenAD_Symbol_189 * OpenAD_Symbol_270) OpenAD_Symbol_294 = (OpenAD_Symbol_190 * OpenAD_Symbol_270) OpenAD_Symbol_295 = (OpenAD_Symbol_180 * OpenAD_Symbol_273) OpenAD_Symbol_297 = (OpenAD_Symbol_181 * OpenAD_Symbol_273) OpenAD_Symbol_298 = (OpenAD_Symbol_168 * OpenAD_Symbol_276) OpenAD_Symbol_299 = (OpenAD_Symbol_169 * OpenAD_Symbol_276) OpenAD_Symbol_300 = (OpenAD_Symbol_159 * OpenAD_Symbol_278) OpenAD_Symbol_302 = (OpenAD_Symbol_160 * OpenAD_Symbol_278) OpenAD_Symbol_303 = (OpenAD_Symbol_147 * OpenAD_Symbol_280) OpenAD_Symbol_305 = (OpenAD_Symbol_148 * OpenAD_Symbol_280) OpenAD_Symbol_306 = (OpenAD_Symbol_138 * OpenAD_Symbol_282) OpenAD_Symbol_308 = (OpenAD_Symbol_139 * OpenAD_Symbol_282) OpenAD_Symbol_309 = (OpenAD_Symbol_126 * OpenAD_Symbol_284) OpenAD_Symbol_311 = (OpenAD_Symbol_127 * OpenAD_Symbol_284) OpenAD_Symbol_312 = (OpenAD_Symbol_117 * OpenAD_Symbol_286) OpenAD_Symbol_314 = (OpenAD_Symbol_118 * OpenAD_Symbol_286) OpenAD_Symbol_315 = (OpenAD_Symbol_105 * OpenAD_Symbol_288) OpenAD_Symbol_316 = (OpenAD_Symbol_106 * OpenAD_Symbol_288) OpenAD_Symbol_317 = (OpenAD_Symbol_96 * OpenAD_Symbol_290) OpenAD_Symbol_319 = (OpenAD_Symbol_97 * OpenAD_Symbol_290) OpenAD_Symbol_320 = (OpenAD_Symbol_191 * OpenAD_Symbol_294) OpenAD_Symbol_322 = (OpenAD_Symbol_192 * OpenAD_Symbol_294) OpenAD_Symbol_323 = (OpenAD_Symbol_182 * OpenAD_Symbol_297) OpenAD_Symbol_325 = (OpenAD_Symbol_183 * OpenAD_Symbol_297) OpenAD_Symbol_326 = (OpenAD_Symbol_170 * OpenAD_Symbol_299) OpenAD_Symbol_328 = (OpenAD_Symbol_171 * OpenAD_Symbol_299) OpenAD_Symbol_329 = (OpenAD_Symbol_161 * OpenAD_Symbol_302) OpenAD_Symbol_330 = (OpenAD_Symbol_162 * OpenAD_Symbol_302) OpenAD_Symbol_331 = (OpenAD_Symbol_149 * OpenAD_Symbol_305) OpenAD_Symbol_333 = (OpenAD_Symbol_150 * OpenAD_Symbol_305) OpenAD_Symbol_334 = (OpenAD_Symbol_140 * OpenAD_Symbol_308) OpenAD_Symbol_336 = (OpenAD_Symbol_141 * OpenAD_Symbol_308) OpenAD_Symbol_337 = (OpenAD_Symbol_128 * OpenAD_Symbol_311) OpenAD_Symbol_338 = (OpenAD_Symbol_129 * OpenAD_Symbol_311) OpenAD_Symbol_339 = (OpenAD_Symbol_119 * OpenAD_Symbol_314) OpenAD_Symbol_340 = (OpenAD_Symbol_120 * OpenAD_Symbol_314) OpenAD_Symbol_341 = (OpenAD_Symbol_107 * OpenAD_Symbol_316) OpenAD_Symbol_342 = (OpenAD_Symbol_108 * OpenAD_Symbol_316) OpenAD_Symbol_343 = (OpenAD_Symbol_98 * OpenAD_Symbol_319) OpenAD_Symbol_344 = (OpenAD_Symbol_99 * OpenAD_Symbol_319) OpenAD_Symbol_345 = (OpenAD_Symbol_193 * OpenAD_Symbol_322) OpenAD_Symbol_346 = (OpenAD_Symbol_345 + OpenAD_Symbol_184 * > OpenAD_Symbol_325) OpenAD_Symbol_347 = (OpenAD_Symbol_172 * OpenAD_Symbol_328) OpenAD_Symbol_348 = (OpenAD_Symbol_347 + OpenAD_Symbol_163 * > OpenAD_Symbol_330) OpenAD_Symbol_349 = (OpenAD_Symbol_151 * OpenAD_Symbol_333) OpenAD_Symbol_350 = (OpenAD_Symbol_349 + OpenAD_Symbol_142 * > OpenAD_Symbol_336) OpenAD_Symbol_351 = (OpenAD_Symbol_130 * OpenAD_Symbol_338) OpenAD_Symbol_352 = (OpenAD_Symbol_351 + OpenAD_Symbol_121 * > OpenAD_Symbol_340) OpenAD_Symbol_353 = (OpenAD_Symbol_109 * OpenAD_Symbol_342) OpenAD_Symbol_354 = (OpenAD_Symbol_353 + OpenAD_Symbol_100 * > OpenAD_Symbol_344) CALL setderiv(OpenAD_Symbol_335,PT) CALL setderiv(OpenAD_Symbol_332,PR) CALL setderiv(OpenAD_Symbol_327,PRR) CALL setderiv(OpenAD_Symbol_324,PTT) CALL setderiv(OpenAD_Symbol_321,PTR) CALL setderiv(OpenAD_Symbol_318,PBB) CALL setderiv(OpenAD_Symbol_313,PBL) CALL setderiv(OpenAD_Symbol_310,PLL) CALL setderiv(OpenAD_Symbol_307,PB) CALL setderiv(OpenAD_Symbol_304,PL) CALL setderiv(OpenAD_Symbol_301,PBR) CALL setderiv(OpenAD_Symbol_296,P) CALL setderiv(OpenAD_Symbol_293,PTL) CALL setderiv(OpenAD_Symbol_275,HY2) CALL setderiv(OpenAD_Symbol_272,HX2) CALL sax(OpenAD_Symbol_271,OpenAD_Symbol_272,PTLAP) CALL saxpy(OpenAD_Symbol_274,OpenAD_Symbol_275,PTLAP) CALL sax(OpenAD_Symbol_277,OpenAD_Symbol_272,PRLAP) CALL saxpy(OpenAD_Symbol_279,OpenAD_Symbol_275,PRLAP) CALL sax(OpenAD_Symbol_281,OpenAD_Symbol_272,PLAP) CALL saxpy(OpenAD_Symbol_283,OpenAD_Symbol_275,PLAP) CALL sax(OpenAD_Symbol_285,OpenAD_Symbol_272,PLLAP) CALL saxpy(OpenAD_Symbol_287,OpenAD_Symbol_275,PLLAP) CALL sax(OpenAD_Symbol_289,OpenAD_Symbol_272,PBLAP) CALL saxpy(OpenAD_Symbol_291,OpenAD_Symbol_275,PBLAP) CALL saxpy(OpenAD_Symbol_292,OpenAD_Symbol_293,PTLAP) CALL saxpy(OpenAD_Symbol_295,OpenAD_Symbol_296,PTLAP) CALL saxpy(OpenAD_Symbol_298,OpenAD_Symbol_296,PRLAP) CALL saxpy(OpenAD_Symbol_300,OpenAD_Symbol_301,PRLAP) CALL saxpy(OpenAD_Symbol_303,OpenAD_Symbol_304,PLAP) CALL saxpy(OpenAD_Symbol_306,OpenAD_Symbol_307,PLAP) CALL saxpy(OpenAD_Symbol_309,OpenAD_Symbol_310,PLLAP) CALL saxpy(OpenAD_Symbol_312,OpenAD_Symbol_313,PLLAP) CALL saxpy(OpenAD_Symbol_315,OpenAD_Symbol_313,PBLAP) CALL saxpy(OpenAD_Symbol_317,OpenAD_Symbol_318,PBLAP) CALL saxpy(OpenAD_Symbol_320,OpenAD_Symbol_321,PTLAP) CALL saxpy(OpenAD_Symbol_323,OpenAD_Symbol_324,PTLAP) CALL saxpy(OpenAD_Symbol_326,OpenAD_Symbol_327,PRLAP) CALL saxpy(OpenAD_Symbol_329,OpenAD_Symbol_321,PRLAP) CALL saxpy(OpenAD_Symbol_331,OpenAD_Symbol_332,PLAP) CALL saxpy(OpenAD_Symbol_334,OpenAD_Symbol_335,PLAP) CALL saxpy(OpenAD_Symbol_337,OpenAD_Symbol_296,PLLAP) CALL saxpy(OpenAD_Symbol_339,OpenAD_Symbol_293,PLLAP) CALL saxpy(OpenAD_Symbol_341,OpenAD_Symbol_301,PBLAP) CALL saxpy(OpenAD_Symbol_343,OpenAD_Symbol_296,PBLAP) CALL saxpy(OpenAD_Symbol_346,OpenAD_Symbol_335,PTLAP) CALL saxpy(OpenAD_Symbol_348,OpenAD_Symbol_332,PRLAP) CALL saxpy(OpenAD_Symbol_350,OpenAD_Symbol_296,PLAP) CALL saxpy(OpenAD_Symbol_352,OpenAD_Symbol_304,PLLAP) CALL saxpy(OpenAD_Symbol_354,OpenAD_Symbol_307,PBLAP) RETURN END SUBROUTINE