Changeset 20500 for issm/trunk/src/m/classes/calvinglevermann.m
- Timestamp:
- 04/12/16 21:32:01 (9 years ago)
- Location:
- issm/trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
- Property svn:ignore
-
old new 1 build-js 2 build-esmf 3 build-gcm 1 4 build-fw 2 5 build-ad
-
- Property svn:mergeinfo changed
/issm/trunk-jpl merged: 19104,19106-19126,19128-19134,19136-19170,19172-19299,19302,19306-19405,19407-19604,19606-19668,19670-20496
- Property svn:ignore
-
issm/trunk/src
- Property svn:mergeinfo changed
-
issm/trunk/src/m/classes/calvinglevermann.m
r19105 r20500 6 6 classdef calvinglevermann 7 7 properties (SetAccess=public) 8 stabilization = 0;9 8 coeff = NaN; 10 9 meltingrate = NaN; … … 35 34 function self = setdefaultparameters(self) % {{{ 36 35 37 %stabilization = 2 by default38 self.stabilization = 2;39 40 36 %Proportionality coefficient in Levermann model 41 37 self.coeff=2e13; … … 43 39 function md = checkconsistency(self,md,solution,analyses) % {{{ 44 40 %Early return 45 if (solution~=TransientSolutionEnum() | md.transient.is calving==0), return; end41 if (solution~=TransientSolutionEnum() | md.transient.ismovingfront==0), return; end 46 42 47 md = checkfield(md,'fieldname','calving.stabilization','values',[0 1 2]);48 43 md = checkfield(md,'fieldname','calving.coeff','>',0,'size',[md.mesh.numberofvertices 1]); 49 md = checkfield(md,'fieldname','calving.meltingrate','NaN',1,' size',[md.mesh.numberofvertices 1],'>=',0);44 md = checkfield(md,'fieldname','calving.meltingrate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1],'>=',0); 50 45 end % }}} 51 46 function disp(self) % {{{ 52 47 disp(sprintf(' Calving Levermann parameters:')); 53 fielddisplay(self,'stabilization','0: no, 1: artificial_diffusivity, 2: streamline upwinding');54 48 fielddisplay(self,'coeff','proportionality coefficient in Levermann model'); 55 49 fielddisplay(self,'meltingrate','melting rate at given location [m/a]'); … … 59 53 yts=365.0*24.0*3600.0; 60 54 WriteData(fid,'enum',CalvingLawEnum(),'data',CalvingLevermannEnum(),'format','Integer'); 61 WriteData(fid,'enum',LevelsetStabilizationEnum(),'data',self.stabilization,'format','Integer');62 55 WriteData(fid,'enum',CalvinglevermannCoeffEnum(),'data',self.coeff,'format','DoubleMat','mattype',1); 63 56 WriteData(fid,'object',self,'fieldname','meltingrate','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'scale',1./yts);
Note:
See TracChangeset
for help on using the changeset viewer.