[CIG-MC] Re: geoid scaling in CitcomS-2.2.1
Eh Tan
tan2 at geodynamics.org
Tue Oct 30 12:02:18 PDT 2007
The radial velocity looks strange to me. The minimal V_r is around -60,
but the maximal V_r is around 6000. Since this is a restarted run, can
you compare the min/max velocity to the previous run? Also, using a more
strict convergence criterion might help.
Eh
Scott King wrote:
>
> These looked pretty good. Even the stresses looked pretty good when
> I plotted them but looking at the minmax I do see some strange
> min/maxes.. No, this was a restart from a previous run, just to get
> the geoid and topography.
>
>
> Stresses:
> (as an example)
> mercl8out.stress.95.0: N = 35937 <-95171/105220>
> <-115840/214550> <-1.7369e+06/397900> <-69303/61590>
> <-358860/278820> <-347420/207560>
> mercl8out.stress.95.1: N = 35937 <-58781/84249>
> <-175000/303600> <-950880/531570> <-52469/48321>
> <-60842/117690> <-345240/173120>
>
>
>
> Velocities
>
> conman:~/Citcom_Mercury sdk$ minmax -H2 mercl8out.velo*
> mercl8out.velo.0.0: N = 35937 <-495.617/460.354>
> <-720.09/699.195> <-40.0003/4330.05> <0.655123/1>
> mercl8out.velo.0.1: N = 35937 <-515.51/415.651>
> <-726.372/692.064> <-45.2284/4258.82> <0.655118/1>
> mercl8out.velo.1.0: N = 35937 <-48.499/131.136>
> <-25.1344/170.102> <-2.63491/3779.79> <0/0.955784>
> mercl8out.velo.1.1: N = 35937 <-44.5856/135.854>
> <-21.8894/155.796> <-3.51104/3808.08> <0/0.955782>
> mercl8out.velo.10.0: N = 35937 <-575.47/559.913>
> <-267.812/304.38> <-42.6444/3384.6> <0.655572/1>
> mercl8out.velo.10.1: N = 35937 <-590.053/556.736>
> <-261.531/311.155> <-50.094/3360.02> <0.655567/1>
> mercl8out.velo.11.0: N = 35937 <-58.3954/72.6513>
> <-34.554/33.267> <-2.9017/2925.71> <0/0.95849>
> mercl8out.velo.11.1: N = 35937 <-52.6087/63.2265>
> <-29.9125/30.8822> <-3.92529/2936.07> <0/0.958488>
> mercl8out.velo.12.0: N = 35937 <-559.914/575.47>
> <-304.381/267.812> <-42.6444/3384.6> <0.655572/1>
> mercl8out.velo.12.1: N = 35937 <-556.736/590.053>
> <-311.155/261.531> <-50.094/3360.02> <0.655567/1>
> mercl8out.velo.13.0: N = 35937 <-72.6512/58.3955>
> <-33.267/34.554> <-2.9017/2925.71> <0/0.95849>
> mercl8out.velo.13.1: N = 35937 <-63.2264/52.6088>
> <-30.8822/29.9126> <-3.92529/2936.07> <0/0.958488>
> mercl8out.velo.14.0: N = 35937 <-579.348/86.1794>
> <-58.9119/297.355> <-42.2866/3384.6> <0.65627/1>
> mercl8out.velo.14.1: N = 35937 <-586.346/135.134>
> <-93.1372/301.098> <-48.4822/3360.02> <0.656265/1>
> mercl8out.velo.15.0: N = 35937 <-0.511314/70.9279>
> <-35.4576/0.291519> <-2.89219/2925.71> <0/0.95849>
> mercl8out.velo.15.1: N = 35937 <-4.00268/62.4625>
> <-31.929/2.06707> <-3.91594/2936.07> <0/0.958488>
> mercl8out.velo.16.0: N = 35937 <-652.384/685.167>
> <-488.217/349.061> <-44.0404/4101.21> <0.655014/1>
> mercl8out.velo.16.1: N = 35937 <-681.07/693.953>
> <-501.859/305.14> <-52.7059/4027.35> <0.655009/1>
> mercl8out.velo.17.0: N = 35937 <-68.5838/110.214>
> <-102.603/27.7676> <-3.17478/3537.83> <0/0.957583>
> mercl8out.velo.17.1: N = 35937 <-61.4193/96.2775>
> <-97.4757/25.9245> <-4.34169/3553.64> <0/0.95758>
> mercl8out.velo.18.0: N = 35937 <-414.956/495.617>
> <-507.847/629.346> <-41.5519/3944.99> <0.654825/1>
> mercl8out.velo.18.1: N = 35937 <-449.233/520.124>
> <-526.764/683.947> <-48.1316/3927.47> <0.65482/1>
> mercl8out.velo.19.0: N = 35937 <-38.5187/100.432>
> <-142.323/38.105> <-2.97485/3543.2> <0/0.959991>
> mercl8out.velo.19.1: N = 35937 <-32.6027/98.6713>
> <-132.425/30.8198> <-4.04284/3595.27> <0/0.959989>
> mercl8out.velo.2.0: N = 35937 <-717.261/729.466>
> <-434.425/457.269> <-43.7896/4658.57> <0.655101/1>
> mercl8out.velo.2.1: N = 35937 <-806.528/788.015>
> <-488.611/462.933> <-49.4163/4660.84> <0.655096/1>
> mercl8out.velo.20.0: N = 35937 <-729.466/717.261>
> <-457.269/434.425> <-43.7897/4658.57> <0.655101/1>
> mercl8out.velo.20.1: N = 35937 <-788.015/806.527>
> <-462.933/488.612> <-49.4163/4660.83> <0.655096/1>
> mercl8out.velo.21.0: N = 35937 <-151.566/70.8971>
> <-66.3183/93.5275> <-2.85027/4290.69> <0/0.957167>
> mercl8out.velo.21.1: N = 35937 <-143.49/54.7005>
> <-65.121/91.8167> <-3.95522/4374.48> <0/0.957165>
> mercl8out.velo.22.0: N = 35937 <-460.354/495.617>
> <-699.194/720.09> <-40.0004/4330.05> <0.655123/1>
> mercl8out.velo.22.1: N = 35937 <-415.651/515.51>
> <-692.064/726.372> <-45.2284/4258.82> <0.655118/1>
> mercl8out.velo.23.0: N = 35937 <-131.136/48.499>
> <-170.102/25.1344> <-2.63491/3779.79> <0/0.955784>
> mercl8out.velo.23.1: N = 35937 <-135.854/44.5856>
> <-155.796/21.8894> <-3.51104/3808.08> <0/0.955782>
> mercl8out.velo.24.0: N = 35937 <-113.515/369.743>
> <-262.708/684.369> <-44.7282/3637.99> <0.652803/1>
> mercl8out.velo.24.1: N = 35937 <-63.0834/416.582>
> <-400.696/792.529> <-56.078/3763.32> <0.652798/1>
> mercl8out.velo.25.0: N = 35937 <-1.59055/134.731>
> <-95.2661/2.46482> <-3.67851/3442.95> <0/0.957604>
> mercl8out.velo.25.1: N = 35937 <-3.78883/135.098>
> <-93.2064/9.96051> <-4.70937/3592.59> <0/0.957603>
> mercl8out.velo.26.0: N = 35937 <-596.206/599.106>
> <-483.107/339.805> <-45.4608/3199.31> <0.653011/1>
> mercl8out.velo.26.1: N = 35937 <-588.448/629.652>
> <-544.088/334.229> <-55.4534/3175.71> <0.653006/1>
> mercl8out.velo.27.0: N = 35937 <-101.722/48.3092>
> <-34.504/88.5194> <-3.67851/2744.67> <0/0.959414>
> mercl8out.velo.27.1: N = 35937 <-98.5019/40.0017>
> <-31.3972/81.9855> <-4.67605/2772.2> <0/0.959413>
> mercl8out.velo.28.0: N = 35937 <-774.532/405.788>
> <-473.06/779.929> <-43.4455/6316.01> <0.657004/1>
> mercl8out.velo.28.1: N = 35937 <-907.767/444.488>
> <-471.411/925.828> <-49.1657/6347.24> <0.656999/1>
> mercl8out.velo.29.0: N = 35937 <-170.508/77.1825>
> <-121.022/57.1721> <-3.03351/6083.03> <0/0.960348>
> mercl8out.velo.29.1: N = 35937 <-173.179/80.6653>
> <-113.717/75.6066> <-3.93105/6204.43> <0/0.960363>
> mercl8out.velo.3.0: N = 35937 <-70.897/151.566>
> <-93.5276/66.3183> <-2.85027/4290.69> <0/0.957167>
> mercl8out.velo.3.1: N = 35937 <-54.7005/143.49>
> <-91.8168/65.1209> <-3.95522/4374.48> <0/0.957165>
> mercl8out.velo.30.0: N = 35937 <-659.25/692.768>
> <-503.557/335.644> <-44.5685/4023.62> <0.655755/1>
> mercl8out.velo.30.1: N = 35937 <-682.874/697.67>
> <-505.873/297.346> <-53.2522/3932.45> <0.65575/1>
> mercl8out.velo.31.0: N = 35937 <-86.268/98.094>
> <-100.427/38.8917> <-3.12328/3480.16> <0/0.958175>
> mercl8out.velo.31.1: N = 35937 <-79.0317/84.4752>
> <-98.8705/35.9701> <-4.12578/3498.93> <0/0.958174>
> mercl8out.velo.32.0: N = 35937 <-615.561/666.04>
> <-363.747/363.692> <-45.004/3531.55> <0.656983/1>
> mercl8out.velo.32.1: N = 35937 <-640.65/688.511>
> <-335.362/335.343> <-55.2954/3517.12> <0.656978/1>
> mercl8out.velo.33.0: N = 35937 <-117.444/63.8421>
> <-79.1824/79.1818> <-3.27651/3137> <0/0.956828>
> mercl8out.velo.33.1: N = 35937 <-107.894/52.7653>
> <-71.6427/71.642> <-4.31339/3175.54> <0/0.956827>
> mercl8out.velo.34.0: N = 35937 <-700.391/645.689>
> <-326.214/493.348> <-42.7479/3981.81> <0.657462/1>
> mercl8out.velo.34.1: N = 35937 <-688.686/657.913>
> <-302.474/484.298> <-46.748/3881.93> <0.657457/1>
> mercl8out.velo.35.0: N = 35937 <-117.444/65.9559>
> <-22.717/129.923> <-2.88399/3367.66> <0/0.95638>
> mercl8out.velo.35.1: N = 35937 <-107.894/54.6108>
> <-22.2502/124.915> <-3.72907/3374.03> <0/0.956378>
> mercl8out.velo.36.0: N = 35937 <-700.393/645.689>
> <-493.352/326.159> <-42.7459/3981.82> <0.657462/1>
> mercl8out.velo.36.1: N = 35937 <-688.688/657.914>
> <-484.301/302.453> <-46.7485/3881.94> <0.657457/1>
> mercl8out.velo.37.0: N = 35937 <-117.444/65.956>
> <-129.924/22.7155> <-2.88399/3367.67> <0/0.95638>
> mercl8out.velo.37.1: N = 35937 <-107.894/54.6109>
> <-124.916/22.2465> <-3.72907/3374.04> <0/0.956378>
> mercl8out.velo.38.0: N = 35937 <-522.218/377.472>
> <-118.007/118.039> <-44.1136/290.844> <0.65611/1>
> mercl8out.velo.38.1: N = 35937 <-569.196/578.38>
> <-158.131/158.104> <-57.8015/346.913> <0.656105/1>
> mercl8out.velo.39.0: N = 35937 <-8.66977/16.7786>
> <-0.946845/0.945456> <-3.07326/95.1445> <0/0.803075>
> mercl8out.velo.39.1: N = 35937 <-5.32647/7.84011>
> <-1.94153/1.94149> <-4.27684/109.725> <0/0.803076>
> mercl8out.velo.4.0: N = 35937 <-495.617/414.957>
> <-629.346/507.847> <-41.5519/3944.99> <0.654825/1>
> mercl8out.velo.4.1: N = 35937 <-520.124/449.234>
> <-683.948/526.764> <-48.1316/3927.47> <0.65482/1>
> mercl8out.velo.40.0: N = 35937 <-652.346/685.154>
> <-348.924/488.081> <-44.0382/4101.07> <0.655011/1>
> mercl8out.velo.40.1: N = 35937 <-680.915/693.901>
> <-305.145/501.705> <-52.7108/4027.27> <0.655006/1>
> mercl8out.velo.41.0: N = 35937 <-68.5946/110.205>
> <-27.7567/102.542> <-3.17441/3538.2> <0/0.957584>
> mercl8out.velo.41.1: N = 35937 <-61.4101/96.2019>
> <-25.9175/97.3781> <-4.34179/3552.81> <0/0.957582>
> mercl8out.velo.42.0: N = 35937 <-729.469/717.263>
> <-434.437/457.282> <-43.7748/4658.52> <0.655101/1>
> mercl8out.velo.42.1: N = 35937 <-787.998/806.527>
> <-488.526/462.945> <-49.4149/4660.82> <0.655096/1>
> mercl8out.velo.43.0: N = 35937 <-151.569/70.8971>
> <-93.5972/66.5373> <-2.85027/4290.69> <0/0.957167>
> mercl8out.velo.43.1: N = 35937 <-143.49/54.7005>
> <-91.8289/65.3082> <-3.95472/4374.45> <0/0.957165>
> mercl8out.velo.44.0: N = 35937 <-414.925/495.62>
> <-629.058/507.625> <-41.5494/3945.97> <0.654826/1>
> mercl8out.velo.44.1: N = 35937 <-449.234/520.177>
> <-683.94/526.4> <-48.1243/3928.82> <0.654821/1>
> mercl8out.velo.45.0: N = 35937 <-38.5317/100.422>
> <-38.1139/142.362> <-2.97415/3545.51> <0/0.959991>
> mercl8out.velo.45.1: N = 35937 <-32.6756/98.6669>
> <-30.8243/132.439> <-4.04151/3597.84> <0/0.959989>
> mercl8out.velo.46.0: N = 35937 <-460.209/495.62>
> <-719.904/699.473> <-40.0031/4330.94> <0.655123/1>
> mercl8out.velo.46.1: N = 35937 <-415.631/515.587>
> <-726.514/692.145> <-45.2375/4260.09> <0.655118/1>
> mercl8out.velo.47.0: N = 35937 <-131.314/48.5557>
> <-25.1227/170.188> <-2.63583/3782.8> <0/0.955778>
> mercl8out.velo.47.1: N = 35937 <-136.049/44.627>
> <-21.9476/155.858> <-3.51206/3811.34> <0/0.955776>
> mercl8out.velo.48.0: N = 35937 <-113.566/369.743>
> <-684.396/262.728> <-44.7259/3637.99> <0.652806/1>
> mercl8out.velo.48.1: N = 35937 <-63.094/416.582>
> <-792.441/400.73> <-56.0853/3763.32> <0.6528/1>
> mercl8out.velo.49.0: N = 35937 <-1.6362/134.731>
> <-2.46058/95.3973> <-3.67888/3442.95> <0/0.957604>
> mercl8out.velo.49.1: N = 35937 <-3.82845/135.098>
> <-9.96575/93.3981> <-4.70919/3592.59> <0/0.957603>
> mercl8out.velo.5.0: N = 35937 <-100.432/38.5186>
> <-38.105/142.323> <-2.97485/3543.2> <0/0.959991>
> mercl8out.velo.5.1: N = 35937 <-98.6712/32.6027>
> <-30.8198/132.425> <-4.04284/3595.27> <0/0.959989>
> mercl8out.velo.50.0: N = 35937 <-774.532/405.788>
> <-779.681/472.419> <-43.4425/6316.01> <0.657004/1>
> mercl8out.velo.50.1: N = 35937 <-907.767/444.488>
> <-925.746/470.738> <-49.1243/6347.24> <0.656999/1>
> mercl8out.velo.51.0: N = 35937 <-170.508/77.1825>
> <-57.1113/121.136> <-3.03189/6083.03> <0/0.960348>
> mercl8out.velo.51.1: N = 35937 <-173.232/80.6653>
> <-75.5528/113.848> <-3.92822/6204.43> <0/0.960363>
> mercl8out.velo.52.0: N = 35937 <-596.433/599.579>
> <-339.783/483.087> <-45.4592/3200.65> <0.653008/1>
> mercl8out.velo.52.1: N = 35937 <-588.601/629.505>
> <-334.142/544.087> <-55.4683/3177.23> <0.653003/1>
> mercl8out.velo.53.0: N = 35937 <-101.672/48.4117>
> <-88.462/34.5135> <-3.67888/2746.69> <0/0.959414>
> mercl8out.velo.53.1: N = 35937 <-98.486/39.9099>
> <-81.9052/31.4169> <-4.67776/2774.35> <0/0.959413>
> mercl8out.velo.54.0: N = 35937 <-659.213/692.843>
> <-335.644/503.293> <-44.5617/4022.31> <0.655755/1>
> mercl8out.velo.54.1: N = 35937 <-682.987/697.677>
> <-297.352/505.589> <-53.2227/3931.01> <0.65575/1>
> mercl8out.velo.55.0: N = 35937 <-86.3097/98.0252>
> <-38.7978/100.383> <-3.12167/3479.55> <0/0.958174>
> mercl8out.velo.55.1: N = 35937 <-78.9729/84.4836>
> <-35.901/98.8526> <-4.12319/3498.01> <0/0.958172>
> mercl8out.velo.56.0: N = 35937 <-86.2151/579.349>
> <-58.9356/297.356> <-42.2867/3384.6> <0.65627/1>
> mercl8out.velo.56.1: N = 35937 <-135.125/586.346>
> <-93.1251/301.099> <-48.4823/3360.02> <0.656265/1>
> mercl8out.velo.57.0: N = 35937 <-70.9279/0.511314>
> <-35.4576/0.29152> <-2.89218/2925.71> <0/0.95849>
> mercl8out.velo.57.1: N = 35937 <-62.4624/4.00268>
> <-31.9289/2.06707> <-3.91594/2936.07> <0/0.958488>
> mercl8out.velo.58.0: N = 35937 <-559.89/575.47>
> <-267.815/304.381> <-42.6443/3384.6> <0.655572/1>
> mercl8out.velo.58.1: N = 35937 <-556.735/590.046>
> <-261.535/311.155> <-50.0941/3360.02> <0.655567/1>
> mercl8out.velo.59.0: N = 35937 <-72.6512/58.3955>
> <-34.5541/33.2669> <-2.9017/2925.71> <0/0.95849>
> mercl8out.velo.59.1: N = 35937 <-63.2264/52.6088>
> <-29.9125/30.882> <-3.92528/2936.07> <0/0.958488>
> mercl8out.velo.6.0: N = 35937 <-685.167/652.384>
> <-349.061/488.217> <-44.0404/4101.21> <0.655014/1>
> mercl8out.velo.6.1: N = 35937 <-693.953/681.071>
> <-305.141/501.859> <-52.7059/4027.35> <0.655009/1>
> mercl8out.velo.60.0: N = 35937 <-575.471/559.89>
> <-304.381/267.815> <-42.6443/3384.6> <0.655572/1>
> mercl8out.velo.60.1: N = 35937 <-590.045/556.735>
> <-311.155/261.535> <-50.0941/3360.02> <0.655567/1>
> mercl8out.velo.61.0: N = 35937 <-58.3956/72.6512>
> <-33.2669/34.554> <-2.9017/2925.71> <0/0.95849>
> mercl8out.velo.61.1: N = 35937 <-52.6089/63.2264>
> <-30.8821/29.9126> <-3.92528/2936.07> <0/0.958488>
> mercl8out.velo.62.0: N = 35937 <-579.348/86.215>
> <-297.355/58.9355> <-42.2866/3384.6> <0.65627/1>
> mercl8out.velo.62.1: N = 35937 <-586.346/135.125>
> <-301.099/93.1249> <-48.4822/3360.02> <0.656265/1>
> mercl8out.velo.63.0: N = 35937 <-0.511314/70.9277>
> <-0.291519/35.4576> <-2.89218/2925.71> <0/0.95849>
> mercl8out.velo.63.1: N = 35937 <-4.00268/62.4623>
> <-2.06707/31.9289> <-3.91594/2936.07> <0/0.958488>
> mercl8out.velo.64.0: N = 35937 <-692.843/659.213>
> <-503.293/335.644> <-44.5617/4022.31> <0.655755/1>
> mercl8out.velo.64.1: N = 35937 <-697.677/682.987>
> <-505.589/297.352> <-53.2227/3931.01> <0.65575/1>
> mercl8out.velo.65.0: N = 35937 <-98.0252/86.3096>
> <-100.383/38.7978> <-3.12168/3479.55> <0/0.958174>
> mercl8out.velo.65.1: N = 35937 <-84.4836/78.9728>
> <-98.8526/35.901> <-4.12319/3498.01> <0/0.958172>
> mercl8out.velo.66.0: N = 35937 <-599.579/596.433>
> <-483.087/339.783> <-45.4592/3200.65> <0.653008/1>
> mercl8out.velo.66.1: N = 35937 <-629.505/588.601>
> <-544.088/334.142> <-55.4683/3177.22> <0.653003/1>
> mercl8out.velo.67.0: N = 35937 <-48.4117/101.672>
> <-34.5134/88.4621> <-3.67888/2746.69> <0/0.959414>
> mercl8out.velo.67.1: N = 35937 <-39.91/98.4859>
> <-31.4169/81.9052> <-4.67776/2774.35> <0/0.959413>
> mercl8out.velo.68.0: N = 35937 <-405.788/774.532>
> <-472.419/779.681> <-43.4425/6316.01> <0.657004/1>
> mercl8out.velo.68.1: N = 35937 <-444.489/907.767>
> <-470.738/925.746> <-49.1243/6347.25> <0.656999/1>
> mercl8out.velo.69.0: N = 35937 <-77.1826/170.508>
> <-121.136/57.1111> <-3.03189/6083.03> <0/0.960348>
> mercl8out.velo.69.1: N = 35937 <-80.6653/173.232>
> <-113.847/75.5526> <-3.92822/6204.43> <0/0.960363>
> mercl8out.velo.7.0: N = 35937 <-110.214/68.5838>
> <-27.7676/102.603> <-3.17478/3537.83> <0/0.957583>
> mercl8out.velo.7.1: N = 35937 <-96.2776/61.4193>
> <-25.9245/97.4758> <-4.34169/3553.64> <0/0.95758>
> mercl8out.velo.70.0: N = 35937 <-369.743/113.566>
> <-262.729/684.396> <-44.7259/3637.99> <0.652806/1>
> mercl8out.velo.70.1: N = 35937 <-416.582/63.0939>
> <-400.73/792.441> <-56.0854/3763.32> <0.6528/1>
> mercl8out.velo.71.0: N = 35937 <-134.731/1.63621>
> <-95.3973/2.46059> <-3.67888/3442.95> <0/0.957604>
> mercl8out.velo.71.1: N = 35937 <-135.098/3.82851>
> <-93.3982/9.96576> <-4.70919/3592.59> <0/0.957603>
> mercl8out.velo.72.0: N = 35937 <-495.62/460.209>
> <-699.473/719.904> <-40.0031/4330.94> <0.655123/1>
> mercl8out.velo.72.1: N = 35937 <-515.587/415.632>
> <-692.145/726.514> <-45.2375/4260.09> <0.655118/1>
> mercl8out.velo.73.0: N = 35937 <-48.5557/131.314>
> <-170.188/25.1227> <-2.63583/3782.8> <0/0.955778>
> mercl8out.velo.73.1: N = 35937 <-44.627/136.049>
> <-155.858/21.9476> <-3.51206/3811.34> <0/0.955776>
> mercl8out.velo.74.0: N = 35937 <-495.62/414.924>
> <-507.625/629.058> <-41.5494/3945.96> <0.654826/1>
> mercl8out.velo.74.1: N = 35937 <-520.177/449.234>
> <-526.4/683.94> <-48.1243/3928.82> <0.654821/1>
> mercl8out.velo.75.0: N = 35937 <-100.422/38.5317>
> <-142.362/38.1139> <-2.97416/3545.51> <0/0.959991>
> mercl8out.velo.75.1: N = 35937 <-98.6671/32.6757>
> <-132.439/30.8243> <-4.04151/3597.84> <0/0.959989>
> mercl8out.velo.76.0: N = 35937 <-717.263/729.469>
> <-457.282/434.438> <-43.7748/4658.52> <0.655101/1>
> mercl8out.velo.76.1: N = 35937 <-806.528/787.999>
> <-462.945/488.526> <-49.4149/4660.82> <0.655096/1>
> mercl8out.velo.77.0: N = 35937 <-70.897/151.569>
> <-66.5373/93.5973> <-2.85027/4290.69> <0/0.957167>
> mercl8out.velo.77.1: N = 35937 <-54.7005/143.49>
> <-65.3082/91.829> <-3.95472/4374.45> <0/0.957165>
> mercl8out.velo.78.0: N = 35937 <-685.153/652.346>
> <-488.081/348.924> <-44.0383/4101.07> <0.655011/1>
> mercl8out.velo.78.1: N = 35937 <-693.9/680.915>
> <-501.705/305.145> <-52.7108/4027.27> <0.655006/1>
> mercl8out.velo.79.0: N = 35937 <-110.205/68.5945>
> <-102.542/27.7567> <-3.17441/3538.2> <0/0.957584>
> mercl8out.velo.79.1: N = 35937 <-96.202/61.4101>
> <-97.3782/25.9174> <-4.34179/3552.81> <0/0.957582>
> mercl8out.velo.8.0: N = 35937 <-86.1792/579.348>
> <-297.355/58.9117> <-42.2866/3384.6> <0.65627/1>
> mercl8out.velo.8.1: N = 35937 <-135.133/586.346>
> <-301.099/93.1371> <-48.4822/3360.02> <0.656265/1>
> mercl8out.velo.80.0: N = 35937 <-377.473/522.218>
> <-118.039/118.007> <-44.1136/290.844> <0.65611/1>
> mercl8out.velo.80.1: N = 35937 <-578.381/569.196>
> <-158.104/158.13> <-57.8015/346.913> <0.656105/1>
> mercl8out.velo.81.0: N = 35937 <-16.7786/8.66981>
> <-0.945463/0.946846> <-3.07326/95.1445> <0/0.803075>
> mercl8out.velo.81.1: N = 35937 <-7.84011/5.32647>
> <-1.94149/1.94153> <-4.27684/109.725> <0/0.803076>
> mercl8out.velo.82.0: N = 35937 <-645.69/700.393>
> <-326.159/493.352> <-42.7459/3981.82> <0.657462/1>
> mercl8out.velo.82.1: N = 35937 <-657.914/688.687>
> <-302.453/484.301> <-46.7485/3881.94> <0.657457/1>
> mercl8out.velo.83.0: N = 35937 <-65.9561/117.444>
> <-22.7155/129.924> <-2.88399/3367.67> <0/0.95638>
> mercl8out.velo.83.1: N = 35937 <-54.611/107.894>
> <-22.2465/124.916> <-3.72907/3374.04> <0/0.956378>
> mercl8out.velo.84.0: N = 35937 <-645.689/700.391>
> <-493.348/326.214> <-42.7479/3981.81> <0.657462/1>
> mercl8out.velo.84.1: N = 35937 <-657.913/688.686>
> <-484.298/302.475> <-46.748/3881.93> <0.657457/1>
> mercl8out.velo.85.0: N = 35937 <-65.956/117.444>
> <-129.923/22.7171> <-2.88399/3367.66> <0/0.95638>
> mercl8out.velo.85.1: N = 35937 <-54.6108/107.894>
> <-124.915/22.2502> <-3.72907/3374.03> <0/0.956378>
> mercl8out.velo.86.0: N = 35937 <-666.04/615.561>
> <-363.692/363.747> <-45.0039/3531.55> <0.656983/1>
> mercl8out.velo.86.1: N = 35937 <-688.511/640.65>
> <-335.343/335.362> <-55.2954/3517.12> <0.656978/1>
> mercl8out.velo.87.0: N = 35937 <-63.8422/117.444>
> <-79.1817/79.1825> <-3.27651/3137> <0/0.956828>
> mercl8out.velo.87.1: N = 35937 <-52.7653/107.894>
> <-71.642/71.6428> <-4.31339/3175.54> <0/0.956827>
> mercl8out.velo.88.0: N = 35937 <-692.768/659.25>
> <-335.644/503.557> <-44.5684/4023.62> <0.655755/1>
> mercl8out.velo.88.1: N = 35937 <-697.67/682.874>
> <-297.346/505.873> <-53.2522/3932.46> <0.65575/1>
> mercl8out.velo.89.0: N = 35937 <-98.094/86.268>
> <-38.8918/100.427> <-3.12328/3480.16> <0/0.958175>
> mercl8out.velo.89.1: N = 35937 <-84.4752/79.0318>
> <-35.9701/98.8706> <-4.12578/3498.93> <0/0.958174>
> mercl8out.velo.9.0: N = 35937 <-70.9278/0.511314>
> <-0.291519/35.4576> <-2.89219/2925.71> <0/0.95849>
> mercl8out.velo.9.1: N = 35937 <-62.4625/4.00268>
> <-2.06707/31.929> <-3.91594/2936.07> <0/0.958488>
> mercl8out.velo.90.0: N = 35937 <-405.788/774.532>
> <-779.929/473.06> <-43.4455/6316.01> <0.657004/1>
> mercl8out.velo.90.1: N = 35937 <-444.489/907.767>
> <-925.827/471.411> <-49.1657/6347.25> <0.656999/1>
> mercl8out.velo.91.0: N = 35937 <-77.1826/170.508>
> <-57.172/121.022> <-3.03351/6083.03> <0/0.960348>
> mercl8out.velo.91.1: N = 35937 <-80.6653/173.179>
> <-75.6066/113.717> <-3.93105/6204.43> <0/0.960363>
> mercl8out.velo.92.0: N = 35937 <-599.106/596.206>
> <-339.806/483.107> <-45.4608/3199.31> <0.653011/1>
> mercl8out.velo.92.1: N = 35937 <-629.651/588.448>
> <-334.229/544.088> <-55.4534/3175.72> <0.653006/1>
> mercl8out.velo.93.0: N = 35937 <-48.3092/101.722>
> <-88.5194/34.504> <-3.67851/2744.67> <0/0.959414>
> mercl8out.velo.93.1: N = 35937 <-40.0018/98.502>
> <-81.9856/31.3971> <-4.67605/2772.21> <0/0.959413>
> mercl8out.velo.94.0: N = 35937 <-369.743/113.515>
> <-684.369/262.708> <-44.7282/3637.99> <0.652803/1>
> mercl8out.velo.94.1: N = 35937 <-416.582/63.0835>
> <-792.529/400.696> <-56.078/3763.32> <0.652798/1>
> mercl8out.velo.95.0: N = 35937 <-134.731/1.59054>
> <-2.46481/95.266> <-3.67851/3442.95> <0/0.957604>
> mercl8out.velo.95.1: N = 35937 <-135.098/3.78884>
> <-9.96052/93.2065> <-4.70937/3592.59> <0/0.957603>
>
>
>
>
> On Oct 30, 2007, at 2:29 PM, Eh Tan wrote:
>
>> Scott King wrote:
>>> initial residue of momentum equation F 1.684117294e+00 9741456
>>> AhatP (000) after 1.88612 seconds with div/v=2.084e-05 dv/v=0.000e+00
>>> and dp/p=0.000e+00 for step 0
>>> bottom heat flux= 26.146639
>>> surface heat flux= 7.580164
>>> initial residue of momentum equation F 1.684157643e+00 9741456
>>> AhatP (000) after 1.87142 seconds with div/v=2.081e-05 dv/v=0.000e+00
>>> and dp/p=0.000e+00 for step 1
>>>
>>> So this is the output from the run. Should I be worried that
>>> dp/p=0.000e+00 ?? How would I tell there is a problem?
>>>
>>
>> Scott,
>>
>> dv/v and dp/p are the relative magnitude of velocity and pressure
>> correction between the Uzawa iterations. They are always 0 for the first
>> iteration. Nothing to worry about them.
>>
>> Your model seems to converge very well at step 0 and 1. Does the
>> problematic geoid occur at these time steps? If yes, what is the min/max
>> of your velocity solution at these time steps? (A convenient GMT command
>> "minmax -H2 *velo*" can query the min/max of all velo files.) If not,
>> can you check the convergence history around the problematic time step?
>>
>>
>> --Eh Tan
>> Staff Scientist
>> Computational Infrastructure for Geodynamics
>> 2750 E. Washington Blvd. Suite 210
>> Pasadena, CA 91107
>> (626) 395-1693
>> http://www.geodynamics.org
>>
>
--
Eh Tan
Staff Scientist
Computational Infrastructure for Geodynamics
2750 E. Washington Blvd. Suite 210
Pasadena, CA 91107
(626) 395-1693
http://www.geodynamics.org
More information about the CIG-MC
mailing list