Changeset 8407

Timestamp:

05/24/11 09:54:40 (14 years ago)

Author:

schlegel

Message:

get forcing length using ength(fieldnames(md.forcings))

File:

• issm/trunk/src/m/model/extrude.m (modified) (1 diff)

issm/trunk/src/m/model/extrude.m

@@ -151 +151 @@
 md.vel_obs_raw=project3d(md,md.vel_obs_raw,'node');
 md.surface_mass_balance=project3d(md,md.surface_mass_balance,'node');
+md.surface_accumulation_rate=project3d(md,md.surface_accumulation_rate,'node');
+md.surface_ablation_rate=project3d(md,md.surface_ablation_rate,'node');
 md.dhdt=project3d(md,md.dhdt,'node');
 md.firn_layer=project3d(md,md.firn_layer,'node',md.numlayers);
 if ~isempty(md.forcings),
         forcingnames=fieldnames(md.forcings);
-        numforcings=length(md.forcings);
+        numforcings=length(fieldnames(md.forcings));
         for i=1:numforcings
                 forcing=md.forcings.(forcingnames{i});