Changeset 18266

Timestamp:

07/18/14 12:59:35 (11 years ago)

Author:

schlegel

Message:

CHG: cleanup gradients

File:

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

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

@@ -72 +72 @@
                 IssmDouble* b_pos       = xNew<IssmDouble>(numvertices); // Hs-SMB relation parameter
                 IssmDouble* b_neg       = xNew<IssmDouble>(numvertices); // Hs-SMB relation paremeter
-                IssmDouble* h           = xNew<IssmDouble>(numvertices); // ice thickness (m)
                 IssmDouble* s           = xNew<IssmDouble>(numvertices); // surface elevation (m)
                 IssmDouble* smb         = xNew<IssmDouble>(numvertices);
@@ -82 +81 @@
                 element->GetInputListOnVertices(b_neg,SurfaceforcingsBNegEnum);
 
-                /*Recover surface elevatio at vertices: */
-                element->GetInputListOnVertices(h,ThicknessEnum);
+                /*Recover surface elevation at vertices: */
                 element->GetInputListOnVertices(s,SurfaceEnum);
 
@@ -98 +96 @@
                                 smb[v]=Smbref[v]+b_neg[v]*(s[v]-Href[v]);
                         }
-                        smb[v]=smb[v]/rho_ice;      // SMB in m/y ice
+                        smb[v]=smb[v]/1000*rho_water/rho_ice;      // SMB in m/y ice
                 }  //end of the loop over the vertices
 
@@ -107 +105 @@
                 xDelete<IssmDouble>(b_pos);
                 xDelete<IssmDouble>(b_neg);
-                xDelete<IssmDouble>(h);
                 xDelete<IssmDouble>(s);
                 xDelete<IssmDouble>(smb);