fig.dat

有限元求解柏松方程。本文采用FORTRAN语言编制程序。程序中大部分变量采用有名公共区存储方式存储

  VARIABLES= "X" "Y" "F" "DIFX" "DIFY" 
 ZONE I=         81J=          81
        .0000       .0000       .0000       .9997       .0000
        .0393       .0000       .0393       .9773      -.0438
        .0785       .0000       .0785       .9753      -.1094
        .1178       .0000       .1175       .9718      -.1749
        .1571       .0000       .1564       .9668      -.2401
        .1963       .0000       .1951       .9603      -.3049
        .2356       .0000       .2334       .9523      -.3692
        .2749       .0000       .2714       .9428      -.4330
        .3142       .0000       .3090       .9319      -.4961
        .3534       .0000       .3461       .9195      -.5584
        .3927       .0000       .3827       .9058      -.6198
        .4320       .0000       .4187       .8906      -.6803
        .4712       .0000       .4540       .8741      -.7398
        .5105       .0000       .4886       .8562      -.7981
        .5498       .0000       .5225       .8370      -.8552
        .5890       .0000       .5556       .8165      -.9110
        .6283       .0000       .5878       .7947      -.9653
        .6676       .0000       .6191       .7717     -1.0182
        .7069       .0000       .6494       .7475     -1.0695
        .7461       .0000       .6788       .7222     -1.1192
        .7854       .0000       .7071       .6958     -1.1671
        .8247       .0000       .7343       .6683     -1.2132
        .8639       .0000       .7604       .6397     -1.2575
        .9032       .0000       .7853       .6102     -1.2999
        .9425       .0000       .8090       .5797     -1.3402
        .9817       .0000       .8315       .5483     -1.3785
       1.0210       .0000       .8526       .5161     -1.4147
       1.0603       .0000       .8725       .4831     -1.4487
       1.0996       .0000       .8910       .4494     -1.4805
       1.1388       .0000       .9081       .4149     -1.5100
       1.1781       .0000       .9239       .3798     -1.5372
       1.2174       .0000       .9382       .3441     -1.5621
       1.2566       .0000       .9511       .3079     -1.5845
       1.2959       .0000       .9625       .2713     -1.6046
       1.3352       .0000       .9724       .2342     -1.6222
       1.3744       .0000       .9808       .1967     -1.6373
       1.4137       .0000       .9877       .1589     -1.6500
       1.4530       .0000       .9931       .1209     -1.6602
       1.4923       .0000       .9969       .0827     -1.6678
       1.5315       .0000       .9992       .0444     -1.6729
       1.5708       .0000      1.0000       .0060     -1.6756
       1.6101       .0000       .9992      -.0325     -1.6757
       1.6493       .0000       .9969      -.0708     -1.6733
       1.6886       .0000       .9931      -.1091     -1.6684
       1.7279       .0000       .9877      -.1472     -1.6610
       1.7671       .0000       .9808      -.1851     -1.6512
       1.8064       .0000       .9724      -.2227     -1.6390
       1.8457       .0000       .9625      -.2599     -1.6244
       1.8850       .0000       .9511      -.2968     -1.6074
       1.9242       .0000       .9382      -.3332     -1.5882
       1.9635       .0000       .9239      -.3691     -1.5667
       2.0028       .0000       .9081      -.4045     -1.5429
       2.0420       .0000       .8910      -.4392     -1.5171
       2.0813       .0000       .8725      -.4733     -1.4892
       2.1206       .0000       .8526      -.5066     -1.4593
       2.1598       .0000       .8315      -.5392     -1.4275
       2.1991       .0000       .8090      -.5709     -1.3938
       2.2384       .0000       .7853      -.6018     -1.3585
       2.2777       .0000       .7604      -.6318     -1.3215
       2.3169       .0000       .7343      -.6608     -1.2830
       2.3562       .0000       .7071      -.6888     -1.2431
       2.3955       .0000       .6788      -.7158     -1.2019
       2.4347       .0000       .6494      -.7417     -1.1596
       2.4740       .0000       .6191      -.7665     -1.1163
       2.5133       .0000       .5878      -.7902     -1.0723
       2.5525       .0000       .5556      -.8127     -1.0275
       2.5918       .0000       .5225      -.8339      -.9824
       2.6311       .0000       .4886      -.8539      -.9370
       2.6704       .0000       .4540      -.8727      -.8916
       2.7096       .0000       .4187      -.8902      -.8464
       2.7489       .0000       .3827      -.9063      -.8018
       2.7882       .0000       .3461      -.9212      -.7580
       2.8274       .0000       .3090      -.9347      -.7154
       2.8667       .0000       .2714      -.9469      -.6743
       2.9060       .0000       .2334      -.9578      -.6353
       2.9452       .0000       .1951      -.9674      -.5988
       2.9845       .0000       .1564      -.9757      -.5655
       3.0238       .0000       .1175      -.9827      -.5361
       3.0631       .0000       .0785      -.9886      -.5116
       3.1023       .0000       .0393      -.9934      -.4935
       3.1416       .0000       .0000      -.9963      -.4847
        .0000       .0393       .0000       .9559      -.0219
        .0393       .0393       .0367       .9340      -.0642
        .0785       .0393       .0733       .9319      -.1277
        .1178       .0393       .1098       .9283      -.1909
        .1571       .0393       .1462       .9232      -.2538
        .1963       .0393       .1823       .9167      -.3163
        .2356       .0393       .2181       .9089      -.3783
        .2749       .0393       .2536       .8996      -.4398
        .3142       .0393       .2887       .8889      -.5005
        .3534       .0393       .3234       .8769      -.5605
        .3927       .0393       .3575       .8635      -.6196
        .4320       .0393       .3912       .8488      -.6778
        .4712       .0393       .4242       .8327      -.7349
        .5105       .0393       .4565       .8154      -.7909
        .5498       .0393       .4882       .7968      -.8457
        .5890       .0393       .5191       .7770      -.8992
        .6283       .0393       .5492       .7560      -.9512
        .6676       .0393       .5784       .7339     -1.0019
        .7069       .0393       .6068       .7105     -1.0509
        .7461       .0393       .6342       .6862     -1.0984
        .7854       .0393       .6607       .6607     -1.1442
        .8247       .0393       .6861       .6342     -1.1882
        .8639       .0393       .7105       .6068     -1.2304
        .9032       .0393       .7337       .5784     -1.2706
        .9425       .0393       .7559       .5491     -1.3090
        .9817       .0393       .7768       .5190     -1.3453
       1.0210       .0393       .7966       .4880     -1.3796
       1.0603       .0393       .8152       .4563     -1.4117
       1.0996       .0393       .8325       .4239     -1.4417
       1.1388       .0393       .8485       .3909     -1.4695
       1.1781       .0393       .8632       .3572     -1.4950
       1.2174       .0393       .8765       .3230     -1.5183
       1.2566       .0393       .8885       .2883     -1.5392
       1.2959       .0393       .8992       .2532     -1.5578
       1.3352       .0393       .9084       .2176     -1.5740
       1.3744       .0393       .9163       .1817     -1.5878
       1.4137       .0393       .9227       .1455     -1.5993
       1.4530       .0393       .9278       .1091     -1.6083
       1.4923       .0393       .9313       .0725     -1.6148
       1.5315       .0393       .9335       .0358     -1.6190
       1.5708       .0393       .9342      -.0009     -1.6207
       1.6101       .0393       .9334      -.0377     -1.6200
       1.6493       .0393       .9313      -.0744     -1.6168
       1.6886       .0393       .9276      -.1110     -1.6113
       1.7279       .0393       .9226      -.1475     -1.6034
       1.7671       .0393       .9161      -.1837     -1.5931
       1.8064       .0393       .9082      -.2197     -1.5805
       1.8457       .0393       .8989      -.2554     -1.5656
       1.8850       .0393       .8882      -.2906     -1.5485
       1.9242       .0393       .8761      -.3255     -1.5292
       1.9635       .0393       .8626      -.3599     -1.5077
       2.0028       .0393       .8479      -.3937     -1.4841
       2.0420       .0393       .8318      -.4270     -1.4585
       2.0813       .0393       .8144      -.4596     -1.4309
       2.1206       .0393       .7957      -.4915     -1.4014
       2.1598       .0393       .7758      -.5228     -1.3701
       2.1991       .0393       .7547      -.5532     -1.3370
       2.2384       .0393       .7324      -.5829     -1.3024
       2.2777       .0393       .7090      -.6117     -1.2662
       2.3169       .0393       .6845      -.6396     -1.2286
       2.3562       .0393       .6588      -.6667     -1.1898
       2.3955       .0393       .6321      -.6927     -1.1497
       2.4347       .0393       .6045      -.7178     -1.1087
       2.4740       .0393       .5758      -.7418     -1.0668
       2.5133       .0393       .5463      -.7648     -1.0241
       2.5525       .0393       .5158      -.7868      -.9810
       2.5918       .0393       .4845      -.8077      -.9375
       2.6311       .0393       .4524      -.8275      -.8939
       2.6704       .0393       .4196      -.8461      -.8504
       2.7096       .0393       .3860      -.8637      -.8072
       2.7489       .0393       .3518      -.8802      -.7647
       2.7882       .0393       .3169      -.8955      -.7232
       2.8274       .0393       .2815      -.9098      -.6830
       2.8667       .0393       .2455      -.9231      -.6445
       2.9060       .0393       .2090      -.9354      -.6082
       2.9452       .0393       .1720      -.9467      -.5746
       2.9845       .0393       .1346      -.9573      -.5444
       3.0238       .0393       .0969      -.9672      -.5183
       3.0631       .0393       .0587      -.9767      -.4975
       3.1023       .0393       .0201      -.9865      -.4832
       3.1416       .0393      -.0189      -.9908      -.4756
        .0000       .0785       .0000       .8931      -.0205
        .0393       .0785       .0343       .8727      -.0600
        .0785       .0785       .0685       .8706      -.1192
        .1178       .0785       .1026       .8673      -.1783
        .1571       .0785       .1366       .8626      -.2371
        .1963       .0785       .1703       .8565      -.2955
        .2356       .0785       .2038       .8492      -.3534
        .2749       .0785       .2369       .8405      -.4108
        .3142       .0785       .2697       .8305      -.4676
        .3534       .0785       .3021       .8193      -.5236
        .3927       .0785       .3341       .8068      -.5788
        .4320       .0785       .3655       .7930      -.6332
        .4712       .0785       .3963       .7780      -.6865
        .5105       .0785       .4265       .7619      -.7388
        .5498       .0785       .4561       .7445      -.7900
        .5890       .0785       .4850       .7260      -.8400
        .6283       .0785       .5131       .7064      -.8886
        .6676       .0785       .5404       .6857      -.9359
        .7069       .0785       .5669       .6639      -.9817
        .7461       .0785       .5926       .6411     -1.0261
        .7854       .0785       .6173       .6173     -1.0688
        .8247       .0785       .6410       .5925     -1.1099
        .8639       .0785       .6638       .5669     -1.1493
        .9032       .0785       .6855       .5404     -1.1870
        .9425       .0785       .7062       .5130     -1.2228
        .9817       .0785       .7258       .4848     -1.2567
       1.0210       .0785       .7443       .4559     -1.2887
       1.0603       .0785       .7616       .4263     -1.3188
       1.0996       .0785       .7778       .3960     -1.3468
       1.1388       .0785       .7927       .3651     -1.3727
       1.1781       .0785       .8065       .3336     -1.3966
       1.2174       .0785       .8189       .3017     -1.4183
       1.2566       .0785       .8302       .2692     -1.4379
       1.2959       .0785       .8401       .2363     -1.4553
       1.3352       .0785       .8487       .2031     -1.4704
       1.3744       .0785       .8561       .1695     -1.4834
       1.4137       .0785       .8621       .1356     -1.4940
       1.4530       .0785       .8667       .1016     -1.5025
       1.4923       .0785       .8701       .0673     -1.5086
       1.5315       .0785       .8720       .0330     -1.5125
       1.5708       .0785       .8727      -.0014     -1.5141
       1.6101       .0785       .8720      -.0359     -1.5135
       1.6493       .0785       .8699      -.0703     -1.5106
       1.6886       .0785       .8665      -.1046     -1.5055
       1.7279       .0785       .8617      -.1387     -1.4981