Changeset 20033 for issm/trunk-jpl/src/m/classes/slr.m

Timestamp:

01/30/16 19:12:44 (9 years ago)

Author:

Eric.Larour

Message:

CHG: parallelized computation of G_elastic in the Sealevelrise UpdateParameters analysis phase.

File:

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

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

@@ -15 +15 @@
                 elastic         = 0;
                 eustatic         = 0;
+                degacc         = 0;
         end
         methods
@@ -38 +39 @@
                 self.elastic=1;
                 self.eustatic=1;
+
+                %numerical discretization accuracy
+                self.degacc=.01;
                 
                 end % }}}
@@ -49 +53 @@
                         md = checkfield(md,'fieldname','slr.abstol','size',[1 1]);
                         md = checkfield(md,'fieldname','slr.maxiter','size',[1 1],'>=',1);
+                        md = checkfield(md,'fieldname','slr.degacc','size',[1 1],'>=',1e-10);
 
                         %check that love numbers are provided at the same level of accuracy: 
@@ -68 +73 @@
                         fielddisplay(self,'elastic','elastic earth graviational potential perturbation');
                         fielddisplay(self,'eustatic','eustatic sea level rise');
+                        fielddisplay(self,'degacc','accuracy (default .01 deg) for numerical discretization of the Green''s functions');
 
                 end % }}}
@@ -80 +86 @@
                         WriteData(fid,'object',self,'class','sealevelrise','fieldname','elastic','format','Boolean');
                         WriteData(fid,'object',self,'class','sealevelrise','fieldname','eustatic','format','Boolean');
+                        WriteData(fid,'object',self,'class','sealevelrise','fieldname','degacc','format','Double');
                 end % }}}
                 function savemodeljs(self,fid,modelname) % {{{
@@ -91 +98 @@
                         writejsdouble(fid,[modelname '.slr.rigid'],self.rigid);
                         writejsdouble(fid,[modelname '.slr.eustatic'],self.eustatic);
+                        writejsdouble(fid,[modelname '.slr.degacc'],self.degacc);
 
                 end % }}}