IMPLICIT REAL *4(A-H,O-Z) DIMENSION CL1(21),CL2(10),PL1(7),PL2(3) C********************************************************************* C* * C* NAME TOPLAG * C* * C* TYPE MAIN * C* * C* PURPOSE TO COMPUTE TOTAL PHASELAG USING APPROXIMATE FORMULA: * C* * C* TPL =C1+C2S+(C3S+C4)EXP(C5S)+C6/(1+C7S+C8S**4)+2.277/S * C* * C* AT 100 KHZ FREQUENCY FOR THE FOLLOWING GROUND PATH * C* AT RESPECTIVE CONDUCTIVITY AND PERMITTIVITY VALUES: * C* * C* GROUND WAVE CONDUCTIVITY PERMITTIVITY * C* PATH (SIEMENS/METRE) ( ESU ) * C* """"""""""" """""""""""" """""""""""" * C* 1.SEA WATER 1.0 TO 5.5 80 * C* 2.FRESH WATER 0.005 80 * C* 3.LAND SOIL 0.005 TO 1.0 ANY VALUE(3 TO 40) * C* 4.LAND SOIL 0.0001 TO 0.005 3,5,10,15,20,25,30 * C* 5.LAND SOIL 0.00001 TO 0.0001 3,4,5 * C* * C* ====> FOR PATH 4,THE SELECTED CONDUCTIVITY VALUES ARE : * C* 0.0001,0.0002,0.0003,0.0004,0.0005,0.0006,0.0007, * C* 0.0008,0.0009,0.001,0.0011,0.0012,0.0013,0.0014, * C* 0.0015,0.00175,0.002,0.0025,0.003,0.004,0.005. * C* ====> FOR PATH 5,THE SELECTED CONDUCTIVITY VALUES ARE : * C* 0.00001,0.00002,0.00003,0.00004,0.00005,0.00006, * C* 0.00007,0.00008,0.00009,0.00010. * C* * C* AUTHOR TUAN BAHAROM MAHMOOD, SEPTEMBER 1986 * C* * C* EXTERNALS COEF00,COEF03,COEF05,COEF10,COEF15,COEF20,COEF25, * C* COEF30,COEF3L,COEF4L,COEF5L,COEFFW * C* NOTE: SUBROUTINE COEF00 COMPUTES THE COEFFICIENTS * C* """" TO BE USED IN THE APPROXIMATE FORMULA FOR * C* PATH 1 AND 3 ABOVE. * C* OTHER SUBROUTINES STORE THESE COEFFICIENTS * C* AT THIER RESPECTIVE PERMITTIVITY VALUES. * C* * C* PARAMETERS * C* COND = GROUND/WATER CONDUCTIVITY (SIEMENS/METRE) * C* PER = GROUND/WATER PERMITTIVITY (ESU) * C* REF = REFRACTIVITY INDEX (N),E.G: N = 338.0 * C* DIST = DISTANCE (METRES) * C* C1 TO C8 = COEFFICIENTS TO BE USED IN THE FORMULA * C* PPL = PRIMARY PHASELAG * C* SPL = SECONDARY PHASELAG * C* TPL = TOTAL PHASELAG * C* S = DISTANCE/100000.0 * C* N = TO DETERMINE THE COEFFICIENTS TO BE USED * C* WITH RESPECTIVE CONDUCTIVITY AND * C* PERMITTIVITY FOR PATH 4 AND 5 * C* * C* DATA INPUT : COND,PER,REF,DIST * C* OUTPUT : PPL,SPL,TPL * C* * C* LANGUAGE FORTRAN * C* * C* REFERENCE MAHMOOD,T.B(1986) "COMPUTATION OF SECONDARY PHASELAG", * C* MASTER OF ENGINEERING REPORT, DEPT. OF SURVEYING * C* ENGINEERING, UNIVERSITY OF NEW BRUNSWICK. * C* * C********************************************************************* DATA CL1/0.0001,0.0002,0.0003,0.0004,0.0005,0.0006,0.0007,0.0008, & 0.0009,0.001,0.0011,0.0012,0.0013,0.0014,0.0015,0.00175, & 0.002,0.0025,0.003,0.004,0.005/ DATA CL2/0.00001,0.00002,0.00003,0.00004,0.00005,0.00006,0.00007, & 0.00008,0.00009,0.0001/ DATA PL1/3.0,5.0,10.0,15.0,20.0,25.0,30.0/ DATA PL2/3.0,4.0,5.0/ WRITE(6,1000) 100 READ(5,*)COND,PER,REF,DIST IF(COND.LE.0.0) GO TO 400 IF(COND.LE.5.5) GO TO 1 WRITE(6,1002) COND GO TO 300 1 IF(COND.LT.0.005) GO TO 101 CALL COEF00(COND,C1,C2,C3,C4,C5,C6,C7,C8) GO TO 200 101 ID1 = 0 ID2 = 0 IF(COND.EQ.0.005.AND.PER.GT.30.0) GO TO 109 IF(COND.LE.0.0001.AND.PER.LE.5.0) GO TO 105 DO 10 I=1,21 IF(COND.NE.CL1(I)) GO TO 10 ID1 = 1 N = I GO TO 102 10 CONTINUE 102 DO 20 I=1,7 IF(PER.NE.PL1(I)) GO TO 20 ID2 = 1 GO TO 103 20 CONTINUE 103 IF(ID1.EQ.1.AND.ID2.EQ.1) GO TO 104 IF(ID1.EQ.0.AND.ID2.EQ.0) WRITE (6,1001) COND,PER IF(ID1.EQ.0.AND.ID2.EQ.1) WRITE (6,1002) COND IF(ID1.EQ.1.AND.ID2.EQ.0) WRITE (6,1003) PER GO TO 300 104 IF(PER.EQ.3.0) CALL COEF03(N,C1,C2,C3,C4,C5,C6,C7,C8) IF(PER.EQ.5.0) CALL COEF05(N,C1,C2,C3,C4,C5,C6,C7,C8) IF(PER.EQ.10.0) CALL COEF10(N,C1,C2,C3,C4,C5,C6,C7,C8) IF(PER.EQ.15.0) CALL COEF15(N,C1,C2,C3,C4,C5,C6,C7,C8) IF(PER.EQ.20.0) CALL COEF20(N,C1,C2,C3,C4,C5,C6,C7,C8) IF(PER.EQ.25.0) CALL COEF25(N,C1,C2,C3,C4,C5,C6,C7,C8) IF(PER.EQ.30.0) CALL COEF30(N,C1,C2,C3,C4,C5,C6,C7,C8) GO TO 200 105 DO 30 I=1,10 IF(COND.NE.CL2(I)) GO TO 30 ID1 = 1 N = I GO TO 106 30 CONTINUE 106 DO 40 I=1,3 IF(PER.NE.PL2(I)) GO TO 40 ID2 = 1 GO TO 107 40 CONTINUE 107 IF(ID1.EQ.1.AND.ID2.EQ.1) GO TO 108 IF(ID1.EQ.0.AND.ID2.EQ.0) WRITE(6,1001) COND,PER IF(ID1.EQ.0.AND.ID2.EQ.1) WRITE(6,1002) COND IF(ID1.EQ.1.AND.ID2.EQ.0) WRITE(6,1003) PER GO TO 300 108 IF(PER.EQ.3.0) CALL COEF3L(N,C1,C2,C3,C4,C5,C6,C7,C8) IF(PER.EQ.4.0) CALL COEF4L(N,C1,C2,C3,C4,C5,C6,C7,C8) IF(PER.EQ.5.0) CALL COEF5L(N,C1,C2,C3,C4,C5,C6,C7,C8) GO TO 200 109 CALL COEFFW(C1,C2,C3,C4,C5,C6,C7,C8) 200 S = DIST/100000.0 TPL = C1+C2*S+(C3*S+C4)*EXP(C5*S)+C6/(1.0+C7*S+C8*S**4)+2.277/S PPL = (DIST/1000.0)*0.338 SPL = TPL-PPL IF(REF.EQ.338.0) GO TO 110 PPL = (DIST/1000.0)*(REF/1000.0) TPL = SPL+PPL 110 WRITE(6,1004) COND,PER,REF WRITE(6,1005)C1,C2,C3,C4,C5,C6,C7,C8 WRITE(6,1006) DIST,PPL,SPL,TPL 300 GO TO 100 400 WRITE(6,1000) 1000 FORMAT('1',' ') 1001 FORMAT(////10X,'***** E R R O R *****'//10X, & 'YOUR CONDUCTIVITY',F10.6,' AND PERMITTIVITY',F5.1, & ' NOT IN THE LIST'/10X,'*********************************', & '********************************') 1002 FORMAT(////10X,'***** E R R O R *****'//10X, & 'YOUR CONDUCTIVITY',F10.6,' NOT IN THE LIST'/10X, & '*******************************************') 1003 FORMAT(////10X,'***** E R R O R *****'//10X, & 'YOUR PERMITTIVITY',F5.1,' NOT IN THE LIST'/10X, & '**************************************') 1004 FORMAT(////10X,'CONDUCTIVITY =',F9.6,2X,'SIEMENS/METRE'/ & 10X,'PERMITTIVITY =',F5.1,7X,'ESU'/ & 10X,'REFRACTIVITY =',F6.1//) 1005 FORMAT(10X,'THE COEFFICIENTS C(I),I = 1 TO 8 ARE:'// & 10X,4F10.3/10X,4F10.3//) 1006 FORMAT(10X,'DISTANCE(M)',5X,'PRI.PHASELAG(M)', & 5X,'SEC.PHASELAG(M)',5X,'TOT.PHASELAG(M)'//11X,F9.1, & 10X,F6.1,15X,F6.1,13X,F6.1) STOP END SUBROUTINE COEF00(COND,C1,C2,C3,C4,C5,C6,C7,C8) IMPLICIT REAL *4(A-H,O-Z) C********************************************************************* C* * C* NAME COEF00 * C* * C* TYPE SUBROUTINE * C* * C* PURPOSE TO COMPUTE THE C COEFFIENTS TO BE USED IN THE * C* COMPUTATION OF TOTAL PHASELAG USING APPROXIMATE * C* FORMULA FOR SEA WATER AND LAND CONDUCTIVITY ABOVE * C* 0.005 SIEMENS/METRE. * C* * C* AUTHOR TUAN BAHAROM MAHMOOD * C* * C* EXTERNALS NONE * C* * C* PARAMETERS ALL EXPLAINED IN THE MAIN * C* INPUT : COND * C* OUTPUT: C1,C2,C3,C4,C5,C6,C7,C8 * C* * C* CALLING CALL COEF00(COND,C1,C2,C3,C4,C5,C6,C7,C8) * C* * C* LANGUAGE FORTRAN * C* * C* REFERENCE MAHMOOD,T.B(1986) "COMPUTATION OF SECONDARY PHSELAG" * C* M.ENG REPORT,DEPT. OF SURVEYING ENG.,UNB * C* * C********************************************************************* CN = COND*10.0 CN12 = CN**0.5 CN14 = CN**0.25 IF(COND.LT.1.0) GO TO 10 C1 = -111.0 C2 = 97.882-11.437/CN+40.006/CN12-13.576/CN14 C3 = -18.496+63.283/CN+5.844/CN12+7.714/CN14 C4 = 112.8 C5 = -0.254 C6 = 0.0 C7 = 0.0 C8 = 0.0 GO TO 20 10 C1 = -128.059+1.085/CN+65.474/CN12+5.902/CN14 C2 = 97.380-0.615/CN+20.238/CN12-4.107/CN14 C3 = 19.048+0.049/CN-4.812/CN12+22.139/CN14 C4 = 126.620-0.901/CN-38.415/CN12-11.674/CN14 C5 = -0.5273-0.0025/CN+0.0387/CN12-0.0927/CN14 C6 = -4.055-0.024/CN-30.266/CN12+23.080/CN14 C7 = 8.421-0.494/CN+5.677/CN12-5.644/CN14 C8 = 315.570+22.365/CN-346.462/CN12+448.310/CN14 20 RETURN END SUBROUTINE COEF03(N,C1,C2,C3,C4,C5,C6,C7,C8) IMPLICIT REAL *4(A-H,O-Z) C********************************************************************* C* * C* THIS SUBROUTINE STORES THE C COEFFICIENTS FOR THE APPROXIMATE * C* FORMULA AT SELECTED CONDUCTIVITIES (SEE MAIN) * C* ......... AND PERMITTIVITY OF 3 ESU ......... * C* N IS DETERMINED IN THE MAIN TO INDICATE THE POSITION * C* OF THE INPUT CONDUCTIVITY WITHIN THE LIST AND C1 TO C8 ARE * C* RETURNED TO THE MAIN. * C* * C********************************************************************* DIMENSION CC1(21),CC2(21),CC3(21),CC4(21),CC5(21) DIMENSION CC6(21),CC7(21),CC8(21) DATA CC1/985.5,1016.8,976.7,933.1,886.0,838.8,790.4, & 743.1,698.4,656.6,618.9,584.4,553.9,527.2, & 503.7,452.6,410.2,347.1,302.1,240.6,199.7/ DATA CC2/187.47,194.08,200.32,204.65,207.70,209.64,210.77, & 211.19,211.03,210.44,209.49,208.32,206.95,205.45, & 203.85,199.87,196.10,189.22,183.31,173.99,167.01/ DATA CC3/160.7,401.8,335.8,271.2,205.2,158.1,124.5, & 100.8,84.3,73.7,65.3,59.7,55.2,52.5, & 49.7,44.3,42.6,42.0,41.7,42.0,43.2/ DATA CC4/229.7,-285.3,-426.2,-492.0,-475.0,-449.1,-415.5, & -381.3,-350.5,-325.1,-300.6,-279.4,-261.6,-253.0, & -246.6,-223.7,-199.1,-166.8,-142.4,-107.8,-85.4/ DATA CC5/-0.514,-0.632,-0.534,-0.478,-0.436,-0.414,-0.402, & -0.398,-0.399,-0.402,-0.406,-0.410,-0.417,-0.434, & -0.457,-0.495,-0.507,-0.534,-0.556,-0.583,-0.596/ DATA CC6/-1040.9,-605.5,-453.1,-362.3,-338.3,-321.5,-310.1, & -300.0,-288.9,-275.2,-264.3,-253.2,-242.5,-226.8, & -211.9,-187.8,-172.7,-146.2,-128.5,-105.8,-90.0/ DATA CC7/6.62,7.110,8.12,9.61,8.85,8.21,7.61, & 7.16,6.87,6.79,6.68,6.65,6.66,6.98, & 7.36,7.85,8.02,8.60,8.97,9.46,9.97/ DATA CC8/30.0,14.0,15.0,54.0,44.0,34.0,23.0, & 17.0,13.0,13.0,12.0,12.0,12.0,16.0, & 21.0,29.0,33.0,53.0,69.0,94.0,142.0/ C1 = CC1(N) C2 = CC2(N) C3 = CC3(N) C4 = CC4(N) C5 = CC5(N) C6 = CC6(N) C7 = CC7(N) C8 = CC8(N) RETURN END SUBROUTINE COEF05(N,C1,C2,C3,C4,C5,C6,C7,C8) IMPLICIT REAL *4(A-H,O-Z) C********************************************************************* C* * C* THIS SUBROUTINE STORES THE C COEFFICIENTS FOR THE APPROXIMATE * C* FORMULA AT THE SELECTED CONDUCTIVITY AND PERMITTIVITY OF 5 ESU * C* """""""""""" """"" * C********************************************************************* DIMENSION CC1(21),CC2(21),CC3(21),CC4(21),CC5(21) DIMENSION CC6(21),CC7(21),CC8(21) DATA CC1/930.9,984.0,956.6,914.2,873.9,829.1,782.6, & 737.2,693.4,652.5,615.5,581.9,552.1,525.5, & 502.3,451.5,409.5,346.7,301.7,240.5,199.4/ DATA CC2/186.62,193.45,199.52,204.08,206.98,208.95,210.12, & 210.56,210.47,209.94,209.04,207.90,206.56,205.10, & 203.53,199.64,195.91,189.09,183.23,173.94,167.00/ DATA CC3/164.8,359.4,320.3,254.5,202.1,156.9,124.2, & 100.6,84.9,74.4,66.0,60.3,55.5,53.1, & 50.1,44.9,42.9,42.2,42.0,42.1,43.9/ DATA CC4/201.7,-233.5,-398.4,-443.9,-459.8,-437.2,-405.9, & -374.1,-345.3,-320.8,-297.0,-276.7,-259.7,-251.6, & -245.4,-222.7,-198.5,-166.6,-142.4,-107.8,-85.4/ DATA CC5/-0.512,-0.609,-0.532,-0.467,-0.437,-0.416,-0.403, & -0.399,-0.399,-0.402,-0.406,-0.410,-0.417,-0.434, & -0.457,-0.495,-0.507,-0.534,-0.556,-0.584,-0.596/ DATA CC6/-959.2,-621.9,-459.1,-387.3,-340.6,-323.1,-311.4, & -301.0,-288.8,-275.2,-264.4,-253.2,-242.4,-226.5, & -211.7,-187.6,-172.4,-145.9,-128.2,-105.6,-89.6/ DATA CC7/6.92,6.76,7.87,8.42,8.70,8.10,7.52, & 7.09,6.85,6.77,6.66,6.64,6.65,6.99, & 7.36,7.85,7.99,8.59,9.01,9.49,10.03/ DATA CC8/30.0,14.0,14.0,25.0,40.0,31.0,22.0, & 16.0,13.0,13.0,12.0,12.0,12.0,16.0, & 21.0,29.0,33.0,54.0,71.0,95.0,145.0/ C1 = CC1(N) C2 = CC2(N) C3 = CC3(N) C4 = CC4(N) C5 = CC5(N) C6 = CC6(N) C7 = CC7(N) C8 = CC8(N) RETURN END SUBROUTINE COEF10(N,C1,C2,C3,C4,C5,C6,C7,C8) IMPLICIT REAL *4(A-H,O-Z) C********************************************************************* C* * C* THIS SUBROUTINE STORES THE C COEFFICIENTS FOR THE APPROXIMATE * C* FORMULA AT THE SELECTED CONDUCTIVITY AND PERMITTIVITY OF 10 ESU * C* """""""""""" """""" * C********************************************************************* DIMENSION CC1(21),CC2(21),CC3(21),CC4(21),CC5(21) DIMENSION CC6(21),CC7(21),CC8(21) DATA CC1/811.6,908.4,906.2,877.9,844.9,805.3,763.2, & 721.4,680.9,642.5,607.6,575.6,546.9,521.2, & 498.9,449.3,408.0,345.7,301.2,240.2,199.1/ DATA CC2/184.68,191.76,197.66,202.15,205.14,207.24,208.52, & 209.08,209.10,208.68,207.88,206.85,205.61,204.24, & 202.74,199.01,195.41,188.76,182.99,173.81,166.92/ DATA CC3/166.1,293.9,284.1,237.0,194.0,153.6,123.4, & 101.6,86.3,75.9,67.5,61.4,56.7,54.3, & 51.1,45.5,43.2,42.7,42.3,42.3,44.1/ DATA CC4/154.7,-161.0,-333.7,-397.9,-425.1,-409.8,-384.1, & -358.0,-333.0,-310.8,-289.7,-270.6,-255.0,-248.3, & -242.5,-220.6,-196.8,-165.9,-142.2,-107.8,-85.5/ DATA CC5/-0.507,-0.574,-0.523,-0.469,-0.442,-0.419,-0.405, & -0.400,-0.399,-0.401,-0.405,-0.409,-0.417,-0.434, & -0.457,-0.495,-0.507,-0.534,-0.556,-0.585,-0.597/ DATA CC6/-798.3,-618.6,-470.2,-394.4,-344.6,-325.6,-313.0, & -300.7,-288.2,-274.8,-263.5,-252.7,-241.7,-225.3, & -211.0,-187.2,-172.5,-145.6,-127.8,-105.4,-89.2/ DATA CC7/7.31,6.64,7.38,7.97,8.38,7.89,7.38, & 7.03,6.82,6.74,6.66,6.62,6.66,7.02, & 7.37,7.83,7.98,8.63,9.05,9.54,10.07/ DATA CC8/28.0,13.0,13.0,20.0,32.0,26.0,19.0, & 15.0,13.0,13.0,12.0,12.0,12.0,16.0, & 21.0,29.0,33.0,55.0,73.0,96.0,149.0/ C1 = CC1(N) C2 = CC2(N) C3 = CC3(N) C4 = CC4(N) C5 = CC5(N) C6 = CC6(N) C7 = CC7(N) C8 = CC8(N) RETURN END SUBROUTINE COEF15(N,C1,C2,C3,C4,C5,C6,C7,C8) IMPLICIT REAL *4(A-H,O-Z) C********************************************************************* C* * C* THIS SUBROUTINE STORES THE C COEFFICIENTS FOR THE APPROXIMATE * C* FORMULA AT THE SELECTED CONDUCTIVITY AND PERMITTIVITY OF 15 ESU * C* """""""""""" """"" * C********************************************************************* DIMENSION CC1(21),CC2(21),CC3(21),CC4(21),CC5(21) DIMENSION CC6(21),CC7(21),CC8(21) DATA CC1/717.0,839.6,857.6,842.2,816.5,782.2,744.4, & 705.9,668.2,632.3,598.8,568.2,541.0,516.8, & 495.2,446.8,406.2,344.7,300.5,239.9,198.9/ DATA CC2/183.06,190.14,195.91,200.32,203.38,205.56,206.94, & 207.62,207.75,207.44,206.78,205.85,204.69,203.38, & 201.96,198.39,194.91,188.43,182.76,173.68,166.84/ DATA CC3/162.6,260.6,257.9,222.3,187.1,150.5,122.0, & 101.7,87.4,77.2,69.5,63.5,58.7,55.3, & 52.0,46.3,43.7,43.1,42.6,42.4,44.2/ DATA CC4/124.5,-135.9,-289.9,-360.3,-396.0,-385.9,-364.1, & -341.6,-320.7,-301.1,-282.2,-264.8,-251.4,-244.4, & -239.2,-218.6,-195.3,-165.1,-141.6,-107.4,-85.3/ DATA CC5/-0.501,-0.555,-0.517,-0.472,-0.447,-0.424,-0.409, & -0.401,-0.400,-0.401,-0.403,-0.407,-0.416,-0.434, & -0.457,-0.495,-0.507,-0.534,-0.556,-0.585,-0.598/ DATA CC6/-681.9,-577.8,-464.9,-394.9,-344.2,-325.6,-313.5, & -300.9,-287.4,-274.0,-262.0,-251.0,-239.4,-224.7, & -210.5,-186.7,-172.2,-145.3,-127.7,-105.3,-89.2/ DATA CC7/7.38,6.91,7.28,7.74,8.21,7.77,7.27, & 6.95,6.79,6.73,6.68,6.66,6.73,7.02, & 7.36,7.84,7.97,8.63,9.04,9.51,10.10/ DATA CC8/24.0,16.0,14.0,18.0,30.0,24.0,18.0, & 14.0,13.0,13.0,13.0,13.0,13.0,17.0, & 21.0,30.0,33.0,55.0,73.0,96.0,149.0/ C1 = CC1(N) C2 = CC2(N) C3 = CC3(N) C4 = CC4(N) C5 = CC5(N) C6 = CC6(N) C7 = CC7(N) C8 = CC8(N) RETURN END SUBROUTINE COEF20(N,C1,C2,C3,C4,C5,C6,C7,C8) IMPLICIT REAL *4(A-H,O-Z) C********************************************************************* C* * C* THIS SUBROUTINE STORES THE C COEFFICIENTS FOR THE APPROXIMATE * C* FORMULA AT THE SELECTED CONDUCTIVITY AND PERMITTIVITY OF 20 ESU * C* """""""""""" """""" * C********************************************************************* DIMENSION CC1(21),CC2(21),CC3(21),CC4(21),CC5(21) DIMENSION CC6(21),CC7(21),CC8(21) DATA CC1/643.3,777.7,810.8,806.9,788.5,758.9,724.9, & 689.9,655.5,622.1,590.6,561.5,535.5,512.4, & 491.4,444.4,404.4,343.6,299.7,239.5,198.7/ DATA CC2/181.72,188.62,194.28,198.61,201.69,203.97,205.45, & 206.24,206.45,206.22,205.66,204.83,203.76,202.52, & 201.18,197.77,194.42,188.10,182.53,173.55,166.76/ DATA CC3/157.5,234.7,234.3,209.5,179.8,146.8,121.1, & 102.3,88.3,78.2,70.5,64.5,59.5,56.1, & 52.9,46.9,44.4,43.4,42.9,42.6,44.3/ DATA CC4/103.0,-116.5,-249.7,-327.6,-367.4,-361.0,-344.2, & -326.4,-308.7,-291.5,-274.6,-258.7,-246.3,-240.3, & -235.9,-216.5,-193.7,-164.1,-140.9,-107.1,-85.2/ DATA CC5/-0.494,-0.540,-0.509,-0.473,-0.451,-0.427,-0.410, & -0.402,-0.400,-0.401,-0.404,-0.407,-0.416,-0.434, & -0.457,-0.495,-0.507,-0.534,-0.556,-0.585,-0.598/ DATA CC6/-596.0,-538.5,-458.0,-391.6,-343.9,-326.2,-313.3, & -299.7,-286.2,-273.0,-261.0,-250.2,-238.8,-224.0, & -209.8,-186.2,-171.7,-145.1,-127.5,-105.2,-89.1/ DATA CC7/7.30,7.09,7.14,7.61,8.04,7.59,7.17, & 6.91,6.77,6.70,6.67,6.66,6.72,7.01, & 7.36,7.84,7.97,8.62,9.03,9.51,10.09/ DATA CC8/21.0,18.0,14.0,18.0,27.0,21.0,17.0, & 14.0,13.0,13.0,13.0,13.0,13.0,17.0, & 21.0,30.0,33.0,55.0,73.0,97.0,149.0/ C1 = CC1(N) C2 = CC2(N) C3 = CC3(N) C4 = CC4(N) C5 = CC5(N) C6 = CC6(N) C7 = CC7(N) C8 = CC8(N) RETURN END SUBROUTINE COEF25(N,C1,C2,C3,C4,C5,C6,C7,C8) IMPLICIT REAL *4(A-H,O-Z) C********************************************************************* C* * C* THIS SUBROUTINE STORES THE C COEFFICIENTS FOR THE APPROXIMATE * C* FORMULA AT THE SELECTED CONDUCTIVITY AND PERMITTIVITY OF 25 ESU * C* """""""""""" """""" * C********************************************************************* DIMENSION CC1(21),CC2(21),CC3(21),CC4(21),CC5(21) DIMENSION CC6(21),CC7(21),CC8(21) DATA CC1/585.9,722.4,767.9,774.4,762.7,738.8,709.0, & 676.2,644.1,613.1,583.5,555.8,530.2,508.1, & 487.6,441.7,402.5,342.5,299.0,239.0,198.4/ DATA CC2/180.58,187.23,192.70,196.90,199.99,202.28,203.84, & 204.77,205.09,204.96,204.49,203.75,202.81,201.63, & 200.40,197.15,193.91,187.76,182.30,173.42,166.68/ DATA CC3/153.6,215.2,218.6,203.6,176.3,146.4,121.9, & 102.4,88.3,78.1,70.4,64.3,59.9,56.6, & 53.6,47.8,45.0,43.8,43.1,42.8,44.4/ DATA CC4/84.7,-102.8,-226.9,-314.0,-352.7,-352.8,-339.5, & -317.4,-299.4,-283.8,-268.2,-253.2,-241.1,-236.6, & -232.3,-214.6,-192.3,-163.0,-140.1,-106.6,-84.7/ DATA CC5/-0.490,-0.528,-0.506,-0.482,-0.461,-0.439,-0.423, & -0.410,-0.405,-0.406,-0.408,-0.412,-0.417,-0.435, & -0.457,-0.495,-0.506,-0.534,-0.556,-0.584,-0.597/ DATA CC6/-529.5,-500.4,-439.5,-373.9,-333.2,-314.9,-302.6, & -294.9,-284.0,-271.5,-260.0,-249.6,-238.4,-223.3, & -209.3,-185.2,-171.2,-144.9,-127.5,-105.2,-89.2/ DATA CC7/7.28,7.23,7.27,7.96,8.25,7.90,7.48, & 7.00,6.78,6.69,6.65,6.62,6.69,7.01, & 7.35,7.86,7.98,8.61,9.02,9.51,10.06/ DATA CC8/20.0,19.0,16.0,25.0,32.0,28.0,22.0, & 15.0,13.0,12.0,12.0,12.0,13.0,17.0, & 21.0,31.0,33.0,54.0,71.0,95.0,145.0/ C1 = CC1(N) C2 = CC2(N) C3 = CC3(N) C4 = CC4(N) C5 = CC5(N) C6 = CC6(N) C7 = CC7(N) C8 = CC8(N) RETURN END SUBROUTINE COEF30(N,C1,C2,C3,C4,C5,C6,C7,C8) IMPLICIT REAL *4(A-H,O-Z) C********************************************************************* C* * C* THIS SUBROUTINE STORES THE C COEFFICIENTS FOR THE APPROXIMATE * C* FORMULA AT THE SELECTED CONDUCTIVITY AND PERMITTIVITY OF 30 ESU * C* """""""""""" """""" * C********************************************************************* DIMENSION CC1(21),CC2(21),CC3(21),CC4(21),CC5(21) DIMENSION CC6(21),CC7(21),CC8(21) DATA CC1/540.3,673.8,727.4,741.8,736.3,717.2,691.6, & 661.6,632.0,602.7,575.1,548.7,524.5,502.9, & 483.4,439.0,400.4,341.4,298.3,238.5,198.2/ DATA CC2/179.57,185.93,191.23,195.35,198.43,200.76,202.36, & 203.39,203.81,203.79,203.40,202.76,201.89,200.83, & 199.64,196.53,193.43,187.42,182.05,173.29,166.60/ DATA CC3/148.7,199.7,205.2,193.9,171.2,144.3,121.9, & 102.8,88.7,78.7,71.0,65.1,60.4,57.9, & 54.5,48.6,45.8,44.0,43.1,43.0,44.5/ DATA CC4/71.4,-92.5,-207.0,-290.3,-332.1,-336.8,-328.7, & -307.4,-289.3,-274.3,-260.3,-246.8,-235.7,-232.7, & -228.8,-212.5,-190.6,-161.9,-139.2,-106.1,-84.3/ DATA CC5/-0.486,-0.519,-0.503,-0.483,-0.466,-0.446,-0.431, & -0.415,-0.408,-0.407,-0.409,-0.412,-0.418,-0.436, & -0.457,-0.495,-0.506,-0.534,-0.556,-0.583,-0.596/ DATA CC6/-479.1,-465.9,-420.5,-365.3,-327.3,-309.3,-296.1, & -290.5,-281.8,-270.5,-259.3,-248.8,-237.9,-221.9, & -208.5,-184.6,-170.6,-144.9,-127.6,-105.1,-89.2/ DATA CC7/7.20,7.33,7.37,7.95,8.26,7.96,7.62, & 7.09,6.79,6.68,6.63,6.61,6.66,7.03, & 7.36,7.88,7.99,8.59,8.98,9.48,10.00/ DATA CC8/18.0,20.0,18.0,26.0,34.0,30.0,24.0, & 17.0,13.0,12.0,12.0,12.0,12.0,17.0, & 21.0,31.0,34.0,53.0,70.0,95.0,143.0/ C1 = CC1(N) C2 = CC2(N) C3 = CC3(N) C4 = CC4(N) C5 = CC5(N) C6 = CC6(N) C7 = CC7(N) C8 = CC8(N) RETURN END SUBROUTINE COEF3L(N,C1,C2,C3,C4,C5,C6,C7,C8) IMPLICIT REAL *4(A-H,O-Z) C********************************************************************* C* * C* THIS SUBROUTINE STORES THE C COEFFICIENTS FOR THE APPROXIMATE * C* FORMULA AT THE SELECTED CONDUCTIVITY AND PERMITTIVITY OF 3 ESU * C* """""""""""" """"" * C********************************************************************* DIMENSION CC1(10),CC2(10),CC3(10),CC4(10),CC5(10) DIMENSION CC6(10),CC7(10),CC8(10) DATA CC1/533.9,697.8,803.9,867.2,901.0,925.9,946.0,961.9, & 975.3,985.9/ DATA CC2/174.44,176.50,178.34,180.10,182.02,183.53,184.73, & 185.75,186.64,187.45/ DATA CC3/135.0,151.2,165.4,164.4,126.5,117.8,120.0,129.1, & 144.1,162.0/ DATA CC4/274.8,260.7,246.8,245.3,276.7,286.4,283.2,271.9, & 252.7,228.2/ DATA CC5/-0.517,-0.538,-0.554,-0.547,-0.492,-0.477,-0.477, & -0.484,-0.500,-0.516/ DATA CC6/-513.8,-707.1,-827.4,-906.1,-979.6,-1019.6,-1041.3, & -1050.6,-1049.3,-1039.7/ DATA CC7/27.00,24.43,19.64,15.68,12.36,10.28,8.88,7.88,7.17,6.63/ DATA CC8/3386.0,2621.0,1485.0,685.0,258.0,136.0,84.0,57.0, & 41.0,30.0/ C1 = CC1(N) C2 = CC2(N) C3 = CC3(N) C4 = CC4(N) C5 = CC5(N) C6 = CC6(N) C7 = CC7(N) C8 = CC8(N) RETURN END SUBROUTINE COEF4L(N,C1,C2,C3,C4,C5,C6,C7,C8) IMPLICIT REAL *4(A-H,O-Z) C********************************************************************* C* * C* THIS SUBROUTINE STORES THE C COEFFICIENTS FOR THE APPROXIMATE * C* FORMULA AT THE SELECTED CONDUCTIVITY AND PERMITTIVITY OF 4 ESU * C* """"""""""""""" * C********************************************************************* DIMENSION CC1(10),CC2(10),CC3(10),CC4(10),CC5(10) DIMENSION CC6(10),CC7(10),CC8(10) DATA CC1/513.8,649.6,749.9,814.4,853.7,885.6,910.5,930.4, & 946.2,959.1/ DATA CC2/174.20,176.09,177.88,179.74,181.64,183.03,184.20, & 185.22,186.13,186.97/ DATA CC3/146.7,160.7,171.4,157.0,126.2,121.4,126.3,137.5, & 151.5,168.4/ DATA CC4/255.9,242.6,231.2,241.4,269.5,273.1,266.6,252.1, & 232.9,208.5/ DATA CC5/-0.531,-0.547,-0.557,-0.533,-0.487,-0.478,-0.481, & -0.492,-0.505,-0.520/ DATA CC6/-494.2,-644.3,-757.3,-847.8,-922.9,-964.7,-988.7, & -999.3,-1000.8,-993.9/ DATA CC7/23.17,21.62,18.44,15.04,12.10,10.26,8.98,8.05,7.35,6.82/ DATA CC8/2197.0,1728.0,1072.0,479.0,203.0,119.0,79.0,56.0, & 41.0,31.0/ C1 = CC1(N) C2 = CC2(N) C3 = CC3(N) C4 = CC4(N) C5 = CC5(N) C6 = CC6(N) C7 = CC7(N) C8 = CC8(N) RETURN END SUBROUTINE COEF5L(N,C1,C2,C3,C4,C5,C6,C7,C8) IMPLICIT REAL *4(A-H,O-Z) C********************************************************************* C* * C* THIS SUBROUTINE STORES THE C COEFFICIENTS FOR THE APPROXIMATE * C* FORMULA AT THE SELECTED CONDUCTIVITY AND PERMITTIVITY OF 5 ESU * C* """"""""""""""""" * C********************************************************************* DIMENSION CC1(10),CC2(10),CC3(10),CC4(10),CC5(10) DIMENSION CC6(10),CC7(10),CC8(10) DATA CC1/495.0,612.0,703.5,757.3,811.2,848.0,877.0,899.8, & 918.2,932.9/ DATA CC2/174.03,175.80,177.61,180.02,181.30,182.58,183.71, & 184.72,185.65,186.51/ DATA CC3/153.3,164.1,163.6,100.3,122.4,124.0,131.6,142.8, & 157.4,172.8/ DATA CC4/241.7,230.2,228.1,285.4,264.9,261.6,251.2,235.8, & 215.2,191.9/ DATA CC5/-0.537,-0.549,-0.543,-0.449,-0.479,-0.478,-0.486, & -0.496,-0.509,-0.522/ DATA CC6/-477.9,-602.2,-710.2,-830.5,-874.5,-915.0,-940.0, & -952.8,-955.8,-951.9/ DATA CC7/20.68,19.55,17.06,13.29,11.71,10.17,9.03,8.17,7.51,6.99/ DATA CC8/1509.0,1165.0,682.0,191.0,157.0,104.0,73.0,54.0, & 41.0,32.0/ C1 = CC1(N) C2 = CC2(N) C3 = CC3(N) C4 = CC4(N) C5 = CC5(N) C6 = CC6(N) C7 = CC7(N) C8 = CC8(N) RETURN END SUBROUTINE COEFFW(C1,C2,C3,C4,C5,C6,C7,C8) IMPLICIT REAL *4(A-H,O-Z) C********************************************************************* C* * C* THIS SUBROUTINE STORES THE C COEFFICIENTS FOR THE APPROXIMATE * C* FORMULA FOR FRESH WATER,CONDUCTIVITY AND PERMITTIVITY OF * C* 0.005 SIEMENS/METRE AND 80 ESU * C* """""""""""""""""""""""""""""" * C********************************************************************* C1 = 195.0 C2 = 165.72 C3 = 46.0 C4 = -82.4 C5 = 0.600 C6 = -87.8 C7 = 10.13 C8 = 150.0 RETURN END