Changeset 19248

Timestamp:

04/03/15 10:59:22 (10 years ago)

Author:

schlegel

Message:

CHG: new pdd class in python

Location:

issm/trunk-jpl/src/m/classes

Files:

• SMBd18opdd.m (modified) (6 diffs)
• SMBd18opdd.py (added)

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

@@ -11 +11 @@
                 rlaps                     = 0;
                 rlapslgm                  = 0; 
-                dpermil                   = 0; 
+                dpermil                   = 0; 
                 Tdiff                     = NaN;
                 sealev                    = NaN;
                 ismungsm                  = 0;
-                isd18opd                  = 0;
+                isd18opd                  = 0;
                 delta18o                  = NaN;
                 delta18o_surface          = NaN;
@@ -46 +46 @@
 
                   self.ismungsm   = 0;
-                  self.isd18opd   = 1;
+                  self.isd18opd   = 1;
                   self.desfac     = 0.5;
                   self.s0p        = 0;
@@ -52 +52 @@
                   self.rlaps      = 6.5;
                   self.rlapslgm   = 6.5;
-                  self.dpermil    = 2.4;
+                  self.dpermil    = 2.4;
                   
                 end % }}}
@@ -63 +63 @@
                                 md = checkfield(md,'fieldname','surfaceforcings.rlaps','>=',0,'numel',1);
                                 md = checkfield(md,'fieldname','surfaceforcings.rlapslgm','>=',0,'numel',1);
-                                if(self.isd18opd==1) 
+                                if(self.isd18opd==1) 
                                         md = checkfield(md,'fieldname','surfaceforcings.temperatures_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1);
                                         md = checkfield(md,'fieldname','surfaceforcings.precipitations_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1);
                                         md = checkfield(md,'fieldname','surfaceforcings.delta18o','NaN',1);
-                                        md = checkfield(md,'fieldname','surfaceforcings.dpermil','>=',0,'numel',1);
+                                        md = checkfield(md,'fieldname','surfaceforcings.dpermil','>=',0,'numel',1);
                                 end
                         end
@@ -103 +103 @@
                         WriteData(fid,'object',self,'class','surfaceforcings','fieldname','rlaps','format','Double');
                         WriteData(fid,'object',self,'class','surfaceforcings','fieldname','rlapslgm','format','Double');
-                        WriteData(fid,'object',self,'class','surfaceforcings','fieldname','Tdiff','format','DoubleMat','mattype',1);
-                        WriteData(fid,'object',self,'class','surfaceforcings','fieldname','sealev','format','DoubleMat','mattype',1);
+                        WriteData(fid,'object',self,'class','surfaceforcings','fieldname','Tdiff','format','DoubleMat','mattype',1);
+                        WriteData(fid,'object',self,'class','surfaceforcings','fieldname','sealev','format','DoubleMat','mattype',1);
 
                         if self.isd18opd
@@ -110 +110 @@
                                 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','precipitations_presentday','format','DoubleMat','mattype',1);
                                 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','delta18o','format','DoubleMat','mattype',1);
-                                WriteData(fid,'object',self,'class','surfaceforcings','fieldname','dpermil','format','Double');
+                                WriteData(fid,'object',self,'class','surfaceforcings','fieldname','dpermil','format','Double');
                         end
                 end % }}}