Index: /issm/trunk-jpl/src/c/analyses/SealevelchangeAnalysis.cpp =================================================================== --- /issm/trunk-jpl/src/c/analyses/SealevelchangeAnalysis.cpp (revision 27096) +++ /issm/trunk-jpl/src/c/analyses/SealevelchangeAnalysis.cpp (revision 27097) @@ -9,4 +9,14 @@ /*Model processing*/ +SealevelchangeAnalysis::~SealevelchangeAnalysis(){/*{{{*/ + + BarystaticContributions* barycontrib=NULL; + GenericParam* barycontribparam=NULL; + + barycontribparam = xDynamicCast*>(femmodel->parameters->FindParamObject(BarystaticContributionsEnum)); + barycontrib=barycontribparam->GetParameterValue(); + + delete barycontrib; +}/*}}}*/ void SealevelchangeAnalysis::CreateConstraints(Constraints* constraints,IoModel* iomodel){/*{{{*/ /*No constraints*/ @@ -285,4 +295,7 @@ parameters->AddObject(new BoolParam(LoveIsTimeEnum,love_istime)); + /*Free allocations:*/ + xDelete(love_timefreq); + // AD performance is sensitive to calls to ensurecontiguous. // // Providing "t" will cause ensurecontiguous to be called. @@ -532,12 +545,17 @@ nt=1; //in elastic, or if we run only selfattraction, we need only one step #ifdef _HAVE_AD_ - G_viscoelastic_interpolated=G_viscoelastic; - U_viscoelastic_interpolated=U_viscoelastic; - if(horiz)H_viscoelastic_interpolated=H_viscoelastic; + G_viscoelastic_interpolated=xNew(M,"t"); + U_viscoelastic_interpolated=xNew(M,"t"); + if(horiz) H_viscoelastic_interpolated=xNew(M,"t"); #else - G_viscoelastic_interpolated=G_viscoelastic; - U_viscoelastic_interpolated=U_viscoelastic; - if(horiz)H_viscoelastic_interpolated=H_viscoelastic; + G_viscoelastic_interpolated=xNew(M); + U_viscoelastic_interpolated=xNew(M); + if(horiz) H_viscoelastic_interpolated=xNew(M); #endif + + + xMemCpy(G_viscoelastic_interpolated,G_viscoelastic,M); + xMemCpy(U_viscoelastic_interpolated,U_viscoelastic,M); + if(horiz) xMemCpy(H_viscoelastic_interpolated,H_viscoelastic,M); if(rotation){ //if this cpu handles degree 2 @@ -558,6 +576,6 @@ #endif - Pmtf_col_interpolated=love_pmtf_colinear; - Pmtf_ortho_interpolated=love_pmtf_ortho; + Pmtf_col_interpolated[0]=love_pmtf_colinear[0]; + Pmtf_ortho_interpolated[0]=love_pmtf_ortho[0]; Pmtf_z_interpolated[0]=1.0+love_k[2]; Love_tk2_interpolated[0]=love_tk[2]; @@ -595,8 +613,11 @@ xDelete(G_viscoelastic); xDelete(G_viscoelastic_local); + xDelete(G_viscoelastic_interpolated); xDelete(U_viscoelastic); + xDelete(U_viscoelastic_interpolated); xDelete(U_viscoelastic_local); if(horiz){ xDelete(H_viscoelastic); + xDelete(H_viscoelastic_interpolated); xDelete(H_viscoelastic_local); } @@ -604,4 +625,13 @@ xDelete(love_pmtf_colinear); xDelete(love_pmtf_ortho); + + xDelete(Love_tk2_interpolated); + xDelete(Love_th2_interpolated); + xDelete(Pmtf_col_interpolated); + xDelete(Pmtf_ortho_interpolated); + xDelete(Pmtf_z_interpolated); + if(horiz){ + xDelete(Love_tl2_interpolated); + } } @@ -632,4 +662,5 @@ iomodel->DeleteData(&requestedoutputs,numoutputs,"md.solidearth.requested_outputs"); /*}}}*/ + }/*}}}*/ Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (revision 27096) +++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (revision 27097) @@ -6583,5 +6583,5 @@ LoveRotRSL = xNewZeroInit(nt); LoveRotU = xNewZeroInit(nt); - LoveRothoriz= xNewZeroInit(nt); + if(horiz)LoveRothoriz= xNewZeroInit(nt); Grot = xNewZeroInit(3*3*nt); //3 polar motion components * 3 vertices * number of time steps GUrot = xNewZeroInit(3*3*nt); @@ -6669,4 +6669,8 @@ } } + /*Free ressources:*/ + xDelete(LoveRotRSL); + xDelete(LoveRotU); + if(horiz)xDelete(LoveRothoriz); } /*}}}*/ @@ -7392,4 +7396,7 @@ /*Keep track of barystatic contributions:*/ barycontrib->Set(this->Sid(),bslcice,bslchydro,bslcbp); + + /*Free ressources:*/ + xDelete(areae); } @@ -7708,4 +7715,7 @@ /*Save viscous stack now that we updated the values:*/ this->inputs->SetArrayInput(viscousenum,this->lid,viscousfield,3*viscousnumsteps); + + /*Free ressources:*/ + xDelete(grdfieldinterp); } @@ -7726,4 +7736,6 @@ for(i=0;i(grdfield); } /*}}}*/ Index: /issm/trunk-jpl/src/c/classes/GrdLoads.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/GrdLoads.cpp (revision 27096) +++ /issm/trunk-jpl/src/c/classes/GrdLoads.cpp (revision 27097) @@ -19,4 +19,5 @@ GrdLoads::GrdLoads(int nel,SealevelGeometry* slgeom){ /*{{{*/ + /*allocate:*/ vloads=new Vector(nel); for (int i=0;i(slgeom->nbar[i]); @@ -26,4 +27,11 @@ vsubsealevelloads=new Vector(slgeom->nbar[SLGEOM_OCEAN]); + + /*make sure all pointers that are not allocated are NULL:*/ + loads=NULL; + for (int i=0;i(loads); + for (int i=0;i(subloads[i]); + } + + /*Serialize:*/ loads=vloads->ToMPISerial(); for (int i=0;i(sealevelloads); + if(subsealevelloads)xDelete(subsealevelloads); + + /*Serialize:*/ sealevelloads=vsealevelloads->ToMPISerial(); subsealevelloads=vsubsealevelloads->ToMPISerial(); } /*}}}*/ -void GrdLoads::SHDegree2Coefficients(IssmDouble* deg2coeff, FemModel* femmodel, SealevelGeometry* slgeom){ +void GrdLoads::SHDegree2Coefficients(IssmDouble* deg2coeff, FemModel* femmodel, SealevelGeometry* slgeom){ /*{{{*/ IssmDouble area,re, S; @@ -119,3 +139,3 @@ ISSM_MPI_Bcast(°2coeff[2],1,ISSM_MPI_DOUBLE,0,IssmComm::GetComm()); -} +} /*}}}*/ Index: /issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp =================================================================== --- /issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp (revision 27096) +++ /issm/trunk-jpl/src/c/cores/sealevelchange_core.cpp (revision 27097) @@ -81,4 +81,8 @@ /*End profiler*/ femmodel->profiler->Stop(SLRCORE); + + /*Free ressources:*/ + delete slgeom; + } /*}}}*/ @@ -364,9 +368,14 @@ //convergence? - if(slcconvergence(loads->vsealevelloads,oldsealevelloads,eps_rel,eps_abs))break; + if(slcconvergence(loads->vsealevelloads,oldsealevelloads,eps_rel,eps_abs)){ + delete oldsealevelloads; break; + } //early return? - if(iterations>=max_nonlinear_iterations)break; + if(iterations>=max_nonlinear_iterations){ + delete oldsealevelloads; break; + } iterations++; + delete oldsealevelloads; } @@ -388,5 +397,5 @@ if(planethasocean){ if (!sealevelloading){ //we haven't done so before, so we need to compute the ocean average and area - loads->sealevelloads=NULL; //needed to trigger the calculation of areas + if(loads->sealevelloads)xDelete(loads->sealevelloads); loads->sealevelloads=NULL; //needed to trigger the calculation of areas /*retrieve sea level average and ocean area:*/ for(Object* & object : femmodel->elements->objects){ @@ -430,4 +439,10 @@ } + /*Free ressources:*/ + delete loads; + delete subelementoceanareas; + delete oceanareas; + xDelete(sealevelpercpu); + } /*}}}*/ @@ -607,4 +622,9 @@ #endif + xDelete(xxe); + xDelete(yye); + xDelete(zze); + xDelete(areae); + return; @@ -672,4 +692,10 @@ femmodel->parameters->AddObject(new DoubleVecParam(ZzeEnum,zze,nel)); femmodel->parameters->AddObject(new DoubleVecParam(AreaeEnum,areae,nel)); + + /*Free ressources:*/ + xDelete(xxe); + xDelete(yye); + xDelete(zze); + xDelete(areae); return slgeom; @@ -870,5 +896,10 @@ xDelete(m3interp); } - } + xDelete(viscoustimes); + xDelete(viscouspolarmotion); + } + xDelete(pmtf_col); + xDelete(pmtf_ortho); + xDelete(pmtf_z); } /*}}}*/ Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp =================================================================== --- /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp (revision 27096) +++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/CreateElementsVerticesAndMaterials.cpp (revision 27097) @@ -223,4 +223,5 @@ { Tria* tria= xDynamicCast(element); + if(tria->hmaterial)delete tria->hmaterial; tria->hmaterial=NULL; }