tAdd standard with specified diffusivity and fix laplacian - cngf-pf - continuu… | |
git clone git://src.adamsgaard.dk/cngf-pf | |
Log | |
Files | |
Refs | |
README | |
LICENSE | |
--- | |
commit d708a8788c52a444e398488a77f8e5776d17370a | |
parent 3d1f6f47161a72d7408db704daa7509a6140d818 | |
Author: Anders Damsgaard <[email protected]> | |
Date: Thu, 9 Jul 2020 19:29:47 +0200 | |
Add standard with specified diffusivity and fix laplacian | |
Diffstat: | |
M fluid.c | 2 +- | |
M test/Makefile | 15 +++++++++++++++ | |
A test/simple_shear_wet_vari_diff.std | 200 +++++++++++++++++++++++++++++… | |
3 files changed, 216 insertions(+), 1 deletion(-) | |
--- | |
diff --git a/fluid.c b/fluid.c | |
t@@ -127,7 +127,7 @@ darcy_pressure_change_1d(const int i, | |
double k_, div_k_grad_p, k_zn, k_zp; | |
if (D > 0.0) | |
- return D*(p_f_ghost_in[i + 2] - 2.0 * p_f_ghost_in[i + 1] - p_… | |
+ return D*(p_f_ghost_in[i + 2] - 2.0 * p_f_ghost_in[i + 1] + p_… | |
else { | |
diff --git a/test/Makefile b/test/Makefile | |
t@@ -10,6 +10,7 @@ TESTS = simple_shear_dry \ | |
simple_shear_wet \ | |
simple_shear_wet_norm \ | |
simple_shear_wet_vari \ | |
+ simple_shear_wet_vari_diff \ | |
simple_shear_wet_vari_pulse \ | |
simple_shear_wet_vari_pulse_vlim \ | |
max_depth_misc\ | |
t@@ -33,6 +34,14 @@ SIMPLE_SHEAR_WET_VARI_OPTS = -L 8.0 \ | |
-q $$( echo "1.0/(3600*24)" | bc -l ) \ | |
-I $$( echo "60*1000" | bc -l ) \ | |
-e $$( echo "60*20" | bc -l ) | |
+SIMPLE_SHEAR_WET_VARI_DIFF_OPTS =-L 8.0 \ | |
+ -n 150e3 \ | |
+ -F -D $$(awk 'BEGIN{print 2e-17/(1.787e-3*(1e… | |
+ -O 50e3 \ | |
+ -a 50e3 \ | |
+ -q $$( echo "1.0/(3600*24)" | bc -l ) \ | |
+ -I $$( echo "60*1000" | bc -l ) \ | |
+ -e $$( echo "60*20" | bc -l ) | |
SIMPLE_SHEAR_WET_VARI_PULSE_OPTS = $(SIMPLE_SHEAR_WET_VARI_OPTS) \ | |
-u 500 | |
SIMPLE_SHEAR_WET_VARI_PULSE_VLIM_OPTS = $(SIMPLE_SHEAR_WET_VARI_OPTS) \ | |
t@@ -70,6 +79,9 @@ simple_shear_wet_norm: $(BIN) | |
simple_shear_wet_vari: $(BIN) | |
./$(BIN) $(SIMPLE_SHEAR_WET_VARI_OPTS) | diff [email protected] - | |
+simple_shear_wet_vari_diff: $(BIN) | |
+ ./$(BIN) $(SIMPLE_SHEAR_WET_VARI_DIFF_OPTS) | diff [email protected] - | |
+ | |
simple_shear_wet_vari_pulse: $(BIN) | |
./$(BIN) $(SIMPLE_SHEAR_WET_VARI_PULSE_OPTS) | diff [email protected] - | |
t@@ -110,6 +122,9 @@ simple_shear_wet_norm.std: $(BIN) | |
simple_shear_wet_vari.std: $(BIN) | |
./$(BIN) $(SIMPLE_SHEAR_WET_VARI_OPTS) > $@ | |
+simple_shear_wet_vari_diff.std: $(BIN) | |
+ ./$(BIN) $(SIMPLE_SHEAR_WET_VARI_DIFF_OPTS) > $@ | |
+ | |
simple_shear_wet_vari_pulse.std: $(BIN) | |
./$(BIN) $(SIMPLE_SHEAR_WET_VARI_PULSE_OPTS) > $@ | |
diff --git a/test/simple_shear_wet_vari_diff.std b/test/simple_shear_wet_vari_d… | |
t@@ -0,0 +1,200 @@ | |
+0 0 233416 69620 0.19278884052507114 3.5752… | |
+0.040201005025125629 4.8050055221904651e-72 232745.56783919601 … | |
+0.080402010050251257 1.9406513266205237e-71 232075.1356783919 … | |
+0.12060301507537688 6.3286427152837201e-71 231404.70351758791 … | |
+0.16080402010050251 1.9483522602725866e-70 230734.27135678395 … | |
+0.20100502512562815 5.8863439985873316e-70 230063.8391959799 … | |
+0.24120603015075376 1.7658964767396653e-69 229393.40703517586 … | |
+0.28140703517587939 5.2805426592428526e-69 228722.9748743719 … | |
+0.32160804020100503 1.5759001489371812e-68 228052.54271356785 … | |
+0.36180904522613067 4.6956158778910409e-68 227382.1105527638 … | |
+0.4020100502512563 1.3970978866884777e-67 226711.67839195981 … | |
+0.44221105527638194 4.1509426934649635e-67 226041.24623115576 … | |
+0.48241206030150752 1.2315580816858545e-66 225370.81407035171 … | |
+0.52261306532663321 3.6487864721758651e-66 224700.38190954775 … | |
+0.56281407035175879 1.0795029658895989e-65 224029.94974874376 … | |
+0.60301507537688448 3.1891584762408827e-65 223359.51758793966 … | |
+0.64321608040201006 9.4080768152498091e-65 222689.08542713569 … | |
+0.68341708542713564 2.7713591961298917e-64 222018.6532663317 … | |
+0.72361809045226133 8.1516849735523802e-64 221348.22110552766 … | |
+0.76381909547738691 2.3941930342519862e-63 220677.78894472361 … | |
+0.8040201005025126 7.0213945689750174e-63 220007.35678391962 … | |
+0.84422110552763818 2.0560540919963838e-62 219336.92462311557 … | |
+0.88442211055276387 6.0115673124570411e-62 218666.49246231152 … | |
+0.92462311557788945 1.7550007663734594e-61 217996.06030150756 … | |
+0.96482412060301503 5.1156196533612182e-61 217325.62814070351 … | |
+1.0050251256281406 1.4888289775083099e-60 216655.19597989946 … | |
+1.0452261306532664 4.3262515699989162e-60 215984.7638190955 … | |
+1.085427135678392 1.2551432437133814e-59 215314.33165829146 … | |
+1.1256281407035176 3.6356659959773821e-59 214643.89949748741 … | |
+1.1658291457286432 1.0514239931109531e-58 213973.46733668342 … | |
+1.206030150753769 3.0357739787100189e-58 213303.03517587937 … | |
+1.2462311557788945 8.7508968583765827e-58 212632.60301507538 … | |
+1.2864321608040201 2.5183815310557172e-57 211962.17085427136 … | |
+1.3266331658291457 7.2355262820756906e-57 211291.73869346737 … | |
+1.3668341708542713 2.0753551534149666e-56 210621.30653266332 … | |
+1.4070351758793971 5.9426768681379838e-56 209950.8743718593 … | |
+1.4472361809045227 1.6987640362005849e-55 209280.44221105531 … | |
+1.4874371859296482 4.8477343301224514e-55 208610.01005025127 … | |
+1.5276381909547738 1.3809979485813292e-54 207939.57788944722 … | |
+1.5678391959798996 3.927255526009741e-54 207269.14572864323 … | |
+1.6080402010050252 1.1148607396955186e-53 206598.71356783918 … | |
+1.6482412060301508 3.1592263073603044e-53 205928.28140703513 … | |
+1.6884422110552764 8.9364019223326389e-53 205257.84924623117 … | |
+1.7286432160804019 2.5232465973849635e-52 204587.41708542712 … | |
+1.7688442211055277 7.111556528471625e-52 203916.98492462307 … | |
+1.8090452261306533 2.0006481168722917e-51 203246.55276381911 … | |
+1.8492462311557789 5.6178540538992398e-51 202576.12060301506 … | |
+1.8894472361809045 1.5745516552389944e-50 201905.68844221107 … | |
+1.9296482412060301 4.4047614607915031e-50 201235.25628140703 … | |
+1.9698492462311559 1.2298717937835865e-49 200564.82412060298 … | |
+2.0100502512562812 3.4273716674460074e-49 199894.39195979899 … | |
+2.050251256281407 9.5327754824579315e-49 199223.95979899497 … | |
+2.0904522613065328 2.6462196246788041e-48 198553.52763819098 … | |
+2.1306532663316582 7.3311754607651057e-48 197883.09547738693 … | |
+2.170854271356784 2.026999677592068e-47 197212.66331658291 … | |
+2.2110552763819098 5.5931814336390175e-47 196542.23115577892 … | |
+2.2512562814070352 1.5402100946785845e-46 195871.79899497487 … | |
+2.291457286432161 4.232623802009668e-46 195201.36683417085 … | |
+2.3316582914572863 1.1607486586421787e-45 194530.93467336684 … | |
+2.3718592964824121 3.1765693982507022e-45 193860.50251256279 … | |
+2.4120603015075379 8.6748061704101453e-45 193190.0703517588 … | |
+2.4522613065326633 2.3639338542249397e-44 192519.63819095478 … | |
+2.4924623115577891 6.4279730177398356e-44 191849.20603015073 … | |
+2.5326633165829144 1.7440896539893465e-43 191178.77386934674 … | |
+2.5728643216080402 4.7218054563094601e-43 190508.34170854272 … | |
+2.613065326633166 1.2755126334172352e-42 189837.9095477387 … | |
+2.6532663316582914 3.4378482919359423e-42 189167.47738693468 … | |
+2.6934673366834172 9.2449974506898511e-42 188497.04522613063 … | |
+2.7336683417085426 2.4804603028666034e-41 187826.61306532664 … | |
+2.7738693467336684 6.639815812208223e-41 187156.1809045226 … | |
+2.8140703517587942 1.7732283595352075e-40 186485.74874371861 … | |
+2.8542713567839195 4.7244475119618528e-40 185815.31658291459 … | |
+2.8944723618090453 1.2557439478213632e-39 185144.88442211054 … | |
+2.9346733668341707 3.3297194569258055e-39 184474.45226130658 … | |
+2.9748743718592965 8.807568414345821e-39 183804.02010050253 … | |
+3.0150753768844223 2.3240149722312978e-38 183133.58793969848 … | |
+3.0552763819095476 6.1170362094943802e-38 182463.15577889449 … | |
+3.0954773869346734 1.6060348283441864e-37 181792.72361809044 … | |
+3.1356783919597992 4.2059562130649428e-37 181122.2914572864 … | |
+3.1758793969849246 1.098658633059263e-36 180451.85929648243 … | |
+3.2160804020100504 2.8624145826865942e-36 179781.42713567839 … | |
+3.2562814070351758 7.4381709155664987e-36 179110.99497487437 … | |
+3.2964824120603016 1.9277314810647299e-35 178440.56281407035 … | |
+3.3366834170854274 4.9827088859692855e-35 177770.13065326633 … | |
+3.3768844221105527 1.2844151587678675e-34 177099.69849246231 … | |
+3.4170854271356785 3.301852941851626e-34 176429.26633165829 … | |
+3.4572864321608039 8.4645567295111818e-34 175758.83417085427 … | |
+3.4974874371859297 2.1638938160800207e-33 175088.40201005025 … | |
+3.5376884422110555 5.5161257689233905e-33 174417.9698492462 … | |
+3.5778894472361809 1.4021336859866479e-32 173747.53768844221 … | |
+3.6180904522613067 3.5537139455124237e-32 173077.1055276382 … | |
+3.658291457286432 8.9805576758732282e-32 172406.67336683415 … | |
+3.6984924623115578 2.2627269788628886e-31 171736.24120603019 … | |
+3.7386934673366836 5.6840569639054206e-31 171065.80904522614 … | |
+3.778894472361809 1.4235142020287359e-30 170395.37688442209 … | |
+3.8190954773869348 3.5541114482580658e-30 169724.9447236181 … | |
+3.8592964824120601 8.8459670430488124e-30 169054.51256281405 … | |
+3.8994974874371859 2.1947883000505131e-29 168384.08040201006 … | |
+3.9396984924623117 5.4281477263964914e-29 167713.64824120604 … | |
+3.9798994974874371 1.3381647892044389e-28 167043.21608040199 … | |
+4.0201005025125625 3.2880966575681042e-28 166372.78391959801 … | |
+4.0603015075376883 8.052726590711344e-28 165702.35175879396 … | |
+4.1005025125628141 1.965541843892534e-27 165031.91959798994 … | |
+4.1407035175879399 4.7813226043655243e-27 164361.48743718595 … | |
+4.1809045226130657 1.1590900765552674e-26 163691.0552763819 … | |
+4.2211055276381906 2.8001042263698857e-26 163020.62311557791 … | |
+4.2613065326633164 6.7405514223191613e-26 162350.19095477386 … | |
+4.3015075376884422 1.6168298846849308e-25 161679.75879396981 … | |
+4.341708542713568 3.8641732300392552e-25 161009.32663316585 … | |
+4.3819095477386938 9.2014139394379743e-25 160338.8944723618 … | |
+4.4221105527638196 2.1828969845228847e-24 159668.46231155776 … | |
+4.4623115577889445 5.1590984650360351e-24 158998.0301507538 … | |
+4.5025125628140703 1.2146481492173996e-23 158327.59798994975 … | |
+4.5427135678391961 2.8486773682998949e-23 157657.16582914573 … | |
+4.5829145728643219 6.6546541854698924e-23 156986.73366834171 … | |
+4.6231155778894477 1.5483740552615012e-22 156316.30150753766 … | |
+4.6633165829145726 3.5881194363679973e-22 155645.86934673367 … | |
+4.7035175879396984 8.2808680549690454e-22 154975.43718592965 … | |
+4.7437185929648242 1.9031554200776046e-21 154305.00502512563 … | |
+4.78391959798995 4.3554997073297867e-21 153634.57286432161 … | |
+4.8241206030150758 9.9251578837658261e-21 152964.14070351759 … | |
+4.8643216080402008 2.2518889392775444e-20 152293.70854271355 … | |
+4.9045226130653266 5.0866961000440195e-20 151623.27638190956 … | |
+4.9447236180904524 1.1438700594873322e-19 150952.84422110554 … | |
+4.9849246231155782 2.5605752298833409e-19 150282.41206030149 … | |
+5.025125628140704 5.7054436227408555e-19 149611.97989949747 … | |
+5.0653266331658289 1.2653153368337802e-18 148941.54773869348 … | |
+5.1055276381909547 2.7927680076201811e-18 148271.11557788946 … | |
+5.1457286432160805 6.1342654600193087e-18 147600.68341708541 … | |
+5.1859296482412063 1.3407542942248825e-17 146930.25125628139 … | |
+5.2261306532663321 2.9158012634584991e-17 146259.8190954774 … | |
+5.266331658291457 6.3089044831858864e-17 145589.38693467336 … | |
+5.3065326633165828 1.3579962760430514e-16 144918.95477386937 … | |
+5.3467336683417086 2.9077346118275755e-16 144248.52261306532 … | |
+5.3869346733668344 6.1927240623565137e-16 143578.0904522613 … | |
+5.4271356783919602 1.3117165683896381e-15 142907.65829145728 … | |
+5.4673366834170851 2.7630359077066654e-15 142237.22613065329 … | |
+5.5075376884422109 5.7873427967340174e-15 141566.79396984924 … | |
+5.5477386934673367 1.2052370236908003e-14 140896.36180904522 … | |
+5.5879396984924625 2.4952918105850065e-14 140225.9296482412 … | |
+5.6281407035175883 5.1354379649382198e-14 139555.49748743715 … | |
+5.6683417085427132 1.050492167116606e-13 138885.06532663316 … | |
+5.708542713567839 2.1355803195132781e-13 138214.63316582915 … | |
+5.7487437185929648 4.3141500968057663e-13 137544.20100502513 … | |
+5.7889447236180906 8.6591647803847289e-13 136873.76884422108 … | |
+5.8291457286432165 1.7266443257985375e-12 136203.33668341709 … | |
+5.8693467336683414 3.4199303290203932e-12 135532.90452261307 … | |
+5.9095477386934672 6.7275923102569735e-12 134862.47236180905 … | |
+5.949748743718593 1.3142174332148021e-11 134192.04020100503 … | |
+5.9899497487437188 2.5490352312778978e-11 133521.60804020101 … | |
+6.0301507537688446 4.9081504298922784e-11 132851.17587939699 … | |
+6.0703517587939704 9.380457662331596e-11 132180.74371859297 … | |
+6.1105527638190953 1.7791858728506302e-10 131510.31155778895 … | |
+6.1507537688442211 3.3483793185881725e-10 130839.87939698491 … | |
+6.1909547738693469 6.2515102096402384e-10 130169.4472361809 … | |
+6.2311557788944727 1.1576853423478197e-09 129499.01507537687 … | |
+6.2713567839195985 2.1260107214288222e-09 128828.58291457285 … | |
+6.3115577889447234 3.8709766529148674e-09 128158.15075376883 … | |
+6.3517587939698492 6.9865372194821131e-09 127487.71859296481 … | |
+6.391959798994975 1.2496604824227431e-08 126817.28643216081 … | |
+6.4321608040201008 2.2146636693232311e-08 126146.85427135677 … | |
+6.4723618090452266 3.8877798826962013e-08 125476.42211055275 … | |
+6.5125628140703515 6.7586539700915413e-08 124805.98994974873 … | |
+6.5527638190954773 1.1632268567222197e-07 124135.55778894473 … | |
+6.5929648241206031 1.9814745364051073e-07 123465.12562814071 … | |
+6.6331658291457289 3.3396428735394409e-07 122794.69346733668 … | |
+6.6733668341708547 5.5674804671559208e-07 122124.26130653266 … | |
+6.7135678391959797 9.1773276199480786e-07 121453.82914572864 … | |
+6.7537688442211055 1.4952542135008859e-06 120783.39698492464 … | |
+6.7939698492462313 2.4070651303889052e-06 120112.9648241206 … | |
+6.8341708542713571 3.8269807241429926e-06 119442.53266331658 … | |
+6.8743718592964829 6.0066196931096725e-06 118772.10050251256 … | |
+6.9145728643216078 9.3026299198293091e-06 118101.66834170854 … | |
+6.9547738693467336 1.4209042807116444e-05 117431.23618090451 … | |
+6.9949748743718594 2.139304037645827e-05 116760.80402010051 … | |
+7.0351758793969852 3.1730411907661513e-05 116090.37185929649 … | |
+7.075376884422111 4.6334087439645592e-05 115419.93969849244 … | |
+7.1155778894472359 6.6565676044677173e-05 114749.50753768845 … | |
+7.1557788944723617 9.4016005621478993e-05 114079.07537688442 … | |
+7.1959798994974875 0.00013043745329413797 113408.64321608041 … | |
+7.2361809045226133 0.00017760924670937584 112738.21105527639 … | |
+7.2763819095477391 0.00023711940423803193 112067.77889447234 … | |
+7.316582914572864 0.00031005509841192872 111397.34673366827 … | |
+7.3567839195979898 0.00039660992690309716 110726.91457286339 … | |
+7.3969849246231156 0.00049785540937986243 110056.4824120531 … | |
+7.4371859296482414 0.00061542341663506329 109386.05025118985 … | |
+7.4773869346733672 0.00075118411031537306 108715.61808993857 … | |
+7.5175879396984921 0.00090699825648039899 108045.18592520704 … | |
+7.557788944723618 0.0010845419925986381 107374.75373618367 … | |
+7.5979899497487438 0.0012851893502205079 106704.32135361977 … | |
+7.6381909547738696 0.0015099293552649426 106033.88762269329 … | |
+7.6783919597989954 0.0017592885915218366 105363.44408836361 … | |
+7.7185929648241203 0.0020332227184863558 104692.93532080733 … | |
+7.7587939698492461 0.0023309275839085283 104021.98960253806 … | |
+7.7989949748743719 0.0026504949176678325 103348.5119453918 … | |
+7.8391959798994977 0.002988293481018256 102660.53002573134 … | |
+7.8793969849246235 0.0033378626127134746 101902.4171755026 … | |
+7.9195979899497493 0.0036878895412303693 100847.20870895487 … | |
+7.9597989949748742 0.0040172168885288431 98782.093561162779 … | |
+8 0.0042720236205609254 94708.202358414579 55291.79764158… |