Changeset 23807

Timestamp:

03/26/19 12:01:25 (6 years ago)

Author:

felicity

Message:

BUG: wrong rate factor conversion from Aoctahedral to Aeffective in BuddJacka

Location:

issm/trunk-jpl/src

Files:

• c/shared/Elements/BuddJacka.cpp (modified) (1 diff)
• m/materials/buddjacka.m (modified) (2 diffs)

issm/trunk-jpl/src/c/shared/Elements/BuddJacka.cpp

@@ -17 +17 @@
 
         if(T<=-40.){
-                B=1e8*(-0.000237326134296*pow(T+50.,3)+ 0.017054655749852*pow(T+50.,2)-0.496435984007500*(T+50.)+7.670967258739796);
+                B=1e9*(-0.000031098521204*pow(T+50.,3)+ 0.002234792114381*pow(T+50.,2)-0.065051516643164*(T+50.)+1.005181071430026);
         }
         else if((-40.<T)  &&  (T<=-35.)){
-         B=1e8*(-0.000237326134296*pow(T+40.,3)+ 0.009934871720961*pow(T+40.,2)-0.226540709299368*(T+40.)+4.174746859353635);
+         B=1e9*(-0.000031098521204*pow(T+40.,3)+ 0.001301836478264*pow(T+40.,2)-0.029685230716715*(T+40.)+0.547046595232583);
         }
         else if((-35.<T)  &&  (T<=-30.)){
-                B=1e8*(-0.000293001369376*pow(T+35.,3)+ 0.006374979706516*pow(T+35.,2)-0.144991452161983*(T+35.)+3.260749339093782);
+                B=1e9*(-0.000038394040864*pow(T+35.,3)+ 0.000835358660205*pow(T+35.,2)-0.018999255024368*(T+35.)+0.427279038455119);
         }
         else if((-30.<T)  &&  (T<=-25.)){
-                B=1e8*(-0.000053702836500*pow(T+30.,3)+ 0.001979959165871*pow(T+30.,2)-0.103216757800049*(T+30.)+2.658541399774723);
+                B=1e9*(-0.000007037062330*pow(T+30.,3)+ 0.000259448047242*pow(T+30.,2)-0.013525221487131*(T+30.)+0.348367474730384);
         }
         else if((-25.<T)  &&  (T<=-20.)){
-                B=1e8*( 0.000006906867543*pow(T+25.,3)+ 0.001174416618375*pow(T+25.,2)-0.087444878878821*(T+25.)+2.185243735358781);
+                B=1e9*( 0.000000905055684*pow(T+25.,3)+ 0.000153892112291*pow(T+25.,2)-0.011458520689465*(T+25.)+0.286347935684521);
         }
         else if((-20.<T)  &&  (T<=-15.)){
-                B=1e8*(-0.000015460250554*pow(T+20.,3)+ 0.001278019631513*pow(T+20.,2)-0.075182697629382*(T+20.)+1.778243114866868);
+                B=1e9*(-0.000002025865930*pow(T+20.,3)+ 0.000167467947546*pow(T+20.,2)-0.009851720390281*(T+20.)+0.233015767004928);
         }
         else if((-15.<T)  &&  (T<=-10.)){
-                B=1e8*(-0.000110386100241*pow(T+15.,3)+ 0.001046115873209*pow(T+15.,2)-0.063562020105772*(T+15.)+1.432347586188582);
+                B=1e9*(-0.000014464671112*pow(T+15.,3)+ 0.000137079958603*pow(T+15.,2)-0.008328980859537*(T+15.)+0.187690630500981);
         }
         else if((-10.<T)  &&  (T<=-5.)){
-                B=1e8*(-0.000108595885218*pow(T+10.,3)+-0.000609675630408*pow(T+10.,2)-0.061379818891767*(T+10.)+1.126892119959808);
+                B=1e9*(-0.000014230086582*pow(T+10.,3)+-0.000079890108083*pow(T+10.,2)-0.008043031606935*(T+10.)+0.147664641279324);
         }
         else if((-5.<T)  &&  (T<=-2.)){
-                B=1e8*( 0.000173187986430*pow(T+5.,3)+-0.002238613908676*pow(T+5.,2)-0.075621266587187*(T+5. )+0.791176649088537);
+                B=1e9*( 0.000022694046251*pow(T+5.,3)+-0.000293341406806*pow(T+5.,2)-0.009909189181377*(T+5. )+0.103673469719891);
         }
         else if((-2.<T)  &&  (T<=-1.)){
-                B=1e8*( 0.000429499435151*pow(T+2.,3)+-0.000679922030808*pow(T+2.,2)-0.084376874405640*(T+2. )+0.548841399782495);
+                B=1e9*( 0.000056280347425*pow(T+2.,3)+-0.000089094990549*pow(T+2.,2)-0.011056498373441*(T+2. )+0.071918568763277);
         }
         else if((-1.<T)){
-                B=1e8*( 0.000429499435151*pow(T+1.,3)+ 0.000608576274646*pow(T+1.,2)-0.084448220161802*(T+1. )+0.464214102781199);
+                B=1e9*( 0.000056280347425*pow(T+1.,3)+ 0.000079746051725*pow(T+1.,2)-0.011065847312265*(T+1. )+0.060829255746712);
         }
         else{
                 /*FIXME: just copying previous case for now.... Felicity?*/
-                B=1e8*( 0.000429499435151*pow(T+1.,3)+ 0.000608576274646*pow(T+1.,2)-0.084448220161802*(T+1. )+0.464214102781199);
+                B=1e9*( 0.000056280347425*pow(T+1.,3)+ 0.000079746051725*pow(T+1.,2)-0.011065847312265*(T+1. )+0.060829255746712);
         }
 

issm/trunk-jpl/src/m/materials/buddjacka.m

