Changeset 19019


Ignore:
Timestamp:
01/20/15 13:21:25 (10 years ago)
Author:
Mathieu Morlighem
Message:

CHG: added some enums

Location:
issm/trunk-jpl/src
Files:
2 added
7 edited

Legend:

Unmodified
Added
Removed

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

    r19002 r19019  
    183183        IssmDouble* eps_xz = xNew<IssmDouble>(numvertices);
    184184        IssmDouble* eps_yz = xNew<IssmDouble>(numvertices);
     185        IssmDouble* eps_ef = xNew<IssmDouble>(numvertices);
    185186
    186187        /* Start looping on the number of vertices: */
    … …  
    200201                        eps_yy[iv]=epsilon[1];
    201202                        eps_xy[iv]=epsilon[2];
     203                        /* eps_eff^2 = 1/2 ( exx^2 + eyy^2 + 2*exy^2 )*/
     204                        eps_ef[iv] = 1./sqrt(2.)*sqrt(epsilon[0]*epsilon[0] + epsilon[1]*epsilon[1] + 2.*epsilon[2]*epsilon[2]);
    202205                }
    203206                else{
    … …  
    209212                        eps_xz[iv]=epsilon[4];
    210213                        eps_yz[iv]=epsilon[5];
     214                        /* eps_eff^2 = exx^2 + eyy^2 + exy^2 + exz^2 + eyz^2 + exx*eyy */
     215                        eps_ef[iv] = sqrt(epsilon[0]*epsilon[0] + epsilon[1]*epsilon[1] + epsilon[3]*epsilon[3] +  epsilon[4]*epsilon[4] + epsilon[5]*epsilon[5] + epsilon[0]*epsilon[1]);
    211216                }
    212217        }
    … …  
    219224        this->AddInput(StrainRateyzEnum,eps_yz,P1Enum);
    220225        this->AddInput(StrainRatezzEnum,eps_zz,P1Enum);
     226        this->AddInput(StrainRateeffectiveEnum,eps_ef,P1Enum);
    221227
    222228        /*Clean up and return*/
    … …  
    13411347                        case StrainRateyzEnum:
    13421348                        case StrainRatezzEnum:
     1349                        case StrainRateeffectiveEnum:
    13431350                                this->ComputeStrainRate();
    13441351                                input=this->inputs->GetInput(output_enum);

  • issm/trunk-jpl/src/c/classes/FemModel.cpp

    r19009 r19019  
    17291729void FemModel::CalvingRateLevermannx(){/*{{{*/
    17301730
    1731            for(int i=0;i<elements->Size();i++){
    1732                               Element* element=dynamic_cast<Element*>(this->elements->GetObjectByOffset(i));
    1733                                               element->CalvingRateLevermann();
    1734                                                            }
     1731        for(int i=0;i<elements->Size();i++){
     1732                Element* element=dynamic_cast<Element*>(this->elements->GetObjectByOffset(i));
     1733                element->CalvingRateLevermann();
     1734        }
    17351735}
    17361736/*}}}*/
    17371737void FemModel::CalvingRatePix(){/*{{{*/
    17381738
    1739            for(int i=0;i<elements->Size();i++){
    1740                               Element* element=dynamic_cast<Element*>(this->elements->GetObjectByOffset(i));
    1741                                               element->CalvingRatePi();
    1742                                                            }
     1739        for(int i=0;i<elements->Size();i++){
     1740                Element* element=dynamic_cast<Element*>(this->elements->GetObjectByOffset(i));
     1741                element->CalvingRatePi();
     1742        }
    17431743}
    17441744/*}}}*/
    17451745void FemModel::StrainRateparallelx(){/*{{{*/
    17461746
    1747            for(int i=0;i<elements->Size();i++){
    1748                               Element* element=dynamic_cast<Element*>(this->elements->GetObjectByOffset(i));
    1749                                               element->StrainRateparallel();
    1750                                                            }
     1747        for(int i=0;i<elements->Size();i++){
     1748                Element* element=dynamic_cast<Element*>(this->elements->GetObjectByOffset(i));
     1749                element->StrainRateparallel();
     1750        }
    17511751}
    17521752/*}}}*/
    17531753void FemModel::StrainRateperpendicularx(){/*{{{*/
    17541754
    1755            for(int i=0;i<elements->Size();i++){
    1756                               Element* element=dynamic_cast<Element*>(this->elements->GetObjectByOffset(i));
    1757                                               element->StrainRateperpendicular();
    1758                                                            }
     1755        for(int i=0;i<elements->Size();i++){
     1756                Element* element=dynamic_cast<Element*>(this->elements->GetObjectByOffset(i));
     1757                element->StrainRateperpendicular();
     1758        }
     1759}
     1760/*}}}*/
     1761void FemModel::DeviatoricStressx(){/*{{{*/
     1762
     1763        for(int i=0;i<elements->Size();i++){
     1764                Element* element=dynamic_cast<Element*>(this->elements->GetObjectByOffset(i));
     1765                element->ComputeDeviatoricStressTensor();
     1766        }
    17591767}
    17601768/*}}}*/

  • issm/trunk-jpl/src/c/classes/FemModel.h

    r19009 r19019  
    8585                void StrainRateparallelx();
    8686                void StrainRateperpendicularx();
     87                void DeviatoricStressx();
    8788                void CalvingRateLevermannx();
    8889                void CalvingRatePix();

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

    r18999 r19019  
    217217        CalvingLevermannEnum,
    218218        CalvingPiEnum,
     219        CalvingDevEnum,
    219220        DefaultCalvingEnum,
    220221        CalvingRequestedOutputsEnum,
    … …  
    229230        StrainRateparallelEnum,
    230231        StrainRateperpendicularEnum,
     232        StrainRateeffectiveEnum,
    231233        MaterialsRhoIceEnum,
    232234        MaterialsRhoSeawaterEnum,

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

    r18999 r19019  
    225225                case CalvingLevermannEnum : return "CalvingLevermann";
    226226                case CalvingPiEnum : return "CalvingPi";
     227                case CalvingDevEnum : return "CalvingDev";
    227228                case DefaultCalvingEnum : return "DefaultCalving";
    228229                case CalvingRequestedOutputsEnum : return "CalvingRequestedOutputs";
    … …  
    237238                case StrainRateparallelEnum : return "StrainRateparallel";
    238239                case StrainRateperpendicularEnum : return "StrainRateperpendicular";
     240                case StrainRateeffectiveEnum : return "StrainRateeffective";
    239241                case MaterialsRhoIceEnum : return "MaterialsRhoIce";
    240242                case MaterialsRhoSeawaterEnum : return "MaterialsRhoSeawater";

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

    r18999 r19019  
    228228              else if (strcmp(name,"CalvingLevermann")==0) return CalvingLevermannEnum;
    229229              else if (strcmp(name,"CalvingPi")==0) return CalvingPiEnum;
     230              else if (strcmp(name,"CalvingDev")==0) return CalvingDevEnum;
    230231              else if (strcmp(name,"DefaultCalving")==0) return DefaultCalvingEnum;
    231232              else if (strcmp(name,"CalvingRequestedOutputs")==0) return CalvingRequestedOutputsEnum;
    … …  
    240241              else if (strcmp(name,"StrainRateparallel")==0) return StrainRateparallelEnum;
    241242              else if (strcmp(name,"StrainRateperpendicular")==0) return StrainRateperpendicularEnum;
     243              else if (strcmp(name,"StrainRateeffective")==0) return StrainRateeffectiveEnum;
    242244              else if (strcmp(name,"MaterialsRhoIce")==0) return MaterialsRhoIceEnum;
    243245              else if (strcmp(name,"MaterialsRhoSeawater")==0) return MaterialsRhoSeawaterEnum;
    … …  
    258260              else if (strcmp(name,"MeshNumberofelements")==0) return MeshNumberofelementsEnum;
    259261              else if (strcmp(name,"MeshNumberoflayers")==0) return MeshNumberoflayersEnum;
    260               else if (strcmp(name,"MeshNumberofvertices2d")==0) return MeshNumberofvertices2dEnum;
    261               else if (strcmp(name,"MeshNumberofvertices")==0) return MeshNumberofverticesEnum;
    262262         else stage=3;
    263263   }
    264264   if(stage==3){
    265               if (strcmp(name,"MeshUpperelements")==0) return MeshUpperelementsEnum;
     265              if (strcmp(name,"MeshNumberofvertices2d")==0) return MeshNumberofvertices2dEnum;
     266              else if (strcmp(name,"MeshNumberofvertices")==0) return MeshNumberofverticesEnum;
     267              else if (strcmp(name,"MeshUpperelements")==0) return MeshUpperelementsEnum;
    266268              else if (strcmp(name,"MeshVertexonbase")==0) return MeshVertexonbaseEnum;
    267269              else if (strcmp(name,"MeshVertexonsurface")==0) return MeshVertexonsurfaceEnum;
    … …  
    381383              else if (strcmp(name,"SurfaceforcingsAccumulation")==0) return SurfaceforcingsAccumulationEnum;
    382384              else if (strcmp(name,"SurfaceforcingsEvaporation")==0) return SurfaceforcingsEvaporationEnum;
    383               else if (strcmp(name,"SurfaceforcingsRunoff")==0) return SurfaceforcingsRunoffEnum;
    384               else if (strcmp(name,"SMBmeltcomponents")==0) return SMBmeltcomponentsEnum;
    385385         else stage=4;
    386386   }
    387387   if(stage==4){
    388               if (strcmp(name,"SurfaceforcingsMelt")==0) return SurfaceforcingsMeltEnum;
     388              if (strcmp(name,"SurfaceforcingsRunoff")==0) return SurfaceforcingsRunoffEnum;
     389              else if (strcmp(name,"SMBmeltcomponents")==0) return SMBmeltcomponentsEnum;
     390              else if (strcmp(name,"SurfaceforcingsMelt")==0) return SurfaceforcingsMeltEnum;
    389391              else if (strcmp(name,"SurfaceforcingsRefreeze")==0) return SurfaceforcingsRefreezeEnum;
    390392              else if (strcmp(name,"SurfaceforcingsIspdd")==0) return SurfaceforcingsIspddEnum;
    … …  
    504506              else if (strcmp(name,"Penpair")==0) return PenpairEnum;
    505507              else if (strcmp(name,"Profiler")==0) return ProfilerEnum;
    506               else if (strcmp(name,"MatrixParam")==0) return MatrixParamEnum;
    507               else if (strcmp(name,"Masscon")==0) return MassconEnum;
    508508         else stage=5;
    509509   }
    510510   if(stage==5){
    511               if (strcmp(name,"MassconName")==0) return MassconNameEnum;
     511              if (strcmp(name,"MatrixParam")==0) return MatrixParamEnum;
     512              else if (strcmp(name,"Masscon")==0) return MassconEnum;
     513              else if (strcmp(name,"MassconName")==0) return MassconNameEnum;
    512514              else if (strcmp(name,"MassconDefinitionenum")==0) return MassconDefinitionenumEnum;
    513515              else if (strcmp(name,"MassconLevelset")==0) return MassconLevelsetEnum;
    … …  
    627629              else if (strcmp(name,"DeviatoricStressxx")==0) return DeviatoricStressxxEnum;
    628630              else if (strcmp(name,"DeviatoricStressxy")==0) return DeviatoricStressxyEnum;
    629               else if (strcmp(name,"DeviatoricStressxz")==0) return DeviatoricStressxzEnum;
    630               else if (strcmp(name,"DeviatoricStressyy")==0) return DeviatoricStressyyEnum;
    631631         else stage=6;
    632632   }
    633633   if(stage==6){
    634               if (strcmp(name,"DeviatoricStressyz")==0) return DeviatoricStressyzEnum;
     634              if (strcmp(name,"DeviatoricStressxz")==0) return DeviatoricStressxzEnum;
     635              else if (strcmp(name,"DeviatoricStressyy")==0) return DeviatoricStressyyEnum;
     636              else if (strcmp(name,"DeviatoricStressyz")==0) return DeviatoricStressyzEnum;
    635637              else if (strcmp(name,"DeviatoricStresszz")==0) return DeviatoricStresszzEnum;
    636638              else if (strcmp(name,"StrainRate")==0) return StrainRateEnum;
    … …  
    750752              else if (strcmp(name,"Outputdefinition70")==0) return Outputdefinition70Enum;
    751753              else if (strcmp(name,"Outputdefinition71")==0) return Outputdefinition71Enum;
    752               else if (strcmp(name,"Outputdefinition72")==0) return Outputdefinition72Enum;
    753               else if (strcmp(name,"Outputdefinition73")==0) return Outputdefinition73Enum;
    754754         else stage=7;
    755755   }
    756756   if(stage==7){
    757               if (strcmp(name,"Outputdefinition74")==0) return Outputdefinition74Enum;
     757              if (strcmp(name,"Outputdefinition72")==0) return Outputdefinition72Enum;
     758              else if (strcmp(name,"Outputdefinition73")==0) return Outputdefinition73Enum;
     759              else if (strcmp(name,"Outputdefinition74")==0) return Outputdefinition74Enum;
    758760              else if (strcmp(name,"Outputdefinition75")==0) return Outputdefinition75Enum;
    759761              else if (strcmp(name,"Outputdefinition76")==0) return Outputdefinition76Enum;
    … …  
    873875              else if (strcmp(name,"Dense")==0) return DenseEnum;
    874876              else if (strcmp(name,"MpiDense")==0) return MpiDenseEnum;
    875               else if (strcmp(name,"MpiSparse")==0) return MpiSparseEnum;
    876               else if (strcmp(name,"Seq")==0) return SeqEnum;
    877877         else stage=8;
    878878   }
    879879   if(stage==8){
    880               if (strcmp(name,"Mpi")==0) return MpiEnum;
     880              if (strcmp(name,"MpiSparse")==0) return MpiSparseEnum;
     881              else if (strcmp(name,"Seq")==0) return SeqEnum;
     882              else if (strcmp(name,"Mpi")==0) return MpiEnum;
    881883              else if (strcmp(name,"Mumps")==0) return MumpsEnum;
    882884              else if (strcmp(name,"Gsl")==0) return GslEnum;

  • issm/trunk-jpl/src/m/enum/EnumDefinitions.py

    r18998 r19019  
    217217def CalvingLevermannEnum(): return StringToEnum("CalvingLevermann")[0]
    218218def CalvingPiEnum(): return StringToEnum("CalvingPi")[0]
     219def CalvingDevEnum(): return StringToEnum("CalvingDev")[0]
    219220def DefaultCalvingEnum(): return StringToEnum("DefaultCalving")[0]
    220221def CalvingRequestedOutputsEnum(): return StringToEnum("CalvingRequestedOutputs")[0]
    … …  
    229230def StrainRateparallelEnum(): return StringToEnum("StrainRateparallel")[0]
    230231def StrainRateperpendicularEnum(): return StringToEnum("StrainRateperpendicular")[0]
     232def StrainRateeffectiveEnum(): return StringToEnum("StrainRateeffective")[0]
    231233def MaterialsRhoIceEnum(): return StringToEnum("MaterialsRhoIce")[0]
    232234def MaterialsRhoSeawaterEnum(): return StringToEnum("MaterialsRhoSeawater")[0]
Note: See TracChangeset for help on using the changeset viewer.