Ignore:
Timestamp:
07/08/10 07:57:49 (15 years ago)
Author:
Mathieu Morlighem
Message:

fixed prognostic2 output result

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk/src/m/utils/Nightly/testsgetfields.m

    r4419 r4455  
    7070elseif (analysis_type==Prognostic2SolutionEnum),
    7171        field_names={'Thickness'};
    72         thickness=zeros(md.numberofgrids,1);
    73         thickness(md.results.Prognostic2Solution.Thickness.index)=md.results.Prognostic2Solution.Thickness.value;
     72        thickness   =zeros(md.numberofgrids,1);
     73        connectivity=sparse(md.results.Prognostic2Solution.Thickness.index(:),1,1);
     74        thickness   =sparse(md.results.Prognostic2Solution.Thickness.index(:),1,md.results.Prognostic2Solution.Thickness.value(:));
     75        thickness=full(thickness./connectivity);
    7476        field_values={thickness};
    7577
Note: See TracChangeset for help on using the changeset viewer.