Changeset 14611

Timestamp:

04/16/13 16:20:56 (12 years ago)

Author:

Mathieu Morlighem

Message:

CHG: netcdf(md) now writes results too

Location:

issm/trunk-jpl/src/m

Files:

• classes/model/model.m (modified) (3 diffs)
• io/structtonc.m (added)

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

@@ -1139 +1139 @@
                 function memory(obj) % {{{
                         
-                disp(sprintf('\nMemory imprint: '));
+                disp(sprintf('\nMemory imprint:\n'));
 
                 fields=properties('model');
@@ -1148 +1148 @@
                         s=whos('field'); 
                         mem=mem+s.bytes/1e6;
-                        disp(sprintf('%19s: %g Mb',fields{i},s.bytes/1e6));
+                        disp(sprintf('%19s: %6.2f Mb',fields{i},s.bytes/1e6));
                 end
                 disp(sprintf('%19s--%10s','--------------','--------------'));
@@ -1166 +1166 @@
                 ncid=netcdf.create(filename,'CLOBBER');
                 netcdf.putAtt(ncid,netcdf.getConstant('NC_GLOBAL'),'Conventions','CF-1.4');
-                netcdf.putAtt(ncid,netcdf.getConstant('NC_GLOBAL'),'Title',['ISSM model' obj.miscellaneous.name]);
+                netcdf.putAtt(ncid,netcdf.getConstant('NC_GLOBAL'),'Title',['ISSM model (' obj.miscellaneous.name ')']);
                 netcdf.putAtt(ncid,netcdf.getConstant('NC_GLOBAL'),'Author',getenv('USER'));
                 netcdf.putAtt(ncid,netcdf.getConstant('NC_GLOBAL'),'Date',datestr(now));
 
-                %2. prepare some simple dimensions
+                %Preallocate variable id, needed to write variables in netcdf file
                 var_id=zeros(300,1);%preallocate
-
-                %2. First round: define variables (step=1), put variable (step=2)
-                fields1=properties('model');
 
                 for step=1:2,
                         counter=0;
-                        for i=1:length(fields1),
-                                fields2=fields(obj.(fields1{i}));
-                                for j=1:length(fields2)
-                                        counter   = counter+1;
-                                        fieldname = [fields1{i} '.' fields2{j}];
-                                        field     = obj.(fields1{i}).(fields2{j});
-                                        if isa(field,'double') & numel(field)==1,
-                                                if step==1,
-                                                        var_id(counter) = netcdf.defVar(ncid,fieldname,'NC_DOUBLE',[]);
-                                                else
-                                                        netcdf.putVar(ncid,var_id(counter),field);
-                                                end
-                                        elseif isa(field,'double') & size(field,2)==1,
-                                                if step==1,
-                                                        dim_id          = netcdf.defDim(ncid,[fieldname '_size1'],size(field,1));
-                                                        var_id(counter) = netcdf.defVar(ncid,fieldname,'NC_DOUBLE',dim_id);
-                                                else
-                                                        netcdf.putVar(ncid,var_id(counter),field);
-                                                end
-                                        elseif isa(field,'double') & size(field,2)>1,
-                                                if step==1,
-                                                        dim1_id         = netcdf.defDim(ncid,[fieldname '_size1'],size(field,1));
-                                                        dim2_id         = netcdf.defDim(ncid,[fieldname '_size2'],size(field,2));
-                                                        var_id(counter) = netcdf.defVar(ncid,fieldname,'NC_DOUBLE',[dim2_id dim1_id]);
-                                                else
-                                                        netcdf.putVar(ncid,var_id(counter),transpose(field));
-                                                end
-                                        elseif isa(field,'char') & ~isempty(field),
-                                                if step==1,
-                                                        dim_id          = netcdf.defDim(ncid,[fieldname '_size1'],numel(field));
-                                                        var_id(counter) = netcdf.defVar(ncid,fieldname,'NC_CHAR',dim_id);
-                                                else
-                                                        netcdf.putVar(ncid,var_id(counter),transpose(field));
-                                                end
-                                        elseif isa(field,'logical') & numel(field)==1,
-                                                if step==1,
-                                                        var_id(counter) = netcdf.defVar(ncid,fieldname,'NC_BYTE',[]);
-                                                else
-                                                        netcdf.putVar(ncid,var_id(counter),int8(field));
-                                                end
-                                        else
-                                                disp(['skipping ' fieldname '...']);
-                                        end
-                                end
-
-                        end
-
+                        var_id=structtonc(ncid,'md',obj,0,var_id,counter,step);
                         if step==1, netcdf.endDef(ncid); end
                 end
+
                 if counter>300,
                         warning('preallocation of var_id need to be updated');