Changeset 18717 for issm/trunk-jpl/src/m/classes/SMBpdd.m

Timestamp:

10/31/14 10:58:45 (10 years ago)

Author:

lemorzad

Message:

NEW: modified lapse rate to PDD scheme

File:

• issm/trunk-jpl/src/m/classes/SMBpdd.m (modified) (5 diffs)

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

@@ -10 +10 @@
                 desfac                    = 0.5;
                 s0p                       = 0;
+                s0t                       = 0;
+                rlaps                     = 6.5;
+                rlapslgm                  = 6.5;                
+                Pfac                      = NaN;
+                Tdiff                     = NaN;
+                sealev                    = NaN;
                 isdelta18o                = 0;
                 delta18o                  = NaN;
@@ -49 +55 @@
                   obj.desfac     = 0.5;
                   obj.s0p        = 0;
+                  obj.s0t        = 0;
+                  obj.rlaps      = 6.5;
+                  obj.rlapslgm   = 6.5;
+                  
                 end % }}}
                 function md = checkconsistency(obj,md,solution,analyses) % {{{
@@ -55 +65 @@
                                 md = checkfield(md,'fieldname','surfaceforcings.desfac','<=',1,'numel',1);
                                 md = checkfield(md,'fieldname','surfaceforcings.s0p','>=',0,'numel',1);
+                                md = checkfield(md,'fieldname','surfaceforcings.s0t','>=',0,'numel',1);
+                                md = checkfield(md,'fieldname','surfaceforcings.rlaps','>=',0,'numel',1);
+                                md = checkfield(md,'fieldname','surfaceforcings.rlapslgm','>=',0,'numel',1);
+                                md = checkfield(md,'fieldname','surfaceforcings.Pfac','NaN',1);
+                                md = checkfield(md,'fieldname','surfaceforcings.Tdiff','NaN',1);
+                                md = checkfield(md,'fieldname','surfaceforcings.sealev','NaN',1);
                                 if(obj.isdelta18o==0)
                                         md = checkfield(md,'fieldname','surfaceforcings.monthlytemperatures','forcing',1,'NaN',1);
@@ -74 +90 @@
                         fielddisplay(obj,'desfac','desertification elevation factor (between 0 and 1, default is 0.5) [m]');
                         fielddisplay(obj,'s0p','should be set to elevation from precip source (between 0 and a few 1000s m, default is 0) [m]');
-                        fielddisplay(obj,'monthlytemperatures','monthly surface temperatures [K], required if pdd is activated and delta18o not activated');
+                        fielddisplay(obj,'s0t','should be set to elevation from temperature source (between 0 and a few 1000s m, default is 0) [m]');
+                        fielddisplay(obj,'rlaps','present day lapse rate [degree/km]');
+                        fielddisplay(obj,'rlapslgm','LGM lapse rate [degree/km]');
+                        fielddisplay(obj,'Pfac','time interpolation parameter for precipitation, 1D(year)');
+                        fielddisplay(obj,'Tdiff','time interpolation parameter for temperature, 1D(year)');
+                        fielddisplay(obj,'sealev','sea level [m], 1D(year)');
+                        fielddisplay(obj,'monthlytemperatures','monthly surface temperatures [K], required if pdd is activated and delta18o not activated');
                         fielddisplay(obj,'precipitation','surface precipitation [m/yr water eq]');
                         fielddisplay(obj,'temperatures_presentday','monthly present day surface temperatures [K], required if pdd is activated and delta18o activated');
@@ -92 +114 @@
                         WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','desfac','format','Double');
                         WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','s0p','format','Double');
+                        WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','s0t','format','Double');
+                        WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','rlaps','format','Double');
+                        WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','rlapslgm','format','Double');
+                        WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','Pfac','format','DoubleMat','mattype',1);
+                        WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','Tdiff','format','DoubleMat','mattype',1);
+                        WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','sealev','format','DoubleMat','mattype',1);
                         WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','isdelta18o','format','Boolean');
                         if obj.isdelta18o