diff options
Diffstat (limited to 'src/Derivatives_8_4_min_err_coeff.F90')
| -rw-r--r-- | src/Derivatives_8_4_min_err_coeff.F90 | 158 |
1 files changed, 79 insertions, 79 deletions
diff --git a/src/Derivatives_8_4_min_err_coeff.F90 b/src/Derivatives_8_4_min_err_coeff.F90 index edff755..0eb43fa 100644 --- a/src/Derivatives_8_4_min_err_coeff.F90 +++ b/src/Derivatives_8_4_min_err_coeff.F90 @@ -39,82 +39,82 @@ subroutine deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, dvar ) q(6,1) = 0.1564962758424078514454013191006671546757_wp q(7,1) = 0.003067629343740613524557473056482141287158_wp q(8,1) = -0.03730827676416971624344442423120269215938_wp - q(9,1) = 0_wp - q(10,1) = 0_wp - q(11,1) = 0_wp - q(12,1) = 0_wp + q(9,1) = 0.0_wp + q(10,1) = 0.0_wp + q(11,1) = 0.0_wp + q(12,1) = 0.0_wp q(1,2) = -0.4338209217401506177055540526837828066711_wp - q(2,2) = 0_wp + q(2,2) = 0.0_wp q(3,2) = 0.06631961192627293470305188523936353462439_wp q(4,2) = 0.6021342012739069461791324177355805603656_wp q(5,2) = -0.2213339967864188453859900127368926875520_wp q(6,2) = -0.001759980003751516130708214237001705688807_wp q(7,2) = -0.03586646217917856587352426166797879385310_wp q(8,2) = 0.02432754750931966421359223835071189877503_wp - q(9,2) = 0_wp - q(10,2) = 0_wp - q(11,2) = 0_wp - q(12,2) = 0_wp + q(9,2) = 0.0_wp + q(10,2) = 0.0_wp + q(11,2) = 0.0_wp + q(12,2) = 0.0_wp q(1,3) = 0.002477771724790106958560398469919805908522_wp q(2,3) = -0.3930241559935857537756812928554282732813_wp - q(3,3) = 0_wp + q(3,3) = 0.0_wp q(4,3) = -0.6080263371305341035810612202824357742461_wp q(5,3) = 2.372533124692864412954923458457782022478_wp q(6,3) = -2.126112217612194556872952028715623246817_wp q(7,3) = 0.9075309435559160631495435310307745215131_wp q(8,3) = -0.1553791292372561688333328461049890555549_wp - q(9,3) = 0_wp - q(10,3) = 0_wp - q(11,3) = 0_wp - q(12,3) = 0_wp + q(9,3) = 0.0_wp + q(10,3) = 0.0_wp + q(11,3) = 0.0_wp + q(12,3) = 0.0_wp q(1,4) = 0.1470043328582935680963279007411074464266_wp q(2,4) = -0.5109179516689331400658365102850902594996_wp q(3,4) = 0.08705685833207777685654755900085198684605_wp - q(4,4) = 0_wp + q(4,4) = 0.0_wp q(5,4) = -0.1306945333157231839535552435788086367322_wp q(6,4) = 0.5981933016362239542598518794410035592720_wp q(7,4) = -0.2031099149801233772210520666270517105660_wp q(8,4) = 0.01246790713818440202771648130798761425319_wp - q(9,4) = 0_wp - q(10,4) = 0_wp - q(11,4) = 0_wp - q(12,4) = 0_wp + q(9,4) = 0.0_wp + q(10,4) = 0.0_wp + q(11,4) = 0.0_wp + q(12,4) = 0.0_wp q(1,5) = -0.1624073990846984662267508265053107632238_wp q(2,5) = 0.8182390368133059538132603839842499379861_wp q(3,5) = -1.480018301574606037840376638130713134769_wp q(4,5) = 0.5694185675497882973361521309100738568707_wp - q(5,5) = 0_wp + q(5,5) = 0.0_wp q(6,5) = 0.02929508293831017534661096502884570435894_wp q(7,5) = 0.2825910005984090924748504462528740123097_wp q(8,5) = -0.04846434332860791750685068920563907354793_wp q(9,5) = -0.008653643911901097396895772334380539984709_wp - q(10,5) = 0_wp - q(11,5) = 0_wp - q(12,5) = 0_wp + q(10,5) = 0.0_wp + q(11,5) = 0.0_wp + q(12,5) = 0.0_wp q(1,6) = -0.03609686950604370828405582087760384548297_wp q(2,6) = 0.002100328577309696555612805397001149641337_wp q(3,6) = 0.4281425817419204360479874008765051769418_wp q(4,6) = -0.8413238238875046736839550701643761633826_wp q(5,6) = -0.009456755727630000971085851751473603919416_wp - q(6,6) = 0_wp + q(6,6) = 0.0_wp q(7,6) = 0.5172389789041733771234394082779034212042_wp q(8,6) = -0.08760813565107717873901807278223424202774_wp q(9,6) = 0.02979718129528502284256573906127239395909_wp q(10,6) = -0.002793485746432970891490538036994286933664_wp - q(11,6) = 0_wp - q(12,6) = 0_wp + q(11,6) = 0.0_wp + q(12,6) = 0.0_wp q(1,7) = -0.0009797678134978722516935350416937185004514_wp q(2,7) = 0.05926834509975463070197111976550139351983_wp q(3,7) = -0.2530570463899371286637621744353665227114_wp q(4,7) = 0.3955555826583941989223787901885502942688_wp q(5,7) = -0.1263166266018192756531792392597193430222_wp q(6,7) = -0.7162192643577544899896537844517283192029_wp - q(7,7) = 0_wp + q(7,7) = 0.0_wp q(8,7) = 0.8209721963136350137518635528607308559584_wp q(9,7) = -0.2166153552278720352907291696202456084323_wp q(10,7) = 0.04126006766245181624585317516576106827282_wp q(11,7) = -0.003868131343354857773048735171790100150577_wp - q(12,7) = 0_wp + q(12,7) = 0.0_wp q(1,8) = 0.01090012273307913185972171872891927027272_wp q(2,8) = -0.03677379943661633440187210478144113487019_wp q(3,8) = 0.03963287609450569641558898819403594845401_wp @@ -122,7 +122,7 @@ subroutine deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, dvar ) q(5,8) = 0.01981665906734246925389947980969432794418_wp q(6,8) = 0.1109698768905366566973791789571759073485_wp q(7,8) = -0.7509903604688148129224205834189298636300_wp - q(8,8) = 0_wp + q(8,8) = 0.0_wp q(9,8) = 0.7926019635554773751160111698352821779736_wp q(10,8) = -0.1981504908888693437790027924588205444934_wp q(11,8) = 0.03774295064549892262457196046834677037969_wp @@ -465,7 +465,7 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva a1(1) = 0.6000000000000000000000000000000000000000_wp a1(2) = -0.1000000000000000000000000000000000000000_wp a1(3) = 0.009523809523809523809523809523809523809524_wp - a1(4) = 0_wp + a1(4) = 0.0_wp q1(1,1) = -1.683432578685956376544499542432311020539_wp q1(2,1) = 2.206194911574014373726323572790957192479_wp @@ -475,10 +475,10 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q1(6,1) = 0.3015760468426855559502835132688146245819_wp q1(7,1) = -0.05430139927462465206917597295687512552949_wp q1(8,1) = -0.02785725195806647467919208850462676788091_wp - q1(9,1) = 0_wp - q1(10,1) = 0_wp - q1(11,1) = 0_wp - q1(12,1) = 0_wp + q1(9,1) = 0.0_wp + q1(10,1) = 0.0_wp + q1(11,1) = 0.0_wp + q1(12,1) = 0.0_wp q1(1,2) = -0.4451389901846271204230470438663315357999_wp q1(2,2) = 0.03979384217746298000464189244717138065595_wp q1(3,2) = 0.04117629055092890361160914968898310705568_wp @@ -487,10 +487,10 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q1(6,2) = -0.1571076966553020594965523100929511511562_wp q1(7,2) = 0.02929412303865500620649355002594213960776_wp q1(8,2) = 0.01325024192729698663163651603043611635147_wp - q1(9,2) = 0_wp - q1(10,2) = 0_wp - q1(11,2) = 0_wp - q1(12,2) = 0_wp + q1(9,2) = 0.0_wp + q1(10,2) = 0.0_wp + q1(11,2) = 0.0_wp + q1(12,2) = 0.0_wp q1(1,3) = 0.1544592597396873173299548977341719367844_wp q1(2,3) = -1.020860666961436408391054234781341175626_wp q1(3,3) = 0.7352247001589518809305910119936419247628_wp @@ -499,10 +499,10 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q1(6,3) = -0.1256847607545176255327746274725689708913_wp q1(7,3) = -0.01218362768224537680217842793138395523727_wp q1(8,3) = 0.01047289768141453270035870286867846966855_wp - q1(9,3) = 0_wp - q1(10,3) = 0_wp - q1(11,3) = 0_wp - q1(12,3) = 0_wp + q1(9,3) = 0.0_wp + q1(10,3) = 0.0_wp + q1(11,3) = 0.0_wp + q1(12,3) = 0.0_wp q1(1,4) = 0.1147228529148158318034738218073197976023_wp q1(2,4) = -0.3538780254601914958115698874206531961307_wp q1(3,4) = -0.1858994785495101740456777336107169678288_wp @@ -511,10 +511,10 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q1(6,4) = 0.3165092134256211709442618752193670178495_wp q1(7,4) = -0.05069213274750316093259254157360772183295_wp q1(8,4) = -0.01753939629032033624782329757829768029059_wp - q1(9,4) = 0_wp - q1(10,4) = 0_wp - q1(11,4) = 0_wp - q1(12,4) = 0_wp + q1(9,4) = 0.0_wp + q1(10,4) = 0.0_wp + q1(11,4) = 0.0_wp + q1(12,4) = 0.0_wp q1(1,5) = -0.01395722952056041916292538026746562805620_wp q1(2,5) = 0.01429909551156929498323180030121632055909_wp q1(3,5) = 0.2332633686362482363893930312775209255800_wp @@ -523,10 +523,10 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q1(6,5) = 0.2825179352075121763892365873707254541218_wp q1(7,5) = 0.03248906901967492013975174489420235738435_wp q1(8,5) = -0.01206707005235509052833799543828477419916_wp - q1(9,5) = 0_wp - q1(10,5) = 0_wp - q1(11,5) = 0_wp - q1(12,5) = 0_wp + q1(9,5) = 0.0_wp + q1(10,5) = 0.0_wp + q1(11,5) = 0.0_wp + q1(12,5) = 0.0_wp q1(1,6) = -0.06956044896551484587496684695381321009857_wp q1(2,6) = 0.1930764794969705526456194423505107491362_wp q1(3,6) = -0.01938619631973248083178041998353003089641_wp @@ -536,9 +536,9 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q1(7,6) = 0.4088327186474708718613662440475307694516_wp q1(8,6) = -0.02316927663483720615835459462174938908831_wp q1(9,6) = 0.007449295323821255710641434765318098489772_wp - q1(10,6) = 0_wp - q1(11,6) = 0_wp - q1(12,6) = 0_wp + q1(10,6) = 0.0_wp + q1(11,6) = 0.0_wp + q1(12,6) = 0.0_wp q1(1,7) = 0.01734328280101122760645584447264624051821_wp q1(2,7) = -0.04840773491893482968357096782164797435370_wp q1(3,7) = 0.01113356043999974094579726795398111697087_wp @@ -549,8 +549,8 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q1(8,7) = 0.6075690064749917706800934212435751620821_wp q1(9,7) = -0.1083076776139360176453645848101228042161_wp q1(10,7) = 0.01031501691561295406146329379144026706820_wp - q1(11,7) = 0_wp - q1(12,7) = 0_wp + q1(11,7) = 0.0_wp + q1(12,7) = 0.0_wp q1(1,8) = 0.008138876723484902470545040666116789759842_wp q1(2,8) = -0.02002921745130282774392822052620050919550_wp q1(3,8) = -0.002671343688148396036480262140129972226295_wp @@ -562,9 +562,9 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q1(9,8) = 0.5944514726666080313370083773764616334802_wp q1(10,8) = -0.09907524544443467188950139622941027224670_wp q1(11,8) = 0.009435737661374730656142990117086692594923_wp - q1(12,8) = 0_wp + q1(12,8) = 0.0_wp - a2(-4) = 0_wp + a2(-4) = 0.0_wp a2(-3) = -0.009523809523809523809523809523809523809524_wp a2(-2) = 0.1000000000000000000000000000000000000000_wp a2(-1) = -0.6000000000000000000000000000000000000000_wp @@ -582,10 +582,10 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q2(6,1) = 0.3015760468426855559502835132688146245819_wp q2(7,1) = -0.05430139927462465206917597295687512552949_wp q2(8,1) = -0.02785725195806647467919208850462676788091_wp - q2(9,1) = 0_wp - q2(10,1) = 0_wp - q2(11,1) = 0_wp - q2(12,1) = 0_wp + q2(9,1) = 0.0_wp + q2(10,1) = 0.0_wp + q2(11,1) = 0.0_wp + q2(12,1) = 0.0_wp q2(1,2) = -0.4264124762187621886561567415382508860795_wp q2(2,2) = -0.03979384217746298000464189244717138065595_wp q2(3,2) = 0.1722618883119834259798412659855476550988_wp @@ -594,10 +594,10 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q2(6,2) = -0.1617893251467682924382748856749713135863_wp q2(7,2) = 0.02929412303865500620649355002594213960776_wp q2(8,2) = 0.01325024192729698663163651603043611635147_wp - q2(9,2) = 0_wp - q2(10,2) = 0_wp - q2(11,2) = 0_wp - q2(12,2) = 0_wp + q2(9,2) = 0.0_wp + q2(10,2) = 0.0_wp + q2(11,2) = 0.0_wp + q2(12,2) = 0.0_wp q2(1,3) = -0.01200671010762254250338646347193566882227_wp q2(2,3) = -0.2440194743406570625021278824861723494614_wp q2(3,3) = -0.7352247001589518809305910119936419247628_wp @@ -606,10 +606,10 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q2(6,3) = 0.09626986570856218757834718746890783658430_wp q2(7,3) = -0.03992795599013035344106865479906855617172_wp q2(8,3) = 0.01047289768141453270035870286867846966855_wp - q2(9,3) = 0_wp - q2(10,3) = 0_wp - q2(11,3) = 0_wp - q2(12,3) = 0_wp + q2(9,3) = 0.0_wp + q2(10,3) = 0.0_wp + q2(11,3) = 0.0_wp + q2(12,3) = 0.0_wp q2(1,4) = 0.1306125202632199301974189512671891155098_wp q2(2,4) = -0.4412711958764140369782680994499344446221_wp q2(3,4) = 0.02066619697974310507560894936758416496903_wp @@ -618,10 +618,10 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q2(6,4) = 0.2052815419867924821866459690002817924969_wp q2(7,4) = -0.01891279805069496414470228265386908601791_wp q2(8,4) = -0.02151181312742136084630957994326500976747_wp - q2(9,4) = 0_wp - q2(10,4) = 0_wp - q2(11,4) = 0_wp - q2(12,4) = 0_wp + q2(9,4) = 0.0_wp + q2(10,4) = 0.0_wp + q2(11,4) = 0.0_wp + q2(12,4) = 0.0_wp q2(1,5) = -0.03126451734436261395671692493622670802562_wp q2(2,5) = 0.1527573981019868533335641576513049603144_wp q2(3,5) = -0.2513406904302132178367702194477893135637_wp @@ -631,9 +631,9 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q2(7,5) = -0.4521149900467865340864115058311078817594_wp q2(8,5) = 0.1263912325380624678219943619118038655562_wp q2(9,5) = -0.01730728782380219479379154466876107996942_wp - q2(10,5) = 0_wp - q2(11,5) = 0_wp - q2(12,5) = 0_wp + q2(10,5) = 0.0_wp + q2(11,5) = 0.0_wp + q2(12,5) = 0.0_wp q2(1,6) = -0.06956044896551484587496684695381321009857_wp q2(2,6) = 0.1874895080041046108626383662765221752689_wp q2(3,6) = 0.02530957562319505343206818860837856004222_wp @@ -644,8 +644,8 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q2(8,6) = -0.1796044784350835760818247246934294573735_wp q2(9,6) = 0.05214506726674878997449004335722668942840_wp q2(10,6) = -0.005586971492865941782981076073988573867329_wp - q2(11,6) = 0_wp - q2(12,6) = 0_wp + q2(11,6) = 0.0_wp + q2(12,6) = 0.0_wp q2(1,7) = 0.01734328280101122760645584447264624051821_wp q2(2,7) = -0.04840773491893482968357096782164797435370_wp q2(3,7) = 0.003397297753290025399699797610400916669720_wp @@ -657,7 +657,7 @@ subroutine up_deriv_gf_8_4_opt ( var, ni, nj, nk, dir, bb, gsize, delta, up, dva q2(9,7) = -0.3249230328418080529360937544303684126484_wp q2(10,7) = 0.07220511840929067843024305654008186947743_wp q2(11,7) = -0.007736262686709715546097470343580200301153_wp - q2(12,7) = 0_wp + q2(12,7) = 0.0_wp q2(1,8) = 0.008138876723484902470545040666116789759842_wp q2(2,8) = -0.02002921745130282774392822052620050919550_wp q2(3,8) = -0.002671343688148396036480262140129972226295_wp |