Changeset 22918 for issm/trunk-jpl/src/m/classes/basalforcingspico.m

Timestamp:

07/13/18 13:01:36 (7 years ago)

Author:

tpelle

Message:

CHG: Adding buoyant plume basal melt-rate parameterization to ISSM, adapted from Lazeroms et al. (2018)

File:

• issm/trunk-jpl/src/m/classes/basalforcingspico.m (modified) (6 diffs)

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

@@ -13 +13 @@
                 farocean_temperature      = NaN;
                 farocean_salinity         = NaN;
+                isplume                   = NaN;
                 geothermalflux            = NaN;
                 groundedice_melting_rate  = NaN;
@@ -38 +39 @@
                                         disp('      no maximum number of boxes set, setting value to 5');
                    end
-
                         if isnan(self.overturning_coeff)
                                 self.overturning_coeff = 1e6; %m^3/s
                                 disp('      no overturning strength set, setting value to 1e6');
                         end
-
                         if isnan(self.gamma_T)
                                 self.gamma_T = 2e-5; %m/s
@@ -56 +55 @@
                 function self = setdefaultparameters(self) % {{{
 
-                        self.maxboxcount = 5;
+                        self.maxboxcount       = 5;
                         self.overturning_coeff = 1e6; %m^3/s
-                        self.gamma_T = 2e-5; %m/s
+                        self.gamma_T           = 2e-5; %m/s
+                        self.isplume           = false;
 
                 end % }}}
@@ -70 +70 @@
                                 md = checkfield(md,'fieldname','basalforcings.farocean_temperature','NaN',1,'Inf',1,'>',0,'size',[md.basalforcings.num_basins+1 NaN]);
                                 md = checkfield(md,'fieldname','basalforcings.farocean_salinity','NaN',1,'Inf',1,'>',0,'size',[md.basalforcings.num_basins+1 NaN]);
+                                md = checkfield(md,'fieldname','basalforcings.isplume','values',[0 1]);
                                 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);
@@ -79 +80 @@
                         fielddisplay(self,'basin_id','basin number assigned to each node [unitless]');
                         fielddisplay(self,'maxboxcount','maximum number of boxes initialized under all ice shelves');
-                        fielddisplay(self,'overturning_coeff','Overturning strength [m^3/s]');
-                        fielddisplay(self,'gamma_T','Turbulent temperature exchange velocity [m/s]');
+                        fielddisplay(self,'overturning_coeff','overturning strength [m^3/s]');
+                        fielddisplay(self,'gamma_T','turbulent temperature exchange velocity [m/s]');
                         fielddisplay(self,'farocean_temperature','depth averaged ocean temperature in front of the ice shelf for basin i [K]');
                         fielddisplay(self,'farocean_salinity','depth averaged ocean salinity in front of the ice shelf for basin i [psu]');
+                        fielddisplay(self,'isplume','boolean to use buoyant plume melt rate parameterization from Lazeroms et al., 2018 (default false)');
                         fielddisplay(self,'geothermalflux','geothermal heat flux [W/m^2]');
                         fielddisplay(self,'groundedice_melting_rate','basal melting rate (positive if melting) [m/yr]');
@@ -98 +100 @@
                         WriteData(fid,prefix,'object',self,'fieldname','farocean_temperature','format','DoubleMat','name','md.basalforcings.farocean_temperature','timeserieslength',md.basalforcings.num_basins+1,'yts',md.constants.yts);
                         WriteData(fid,prefix,'object',self,'fieldname','farocean_salinity','format','DoubleMat','name','md.basalforcings.farocean_salinity','timeserieslength',md.basalforcings.num_basins+1,'yts',md.constants.yts);
-                        %WriteData(fid,prefix,'object',self,'fieldname','basin_id','format','DoubleMat','name','md.basalforcings.basin_id','mattype',2);
                         WriteData(fid,prefix,'object',self,'fieldname','basin_id','data',self.basin_id-1,'name','md.basalforcings.basin_id','format','IntMat','mattype',2);   %Change to 0-indexing
+                        WriteData(fid,prefix,'object',self,'fieldname','isplume','format','Boolean');
                         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);