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

Timestamp:

12/05/14 13:24:05 (10 years ago)

Author:

lemorzad

Message:

minor changes to PDD method

File:

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

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

@@ -11 +11 @@
                 s0p                       = 0;
                 s0t                       = 0;
-      rlaps                     = 0;
-      rlapslgm                  = 0;                
-      Pfac                      = NaN;
-      Tdiff                     = NaN;
-      sealev                    = NaN;
+                rlaps                     = 0;
+                rlapslgm                  = 0;                
+                Pfac                      = NaN;
+                Tdiff                     = NaN;
+                sealev                    = NaN;
                 isdelta18o                = 0;
                 delta18o                  = NaN;
@@ -22 +22 @@
                 temperatures_lgm          = NaN;
                 precipitations_presentday = NaN;
+                precipitations_lgm        = NaN;
         end
         methods
@@ -39 +40 @@
                         if(self.isdelta18o),self.temperatures_presentday=project3d(md,'vector',self.temperatures_presentday,'type','node');end
                         if(self.isdelta18o),self.precipitations_presentday=project3d(md,'vector',self.precipitations_presentday,'type','node');end
+                        if(self.isdelta18o),self.precipitations_lgm=project3d(md,'vector',self.precipitations_lgm,'type','node');end
 
 
@@ -80 +82 @@
                                         md = checkfield(md,'fieldname','surfaceforcings.temperatures_lgm','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.precipitations_lgm','size',[md.mesh.numberofvertices+1 12],'NaN',1);
                                 end
                         end
@@ -96 +99 @@
                         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,'monthlytemperatures',['CURRENTLY NOT USED monthly surface temperatures [K], required if pdd is activated and delta18o not activated']);
+                        fielddisplay(obj,'precipitation',['CURRENTLY NOT USED monthly surface precipitation [m/yr water eq]']);
                         fielddisplay(obj,'temperatures_presentday','monthly present day surface temperatures [K], required if pdd is activated and delta18o activated');
                         fielddisplay(obj,'temperatures_lgm','monthly LGM surface temperatures [K], required if pdd is activated and delta18o activated');
                         fielddisplay(obj,'precipitations_presentday','monthly surface precipitation [m/yr water eq], required if pdd is activated and delta18o activated');
+                        fielddisplay(obj,'precipitations_lgm','monthly surface precipitation [m/yr water eq], required if pdd is activated and delta18o activated');
                         fielddisplay(obj,'delta18o','delta18o, required if pdd is activated and delta18o activated');
                         fielddisplay(obj,'delta18o_surface','surface elevation of the delta18o site, required if pdd is activated and delta18o activated');
@@ -125 +129 @@
                                 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','temperatures_lgm','format','DoubleMat','mattype',1);
                                 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','precipitations_presentday','format','DoubleMat','mattype',1);
+                                WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','precipitations_lgm','format','DoubleMat','mattype',1);
                                 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','delta18o_surface','format','DoubleMat','mattype',1);
                                 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','delta18o','format','DoubleMat','mattype',1);