Changeset 22982

Timestamp:

07/19/18 16:20:56 (7 years ago)

Author:

Eric.Larour

Message:

CHG: returning inertia tensor components for test 2010

File:

• issm/trunk-jpl/src/c/cores/sealevelrise_core.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/cores/sealevelrise_core.cpp

@@ -371 +371 @@
         IssmDouble          *radius    = NULL;
         IssmDouble           eustatic;
+        IssmDouble                      Ixz, Iyz, Izz; 
         int                      loop;
 
@@ -417 +418 @@
 
                 if(rotation){
+
                         /*call rotational feedback  module: */
                         RSLgo_rot = new Vector<IssmDouble>(gsize); RSLgo_rot->Assemble();
-                        femmodel->SealevelriseRotationalFeedback(RSLgo_rot,RSLg_old,NULL,NULL,NULL,latitude,longitude,radius); 
+                        femmodel->SealevelriseRotationalFeedback(RSLgo_rot,RSLg_old,&Ixz,&Iyz,&Izz,latitude,longitude,radius); 
                         RSLgo_rot->Assemble(); 
+                        
+                        /*save changes in inertia tensor as results: */
+                        femmodel->results->AddResult(new GenericExternalResult<IssmDouble>(femmodel->results->Size()+1,SealevelInertiaTensorXZEnum,Ixz));
+                        femmodel->results->AddResult(new GenericExternalResult<IssmDouble>(femmodel->results->Size()+1,SealevelInertiaTensorYZEnum,Iyz));
+                        femmodel->results->AddResult(new GenericExternalResult<IssmDouble>(femmodel->results->Size()+1,SealevelInertiaTensorZZEnum,Izz));
 
                         RSLgo->AXPY(RSLgo_rot,1);