Changeset 16181


Ignore:
Timestamp:
09/19/13 11:02:41 (12 years ago)
Author:
Eric.Larour
Message:

NEW: hooked up damage core into ISSM, created 271 solution to start running it. More remains to be done inside the damage_core to hook up to the
Element Tria CreateKMatrix and CreatePVector.

Location:
issm/trunk-jpl
Files:
45 added
1 deleted
19 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/m4/issm_options.m4

    r16173 r16181  
    17111711        AC_MSG_RESULT($HAVE_BALANCED)
    17121712        dnl }}}
     1713        dnl with-damage{{{
     1714        AC_ARG_WITH([damage],
     1715                AS_HELP_STRING([--with-damage = YES], [compile with damage capabilities (default is yes)]),
     1716                [DAMAGE=$withval],[DAMAGE=yes])
     1717        AC_MSG_CHECKING(for damage capability compilation)
     1718
     1719        if test "x$DAMAGE" = "xyes"; then
     1720
     1721                dnl defaults
     1722                HAVE_DAMAGE=yes
     1723
     1724                AC_DEFINE([_HAVE_DAMAGE_],[1],[with damage capability])
     1725        else
     1726                HAVE_DAMAGE=no
     1727        fi
     1728        AM_CONDITIONAL([DAMAGE], [test x$HAVE_DAMAGE = xyes])
     1729        AC_MSG_RESULT($HAVE_DAMAGE)
     1730        dnl }}}
    17131731        dnl with-responses{{{
    17141732        AC_ARG_WITH([responses],
  • issm/trunk-jpl/src/c/Makefile.am

    r16167 r16181  
    546546
    547547#}}}
     548#Damage sources  {{{
     549damage_sources =  ./analyses/damage_core.cpp\
     550                                        ./modules/ModelProcessorx/Damage/UpdateElementsDamage.cpp\
     551                                        ./modules/ModelProcessorx/Damage/CreateNodesDamage.cpp \
     552                                        ./modules/ModelProcessorx/Damage/CreateConstraintsDamage.cpp\
     553                                        ./modules/ModelProcessorx/Damage/CreateParametersDamage.cpp\
     554                                        ./modules/ModelProcessorx/Damage/CreateLoadsDamage.cpp
     555
     556#}}}
    548557#Groundingline sources  {{{
    549558groundingline_sources= ./modules/GroundinglineMigrationx/GroundinglineMigrationx.cpp\
     
    950959issm_sources +=  $(pkriging_sources)
    951960endif
     961
     962if DAMAGE
     963issm_sources +=  $(damage_sources)
     964endif
     965
    952966#}}}
    953967#Library flags and sources {{{
  • issm/trunk-jpl/src/c/analyses/AnalysisConfiguration.cpp

    r16007 r16181  
    113113                        analyses[0]=GiaAnalysisEnum;
    114114                        break;
     115               
     116                case DamageEvolutionSolutionEnum:
     117                        numanalyses=1;
     118                        analyses=xNew<int>(numanalyses);
     119                        analyses[0]=DamageEvolutionAnalysisEnum;
     120                        break;
    115121
    116122                case TransientSolutionEnum:
  • issm/trunk-jpl/src/c/analyses/CorePointerFromSolutionEnum.cpp

    r16007 r16181  
    115115                        #endif
    116116                        break;
     117                case DamageEvolutionSolutionEnum:
     118                        #ifdef _HAVE_DAMAGE_
     119                        solutioncore=&damage_core;
     120                        #else
     121                        _error_("ISSM was not compiled with damage evolution capabilities. Exiting");
     122                        #endif
     123                        break;
    117124
    118125                default:
  • issm/trunk-jpl/src/c/analyses/analyses.h

    r16007 r16181  
    3838void dummy_core(FemModel* femmodel);
    3939void gia_core(FemModel* femmodel);
     40void damage_core(FemModel* femmodel);
    4041IssmDouble objectivefunction(IssmDouble search_scalar,OptArgs* optargs);
    4142
  • issm/trunk-jpl/src/c/classes/Node.cpp

    r16164 r16181  
    9292                                analysis_type==BalancethicknessAnalysisEnum ||
    9393                                analysis_type==HydrologyDCInefficientAnalysisEnum ||
     94                                analysis_type==DamageEvolutionAnalysisEnum ||
    9495                                analysis_type==HydrologyDCEfficientAnalysisEnum
    9596                                ){
  • issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateDataSets.cpp

    r16142 r16181  
    130130                        break;
    131131                #endif
     132
     133                #ifdef _HAVE_DAMAGE_
     134                case DamageEvolutionAnalysisEnum:
     135                        CreateNodesDamage(pnodes, iomodel);
     136                        CreateConstraintsDamage(pconstraints,iomodel);
     137                        CreateLoadsDamage(ploads,iomodel);
     138                        UpdateElementsDamage(elements,iomodel,analysis_counter,analysis_type);
     139                        break;
     140                #endif
     141
    132142
    133143                #ifdef _HAVE_SLOPE_
  • issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp

    r16164 r16181  
    109109        parameters->AddObject(iomodel->CopyConstantObject(SurfaceforcingsIssmbgradientsEnum));
    110110        parameters->AddObject(iomodel->CopyConstantObject(GiaCrossSectionShapeEnum));
    111 
     111       
    112112        iomodel->Constant(&ispdd,SurfaceforcingsIspddEnum);
    113113        if(ispdd){
     
    237237        CreateParametersHydrologyDCEfficient(&parameters,iomodel,solution_type,analysis_type);
    238238        #endif
     239       
     240        #ifdef _HAVE_DAMAGE_
     241        CreateParametersDamage(&parameters,iomodel,solution_type,analysis_type);
     242        #endif
    239243
    240244        /*Before returning, create parameters in case we are running Qmu or control types runs: */
  • issm/trunk-jpl/src/c/modules/ModelProcessorx/DistributeNumDofs.cpp

    r16007 r16181  
    112112                        numdofs=1;
    113113                        break;
     114                case DamageEvolutionAnalysisEnum:
     115                        numdofs=1;
     116                        break;
    114117                case BalancethicknessAnalysisEnum:
    115118                        numdofs=1;
  • issm/trunk-jpl/src/c/modules/ModelProcessorx/ModelProcessorx.h

    r16142 r16181  
    5151void CreateLoadsGia(Loads** ploads, IoModel* iomodel);
    5252void UpdateElementsGia(Elements* elements,IoModel* iomodel,int analysis_counter,int analysis_type);
     53#endif
     54
     55#ifdef _HAVE_DAMAGE_
     56/*damage*/
     57void CreateNodesDamage(Nodes** pnodes,IoModel* iomodel);
     58void CreateConstraintsDamage(Constraints** pconstraints,IoModel* iomodel);
     59void CreateLoadsDamage(Loads** ploads, IoModel* iomodel);
     60void UpdateElementsDamage(Elements* elements,IoModel* iomodel,int analysis_counter,int analysis_type);
     61void CreateParametersDamage(Parameters** pparameters,IoModel* iomodel,int solution_type,int analysis_type);
    5362#endif
    5463
  • issm/trunk-jpl/src/c/shared/Enum/EnumDefinitions.h

    r16167 r16181  
    153153        MaterialsRheologyNEnum,
    154154        DamageDEnum,
     155        QmuDamageDEnum,
    155156        DamageDbarEnum,
    156         DamageIsdamageEnum,
    157157        DamageLawEnum,
    158158        DamageC1Enum,
     
    160160        DamageC3Enum,
    161161        DamageC4Enum,
     162        DamageHealingEnum,
    162163        DamageStressThresholdEnum,
     164        DamageStabilizationEnum,
     165        DamageSpcdamageEnum,
    163166        MaterialsRhoIceEnum,
    164167        MaterialsRhoWaterEnum,
     
    264267        TransientIsthermalEnum,
    265268        TransientIsgiaEnum,
     269        TransientIsdamageEnum,
    266270        TransientNumRequestedOutputsEnum,
    267271        TransientRequestedOutputsEnum,
     
    287291        BedSlopeXAnalysisEnum,
    288292        BedSlopeYAnalysisEnum,
     293        DamageEvolutionSolutionEnum,
     294        DamageEvolutionAnalysisEnum,
    289295        StressbalanceAnalysisEnum,
    290296        StressbalanceSIAAnalysisEnum,
  • issm/trunk-jpl/src/c/shared/Enum/EnumToStringx.cpp

    r16167 r16181  
    161161                case MaterialsRheologyNEnum : return "MaterialsRheologyN";
    162162                case DamageDEnum : return "DamageD";
     163                case QmuDamageDEnum : return "QmuDamageD";
    163164                case DamageDbarEnum : return "DamageDbar";
    164                 case DamageIsdamageEnum : return "DamageIsdamage";
    165165                case DamageLawEnum : return "DamageLaw";
    166166                case DamageC1Enum : return "DamageC1";
     
    168168                case DamageC3Enum : return "DamageC3";
    169169                case DamageC4Enum : return "DamageC4";
     170                case DamageHealingEnum : return "DamageHealing";
    170171                case DamageStressThresholdEnum : return "DamageStressThreshold";
     172                case DamageStabilizationEnum : return "DamageStabilization";
     173                case DamageSpcdamageEnum : return "DamageSpcdamage";
    171174                case MaterialsRhoIceEnum : return "MaterialsRhoIce";
    172175                case MaterialsRhoWaterEnum : return "MaterialsRhoWater";
     
    272275                case TransientIsthermalEnum : return "TransientIsthermal";
    273276                case TransientIsgiaEnum : return "TransientIsgia";
     277                case TransientIsdamageEnum : return "TransientIsdamage";
    274278                case TransientNumRequestedOutputsEnum : return "TransientNumRequestedOutputs";
    275279                case TransientRequestedOutputsEnum : return "TransientRequestedOutputs";
     
    293297                case BedSlopeXAnalysisEnum : return "BedSlopeXAnalysis";
    294298                case BedSlopeYAnalysisEnum : return "BedSlopeYAnalysis";
     299                case DamageEvolutionSolutionEnum : return "DamageEvolutionSolution";
     300                case DamageEvolutionAnalysisEnum : return "DamageEvolutionAnalysis";
    295301                case StressbalanceAnalysisEnum : return "StressbalanceAnalysis";
    296302                case StressbalanceSIAAnalysisEnum : return "StressbalanceSIAAnalysis";
  • issm/trunk-jpl/src/c/shared/Enum/StringToEnumx.cpp

    r16167 r16181  
    164164              else if (strcmp(name,"MaterialsRheologyN")==0) return MaterialsRheologyNEnum;
    165165              else if (strcmp(name,"DamageD")==0) return DamageDEnum;
     166              else if (strcmp(name,"QmuDamageD")==0) return QmuDamageDEnum;
    166167              else if (strcmp(name,"DamageDbar")==0) return DamageDbarEnum;
    167               else if (strcmp(name,"DamageIsdamage")==0) return DamageIsdamageEnum;
    168168              else if (strcmp(name,"DamageLaw")==0) return DamageLawEnum;
    169169              else if (strcmp(name,"DamageC1")==0) return DamageC1Enum;
     
    171171              else if (strcmp(name,"DamageC3")==0) return DamageC3Enum;
    172172              else if (strcmp(name,"DamageC4")==0) return DamageC4Enum;
     173              else if (strcmp(name,"DamageHealing")==0) return DamageHealingEnum;
    173174              else if (strcmp(name,"DamageStressThreshold")==0) return DamageStressThresholdEnum;
     175              else if (strcmp(name,"DamageStabilization")==0) return DamageStabilizationEnum;
     176              else if (strcmp(name,"DamageSpcdamage")==0) return DamageSpcdamageEnum;
    174177              else if (strcmp(name,"MaterialsRhoIce")==0) return MaterialsRhoIceEnum;
    175178              else if (strcmp(name,"MaterialsRhoWater")==0) return MaterialsRhoWaterEnum;
     
    257260              else if (strcmp(name,"ThermalIsdynamicbasalspc")==0) return ThermalIsdynamicbasalspcEnum;
    258261              else if (strcmp(name,"ThermalMaxiter")==0) return ThermalMaxiterEnum;
    259               else if (strcmp(name,"ThermalPenaltyFactor")==0) return ThermalPenaltyFactorEnum;
    260               else if (strcmp(name,"ThermalPenaltyLock")==0) return ThermalPenaltyLockEnum;
    261               else if (strcmp(name,"ThermalPenaltyThreshold")==0) return ThermalPenaltyThresholdEnum;
    262262         else stage=3;
    263263   }
    264264   if(stage==3){
    265               if (strcmp(name,"ThermalSpctemperature")==0) return ThermalSpctemperatureEnum;
     265              if (strcmp(name,"ThermalPenaltyFactor")==0) return ThermalPenaltyFactorEnum;
     266              else if (strcmp(name,"ThermalPenaltyLock")==0) return ThermalPenaltyLockEnum;
     267              else if (strcmp(name,"ThermalPenaltyThreshold")==0) return ThermalPenaltyThresholdEnum;
     268              else if (strcmp(name,"ThermalSpctemperature")==0) return ThermalSpctemperatureEnum;
    266269              else if (strcmp(name,"ThermalStabilization")==0) return ThermalStabilizationEnum;
    267270              else if (strcmp(name,"GiaMantleViscosity")==0) return GiaMantleViscosityEnum;
     
    278281              else if (strcmp(name,"TransientIsthermal")==0) return TransientIsthermalEnum;
    279282              else if (strcmp(name,"TransientIsgia")==0) return TransientIsgiaEnum;
     283              else if (strcmp(name,"TransientIsdamage")==0) return TransientIsdamageEnum;
    280284              else if (strcmp(name,"TransientNumRequestedOutputs")==0) return TransientNumRequestedOutputsEnum;
    281285              else if (strcmp(name,"TransientRequestedOutputs")==0) return TransientRequestedOutputsEnum;
     
    299303              else if (strcmp(name,"BedSlopeXAnalysis")==0) return BedSlopeXAnalysisEnum;
    300304              else if (strcmp(name,"BedSlopeYAnalysis")==0) return BedSlopeYAnalysisEnum;
     305              else if (strcmp(name,"DamageEvolutionSolution")==0) return DamageEvolutionSolutionEnum;
     306              else if (strcmp(name,"DamageEvolutionAnalysis")==0) return DamageEvolutionAnalysisEnum;
    301307              else if (strcmp(name,"StressbalanceAnalysis")==0) return StressbalanceAnalysisEnum;
    302308              else if (strcmp(name,"StressbalanceSIAAnalysis")==0) return StressbalanceSIAAnalysisEnum;
     
    377383              else if (strcmp(name,"SSA3dIceFront")==0) return SSA3dIceFrontEnum;
    378384              else if (strcmp(name,"Matice")==0) return MaticeEnum;
    379               else if (strcmp(name,"Matpar")==0) return MatparEnum;
     385         else stage=4;
     386   }
     387   if(stage==4){
     388              if (strcmp(name,"Matpar")==0) return MatparEnum;
    380389              else if (strcmp(name,"Node")==0) return NodeEnum;
    381390              else if (strcmp(name,"Numericalflux")==0) return NumericalfluxEnum;
     
    383392              else if (strcmp(name,"Param")==0) return ParamEnum;
    384393              else if (strcmp(name,"L1L2IceFront")==0) return L1L2IceFrontEnum;
    385          else stage=4;
    386    }
    387    if(stage==4){
    388               if (strcmp(name,"HOIceFront")==0) return HOIceFrontEnum;
     394              else if (strcmp(name,"HOIceFront")==0) return HOIceFrontEnum;
    389395              else if (strcmp(name,"Pengrid")==0) return PengridEnum;
    390396              else if (strcmp(name,"Penpair")==0) return PenpairEnum;
     
    500506              else if (strcmp(name,"StressTensor")==0) return StressTensorEnum;
    501507              else if (strcmp(name,"StressTensorxx")==0) return StressTensorxxEnum;
    502               else if (strcmp(name,"StressTensorxy")==0) return StressTensorxyEnum;
     508         else stage=5;
     509   }
     510   if(stage==5){
     511              if (strcmp(name,"StressTensorxy")==0) return StressTensorxyEnum;
    503512              else if (strcmp(name,"StressTensorxz")==0) return StressTensorxzEnum;
    504513              else if (strcmp(name,"StressTensoryy")==0) return StressTensoryyEnum;
     
    506515              else if (strcmp(name,"StressTensorzz")==0) return StressTensorzzEnum;
    507516              else if (strcmp(name,"GiaCrossSectionShape")==0) return GiaCrossSectionShapeEnum;
    508          else stage=5;
    509    }
    510    if(stage==5){
    511               if (strcmp(name,"GiadWdt")==0) return GiadWdtEnum;
     517              else if (strcmp(name,"GiadWdt")==0) return GiadWdtEnum;
    512518              else if (strcmp(name,"GiaW")==0) return GiaWEnum;
    513519              else if (strcmp(name,"P0")==0) return P0Enum;
  • issm/trunk-jpl/src/m/classes/damage.m

    r16179 r16181  
    88        % {{{
    99                %damage
    10                 isdamage            = NaN;
    1110                D                   = NaN;
    1211                law                 = '';
     12                spcdamage           = NaN;
     13                stabilization       = NaN;
    1314               
    1415                %parameters for law 'initial':
     
    1819                c3                  = NaN;
    1920                c4                  = NaN;
     21                healing             = NaN;
    2022
    2123        end % }}}
     
    4244
    4345                        %damage parameters:
    44                         obj.isdamage=0;
    4546                        obj.D=0;
    4647                        obj.law='undamaged';
    47                        
     48               
     49                        %Type of stabilization used
     50                        obj.stabilization=2;
     51       
    4852                        %damage parameters for 'initial' law of damage propagation
    4953                        obj.stress_threshold=0;
     54                        obj.healing=0;
    5055                        obj.c1=0;
    5156                        obj.c2=0;
     
    5762                       
    5863                        md = checkfield(md,'damage.D','>=0',0,'size',[md.mesh.numberofvertices 1]);
    59                         if obj.isdamage,
    60                                 md = checkfield(md,'damage.law','values',{'undamaged','pralong'});
    61                         end
     64                        md = checkfield(md,'damage.law','values',{'undamaged','pralong'});
     65                        md = checkfield(md,'damage.spcdamage','forcing',1);
     66                        md = checkfield(md,'damage.stabilization','numel',[1],'values',[0 1 2]);
    6267                        if strcmpi(obj.law,'pralong'),
     68                                md = checkfield(md,'damage.healing','>=',0);
    6369                                md = checkfield(md,'damage.c1','>=',0);
    6470                                md = checkfield(md,'damage.c2','>=',0);
     
    7480                        fielddisplay(obj,'D','damage tensor (scalar)');
    7581                        fielddisplay(obj,'law','damage law (string) from {''undamaged'',''pralong''}');
     82                        fielddisplay(obj,'spcdamage','damage constraints (NaN means no constraint) [m]');
     83                        fielddisplay(obj,'stabilization','0: no, 1: artificial_diffusivity, 2: SUPG');
    7684                        if strcmpi(obj.law,'pralong'),
    7785                                fielddisplay(obj,'c1','damage parameter 1');
     
    7987                                fielddisplay(obj,'c3','damage parameter 3');
    8088                                fielddisplay(obj,'c4','damage parameter 4');
     89                                fielddisplay(obj,'healing','damage healing parameter 1');
    8190                                fielddisplay(obj,'stress_threshold','damage stress threshold [Pa]');
    8291                        end
     
    8594                function marshall(obj,md,fid) % {{{
    8695               
    87                         WriteData(fid,'object',obj,'fieldname','isdamage','format','DoubleMat','mattype',1);
    8896                        WriteData(fid,'object',obj,'fieldname','D','format','DoubleMat','mattype',1);
    89                         if obj.isdamage,
    90                                 WriteData(fid,'object',obj,'fieldname','law','format','String');
    91                         end
     97                        WriteData(fid,'object',obj,'fieldname','law','format','String');
     98                        WriteData(fid,'object',obj,'fieldname','stabilization','format','Integer');
     99                        WriteData(fid,'object',obj,'fieldname','spcdamage','format','DoubleMat','mattype',1,'forcinglength',md.mesh.numberofvertices+1);
    92100                        if strcmpi(obj.law,'pralong'),
    93101                                WriteData(fid,'object',obj,'fieldname','c1','format','Double');
     
    96104                                WriteData(fid,'object',obj,'fieldname','c4','format','Double');
    97105                                WriteData(fid,'object',obj,'fieldname','stress_threshold','format','Double');
     106                                WriteData(fid,'object',obj,'fieldname','healing','format','Double');
    98107                        end
    99108
  • issm/trunk-jpl/src/m/classes/transient.m

    r15771 r16181  
    1111                isgroundingline   = 0;
    1212                isgia             = 0;
     13                isdamage             = 0;
    1314                requested_outputs = NaN;
    1415        end
     
    3031                        obj.isgroundingline=0;
    3132                        obj.isgia=0;
     33                        obj.isdamage=0;
    3234
    3335                end % }}}
     
    4244                        md = checkfield(md,'transient.isgroundingline','numel',[1],'values',[0 1]);
    4345                        md = checkfield(md,'transient.isgia','numel',[1],'values',[0 1]);
     46                        md = checkfield(md,'transient.isdamage','numel',[1],'values',[0 1]);
    4447                        md = checkfield(md,'transient.requested_outputs','size',[NaN 1]);
    4548
     
    4851                        disp(sprintf('   transient solution parameters:'));
    4952
    50                         fielddisplay(obj,'ismasstransport','indicates if a masstransport solution is used in the transient');
    51                         fielddisplay(obj,'isstressbalance','indicates if a stressbalance solution is used in the transient');
    52                         fielddisplay(obj,'isthermal','indicates if a thermal solution is used in the transient');
    53                         fielddisplay(obj,'isgroundingline','indicates if a groundingline migration is used in the transient');
    54                         fielddisplay(obj,'isgia','indicates if a postglacial rebound model is used in the transient');
     53                        fielddisplay(obj,'ismasstransport','indicates whether a masstransport solution is used in the transient');
     54                        fielddisplay(obj,'isstressbalance','indicates whether a stressbalance solution is used in the transient');
     55                        fielddisplay(obj,'isthermal','indicates whether a thermal solution is used in the transient');
     56                        fielddisplay(obj,'isgroundingline','indicates whether a groundingline migration is used in the transient');
     57                        fielddisplay(obj,'isgia','indicates whether a postglacial rebound model is used in the transient');
     58                        fielddisplay(obj,'isdamage','indicates whether damage evolution is used in the transient');
    5559                        fielddisplay(obj,'requested_outputs','list of additional outputs requested');
    5660
     
    6266                        WriteData(fid,'object',obj,'fieldname','isgroundingline','format','Boolean');
    6367                        WriteData(fid,'object',obj,'fieldname','isgia','format','Boolean');
     68                        WriteData(fid,'object',obj,'fieldname','isdamage','format','Boolean');
    6469                        WriteData(fid,'object',obj,'fieldname','requested_outputs','format','DoubleMat','mattype',3);
    6570                end % }}}
  • issm/trunk-jpl/src/m/consistency/ismodelselfconsistent.m

    r16008 r16181  
    100100                numanalyses=3;
    101101                analyses=[BedSlopeAnalysisEnum();SurfaceSlopeAnalysisEnum();HydrologyShreveAnalysisEnum();HydrologyDCInefficientAnalysisEnum();HydrologyDCEfficientAnalysisEnum()];
     102       
     103        case DamageEvolutionSolutionEnum(),
     104                numanalyses=1;
     105                analyses=[DamageEvolutionAnalysisEnum()];
    102106
    103107        otherwise
  • issm/trunk-jpl/src/m/enum/EnumDefinitions.py

    r16167 r16181  
    153153def MaterialsRheologyNEnum(): return StringToEnum("MaterialsRheologyN")[0]
    154154def DamageDEnum(): return StringToEnum("DamageD")[0]
     155def QmuDamageDEnum(): return StringToEnum("QmuDamageD")[0]
    155156def DamageDbarEnum(): return StringToEnum("DamageDbar")[0]
    156 def DamageIsdamageEnum(): return StringToEnum("DamageIsdamage")[0]
    157157def DamageLawEnum(): return StringToEnum("DamageLaw")[0]
    158158def DamageC1Enum(): return StringToEnum("DamageC1")[0]
     
    160160def DamageC3Enum(): return StringToEnum("DamageC3")[0]
    161161def DamageC4Enum(): return StringToEnum("DamageC4")[0]
     162def DamageHealingEnum(): return StringToEnum("DamageHealing")[0]
    162163def DamageStressThresholdEnum(): return StringToEnum("DamageStressThreshold")[0]
     164def DamageStabilizationEnum(): return StringToEnum("DamageStabilization")[0]
     165def DamageSpcdamageEnum(): return StringToEnum("DamageSpcdamage")[0]
    163166def MaterialsRhoIceEnum(): return StringToEnum("MaterialsRhoIce")[0]
    164167def MaterialsRhoWaterEnum(): return StringToEnum("MaterialsRhoWater")[0]
     
    264267def TransientIsthermalEnum(): return StringToEnum("TransientIsthermal")[0]
    265268def TransientIsgiaEnum(): return StringToEnum("TransientIsgia")[0]
     269def TransientIsdamageEnum(): return StringToEnum("TransientIsdamage")[0]
    266270def TransientNumRequestedOutputsEnum(): return StringToEnum("TransientNumRequestedOutputs")[0]
    267271def TransientRequestedOutputsEnum(): return StringToEnum("TransientRequestedOutputs")[0]
     
    285289def BedSlopeXAnalysisEnum(): return StringToEnum("BedSlopeXAnalysis")[0]
    286290def BedSlopeYAnalysisEnum(): return StringToEnum("BedSlopeYAnalysis")[0]
     291def DamageEvolutionSolutionEnum(): return StringToEnum("DamageEvolutionSolution")[0]
     292def DamageEvolutionAnalysisEnum(): return StringToEnum("DamageEvolutionAnalysis")[0]
    287293def StressbalanceAnalysisEnum(): return StringToEnum("StressbalanceAnalysis")[0]
    288294def StressbalanceSIAAnalysisEnum(): return StringToEnum("StressbalanceSIAAnalysis")[0]
  • issm/trunk-jpl/src/m/solve/process_solve_options.m

    r16008 r16181  
    1313                BalancethicknessSolutionEnum(),BalancethicknessSoftSolutionEnum(),...
    1414                BalancevelocitySolutionEnum(),BedSlopeSolutionEnum(),...
    15                 SurfaceSlopeSolutionEnum(),HydrologySolutionEnum(),FlaimSolutionEnum(),GiaSolutionEnum()]),
     15                SurfaceSlopeSolutionEnum(),HydrologySolutionEnum(),FlaimSolutionEnum(),GiaSolutionEnum(),DamageEvolutionSolutionEnum()]),
    1616        error(['process_solve_options error message: solution_type ' EnumToString(solution_type) ' not supported yet!']);
    1717end
  • issm/trunk-jpl/test/NightlyRun/IdToName.m

    r16044 r16181  
    6666        case 237, name='SquareShelfTranIspddIsdeltaSSA3d';
    6767        case 270, name='SquareShelfStressSSA2dDamage';
     68        case 271, name='SquareShelfStressSSA2dDamageTran';
    6869        case 272, name='SquareShelfCMZSSA2dDamage';
    6970        case 274, name='SquareShelfStressSSA2dDamageRift';
Note: See TracChangeset for help on using the changeset viewer.