Changeset 13109

Timestamp:

08/21/12 13:15:42 (13 years ago)

Author:

seroussi

Message:

CHG: now use .nc archives for Nightly Runs

Location:

issm/trunk-jpl/test/NightlyRun

Files:

• runme.m (modified) (2 diffs)
• runme2.m (deleted)

issm/trunk-jpl/test/NightlyRun/runme.m

@@ -124 +124 @@
                 archive_name=['Archive' num2str(id) ];
                 if strcmpi(procedure,'update'),
-
-                        if ~strcmp(oshostname(),'larsen');
-                                error(['Nightly run archives must be saved on "larsen" (hostname is "' oshostname() '")']);
-                        end
+                        delete(['../Archives/' archive_name '.nc'])
                         for k=1:length(field_names),
                                 field=field_values{k};
-                                eval([ archive_name '_field' num2str(k) ' =  field ;']);
+                                % matlab writes the dimensions reversed and matrices transposed into netcdf, so compensate for that
+                                nccreate(['../Archives/' archive_name '.nc'],[archive_name '_field' num2str(k)],...
+                                         'Dimensions',{[archive_name '_field' num2str(k) '_2'] size(field,2) [archive_name '_field' num2str(k) '_1'] size(field,1)},...
+                                         'Format','classic');
+                                ncwrite(['../Archives/' archive_name '.nc'],[archive_name '_field' num2str(k)],transpose(field));
                         end
-                        eval(['save ../Archives/' archive_name ' ' archive_name '_field*']);
-                        disp(sprintf(['File ./../Archives/' archive_name ' saved\n']));
+                        disp(sprintf(['File ./../Archives/' archive_name '.nc saved\n']));
 
                 %ELSE: CHECK TEST
                 else,
-
-                        %load archive
-                        load(['../Archives/' archive_name ]);
-
                         for k=1:length(field_names),
 
@@ -150 +146 @@
 
                                         %compare to archive
-                                        eval(['archive=' archive_name '_field' num2str(k) ';']);
+                                        % matlab reads the dimensions reversed and matrices transposed from netcdf, so compensate for that
+                                        archive=transpose(ncread(['../Archives/' archive_name '.nc'],[archive_name '_field' num2str(k)]));
                                         error_diff=full(max(abs(archive(:)-field(:)))/(max(abs(archive))+eps));