[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