Changeset 26481

Timestamp:

10/13/21 17:34:00 (3 years ago)

Author:

Mathieu Morlighem

Message:

CHG: cosmetics

File:

• issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp (modified) (2 diffs)

issm/trunk-jpl/src/c/modules/SurfaceMassBalancex/SurfaceMassBalancex.cpp

@@ -208 +208 @@
         femmodel->parameters->FindParam(&arorder,SmbAutoregressiveOrderEnum);
         IssmDouble tinit_ar;
-        IssmDouble* beta0                = xNew<IssmDouble>(numbasins);
-        IssmDouble* beta1                = xNew<IssmDouble>(numbasins);
-        IssmDouble* phi                  = xNew<IssmDouble>(numbasins*arorder);
-        IssmDouble* noiseterms           = xNew<IssmDouble>(numbasins);
-        IssmDouble* covmat               = xNew<IssmDouble>(numbasins*numbasins);
+        IssmDouble* beta0      = xNew<IssmDouble>(numbasins);
+        IssmDouble* beta1      = xNew<IssmDouble>(numbasins);
+        IssmDouble* phi        = xNew<IssmDouble>(numbasins*arorder);
+        IssmDouble* noiseterms = xNew<IssmDouble>(numbasins);
+        IssmDouble* covmat     = xNew<IssmDouble>(numbasins*numbasins);
         femmodel->parameters->FindParam(&tinit_ar,SmbAutoregressionInitialTimeEnum);
         femmodel->parameters->FindParam(&beta0,&M,SmbBeta0Enum);    _assert_(M==numbasins);
@@ -218 +218 @@
         femmodel->parameters->FindParam(&phi,&M,&Nphi,SmbPhiEnum);  _assert_(M==numbasins); _assert_(Nphi==arorder);
         femmodel->parameters->FindParam(&covmat,&M,&N,SmbCovmatEnum); _assert_(M==numbasins); _assert_(N==numbasins);
-        IssmDouble telapsed_ar = time-tinit_ar; //time elapsed with respect to AR model initial time
+
+        /*time elapsed with respect to AR model initial time*/
+        IssmDouble telapsed_ar = time-tinit_ar; 
+
         /*Before looping through elements: compute noise term specific to each basin from covmat*/
         multivariateNormal(&noiseterms,numbasins,0.0,covmat);
+
         /*Loop over each element to compute SMB at vertices*/
         for(Object* &object:femmodel->elements->objects){
-                Element* element      = xDynamicCast<Element*>(object); //generate element object
+                Element* element = xDynamicCast<Element*>(object);
                 element->Smbautoregression(isstepforar,arorder,telapsed_ar,beta0,beta1,phi,noiseterms);
         }
+
         /*Cleanup*/
         xDelete<IssmDouble>(beta0);