Changeset 19087
- Timestamp:
- 02/04/15 09:13:18 (10 years ago)
- Location:
- issm/trunk-jpl/src
- Files:
-
- 51 deleted
- 15 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/Makefile.am
r19047 r19087 55 55 ./classes/Materials/Materials.cpp\ 56 56 ./classes/Materials/Matice.cpp\ 57 ./classes/Materials/Matseaice.cpp\58 57 ./classes/Materials/Matpar.cpp\ 59 58 ./classes/Constraints/Constraints.cpp\ … … 410 409 issm_sources += ./analyses/DepthAverageAnalysis.cpp 411 410 endif 412 if SEAICE413 issm_sources += ./cores/seaice_core.cpp414 issm_sources += ./analyses/SeaiceAnalysis.cpp415 endif416 411 if THERMAL 417 412 issm_sources += ./analyses/ThermalAnalysis.cpp -
issm/trunk-jpl/src/c/analyses/EnumToAnalysis.cpp
r18593 r19087 98 98 case GiaAnalysisEnum : return new GiaAnalysis(); 99 99 #endif 100 #ifdef _HAVE_SEAICE_101 case SeaiceAnalysisEnum : return new SeaiceAnalysis();102 #endif103 100 #ifdef _HAVE_MESHDEFORMATION_ 104 101 case MeshdeformationAnalysisEnum : return new MeshdeformationAnalysis(); -
issm/trunk-jpl/src/c/analyses/analyses.h
r18593 r19087 33 33 #include "./MeshdeformationAnalysis.h" 34 34 #include "./SmoothAnalysis.h" 35 #include "./SeaiceAnalysis.h"36 35 #include "./StressbalanceAnalysis.h" 37 36 #include "./StressbalanceSIAAnalysis.h" -
issm/trunk-jpl/src/c/classes/FemModel.cpp
r19062 r19087 586 586 analyses_temp[numanalyses++]=L2ProjectionBaseAnalysisEnum; 587 587 } 588 break;589 590 case SeaiceSolutionEnum:591 analyses_temp[numanalyses++]=SeaiceAnalysisEnum;592 588 break; 593 589 -
issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp
r18946 r19087 99 99 iomodel->Constant(&this->mantle_shear_modulus,MaterialsMantleShearModulusEnum); 100 100 iomodel->Constant(&this->mantle_density,MaterialsMantleDensityEnum); 101 break;102 case MatseaiceEnum:103 iomodel->Constant(&this->rho_ice,MaterialsRhoIceEnum);104 iomodel->Constant(&this->poisson,MaterialsPoissonEnum);105 iomodel->Constant(&this->young_modulus,MaterialsYoungModulusEnum);106 iomodel->Constant(&this->ridging_exponent,MaterialsRidgingExponentEnum);107 iomodel->Constant(&this->cohesion,MaterialsCohesionEnum);108 iomodel->Constant(&this->internal_friction_coef,MaterialsInternalFrictionCoefEnum);109 iomodel->Constant(&this->compression_coef,MaterialsCompressionCoefEnum);110 iomodel->Constant(&this->traction_coef,MaterialsTractionCoefEnum);111 iomodel->Constant(&this->g,ConstantsGEnum);112 iomodel->Constant(&this->omega,ConstantsOmegaEnum);113 iomodel->Constant(&this->time_relaxation_stress,MaterialsTimeRelaxationStressEnum);114 iomodel->Constant(&this->time_relaxation_damage,MaterialsTimeRelaxationDamageEnum);115 101 break; 116 102 default: … … 409 395 case HydrologydcWaterCompressibilityEnum: return this->water_compressibility; 410 396 case ConstantsGEnum: return this->g; 411 case MaterialsPoissonEnum: return this->poisson;412 case MaterialsYoungModulusEnum: return this->young_modulus;413 case MaterialsRidgingExponentEnum: return this->ridging_exponent;414 case MaterialsCohesionEnum: return this->cohesion;415 case MaterialsInternalFrictionCoefEnum: return this->internal_friction_coef;416 case MaterialsCompressionCoefEnum: return this->compression_coef;417 case MaterialsTractionCoefEnum: return this->traction_coef;418 397 case SurfaceforcingsDesfacEnum: return this->desfac; 419 398 case SurfaceforcingsS0pEnum: return this->s0p; … … 425 404 case MaterialsMantleDensityEnum: return this->mantle_density; 426 405 case MaterialsMantleShearModulusEnum: return this->mantle_shear_modulus; 427 case ConstantsOmegaEnum: return this->omega;428 case MaterialsTimeRelaxationStressEnum: return this->time_relaxation_stress;429 case MaterialsTimeRelaxationDamageEnum: return this->time_relaxation_damage;430 406 default: _error_("Enum "<<EnumToStringx(enum_in)<<" not supported yet"); 431 407 } -
issm/trunk-jpl/src/c/classes/classes.h
r18812 r19087 80 80 #include "./Materials/Material.h" 81 81 #include "./Materials/Matice.h" 82 #include "./Materials/Matseaice.h"83 82 #include "./Materials/Matpar.h" 84 83 -
issm/trunk-jpl/src/c/cores/CorePointerFromSolutionEnum.cpp
r18492 r19087 69 69 solutioncore=&damage_core; 70 70 break; 71 case SeaiceSolutionEnum:72 solutioncore=&seaice_core;73 break;74 71 default: 75 72 _error_("solution type: " << EnumToStringx(solutiontype) << " not supported yet!"); -
issm/trunk-jpl/src/c/cores/cores.h
r18896 r19087 47 47 void gia_core(FemModel* femmodel); 48 48 void damage_core(FemModel* femmodel); 49 void seaice_core(FemModel* femmodel);50 49 IssmDouble objectivefunction(IssmDouble search_scalar,FemModel* femmodel); 51 50 -
issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp
r18504 r19087 79 79 } 80 80 break; 81 case MatseaiceEnum:82 iomodel->FetchDataToInput(elements,DamageEnum);83 for (i=0;i<iomodel->numberofelements;i++) if(iomodel->my_elements[i]) materials->AddObject(new Matseaice(i+1,i,iomodel));84 break;85 81 default: 86 82 _error_("Materials "<<EnumToStringx(materials_type)<<" not supported"); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp
r18968 r19087 68 68 parameters->AddObject(iomodel->CopyConstantObject(InversionTypeEnum)); 69 69 parameters->AddObject(iomodel->CopyConstantObject(CalvingLawEnum)); 70 if(solution_type==SeaiceSolutionEnum){ 71 } 72 else{ 70 {/*This is specific to ice...*/ 73 71 parameters->AddObject(iomodel->CopyConstantObject(MeshElementtypeEnum)); 74 72 parameters->AddObject(iomodel->CopyConstantObject(SteadystateReltolEnum)); … … 160 158 /*Nothing to add to parameters*/ 161 159 break; 162 case SeaiceatmEnum:163 parameters->AddObject(iomodel->CopyConstantObject(SurfaceforcingsRhoAirEnum));164 parameters->AddObject(iomodel->CopyConstantObject(SurfaceforcingsAirCoefEnum));165 parameters->AddObject(iomodel->CopyConstantObject(SurfaceforcingsAirLinDragCoefEnum));166 parameters->AddObject(iomodel->CopyConstantObject(SurfaceforcingsAirQuadDragCoefEnum));167 break;168 160 default: 169 161 _error_("Surface mass balance model "<<EnumToStringx(smb_model)<<" not supported yet"); … … 181 173 parameters->AddObject(iomodel->CopyConstantObject(BasalforcingsDeepwaterElevationEnum)); 182 174 parameters->AddObject(iomodel->CopyConstantObject(BasalforcingsUpperwaterElevationEnum)); 183 break;184 case SeaiceoceanEnum:185 parameters->AddObject(iomodel->CopyConstantObject(BasalforcingsRhoOceanEnum));186 parameters->AddObject(iomodel->CopyConstantObject(BasalforcingsOceanCoefEnum));187 parameters->AddObject(iomodel->CopyConstantObject(BasalforcingsOceanLinDragCoefEnum));188 parameters->AddObject(iomodel->CopyConstantObject(BasalforcingsOceanQuadDragCoefEnum));189 parameters->AddObject(iomodel->CopyConstantObject(BasalforcingsOceanTurningAngleEnum));190 175 break; 191 176 default: -
issm/trunk-jpl/src/c/shared/Enum/EnumDefinitions.h
r19066 r19087 430 430 GiaSolutionEnum, 431 431 GiaAnalysisEnum, 432 SeaiceSolutionEnum,433 SeaiceAnalysisEnum,434 432 MeshdeformationSolutionEnum, 435 433 MeshdeformationAnalysisEnum, … … 918 916 LevelsetfunctionPicardEnum, 919 917 /*}}}*/ 920 /*Sea Ice {{{*/921 SeaiceatmEnum,922 SeaiceoceanEnum,923 SeaiceThicknessEnum,924 SeaiceConcentrationEnum,925 SeaiceMinConcentrationEnum,926 SeaiceMinThicknessEnum,927 SeaiceMaxThicknessEnum,928 SeaiceSpcvxEnum,929 SeaiceSpcvyEnum,930 SeaiceCoriolisFactorEnum,931 BasalforcingsRhoOceanEnum,932 BasalforcingsOceanCoefEnum,933 BasalforcingsOceanLinDragCoefEnum,934 BasalforcingsOceanQuadDragCoefEnum,935 BasalforcingsOceanTurningAngleEnum,936 BasalforcingsOceanSshEnum,937 BasalforcingsOceanVxEnum,938 BasalforcingsOceanVyEnum,939 SurfaceforcingsRhoAirEnum,940 SurfaceforcingsAirCoefEnum,941 SurfaceforcingsAirLinDragCoefEnum,942 SurfaceforcingsAirQuadDragCoefEnum,943 SurfaceforcingsWindVxEnum,944 SurfaceforcingsWindVyEnum,945 MatseaiceEnum,946 MaterialsPoissonEnum,947 MaterialsYoungModulusEnum,948 MaterialsTimeRelaxationStressEnum,949 MaterialsTimeRelaxationDamageEnum,950 MaterialsRidgingExponentEnum,951 MaterialsCohesionEnum,952 MaterialsInternalFrictionCoefEnum,953 MaterialsCompressionCoefEnum,954 MaterialsTractionCoefEnum,955 VxStarEnum,956 VyStarEnum,957 StressTensorPredictorxxEnum,958 StressTensorPredictoryyEnum,959 StressTensorPredictorxyEnum,960 /*}}}*/961 918 MaximumNumberOfDefinitionsEnum 962 919 }; -
issm/trunk-jpl/src/c/shared/Enum/EnumToStringx.cpp
r19066 r19087 433 433 case GiaSolutionEnum : return "GiaSolution"; 434 434 case GiaAnalysisEnum : return "GiaAnalysis"; 435 case SeaiceSolutionEnum : return "SeaiceSolution";436 case SeaiceAnalysisEnum : return "SeaiceAnalysis";437 435 case MeshdeformationSolutionEnum : return "MeshdeformationSolution"; 438 436 case MeshdeformationAnalysisEnum : return "MeshdeformationAnalysis"; … … 878 876 case LevelsetfunctionSlopeYEnum : return "LevelsetfunctionSlopeY"; 879 877 case LevelsetfunctionPicardEnum : return "LevelsetfunctionPicard"; 880 case SeaiceatmEnum : return "Seaiceatm";881 case SeaiceoceanEnum : return "Seaiceocean";882 case SeaiceThicknessEnum : return "SeaiceThickness";883 case SeaiceConcentrationEnum : return "SeaiceConcentration";884 case SeaiceMinConcentrationEnum : return "SeaiceMinConcentration";885 case SeaiceMinThicknessEnum : return "SeaiceMinThickness";886 case SeaiceMaxThicknessEnum : return "SeaiceMaxThickness";887 case SeaiceSpcvxEnum : return "SeaiceSpcvx";888 case SeaiceSpcvyEnum : return "SeaiceSpcvy";889 case SeaiceCoriolisFactorEnum : return "SeaiceCoriolisFactor";890 case BasalforcingsRhoOceanEnum : return "BasalforcingsRhoOcean";891 case BasalforcingsOceanCoefEnum : return "BasalforcingsOceanCoef";892 case BasalforcingsOceanLinDragCoefEnum : return "BasalforcingsOceanLinDragCoef";893 case BasalforcingsOceanQuadDragCoefEnum : return "BasalforcingsOceanQuadDragCoef";894 case BasalforcingsOceanTurningAngleEnum : return "BasalforcingsOceanTurningAngle";895 case BasalforcingsOceanSshEnum : return "BasalforcingsOceanSsh";896 case BasalforcingsOceanVxEnum : return "BasalforcingsOceanVx";897 case BasalforcingsOceanVyEnum : return "BasalforcingsOceanVy";898 case SurfaceforcingsRhoAirEnum : return "SurfaceforcingsRhoAir";899 case SurfaceforcingsAirCoefEnum : return "SurfaceforcingsAirCoef";900 case SurfaceforcingsAirLinDragCoefEnum : return "SurfaceforcingsAirLinDragCoef";901 case SurfaceforcingsAirQuadDragCoefEnum : return "SurfaceforcingsAirQuadDragCoef";902 case SurfaceforcingsWindVxEnum : return "SurfaceforcingsWindVx";903 case SurfaceforcingsWindVyEnum : return "SurfaceforcingsWindVy";904 case MatseaiceEnum : return "Matseaice";905 case MaterialsPoissonEnum : return "MaterialsPoisson";906 case MaterialsYoungModulusEnum : return "MaterialsYoungModulus";907 case MaterialsTimeRelaxationStressEnum : return "MaterialsTimeRelaxationStress";908 case MaterialsTimeRelaxationDamageEnum : return "MaterialsTimeRelaxationDamage";909 case MaterialsRidgingExponentEnum : return "MaterialsRidgingExponent";910 case MaterialsCohesionEnum : return "MaterialsCohesion";911 case MaterialsInternalFrictionCoefEnum : return "MaterialsInternalFrictionCoef";912 case MaterialsCompressionCoefEnum : return "MaterialsCompressionCoef";913 case MaterialsTractionCoefEnum : return "MaterialsTractionCoef";914 case VxStarEnum : return "VxStar";915 case VyStarEnum : return "VyStar";916 case StressTensorPredictorxxEnum : return "StressTensorPredictorxx";917 case StressTensorPredictoryyEnum : return "StressTensorPredictoryy";918 case StressTensorPredictorxyEnum : return "StressTensorPredictorxy";919 878 case MaximumNumberOfDefinitionsEnum : return "MaximumNumberOfDefinitions"; 920 879 default : return "unknown"; -
issm/trunk-jpl/src/c/shared/Enum/StringToEnumx.cpp
r19066 r19087 442 442 else if (strcmp(name,"GiaSolution")==0) return GiaSolutionEnum; 443 443 else if (strcmp(name,"GiaAnalysis")==0) return GiaAnalysisEnum; 444 else if (strcmp(name,"SeaiceSolution")==0) return SeaiceSolutionEnum;445 else if (strcmp(name,"SeaiceAnalysis")==0) return SeaiceAnalysisEnum;446 444 else if (strcmp(name,"MeshdeformationSolution")==0) return MeshdeformationSolutionEnum; 447 445 else if (strcmp(name,"MeshdeformationAnalysis")==0) return MeshdeformationAnalysisEnum; … … 506 504 else if (strcmp(name,"Penpair")==0) return PenpairEnum; 507 505 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; 508 508 else stage=5; 509 509 } 510 510 if(stage==5){ 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; 511 if (strcmp(name,"MassconName")==0) return MassconNameEnum; 514 512 else if (strcmp(name,"MassconDefinitionenum")==0) return MassconDefinitionenumEnum; 515 513 else if (strcmp(name,"MassconLevelset")==0) return MassconLevelsetEnum; … … 629 627 else if (strcmp(name,"DeviatoricStressxx")==0) return DeviatoricStressxxEnum; 630 628 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; 631 631 else stage=6; 632 632 } 633 633 if(stage==6){ 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; 634 if (strcmp(name,"DeviatoricStressyz")==0) return DeviatoricStressyzEnum; 637 635 else if (strcmp(name,"DeviatoricStresszz")==0) return DeviatoricStresszzEnum; 638 636 else if (strcmp(name,"StrainRate")==0) return StrainRateEnum; … … 752 750 else if (strcmp(name,"Outputdefinition70")==0) return Outputdefinition70Enum; 753 751 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; 754 754 else stage=7; 755 755 } 756 756 if(stage==7){ 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; 757 if (strcmp(name,"Outputdefinition74")==0) return Outputdefinition74Enum; 760 758 else if (strcmp(name,"Outputdefinition75")==0) return Outputdefinition75Enum; 761 759 else if (strcmp(name,"Outputdefinition76")==0) return Outputdefinition76Enum; … … 875 873 else if (strcmp(name,"XY")==0) return XYEnum; 876 874 else if (strcmp(name,"XYZ")==0) return XYZEnum; 875 else if (strcmp(name,"Dense")==0) return DenseEnum; 876 else if (strcmp(name,"MpiDense")==0) return MpiDenseEnum; 877 877 else stage=8; 878 878 } 879 879 if(stage==8){ 880 if (strcmp(name,"Dense")==0) return DenseEnum; 881 else if (strcmp(name,"MpiDense")==0) return MpiDenseEnum; 882 else if (strcmp(name,"MpiSparse")==0) return MpiSparseEnum; 880 if (strcmp(name,"MpiSparse")==0) return MpiSparseEnum; 883 881 else if (strcmp(name,"Seq")==0) return SeqEnum; 884 882 else if (strcmp(name,"Mpi")==0) return MpiEnum; … … 899 897 else if (strcmp(name,"LevelsetfunctionSlopeY")==0) return LevelsetfunctionSlopeYEnum; 900 898 else if (strcmp(name,"LevelsetfunctionPicard")==0) return LevelsetfunctionPicardEnum; 901 else if (strcmp(name,"Seaiceatm")==0) return SeaiceatmEnum;902 else if (strcmp(name,"Seaiceocean")==0) return SeaiceoceanEnum;903 else if (strcmp(name,"SeaiceThickness")==0) return SeaiceThicknessEnum;904 else if (strcmp(name,"SeaiceConcentration")==0) return SeaiceConcentrationEnum;905 else if (strcmp(name,"SeaiceMinConcentration")==0) return SeaiceMinConcentrationEnum;906 else if (strcmp(name,"SeaiceMinThickness")==0) return SeaiceMinThicknessEnum;907 else if (strcmp(name,"SeaiceMaxThickness")==0) return SeaiceMaxThicknessEnum;908 else if (strcmp(name,"SeaiceSpcvx")==0) return SeaiceSpcvxEnum;909 else if (strcmp(name,"SeaiceSpcvy")==0) return SeaiceSpcvyEnum;910 else if (strcmp(name,"SeaiceCoriolisFactor")==0) return SeaiceCoriolisFactorEnum;911 else if (strcmp(name,"BasalforcingsRhoOcean")==0) return BasalforcingsRhoOceanEnum;912 else if (strcmp(name,"BasalforcingsOceanCoef")==0) return BasalforcingsOceanCoefEnum;913 else if (strcmp(name,"BasalforcingsOceanLinDragCoef")==0) return BasalforcingsOceanLinDragCoefEnum;914 else if (strcmp(name,"BasalforcingsOceanQuadDragCoef")==0) return BasalforcingsOceanQuadDragCoefEnum;915 else if (strcmp(name,"BasalforcingsOceanTurningAngle")==0) return BasalforcingsOceanTurningAngleEnum;916 else if (strcmp(name,"BasalforcingsOceanSsh")==0) return BasalforcingsOceanSshEnum;917 else if (strcmp(name,"BasalforcingsOceanVx")==0) return BasalforcingsOceanVxEnum;918 else if (strcmp(name,"BasalforcingsOceanVy")==0) return BasalforcingsOceanVyEnum;919 else if (strcmp(name,"SurfaceforcingsRhoAir")==0) return SurfaceforcingsRhoAirEnum;920 else if (strcmp(name,"SurfaceforcingsAirCoef")==0) return SurfaceforcingsAirCoefEnum;921 else if (strcmp(name,"SurfaceforcingsAirLinDragCoef")==0) return SurfaceforcingsAirLinDragCoefEnum;922 else if (strcmp(name,"SurfaceforcingsAirQuadDragCoef")==0) return SurfaceforcingsAirQuadDragCoefEnum;923 else if (strcmp(name,"SurfaceforcingsWindVx")==0) return SurfaceforcingsWindVxEnum;924 else if (strcmp(name,"SurfaceforcingsWindVy")==0) return SurfaceforcingsWindVyEnum;925 else if (strcmp(name,"Matseaice")==0) return MatseaiceEnum;926 else if (strcmp(name,"MaterialsPoisson")==0) return MaterialsPoissonEnum;927 else if (strcmp(name,"MaterialsYoungModulus")==0) return MaterialsYoungModulusEnum;928 else if (strcmp(name,"MaterialsTimeRelaxationStress")==0) return MaterialsTimeRelaxationStressEnum;929 else if (strcmp(name,"MaterialsTimeRelaxationDamage")==0) return MaterialsTimeRelaxationDamageEnum;930 else if (strcmp(name,"MaterialsRidgingExponent")==0) return MaterialsRidgingExponentEnum;931 else if (strcmp(name,"MaterialsCohesion")==0) return MaterialsCohesionEnum;932 else if (strcmp(name,"MaterialsInternalFrictionCoef")==0) return MaterialsInternalFrictionCoefEnum;933 else if (strcmp(name,"MaterialsCompressionCoef")==0) return MaterialsCompressionCoefEnum;934 else if (strcmp(name,"MaterialsTractionCoef")==0) return MaterialsTractionCoefEnum;935 else if (strcmp(name,"VxStar")==0) return VxStarEnum;936 else if (strcmp(name,"VyStar")==0) return VyStarEnum;937 else if (strcmp(name,"StressTensorPredictorxx")==0) return StressTensorPredictorxxEnum;938 else if (strcmp(name,"StressTensorPredictoryy")==0) return StressTensorPredictoryyEnum;939 else if (strcmp(name,"StressTensorPredictorxy")==0) return StressTensorPredictorxyEnum;940 899 else if (strcmp(name,"MaximumNumberOfDefinitions")==0) return MaximumNumberOfDefinitionsEnum; 941 900 else stage=9; -
issm/trunk-jpl/src/m/classes/model.m
r19050 r19087 40 40 gia = 0; 41 41 42 seaice = 0;43 44 42 autodiff = 0; 45 43 flaim = 0; … … 127 125 end 128 126 end 129 %2014 September 11130 if isa(md.seaice,'double');131 md.seaice=seaice();132 end133 127 %2014 November 12 134 128 if isa(md.calving,'double'); … … 144 138 md=setdefaultparameters(md); 145 139 case 1 146 if strcmpi(varargin{1},'seaice'), 147 md=setdefaultparameters(md); 148 149 %Specific subject for sea ice model 150 md.materials = matseaice(); 151 md.surfaceforcings = seaiceatm(); 152 md.basalforcings = seaiceocean(); 153 md.initialization = seaiceinitialization(); 154 155 %Change some of the defauls 156 md.timestepping.in_years = false; 157 md.timestepping.time_step = 600; 158 md.timestepping.final_time = 600; 159 md.constants.g = 9.80616; %Same as TOPAZ 160 else 161 error('model constructor not supported yet'); 162 end 140 error('model constructor not supported yet'); 163 141 164 142 otherwise … … 1092 1070 md.calving = calving(); 1093 1071 md.gia = gia(); 1094 md.seaice = seaice();1095 1072 md.autodiff = autodiff(); 1096 1073 md.flaim = flaim(); … … 1271 1248 disp(sprintf('%19s: %-22s -- %s','radaroverlay' ,['[1x1 ' class(self.radaroverlay) ']'],'radar image for plot overlay')); 1272 1249 disp(sprintf('%19s: %-22s -- %s','miscellaneous' ,['[1x1 ' class(self.miscellaneous) ']'],'miscellaneous fields')); 1273 disp(sprintf('%19s: %-22s -- %s','seaice' ,['[1x1 ' class(self.seaice) ']'],'parameters for Sea Ice solution'));1274 1250 end % }}} 1275 1251 function memory(self) % {{{ -
issm/trunk-jpl/src/m/enum/EnumDefinitions.py
r19065 r19087 425 425 def GiaSolutionEnum(): return StringToEnum("GiaSolution")[0] 426 426 def GiaAnalysisEnum(): return StringToEnum("GiaAnalysis")[0] 427 def SeaiceSolutionEnum(): return StringToEnum("SeaiceSolution")[0]428 def SeaiceAnalysisEnum(): return StringToEnum("SeaiceAnalysis")[0]429 427 def MeshdeformationSolutionEnum(): return StringToEnum("MeshdeformationSolution")[0] 430 428 def MeshdeformationAnalysisEnum(): return StringToEnum("MeshdeformationAnalysis")[0] … … 870 868 def LevelsetfunctionSlopeYEnum(): return StringToEnum("LevelsetfunctionSlopeY")[0] 871 869 def LevelsetfunctionPicardEnum(): return StringToEnum("LevelsetfunctionPicard")[0] 872 def SeaiceatmEnum(): return StringToEnum("Seaiceatm")[0]873 def SeaiceoceanEnum(): return StringToEnum("Seaiceocean")[0]874 def SeaiceThicknessEnum(): return StringToEnum("SeaiceThickness")[0]875 def SeaiceConcentrationEnum(): return StringToEnum("SeaiceConcentration")[0]876 def SeaiceMinConcentrationEnum(): return StringToEnum("SeaiceMinConcentration")[0]877 def SeaiceMinThicknessEnum(): return StringToEnum("SeaiceMinThickness")[0]878 def SeaiceMaxThicknessEnum(): return StringToEnum("SeaiceMaxThickness")[0]879 def SeaiceSpcvxEnum(): return StringToEnum("SeaiceSpcvx")[0]880 def SeaiceSpcvyEnum(): return StringToEnum("SeaiceSpcvy")[0]881 def SeaiceCoriolisFactorEnum(): return StringToEnum("SeaiceCoriolisFactor")[0]882 def BasalforcingsRhoOceanEnum(): return StringToEnum("BasalforcingsRhoOcean")[0]883 def BasalforcingsOceanCoefEnum(): return StringToEnum("BasalforcingsOceanCoef")[0]884 def BasalforcingsOceanLinDragCoefEnum(): return StringToEnum("BasalforcingsOceanLinDragCoef")[0]885 def BasalforcingsOceanQuadDragCoefEnum(): return StringToEnum("BasalforcingsOceanQuadDragCoef")[0]886 def BasalforcingsOceanTurningAngleEnum(): return StringToEnum("BasalforcingsOceanTurningAngle")[0]887 def BasalforcingsOceanSshEnum(): return StringToEnum("BasalforcingsOceanSsh")[0]888 def BasalforcingsOceanVxEnum(): return StringToEnum("BasalforcingsOceanVx")[0]889 def BasalforcingsOceanVyEnum(): return StringToEnum("BasalforcingsOceanVy")[0]890 def SurfaceforcingsRhoAirEnum(): return StringToEnum("SurfaceforcingsRhoAir")[0]891 def SurfaceforcingsAirCoefEnum(): return StringToEnum("SurfaceforcingsAirCoef")[0]892 def SurfaceforcingsAirLinDragCoefEnum(): return StringToEnum("SurfaceforcingsAirLinDragCoef")[0]893 def SurfaceforcingsAirQuadDragCoefEnum(): return StringToEnum("SurfaceforcingsAirQuadDragCoef")[0]894 def SurfaceforcingsWindVxEnum(): return StringToEnum("SurfaceforcingsWindVx")[0]895 def SurfaceforcingsWindVyEnum(): return StringToEnum("SurfaceforcingsWindVy")[0]896 def MatseaiceEnum(): return StringToEnum("Matseaice")[0]897 def MaterialsPoissonEnum(): return StringToEnum("MaterialsPoisson")[0]898 def MaterialsYoungModulusEnum(): return StringToEnum("MaterialsYoungModulus")[0]899 def MaterialsTimeRelaxationStressEnum(): return StringToEnum("MaterialsTimeRelaxationStress")[0]900 def MaterialsTimeRelaxationDamageEnum(): return StringToEnum("MaterialsTimeRelaxationDamage")[0]901 def MaterialsRidgingExponentEnum(): return StringToEnum("MaterialsRidgingExponent")[0]902 def MaterialsCohesionEnum(): return StringToEnum("MaterialsCohesion")[0]903 def MaterialsInternalFrictionCoefEnum(): return StringToEnum("MaterialsInternalFrictionCoef")[0]904 def MaterialsCompressionCoefEnum(): return StringToEnum("MaterialsCompressionCoef")[0]905 def MaterialsTractionCoefEnum(): return StringToEnum("MaterialsTractionCoef")[0]906 def VxStarEnum(): return StringToEnum("VxStar")[0]907 def VyStarEnum(): return StringToEnum("VyStar")[0]908 def StressTensorPredictorxxEnum(): return StringToEnum("StressTensorPredictorxx")[0]909 def StressTensorPredictoryyEnum(): return StringToEnum("StressTensorPredictoryy")[0]910 def StressTensorPredictorxyEnum(): return StringToEnum("StressTensorPredictorxy")[0]911 870 def MaximumNumberOfDefinitionsEnum(): return StringToEnum("MaximumNumberOfDefinitions")[0]
Note:
See TracChangeset
for help on using the changeset viewer.