Changeset 24044 for issm/trunk-jpl/src/m/classes/basalforcingsismip6.m

Timestamp:

06/24/19 07:44:12 (6 years ago)

Author:

seroussi

Message:

NEW: added anomaly to ISMIP6 forcing

File:

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

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

@@ -15 +15 @@
                 geothermalflux            = NaN;
                 groundedice_melting_rate  = NaN;
+                melt_anomaly              = NaN;
         end
         methods
@@ -24 +25 @@
                         self.geothermalflux=project3d(md,'vector',self.geothermalflux,'type','element','layer',1); %bedrock only gets geothermal flux
                         self.groundedice_melting_rate=project3d(md,'vector',self.groundedice_melting_rate,'type','node','layer',1);
+                        self.melt_anomaly=project3d(md,'vector',self.melt_anomaly,'type','element','layer',1); %bedrock only gets geothermal flux
                 end % }}}
                 function self = basalforcingsismip6(varargin) % {{{
@@ -61 +63 @@
                         md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'Inf',1,'>=',0,'timeseries',1);
                         md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1);
+                        if length(md.basalforcings.melt_anomaly)>1,
+                                md = checkfield(md,'fieldname','basalforcings.melt_anomaly','NaN',1,'Inf',1,'timeseries',1);
+                        end
 
                         md = checkfield(md,'fieldname','basalforcings.tf','size',[1,1,numel(md.basalforcings.tf_depths)]);
@@ -79 +84 @@
                         fielddisplay(self,'geothermalflux','geothermal heat flux (W/m^2)');
                         fielddisplay(self,'groundedice_melting_rate','basal melting rate (positive if melting) (m/yr)');
+                        fielddisplay(self,'melt_anomaly','floating ice basal melt anomaly (m/yr)');
 
                 end % }}}
@@ -95 +101 @@
                         WriteData(fid,prefix,'object',self,'fieldname','geothermalflux','format','DoubleMat','name','md.basalforcings.geothermalflux','mattype',1,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts);
                         WriteData(fid,prefix,'object',self,'fieldname','groundedice_melting_rate','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
+                        WriteData(fid,prefix,'object',self,'fieldname','melt_anomaly','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);
 
                 end % }}}