Changeset 22576


Ignore:
Timestamp:
Mar 20, 2018, 11:33:59 AM (9 years ago)
Author:
Mathieu Morlighem
Message:

CHG: no more viscosity_overhoot, makes models crash and they don't go facter

Location:
issm/trunk-jpl/src
Files:
12 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/c/analyses/StressbalanceSIAAnalysis.cpp

    r22307 r22576  
    621621        }
    622622
    623         /*Now, we have to move the previous Vx and Vy inputs  to old
    624          * status, otherwise, we'll wipe them off: */
    625         element->InputChangeName(VxEnum,VxPicardEnum);
    626         element->InputChangeName(VyEnum,VyPicardEnum);
    627         element->InputChangeName(PressureEnum,PressurePicardEnum);
    628 
    629623        /*Add vx and vy as inputs to the tria element: */
    630624        element->AddInput(VxEnum,vx,P1Enum);

  • issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.cpp

    r22535 r22576  
    567567                for(i=0;i<numnodes;i++) pressure[i]=rho_ice*g*(surface[i]-xyz_list[i*3+2]);
    568568        }
    569 
    570         /*Now, we have to move the previous Vz inputs to old
    571          * status, otherwise, we'll wipe them off and add the new inputs: */
    572         element->InputChangeName(VzEnum,VzPicardEnum);
    573 
    574569        if(approximation!=HOFSApproximationEnum && approximation!=SSAFSApproximationEnum){
    575                 element->InputChangeName(PressureEnum,PressurePicardEnum);
    576570                element->AddInput(PressureEnum,pressure,element->GetElementType());
    577571        }

  • issm/trunk-jpl/src/c/classes/Elements/Element.cpp

    r22539 r22576  
    19951995                                name==LoadingforceYEnum ||
    19961996                                name==VelEnum ||
    1997                                 name==VxPicardEnum ||
    1998                                 name==VyPicardEnum ||
    19991997                                name==DeviatoricStressxxEnum ||
    20001998                                name==DeviatoricStressyyEnum ||

  • issm/trunk-jpl/src/c/shared/Enum/EnumDefinitions.h

    r22558 r22576  
    550550        MisfitEnum,
    551551        PressureEnum,
    552         PressurePicardEnum,
    553552        AndroidFrictionCoefficientEnum,
    554553        ResetPenaltiesEnum,
    … …  
    573572        VxAverageEnum,
    574573        VxEnum,
    575         VxPicardEnum,
    576574        VyAverageEnum,
    577575        VyEnum,
    578         VyPicardEnum,
    579576        VzEnum,
    580577        VzSSAEnum,
    581578        VzHOEnum,
    582         VzPicardEnum,
    583579        VzFSEnum,
    584580        VxMeshEnum,

  • issm/trunk-jpl/src/c/shared/Enum/EnumToStringx.cpp

    r22558 r22576  
    548548                case MisfitEnum : return "Misfit";
    549549                case PressureEnum : return "Pressure";
    550                 case PressurePicardEnum : return "PressurePicard";
    551550                case AndroidFrictionCoefficientEnum : return "AndroidFrictionCoefficient";
    552551                case ResetPenaltiesEnum : return "ResetPenalties";
    … …  
    571570                case VxAverageEnum : return "VxAverage";
    572571                case VxEnum : return "Vx";
    573                 case VxPicardEnum : return "VxPicard";
    574572                case VyAverageEnum : return "VyAverage";
    575573                case VyEnum : return "Vy";
    576                 case VyPicardEnum : return "VyPicard";
    577574                case VzEnum : return "Vz";
    578575                case VzSSAEnum : return "VzSSA";
    579576                case VzHOEnum : return "VzHO";
    580                 case VzPicardEnum : return "VzPicard";
    581577                case VzFSEnum : return "VzFS";
    582578                case VxMeshEnum : return "VxMesh";

  • issm/trunk-jpl/src/c/shared/Enum/StringToEnumx.cpp

    r22558 r22576  
    560560              else if (strcmp(name,"Misfit")==0) return MisfitEnum;
    561561              else if (strcmp(name,"Pressure")==0) return PressureEnum;
    562               else if (strcmp(name,"PressurePicard")==0) return PressurePicardEnum;
    563562              else if (strcmp(name,"AndroidFrictionCoefficient")==0) return AndroidFrictionCoefficientEnum;
    564563              else if (strcmp(name,"ResetPenalties")==0) return ResetPenaltiesEnum;
    … …  
    583582              else if (strcmp(name,"VxAverage")==0) return VxAverageEnum;
    584583              else if (strcmp(name,"Vx")==0) return VxEnum;
    585               else if (strcmp(name,"VxPicard")==0) return VxPicardEnum;
    586584              else if (strcmp(name,"VyAverage")==0) return VyAverageEnum;
    587585              else if (strcmp(name,"Vy")==0) return VyEnum;
    588               else if (strcmp(name,"VyPicard")==0) return VyPicardEnum;
    589586              else if (strcmp(name,"Vz")==0) return VzEnum;
    590587              else if (strcmp(name,"VzSSA")==0) return VzSSAEnum;
    591588              else if (strcmp(name,"VzHO")==0) return VzHOEnum;
    592               else if (strcmp(name,"VzPicard")==0) return VzPicardEnum;
    593589              else if (strcmp(name,"VzFS")==0) return VzFSEnum;
    594590              else if (strcmp(name,"VxMesh")==0) return VxMeshEnum;
    … …  
    629625              else if (strcmp(name,"DeviatoricStressxz")==0) return DeviatoricStressxzEnum;
    630626              else if (strcmp(name,"DeviatoricStressyy")==0) return DeviatoricStressyyEnum;
     627              else if (strcmp(name,"DeviatoricStressyz")==0) return DeviatoricStressyzEnum;
     628              else if (strcmp(name,"DeviatoricStresszz")==0) return DeviatoricStresszzEnum;
     629              else if (strcmp(name,"DeviatoricStresseffective")==0) return DeviatoricStresseffectiveEnum;
     630              else if (strcmp(name,"LambdaS")==0) return LambdaSEnum;
    631631         else stage=6;
    632632   }
    633633   if(stage==6){
    634               if (strcmp(name,"DeviatoricStressyz")==0) return DeviatoricStressyzEnum;
    635               else if (strcmp(name,"DeviatoricStresszz")==0) return DeviatoricStresszzEnum;
    636               else if (strcmp(name,"DeviatoricStresseffective")==0) return DeviatoricStresseffectiveEnum;
    637               else if (strcmp(name,"LambdaS")==0) return LambdaSEnum;
    638               else if (strcmp(name,"StrainRate")==0) return StrainRateEnum;
     634              if (strcmp(name,"StrainRate")==0) return StrainRateEnum;
    639635              else if (strcmp(name,"StrainRatexx")==0) return StrainRatexxEnum;
    640636              else if (strcmp(name,"StrainRatexy")==0) return StrainRatexyEnum;
    … …  
    752748              else if (strcmp(name,"Outputdefinition89")==0) return Outputdefinition89Enum;
    753749              else if (strcmp(name,"Outputdefinition90")==0) return Outputdefinition90Enum;
     750              else if (strcmp(name,"Outputdefinition91")==0) return Outputdefinition91Enum;
     751              else if (strcmp(name,"Outputdefinition92")==0) return Outputdefinition92Enum;
     752              else if (strcmp(name,"Outputdefinition93")==0) return Outputdefinition93Enum;
     753              else if (strcmp(name,"Outputdefinition94")==0) return Outputdefinition94Enum;
    754754         else stage=7;
    755755   }
    756756   if(stage==7){
    757               if (strcmp(name,"Outputdefinition91")==0) return Outputdefinition91Enum;
    758               else if (strcmp(name,"Outputdefinition92")==0) return Outputdefinition92Enum;
    759               else if (strcmp(name,"Outputdefinition93")==0) return Outputdefinition93Enum;
    760               else if (strcmp(name,"Outputdefinition94")==0) return Outputdefinition94Enum;
    761               else if (strcmp(name,"Outputdefinition95")==0) return Outputdefinition95Enum;
     757              if (strcmp(name,"Outputdefinition95")==0) return Outputdefinition95Enum;
    762758              else if (strcmp(name,"Outputdefinition96")==0) return Outputdefinition96Enum;
    763759              else if (strcmp(name,"Outputdefinition97")==0) return Outputdefinition97Enum;
    … …  
    875871              else if (strcmp(name,"SealevelriseTransitions")==0) return SealevelriseTransitionsEnum;
    876872              else if (strcmp(name,"SealevelriseRequestedOutputs")==0) return SealevelriseRequestedOutputsEnum;
     873              else if (strcmp(name,"SealevelriseNumRequestedOutputs")==0) return SealevelriseNumRequestedOutputsEnum;
     874              else if (strcmp(name,"LoveNfreq")==0) return LoveNfreqEnum;
     875              else if (strcmp(name,"LoveFrequencies")==0) return LoveFrequenciesEnum;
     876              else if (strcmp(name,"LoveShNmax")==0) return LoveShNmaxEnum;
    877877         else stage=8;
    878878   }
    879879   if(stage==8){
    880               if (strcmp(name,"SealevelriseNumRequestedOutputs")==0) return SealevelriseNumRequestedOutputsEnum;
    881               else if (strcmp(name,"LoveNfreq")==0) return LoveNfreqEnum;
    882               else if (strcmp(name,"LoveFrequencies")==0) return LoveFrequenciesEnum;
    883               else if (strcmp(name,"LoveShNmax")==0) return LoveShNmaxEnum;
    884               else if (strcmp(name,"LoveShNmin")==0) return LoveShNminEnum;
     880              if (strcmp(name,"LoveShNmin")==0) return LoveShNminEnum;
    885881              else if (strcmp(name,"LoveG0")==0) return LoveG0Enum;
    886882              else if (strcmp(name,"LoveR0")==0) return LoveR0Enum;
    … …  
    998994              else if (strcmp(name,"SegInput")==0) return SegInputEnum;
    999995              else if (strcmp(name,"Tria")==0) return TriaEnum;
     996              else if (strcmp(name,"TriaInput")==0) return TriaInputEnum;
     997              else if (strcmp(name,"Tetra")==0) return TetraEnum;
     998              else if (strcmp(name,"TetraInput")==0) return TetraInputEnum;
     999              else if (strcmp(name,"Penta")==0) return PentaEnum;
    10001000         else stage=9;
    10011001   }
    10021002   if(stage==9){
    1003               if (strcmp(name,"TriaInput")==0) return TriaInputEnum;
    1004               else if (strcmp(name,"Tetra")==0) return TetraEnum;
    1005               else if (strcmp(name,"TetraInput")==0) return TetraInputEnum;
    1006               else if (strcmp(name,"Penta")==0) return PentaEnum;
    1007               else if (strcmp(name,"PentaInput")==0) return PentaInputEnum;
     1003              if (strcmp(name,"PentaInput")==0) return PentaInputEnum;
    10081004              else if (strcmp(name,"Vertex")==0) return VertexEnum;
    10091005              else if (strcmp(name,"VertexPId")==0) return VertexPIdEnum;
    … …  
    11211117              else if (strcmp(name,"P1xP2")==0) return P1xP2Enum;
    11221118              else if (strcmp(name,"P1xP3")==0) return P1xP3Enum;
     1119              else if (strcmp(name,"P1xP4")==0) return P1xP4Enum;
     1120              else if (strcmp(name,"P2xP4")==0) return P2xP4Enum;
     1121              else if (strcmp(name,"P1P1")==0) return P1P1Enum;
     1122              else if (strcmp(name,"P1P1GLS")==0) return P1P1GLSEnum;
    11231123         else stage=10;
    11241124   }
    11251125   if(stage==10){
    1126               if (strcmp(name,"P1xP4")==0) return P1xP4Enum;
    1127               else if (strcmp(name,"P2xP4")==0) return P2xP4Enum;
    1128               else if (strcmp(name,"P1P1")==0) return P1P1Enum;
    1129               else if (strcmp(name,"P1P1GLS")==0) return P1P1GLSEnum;
    1130               else if (strcmp(name,"MINI")==0) return MINIEnum;
     1126              if (strcmp(name,"MINI")==0) return MINIEnum;
    11311127              else if (strcmp(name,"MINIcondensed")==0) return MINIcondensedEnum;
    11321128              else if (strcmp(name,"TaylorHood")==0) return TaylorHoodEnum;

  • issm/trunk-jpl/src/m/classes/snowpack.m

    r20902 r22576  
    463463                        WriteData(fid,prefix,'object',self,'class','snowpack','fieldname','isnewton','format','Integer');
    464464                        WriteData(fid,prefix,'object',self,'class','snowpack','fieldname','FSreconditioning','format','Double');
    465                         WriteData(fid,prefix,'object',self,'class','snowpack','fieldname','viscosity_overshoot','format','Double');
    466465                        WriteData(fid,prefix,'object',self,'class','snowpack','fieldname','maxiter','format','Integer');
    467466                        WriteData(fid,prefix,'object',self,'class','snowpack','fieldname','shelf_dampening','format','Integer');

  • issm/trunk-jpl/src/m/classes/stressbalance.js

    r21065 r22576  
    2121                //Penalty factor applied kappa=max(stiffness matrix)*10^penalty_factor
    2222                this.penalty_factor=3;
    23 
    24                 //coefficient to update the viscosity between each iteration of
    25                 //a stressbalance according to the following formula
    26                 //viscosity(n)=viscosity(n)+viscosity_overshoot(viscosity(n)-viscosity(n-1))
    27                 this.viscosity_overshoot=0;
    2823
    2924                //Stop the iterations of rift if below a threshold
    … …  
    4843                fielddisplay(this,'isnewton',"0: Picard's fixed point, 1: Newton's method, 2: hybrid");
    4944                fielddisplay(this,'maxiter','maximum number of nonlinear iterations');
    50                 fielddisplay(this,'viscosity_overshoot','over-shooting constant new=new+C*(new-old)');
    5145
    5246                console.log(sprintf('\n      %s','boundary conditions:'));
    … …  
    9791                checkfield(md,'fieldname','stressbalance.isnewton','numel',[1],'values',[0, 1, 2]);
    9892                checkfield(md,'fieldname','stressbalance.FSreconditioning','size',[1, 1],'NaN',1,'Inf',1);
    99                 checkfield(md,'fieldname','stressbalance.viscosity_overshoot','size',[1, 1],'NaN',1,'Inf',1);
    10093                checkfield(md,'fieldname','stressbalance.maxiter','size',[1, 1],'>=',1);
    10194                checkfield(md,'fieldname','stressbalance.referential','size',[md.mesh.numberofvertices, 6]);
    … …  
    158151                WriteData(fid,prefix,'object',this,'class','stressbalance','fieldname','isnewton','format','Integer');
    159152                WriteData(fid,prefix,'object',this,'class','stressbalance','fieldname','FSreconditioning','format','Double');
    160                 WriteData(fid,prefix,'object',this,'class','stressbalance','fieldname','viscosity_overshoot','format','Double');
    161153                WriteData(fid,prefix,'object',this,'class','stressbalance','fieldname','maxiter','format','Integer');
    162154                WriteData(fid,prefix,'object',this,'class','stressbalance','fieldname','shelf_dampening','format','Integer');
    … …  
    214206        this.isnewton               = 0;
    215207        this.FSreconditioning       = 0;
    216         this.viscosity_overshoot    = 0;
    217208        this.maxiter                = 0;
    218209        this.shelf_dampening        = 0;

  • issm/trunk-jpl/src/m/classes/stressbalance.m

    r21049 r22576  
    1414                isnewton               = 0;
    1515                FSreconditioning       = 0;
    16                 viscosity_overshoot    = 0;
    1716                maxiter                = 0;
    1817                shelf_dampening        = 0;
    … …  
    6867                         self.penalty_factor=3;
    6968
    70                          %coefficient to update the viscosity between each iteration of
    71                          %a stressbalance according to the following formula
    72                          %viscosity(n)=viscosity(n)+viscosity_overshoot(viscosity(n)-viscosity(n-1))
    73                          self.viscosity_overshoot=0;
    74 
    7569                         %Stop the iterations of rift if below a threshold
    7670                         self.rift_penalty_threshold=0;
    … …  
    9993                        md = checkfield(md,'fieldname','stressbalance.isnewton','numel',[1],'values',[0 1 2]);
    10094                        md = checkfield(md,'fieldname','stressbalance.FSreconditioning','size',[1 1],'NaN',1,'Inf',1);
    101                         md = checkfield(md,'fieldname','stressbalance.viscosity_overshoot','size',[1 1],'NaN',1,'Inf',1);
    10295                        md = checkfield(md,'fieldname','stressbalance.maxiter','size',[1 1],'>=',1);
    10396                        md = checkfield(md,'fieldname','stressbalance.referential','size',[md.mesh.numberofvertices 6]);
    … …  
    150143                        fielddisplay(self,'isnewton','0: Picard''s fixed point, 1: Newton''s method, 2: hybrid');
    151144                        fielddisplay(self,'maxiter','maximum number of nonlinear iterations');
    152                         fielddisplay(self,'viscosity_overshoot','over-shooting constant new=new+C*(new-old)');
    153145
    154146                        disp(sprintf('\n      %s','boundary conditions:'));
    … …  
    187179                        WriteData(fid,prefix,'object',self,'class','stressbalance','fieldname','isnewton','format','Integer');
    188180                        WriteData(fid,prefix,'object',self,'class','stressbalance','fieldname','FSreconditioning','format','Double');
    189                         WriteData(fid,prefix,'object',self,'class','stressbalance','fieldname','viscosity_overshoot','format','Double');
    190181                        WriteData(fid,prefix,'object',self,'class','stressbalance','fieldname','maxiter','format','Integer');
    191182                        WriteData(fid,prefix,'object',self,'class','stressbalance','fieldname','shelf_dampening','format','Integer');
    … …  
    220211                        writejsdouble(fid,[modelname '.stressbalance.isnewton'],self.isnewton);
    221212                        writejsdouble(fid,[modelname '.stressbalance.FSreconditioning'],self.FSreconditioning);
    222                         writejsdouble(fid,[modelname '.stressbalance.viscosity_overshoot'],self.viscosity_overshoot);
    223213                        writejsdouble(fid,[modelname '.stressbalance.maxiter'],self.maxiter);
    224214                        writejsdouble(fid,[modelname '.stressbalance.shelf_dampening'],self.shelf_dampening);

  • issm/trunk-jpl/src/m/classes/stressbalance.py

    r21303 r22576  
    1717
    1818        def __init__(self): # {{{
    19                 self.spcvx                    = float('NaN')
    20                 self.spcvy                    = float('NaN')
    21                 self.spcvz                    = float('NaN')
    22                 self.restol                   = 0
    23                 self.reltol                   = 0
    24                 self.abstol                   = 0
    25                 self.isnewton                 = 0
    26                 self.FSreconditioning     = 0
    27                 self.viscosity_overshoot      = 0
    28                 self.icefront                 = float('NaN')
    29                 self.maxiter                  = 0
    30                 self.shelf_dampening          = 0
    31                 self.vertex_pairing           = float('NaN')
    32                 self.penalty_factor           = float('NaN')
    33                 self.rift_penalty_lock        = float('NaN')
    34                 self.rift_penalty_threshold   = 0
    35                 self.referential              = float('NaN')
    36                 self.loadingforce             = float('NaN')
    37                 self.requested_outputs        = []
     19                self.spcvx                  = float('NaN')
     20                self.spcvy                  = float('NaN')
     21                self.spcvz                  = float('NaN')
     22                self.restol                 = 0
     23                self.reltol                 = 0
     24                self.abstol                 = 0
     25                self.isnewton               = 0
     26                self.FSreconditioning       = 0
     27                self.icefront               = float('NaN')
     28                self.maxiter                = 0
     29                self.shelf_dampening        = 0
     30                self.vertex_pairing         = float('NaN')
     31                self.penalty_factor         = float('NaN')
     32                self.rift_penalty_lock      = float('NaN')
     33                self.rift_penalty_threshold = 0
     34                self.referential            = float('NaN')
     35                self.loadingforce           = float('NaN')
     36                self.requested_outputs      = []
    3837
    3938                #set defaults
    … …  
    5049                string="%s\n%s"%(string,fielddisplay(self,'isnewton',"0: Picard's fixed point, 1: Newton's method, 2: hybrid"))
    5150                string="%s\n%s"%(string,fielddisplay(self,'maxiter','maximum number of nonlinear iterations'))
    52                 string="%s\n%s"%(string,fielddisplay(self,'viscosity_overshoot','over-shooting constant new=new+C*(new-old)'))
    5351
    5452                string="%s\n%s"%(string,'\n      boundary conditions:')
    … …  
    9896                #Penalty factor applied kappa=max(stiffness matrix)*10^penalty_factor
    9997                self.penalty_factor=3
    100 
    101                 #coefficient to update the viscosity between each iteration of
    102                 #a stressbalance according to the following formula
    103                 #viscosity(n)=viscosity(n)+viscosity_overshoot(viscosity(n)-viscosity(n-1))
    104                 self.viscosity_overshoot=0
    10598
    10699                #Stop the iterations of rift if below a threshold
    … …  
    141134                md = checkfield(md,'fieldname','stressbalance.isnewton','numel',[1],'values',[0,1,2])
    142135                md = checkfield(md,'fieldname','stressbalance.FSreconditioning','size',[1],'NaN',1,'Inf',1)
    143                 md = checkfield(md,'fieldname','stressbalance.viscosity_overshoot','size',[1],'NaN',1,'Inf',1)
    144136                md = checkfield(md,'fieldname','stressbalance.maxiter','size',[1],'>=',1)
    145137                md = checkfield(md,'fieldname','stressbalance.referential','size',[md.mesh.numberofvertices,6])
    … …  
    187179                WriteData(fid,prefix,'object',self,'class','stressbalance','fieldname','isnewton','format','Integer')
    188180                WriteData(fid,prefix,'object',self,'class','stressbalance','fieldname','FSreconditioning','format','Double')
    189                 WriteData(fid,prefix,'object',self,'class','stressbalance','fieldname','viscosity_overshoot','format','Double')
    190181                WriteData(fid,prefix,'object',self,'class','stressbalance','fieldname','maxiter','format','Integer')
    191182                WriteData(fid,prefix,'object',self,'class','stressbalance','fieldname','shelf_dampening','format','Integer')

  • issm/trunk-jpl/src/py3/classes/stressbalance.py

    r19903 r22576  
    1818
    1919        def __init__(self): # {{{
    20                 self.spcvx                    = float('NaN')
    21                 self.spcvy                    = float('NaN')
    22                 self.spcvz                    = float('NaN')
    23                 self.restol                   = 0
    24                 self.reltol                   = 0
    25                 self.abstol                   = 0
    26                 self.isnewton                 = 0
    27                 self.FSreconditioning     = 0
    28                 self.viscosity_overshoot      = 0
    29                 self.icefront                 = float('NaN')
    30                 self.maxiter                  = 0
    31                 self.shelf_dampening          = 0
    32                 self.vertex_pairing           = float('NaN')
    33                 self.penalty_factor           = float('NaN')
    34                 self.rift_penalty_lock        = float('NaN')
    35                 self.rift_penalty_threshold   = 0
    36                 self.referential              = float('NaN')
    37                 self.loadingforce             = float('NaN')
    38                 self.requested_outputs        = []
     20                self.spcvx                  = float('NaN')
     21                self.spcvy                  = float('NaN')
     22                self.spcvz                  = float('NaN')
     23                self.restol                 = 0
     24                self.reltol                 = 0
     25                self.abstol                 = 0
     26                self.isnewton               = 0
     27                self.FSreconditioning       = 0
     28                self.icefront               = float('NaN')
     29                self.maxiter                = 0
     30                self.shelf_dampening        = 0
     31                self.vertex_pairing         = float('NaN')
     32                self.penalty_factor         = float('NaN')
     33                self.rift_penalty_lock      = float('NaN')
     34                self.rift_penalty_threshold = 0
     35                self.referential            = float('NaN')
     36                self.loadingforce           = float('NaN')
     37                self.requested_outputs      = []
    3938
    4039                #set defaults
    … …  
    5150                string="%s\n%s"%(string,fielddisplay(self,'isnewton',"0: Picard's fixed point, 1: Newton's method, 2: hybrid"))
    5251                string="%s\n%s"%(string,fielddisplay(self,'maxiter','maximum number of nonlinear iterations'))
    53                 string="%s\n%s"%(string,fielddisplay(self,'viscosity_overshoot','over-shooting constant new=new+C*(new-old)'))
    5452
    5553                string="%s\n%s"%(string,'\n      boundary conditions:')
    … …  
    10098                self.penalty_factor=3
    10199
    102                 #coefficient to update the viscosity between each iteration of
    103                 #a stressbalance according to the following formula
    104                 #viscosity(n)=viscosity(n)+viscosity_overshoot(viscosity(n)-viscosity(n-1))
    105                 self.viscosity_overshoot=0
    106 
    107100                #Stop the iterations of rift if below a threshold
    108101                self.rift_penalty_threshold=0
    … …  
    142135                md = checkfield(md,'fieldname','stressbalance.isnewton','numel',[1],'values',[0,1,2])
    143136                md = checkfield(md,'fieldname','stressbalance.FSreconditioning','size',[1],'NaN',1,'Inf',1)
    144                 md = checkfield(md,'fieldname','stressbalance.viscosity_overshoot','size',[1],'NaN',1,'Inf',1)
    145137                md = checkfield(md,'fieldname','stressbalance.maxiter','size',[1],'>=',1)
    146138                md = checkfield(md,'fieldname','stressbalance.referential','size',[md.mesh.numberofvertices,6])
    … …  
    186178                WriteData(fid,'object',self,'class','stressbalance','fieldname','isnewton','format','Integer')
    187179                WriteData(fid,'object',self,'class','stressbalance','fieldname','FSreconditioning','format','Double')
    188                 WriteData(fid,'object',self,'class','stressbalance','fieldname','viscosity_overshoot','format','Double')
    189180                WriteData(fid,'object',self,'class','stressbalance','fieldname','maxiter','format','Integer')
    190181                WriteData(fid,'object',self,'class','stressbalance','fieldname','shelf_dampening','format','Integer')

  • issm/trunk-jpl/src/py3/modifier.sh

    r19895 r22576  
    260260sed -i 's/(md,'fieldname','stressbalance.isnewton','numel',[1],'values',[0,1,2])/(md,fieldname='stressbalance.isnewton',numel=[1],values=[0,1,2])/g' ./classes/stressbalance.py:               
    261261sed -i 's/(md,'fieldname','stressbalance.FSreconditioning',size=[1],'NaN',1)/(md,fieldname='stressbalance.FSreconditioning','size',[1],NaN=1)/g' ./classes/stressbalance.py:           
    262 sed -i 's/(md,'fieldname','stressbalance.viscosity_overshoot',size=[1],'NaN',1)/(md,fieldname='stressbalance.viscosity_overshoot','size',[1],NaN=1)/g' ./classes/stressbalance.py:             
    263262sed -i 's/(md,'fieldname','stressbalance.maxiter',size=[1],'>=',1)/(md,fieldname='stressbalance.maxiter','size',[1],ge=1)/g' ./classes/stressbalance.py:               
    264263sed -i 's/(md,'fieldname','stressbalance.referential',size=[md.mesh.numberofvertices,6])/(md,fieldname='stressbalance.referential','size',[md.mesh.numberofvertices,6])/g' ./classes/stressbalance.py:         
Note: See TracChangeset for help on using the changeset viewer.