Changeset 21326 for issm/branches/trunk-larour-NatGeoScience2016/src/m/classes/slr.m

Timestamp:

10/29/16 20:14:32 (8 years ago)

Author:

Eric.Larour

Message:

CHG: new steric rate formualtion.

File:

• issm/branches/trunk-larour-NatGeoScience2016/src/m/classes/slr.m (modified) (6 diffs)

issm/branches/trunk-larour-NatGeoScience2016/src/m/classes/slr.m

@@ -21 +21 @@
                 rotation       = 0;
                 ocean_area_scaling = 0;
+                steric_rate    = 0; %rate of ocean expansion from steric effects. 
                 degacc         = 0;
                 requested_outputs      = {};
@@ -55 +56 @@
                 %numerical discretization accuracy
                 self.degacc=.01;
+
+                %steric: 
+                self.steric_rate=0;
                 
                 %output default:
@@ -77 +81 @@
                         md = checkfield(md,'fieldname','slr.abstol','size',[1 1]);
                         md = checkfield(md,'fieldname','slr.maxiter','size',[1 1],'>=',1);
+                        md = checkfield(md,'fieldname','slr.steric_rate','size',[1 1]);
                         md = checkfield(md,'fieldname','slr.degacc','size',[1 1],'>=',1e-10);
                         md = checkfield(md,'fieldname','slr.requested_outputs','stringrow',1);
@@ -113 +118 @@
                         fielddisplay(self,'rotation','earth rotational potential perturbation');
                         fielddisplay(self,'ocean_area_scaling','correction for model representation of ocean area [default: No correction]'); 
+                        fielddisplay(self,'steric_rate','rate of steric ocean expansion (in mm/yr)'); 
                         fielddisplay(self,'rigid','rigid earth graviational potential perturbation');
                         fielddisplay(self,'elastic','elastic earth graviational potential perturbation');
@@ -136 +142 @@
                         WriteData(fid,prefix,'object',self,'fieldname','rotation','format','Boolean');
                         WriteData(fid,prefix,'object',self,'fieldname','ocean_area_scaling','format','Boolean');
+                        WriteData(fid,prefix,'object',self,'fieldname','steric_rate','format','Double','scale',1e-3/yts);
                         WriteData(fid,prefix,'object',self,'fieldname','degacc','format','Double');
                         WriteData(fid,prefix,'object',self,'fieldname','transitions','format','MatArray');
@@ -165 +172 @@
                         writejsdouble(fid,[modelname '.slr.rotation'],self.rotation);
                         writejsdouble(fid,[modelname '.slr.ocean_area_scaling'],self.ocean_area_scaling);
+                        writejsdouble(fid,[modelname '.slr.steric_rate'],self.steric_rate);
                         writejsdouble(fid,[modelname '.slr.elastic'],self.elastic);
                         writejsdouble(fid,[modelname '.slr.degacc'],self.degacc);