@@ -29 +29 @@
 %Temp=BJtable(:,1);
 %Ao=BJtable(:,2)*1e-18; %conversion from MPa^-3 to Pa^-3
-%Ae=Ao*(2/3)^((1-3)/2);
+%Ae=Ao*(2/3)^((3-1)/2);
 %B=Ae.^(-1/3);
 %fittedmodel=fit(Temp,B,'cubicspline');
@@ -43 +43 @@
 %semilogy(Temp,Bcall,'--k',Temp,Bbjall,'k')
 
+
 rigidity=zeros(size(T));
 pos=find(T<=-40);
-rigidity(pos)=1e8*(-0.000237326134296*(T(pos)+50).^3+ 0.017054655749852*(T(pos)+50).^2-0.496435984007500*(T(pos)+50)+7.670967258739796);
+rigidity(pos)=1e9*(-0.000031098521204*(T(pos)+50).^3+ 0.002234792114381*(T(pos)+50).^2-0.065051516643164*(T(pos)+50)+1.005181071430026);
 pos=find(-40<T & T<=-35);
-rigidity(pos)=1e8*(-0.000237326134296*(T(pos)+40).^3+ 0.009934871720961*(T(pos)+40).^2-0.226540709299368*(T(pos)+40)+4.174746859353635);
+rigidity(pos)=1e9*(-0.000031098521204*(T(pos)+40).^3+ 0.001301836478264*(T(pos)+40).^2-0.029685230716715*(T(pos)+40)+0.547046595232583);
 pos=find(-35<T & T<=-30);
-rigidity(pos)=1e8*(-0.000293001369376*(T(pos)+35).^3+ 0.006374979706516*(T(pos)+35).^2-0.144991452161983*(T(pos)+35)+3.260749339093782);
+rigidity(pos)=1e9*(-0.000038394040864*(T(pos)+35).^3+ 0.000835358660205*(T(pos)+35).^2-0.018999255024368*(T(pos)+35)+0.427279038455119);
 pos=find(-30<T & T<=-25);
-rigidity(pos)=1e8*(-0.000053702836500*(T(pos)+30).^3+ 0.001979959165871*(T(pos)+30).^2-0.103216757800049*(T(pos)+30)+2.658541399774723);
+rigidity(pos)=1e9*(-0.000007037062330*(T(pos)+30).^3+ 0.000259448047242*(T(pos)+30).^2-0.013525221487131*(T(pos)+30)+0.348367474730384);
 pos=find(-25<T & T<=-20);
-rigidity(pos)=1e8*( 0.000006906867543*(T(pos)+25).^3+ 0.001174416618375*(T(pos)+25).^2-0.087444878878821*(T(pos)+25)+2.185243735358781);
+rigidity(pos)=1e9*( 0.000000905055684*(T(pos)+25).^3+ 0.000153892112291*(T(pos)+25).^2-0.011458520689465*(T(pos)+25)+0.286347935684521);
 pos=find(-20<T & T<=-15);
-rigidity(pos)=1e8*(-0.000015460250554*(T(pos)+20).^3+ 0.001278019631513*(T(pos)+20).^2-0.075182697629382*(T(pos)+20)+1.778243114866868);
+rigidity(pos)=1e9*(-0.000002025865930*(T(pos)+20).^3+ 0.000167467947546*(T(pos)+20).^2-0.009851720390281*(T(pos)+20)+0.233015767004928);
 pos=find(-15<T & T<=-10);
-rigidity(pos)=1e8*(-0.000110386100241*(T(pos)+15).^3+ 0.001046115873209*(T(pos)+15).^2-0.063562020105772*(T(pos)+15)+1.432347586188582);
+rigidity(pos)=1e9*(-0.000014464671112*(T(pos)+15).^3+ 0.000137079958603*(T(pos)+15).^2-0.008328980859537*(T(pos)+15)+0.187690630500981);
 pos=find(-10<T & T<=-5);
-rigidity(pos)=1e8*(-0.000108595885218*(T(pos)+10).^3+-0.000609675630408*(T(pos)+10).^2-0.061379818891767*(T(pos)+10)+1.126892119959808);
+rigidity(pos)=1e9*(-0.000014230086582*(T(pos)+10).^3+-0.000079890108083*(T(pos)+10).^2-0.008043031606935*(T(pos)+10)+0.147664641279324);
 pos=find(-5<T & T<=-2);
-rigidity(pos)=1e8*( 0.000173187986430*(T(pos)+5).^3+-0.002238613908676*(T(pos)+5).^2-0.075621266587187*(T(pos)+5)+0.791176649088537);
+rigidity(pos)=1e9*( 0.000022694046251*(T(pos)+5).^3+-0.000293341406806*(T(pos)+5).^2-0.009909189181377*(T(pos)+5)+0.103673469719891);
 pos=find(-2<T & T<=-1);
-rigidity(pos)=1e8*( 0.000429499435151*(T(pos)+2).^3+-0.000679922030808*(T(pos)+2).^2-0.084376874405640*(T(pos)+2)+0.548841399782495);
+rigidity(pos)=1e9*( 0.000056280347425*(T(pos)+2).^3+-0.000089094990549*(T(pos)+2).^2-0.011056498373441*(T(pos)+2)+0.071918568763277);
 pos=find(-1<T);
-rigidity(pos)=1e8*( 0.000429499435151*(T(pos)+1).^3+ 0.000608576274646*(T(pos)+1).^2-0.084448220161802*(T(pos)+1)+0.464214102781199);
+rigidity(pos)=1e9*( 0.000056280347425*(T(pos)+1).^3+ 0.000079746051725*(T(pos)+1).^2-0.011065847312265*(T(pos)+1)+0.060829255746712);
 
 %Now make sure that rigidity is positive