Changeset 25037

Timestamp:

06/15/20 21:07:29 (5 years ago)

Author:

Eric.Larour

Message:

CHG: ismip6 class.

File:

• issm/trunk-jpl/src/m/contrib/larour/ismip6.m (modified) (4 diffs)

issm/trunk-jpl/src/m/contrib/larour/ismip6.m

@@ -13 +13 @@
                 directories   = {};   %directories where the files are
                 experiments   = {};   %names of experiments
-
+                thickness     = {};   %placeholder for thicknesses
+                deltathickness     = {};   %placeholder for delta thicknesses
+                icemask     = {};   %placeholder for ice masks
+                oceanmask     = {};   %placeholder for ocean masks
+                time     = {};   %placeholder for times
         end
         methods
@@ -46 +50 @@
                                 end
 
+                                %initialize fields: 
+                                self.thickness=cell(self.n,1);
+                                self.deltathickness=cell(self.n,1);
+                                self.icemask=cell(self.n,1);
+                                self.oceanmask=cell(self.n,1);
+
                         end
                 end
@@ -55 +65 @@
                         disp('   CMIP5 Ocean MIP:');
 
-                                fielddisplay(self,'n','number of files');
-                                fielddisplay(self,'root','where are the files for ISMIP6');
-                                fielddisplay(self,'directories','IMSIP6 directories');
+                        fielddisplay(self,'n','number of files');
+                        fielddisplay(self,'root','where are the files for ISMIP6');
+                        fielddisplay(self,'directories','IMSIP6 directories');
                 end % }}}
                 function listexp(self) % {{{
@@ -66 +76 @@
 
                 end % }}}
-        end 
+                function [output,time]=read(self,experiment,field) % {{{
+
+                        %go through list of experiments and find the right one: 
+                        if strcmpi(class(experiment),'double'),
+                                ind=experiment;
+                        elseif strcmpi(class(experiment),'char'),
+                                for i=1:self.n,
+                                        if strcmpi(experiment,self.experiments{i}),
+                                                ind=i;
+                                                break;
+                                        end
+                                end
+                        else 
+                                error(['ismip6 read error message: experiment should be a string or index']);
+                        end
+
+                        if ind==0 
+                                error(['ismip6 read error message: experiment ' experiment ' could not be found!']);
+                        end;
+
+                        %figure out the files in this directory: 
+                        dir=self.directories{ind};
+                        currentdir=pwd;
+                        cd([self.root '/' dir]); 
+                        list=listfiles;
+                        cd(currentdir);
+
+                        %go through list of files and figure out which one starts with the field: 
+                        for i=1:length(list),
+                                file=list{i};
+                                ind=findstr(file,'_');
+                                file_field=file(1:ind-1);
+                                if strcmpi(file_field,field),
+                                        break;
+                                end
+                        end
+
+                        %read file: 
+                        %output=ncread([self.root '/' dir '/' file],'file_field');
+                        time=ncread([self.root '/' dir '/' file],'time');
+                        output=ncread([self.root '/' dir '/' file ],field);
+
+                end % }}}
+                function interpolate(self,field,ismip2mesh,ismip2mesh_correction) % {{{
+
+                        for i=1:self.n,
+                                disp(['reading and interpolating field ' field ' for model ' self.experiments{i}]);
+
+                                %read field from disk: 
+                                [h,t]=self.read(i,field); nt=length(t);
+
+                                %map onto 1 dimension field: 
+                                ht=zeros(size(h,1)*size(h,2),nt);
+                                for j=1:size(h,3),
+                                        hj= h(:,:,j)'; hj=hj(:); ht(:,j)=double(hj);
+                                end
+
+                                %map onto mesh: 
+                                hg=ismip2mesh_correction.*(ismip2mesh*ht); 
+                                pos=find(isnan(hg)); hg(pos)=0;
+
+                                %keep field:
+                                if strcmpi(field,'lithk'),
+                                        self.thickness{i}=hg; 
+                                end
+                                self.time{i}=t;
+                        end
+                end  % }}}
+        end
 end