Changeset 16144
- Timestamp:
- 09/17/13 12:21:42 (12 years ago)
- Location:
- issm/trunk-jpl/src
- Files:
-
- 21 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/c/analyses/ad_core.cpp
r16141 r16144 34 34 double *xp = NULL; 35 35 int my_rank=IssmComm::GetRank(); 36 37 36 38 37 /*AD mode on?: */ -
issm/trunk-jpl/src/c/analyses/balancevelocity_core.cpp
r16018 r16144 14 14 /*parameters: */ 15 15 bool save_results; 16 17 16 18 17 /*recover parameters: */ -
issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
r16125 r16144 5417 5417 /*FUNCTION Penta::DrainWaterfraction{{{*/ 5418 5418 void Penta::DrainWaterfraction(void){ 5419 5419 5420 5420 /*Intermediaries*/ 5421 5421 bool isenthalpy; … … 5754 5754 /*FUNCTION Penta::CreatePVectorAdjointHO{{{*/ 5755 5755 ElementVector* Penta::CreatePVectorAdjointHO(void){ 5756 5757 5756 5758 5757 /*Nothing to be done if not on surface*/ … … 7592 7591 ElementMatrix* Ke=new ElementMatrix(Ke1,Ke2); 7593 7592 delete Ke1; delete Ke2; 7594 7593 7595 7594 /*Compute HO Matrix with P1 element type\n");*/ 7596 7595 this->element_type=P1Enum; … … 8282 8281 D_scalar=gauss->weight*Jdet; 8283 8282 8284 8285 8283 /*Add stabilization*/ 8286 8284 GetNodalFunctionsP1Derivatives(&dbasis[0][0],&xyz_list[0][0],gauss); -
issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
r16125 r16144 1766 1766 /*FUNCTION Tria::InputUpdateFromVector(IssmDouble* vector, int name, int type);{{{*/ 1767 1767 void Tria::InputUpdateFromVector(IssmDouble* vector, int name, int type){ 1768 1769 1768 1770 1769 /*Check that name is an element input*/ -
issm/trunk-jpl/src/c/classes/FemModel.cpp
r16143 r16144 189 189 /*Open input file on cpu 0: */ 190 190 if(my_rank==0) IOMODEL = pfopen0(inputfilename ,"rb"); 191 191 192 192 /*Open toolkits file: */ 193 193 toolkitsoptionsfid=pfopen(toolkitsfilename,"r"); -
issm/trunk-jpl/src/c/classes/IndependentObject.cpp
r16116 r16144 106 106 if(my_rank==0){ 107 107 if(fread(&pscalar,sizeof(IssmPDouble),1,fid)!=1)_error_("could not read scalar "); 108 108 109 109 /*Now, before we even broadcast this to other nodes, declare the scalar as an independent variable!: */ 110 110 scalar<<=pscalar; 111 111 } 112 112 113 114 113 ISSM_MPI_Bcast(&scalar,1,ISSM_MPI_DOUBLE,0,IssmComm::GetComm()); 115 114 -
issm/trunk-jpl/src/c/modules/AllocateSystemMatricesx/AllocateSystemMatricesx.h
r16126 r16144 11 11 void MatrixNonzeros(int** pd_nnz,int** po_nnz,FemModel* femmodel,int set1enum,int set2enum); 12 12 13 14 13 #endif -
issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateParameters.cpp
r16142 r16144 245 245 CreateParametersDakota(¶meters,iomodel,rootpath,solution_type,analysis_type); 246 246 #endif 247 247 248 248 /*Now, deal with toolkits options, which need to be put into the parameters dataset: */ 249 249 ParseToolkitsOptionsx(parameters,toolkitsoptionsfid); -
issm/trunk-jpl/src/c/modules/ModelProcessorx/Stressbalance/UpdateElementsStressbalance.cpp
r15986 r16144 38 38 iomodel->FetchData(1,FlowequationElementEquationEnum); 39 39 finiteelement_list=xNewZeroInit<int>(iomodel->numberofelements); 40 41 40 42 41 /*Do we have coupling*/ -
issm/trunk-jpl/src/c/modules/ResetConstraintsx/ResetConstraintsx.cpp
r15849 r16144 15 15 16 16 void ResetConstraintsx(FemModel* femmodel){ 17 18 17 19 18 /*Display message*/ -
issm/trunk-jpl/src/c/shared/Elements/DrainageFunctionWaterfraction.cpp
r16026 r16144 27 27 else 28 28 Dret=D2; 29 29 30 30 /*check if dt*Dret>waterfraction. If so, drain whole waterfraction*/ 31 31 if(dt==0.){ -
issm/trunk-jpl/src/c/shared/Elements/LliboutryDuval.cpp
r15843 r16144 6 6 #include "../Numerics/types.h" 7 7 #include "../Exceptions/exceptions.h" 8 9 8 10 9 /* get ice stiffness B from enthalpy, pressure and flow law exponent*/ … … 35 34 * 36 35 * Convert A to B : B = A^(-1/n) */ 37 36 38 37 /*Some physical constants (Aschwanden 2012)*/ 39 38 /*TODO: get those constants from model*/ … … 46 45 /*Intermediaries*/ 47 46 IssmDouble A,B,Tstar,Tpmp,H_sp,waterfraction; 48 47 49 48 _assert_(pressure>0); 50 49 _assert_(enthalpy>0); … … 81 80 * get enthalpy from temperature and water fraction, 82 81 * and use LliboutryDuval(IssmDouble enthalpy, IssmDouble pressure,IssmDouble n) */ 83 82 84 83 IssmDouble rho_ice=910; // kg/m^3 85 84 IssmDouble g=9.81; //kg*m/s^2 … … 103 102 enthalpy=heatcapacity*(273.15 - Tref) + waterfraction*latentheat; 104 103 } 105 104 106 105 B=LliboutryDuval(enthalpy, pressure, n); 107 106 108 107 return B; 109 108 } 110 -
issm/trunk-jpl/src/c/shared/MemOps/MemOps.cpp
r16113 r16144 5 5 * Author: utke 6 6 */ 7 8 7 9 8 #ifdef HAVE_CONFIG_H -
issm/trunk-jpl/src/c/toolkits/issm/IssmMpiDenseMat.h
r16142 r16144 513 513 pf=(IssmMpiVec<IssmDouble>*)pfin; 514 514 515 516 515 switch(IssmSolverTypeFromToolkitOptions()){ 517 516 case MumpsEnum: { -
issm/trunk-jpl/src/c/toolkits/issm/IssmToolkitUtils.cpp
r16142 r16144 65 65 *we try and stick with the Petsc vector types: */ 66 66 vec_type=ToolkitOptions::GetToolkitOptionValue("vec_type"); 67 67 68 68 if (strcmp(vec_type,"mpi")==0){ 69 69 vec_type_enum=MpiEnum; -
issm/trunk-jpl/src/c/toolkits/mpi/issmmpi.cpp
r16085 r16144 110 110 }/*}}}*/ 111 111 int ISSM_MPI_Allreduce(void *sendbuf, void *recvbuf, int count, ISSM_MPI_Datatype datatype, ISSM_MPI_Op op, ISSM_MPI_Comm comm){/*{{{*/ 112 112 113 113 int rc=0; 114 114 #ifdef _HAVE_MPI_ … … 178 178 }/*}}}*/ 179 179 int ISSM_MPI_Comm_free(ISSM_MPI_Comm *comm){ /*{{{*/ 180 180 181 181 int rc=0; 182 182 #ifdef _HAVE_MPI_ … … 203 203 }/*}}}*/ 204 204 int ISSM_MPI_Comm_size( ISSM_MPI_Comm comm, int *size){ /*{{{*/ 205 205 206 206 int rc=0; 207 207 #ifdef _HAVE_MPI_ … … 263 263 }/*}}}*/ 264 264 int ISSM_MPI_Gatherv(void *sendbuf, int sendcnt, ISSM_MPI_Datatype sendtype, void *recvbuf, int *recvcnts, int *displs, ISSM_MPI_Datatype recvtype, int root, ISSM_MPI_Comm comm){/*{{{*/ 265 265 266 266 int rc=0; 267 267 assert(sendtype==recvtype); // we handle only identical representations -
issm/trunk-jpl/src/c/toolkits/mpi/issmmpi.h
r16090 r16144 51 51 #define ISSM_MPI_ANY_TAG MPI_ANY_TAG 52 52 #define ISSM_MPI_ANY_SOURCE MPI_ANY_SOURCE 53 53 54 54 /*other include files: */ 55 55 #include "./commops/commops.h" -
issm/trunk-jpl/src/m/classes/snowpack.m
r15931 r16144 124 124 end % }}} 125 125 function obj = setdefaultparameters(obj) % {{{ 126 126 127 127 %snowpack: %{{{ 128 128 obj.snowpack_meas_tss = 1; … … 279 279 %filters {{{ 280 280 filter_values={'MIN_MAX','RATE_FILTER1','RATE_FILTER2','UNHEATED_RAIN_GAUGE_FILTER','WMO_UNDERCATCH_FILTER','WMO_UNDERCATCH_FILTER-SIMPLIFIED','UNVENTILLATED_TEMPERATURE_SENSOR','ADD_AN_OFFSET'}; 281 282 281 283 282 md=checkfield(md,'snowpack.filters_ta_filter1','values',{filter_values}); -
issm/trunk-jpl/src/m/classes/thermal.m
r16026 r16144 40 40 %Should we use cold ice (default) or enthalpy formulation 41 41 obj.isenthalpy=0; 42 42 43 43 %will basal boundary conditions be set dynamically 44 44 obj.isdynamicbasalspc=0; … … 69 69 fielddisplay(obj,'isenthalpy','use an enthalpy formulation to include temperate ice (default is 0)'); 70 70 fielddisplay(obj,'isdynamicbasalspc',['enable dynamic setting of basal forcing. required for enthalpy formulation (default is 0)']); 71 71 72 72 end % }}} 73 73 function marshall(obj,md,fid) % {{{ -
issm/trunk-jpl/src/m/exp/Shp2Exp.m
r16049 r16144 9 9 % 10 10 % See also EXPMASTER, EXPDOC 11 12 11 13 12 if ~exist(shapefilename,'file'), … … 22 21 ids=find(isnan(x)); 23 22 x(ids)=[]; y(ids)=[]; 24 23 25 24 exp(end+1).x=x; 26 25 exp(end).y=y; … … 33 32 34 33 expwrite(exp,expfilename); 35 -
issm/trunk-jpl/src/wrappers/python/include/wrapper_macros.h
r15107 r16144 43 43 /* WRAPPER 3.2 {{{*/ 44 44 #define WRAPPER(modulename,...) \ 45 \46 45 static PyObject* modulename(PyObject* self,PyObject* args);\ 47 46 static PyMethodDef modulename##_funcs[] = {\ … … 49 48 {NULL,NULL,0,NULL}\ 50 49 };\ 51 \52 50 static struct PyModuleDef modulename##module= {\ 53 51 PyModuleDef_HEAD_INIT,\ … … 58 56 modulename##_funcs\ 59 57 };\ 60 \61 58 PyMODINIT_FUNC PyInit_##modulename(void){\ 62 \63 59 import_array();\ 64 60 return PyModule_Create(&modulename##module);\ 65 61 }\ 66 \67 62 static PyObject* modulename(PyObject* self,PyObject* args) 68 63 /*}}}*/ … … 70 65 /* WRAPPER 2.7 {{{*/ 71 66 #define WRAPPER(modulename,...) \ 72 \73 67 static PyObject* modulename(PyObject* self,PyObject* args);\ 74 68 static PyMethodDef modulename##_funcs[] = {\ … … 76 70 {NULL,NULL,0,NULL}\ 77 71 };\ 78 \79 72 PyMODINIT_FUNC init##modulename(void){\ 80 \81 73 import_array();\ 82 74 (void) Py_InitModule(#modulename, modulename##_funcs);\ 83 75 }\ 84 \85 76 static PyObject* modulename(PyObject* self,PyObject* args) 86 77 /*}}}*/
Note:
See TracChangeset
for help on using the changeset viewer.