Changeset 15621 for issm/trunk-jpl/src/m/classes/stressbalance.py

Timestamp:

07/25/13 14:54:51 (12 years ago)

Author:

Mathieu Morlighem

Message:

CHG: spcvx vy and vy need to be scaled before they are marshalled, removed unit conversion in CreateConstraintsDiagnosticHoriz and preparing HO with
quadratic elements

File:

• issm/trunk-jpl/src/m/classes/stressbalance.py (modified) (1 diff)

issm/trunk-jpl/src/m/classes/stressbalance.py

@@ -154 +154 @@
         # }}}
         def marshall(self,md,fid):    # {{{
-                WriteData(fid,'object',self,'class','diagnostic','fieldname','spcvx','format','DoubleMat','mattype',1,'forcinglength',md.mesh.numberofvertices+1)
-                WriteData(fid,'object',self,'class','diagnostic','fieldname','spcvy','format','DoubleMat','mattype',1,'forcinglength',md.mesh.numberofvertices+1)
-                WriteData(fid,'object',self,'class','diagnostic','fieldname','spcvz','format','DoubleMat','mattype',1,'forcinglength',md.mesh.numberofvertices+1)
+
+                yts=365.0*24.0*3600.0
+
+                WriteData(fid,'object',self,'class','diagnostic','fieldname','spcvx','format','DoubleMat','mattype',1,'scale',1./yts,'forcinglength',md.mesh.numberofvertices+1)
+                WriteData(fid,'object',self,'class','diagnostic','fieldname','spcvy','format','DoubleMat','mattype',1,'scale',1./yts,'forcinglength',md.mesh.numberofvertices+1)
+                WriteData(fid,'object',self,'class','diagnostic','fieldname','spcvz','format','DoubleMat','mattype',1,'scale',1./yts,'forcinglength',md.mesh.numberofvertices+1)
                 WriteData(fid,'object',self,'class','diagnostic','fieldname','restol','format','Double')
                 WriteData(fid,'object',self,'class','diagnostic','fieldname','reltol','format','Double')