Changeset 27576

Timestamp:

02/15/23 12:06:08 (2 years ago)

Author:

Mathieu Morlighem

Message:

CHG: minor

File:

• issm/trunk-jpl/src/m/classes/basalforcings.m (modified) (5 diffs)

issm/trunk-jpl/src/m/classes/basalforcings.m

@@ -8 +8 @@
                 groundedice_melting_rate  = NaN;
                 floatingice_melting_rate  = NaN;
+                perturbation_melting_rate = NaN;
                 geothermalflux            = NaN;
         end
@@ -26 +27 @@
                         fielddisplay(self,'groundedice_melting_rate','basal melting rate (positive if melting) [m/yr]');
                         fielddisplay(self,'floatingice_melting_rate','basal melting rate (positive if melting) [m/yr]');
+                        fielddisplay(self,'perturbation_melting_rate','(optional) perturbation in basal melting rate under floating ice [m/yr]');
                         fielddisplay(self,'geothermalflux','geothermal heat flux [W/m^2]');
 
@@ -32 +34 @@
                         self.groundedice_melting_rate=project3d(md,'vector',self.groundedice_melting_rate,'type','node','layer',1); 
                         self.floatingice_melting_rate=project3d(md,'vector',self.floatingice_melting_rate,'type','node','layer',1); 
+                        self.perturbation_melting_rate=project3d(md,'vector',self.perturbation_melting_rate,'type','node','layer',1); 
                         self.geothermalflux=project3d(md,'vector',self.geothermalflux,'type','node','layer',1); %bedrock only gets geothermal flux
                 end % }}}
@@ -74 +77 @@
                         WriteData(fid,prefix,'object',self,'fieldname','floatingice_melting_rate','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',yts)
                         WriteData(fid,prefix,'object',self,'fieldname','geothermalflux','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'yts',yts);
+                        WriteData(fid,prefix,'object',self,'fieldname','perturbation_melting_rate','format','DoubleMat','name','md.basalforcings.perturbation_melting_rate','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts)
                 end % }}}
                 function savemodeljs(self,fid,modelname) % {{{
@@ -79 +83 @@
                         writejs1Darray(fid,[modelname '.basalforcings.groundedice_melting_rate'],self.groundedice_melting_rate);
                         writejs1Darray(fid,[modelname '.basalforcings.floatingice_melting_rate'],self.floatingice_melting_rate);
+                        writejs1Darray(fid,[modelname '.basalforcings.perturbation_melting_rate'],self.perturbation_melting_rate);
                         writejs1Darray(fid,[modelname '.basalforcings.geothermalflux'],self.geothermalflux);