Skip to content

Commit

Permalink
BUG: ocean-averaged sea level loads
Browse files Browse the repository at this point in the history
  • Loading branch information
AdhikariJPL committed Oct 1, 2024
1 parent ee8d2e1 commit 40790d0
Showing 1 changed file with 2 additions and 7 deletions.
9 changes: 2 additions & 7 deletions src/c/cores/sealevelchange_core.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -840,13 +840,8 @@ IssmDouble SealevelloadsOceanAverage(GrdLoads* loads, Vector<IssmDouble>* ocean
IssmDouble sealevelloadsaverage;
IssmDouble subsealevelloadsaverage;

Vector<IssmDouble>* vsealevelloadsvolume=loads->vsealevelloads->Duplicate();
Vector<IssmDouble>* vsubsealevelloadsvolume=loads->vsubsealevelloads->Duplicate();

vsealevelloadsvolume->Sum(&sealevelloadsaverage);
vsubsealevelloadsvolume->Sum(&subsealevelloadsaverage);
delete vsealevelloadsvolume;
delete vsubsealevelloadsvolume;
loads->vsealevelloads->Sum(&sealevelloadsaverage);
   loads->vsubsealevelloads->Sum(&subsealevelloadsaverage);

return (sealevelloadsaverage+subsealevelloadsaverage)/totaloceanarea;
} /*}}}*/
Expand Down

0 comments on commit 40790d0

Please sign in to comment.