/[MITgcm]/MITgcm_contrib/gael/profilesMatlabProcessing/profiles_misc/model_load.m
ViewVC logotype

Diff of /MITgcm_contrib/gael/profilesMatlabProcessing/profiles_misc/model_load.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by roquet, Wed Apr 13 21:16:01 2011 UTC revision 1.3 by gforget, Mon May 5 13:00:40 2014 UTC
# Line 60  switch model Line 60  switch model
60          fileS=[dir_in 'S_OWPv1_M_eccollc_90x50.bin'];          fileS=[dir_in 'S_OWPv1_M_eccollc_90x50.bin'];
61                    
62          %load and save grid :          %load and save grid :
63          global fileFormat; fileFormat='compact';          grid_load(dirGrid,5,'compact'); gcmfaces_bindata;
         grid_load(dirGrid,5); gcmfaces_bindata;  
64          mygrid=rmfield(mygrid,{'XG','YG','RAC','RAZ','DXC','DYC','DXG','DYG'});          mygrid=rmfield(mygrid,{'XG','YG','RAC','RAZ','DXC','DYC','DXG','DYG'});
65          mygrid=rmfield(mygrid,{'hFacC','hFacW','hFacS','Depth','AngleCS','AngleSN'});          mygrid=rmfield(mygrid,{'hFacC','hFacW','hFacS','Depth','AngleCS','AngleSN'});
66          mygrid=rmfield(mygrid,{'hFacCsurf','mskW','mskS','DRC','DRF','RF'});          mygrid=rmfield(mygrid,{'hFacCsurf','mskW','mskS','DRC','DRF','RF'});
67          % list_param={'XC','YC','RAC','RC'};          % list_param={'XC','YC','RAC','RC'};
68          % grid_load(dirGrid,5,list_param);          % grid_load(dirGrid,5,list_param);
69          % gcmfaces_bindata;          % gcmfaces_bindata;
70          MYBASININDEX=convert2array(v4_read_bin(fileBasin,1,0));          MYBASININDEX=convert2array(read_bin(fileBasin,1,0));
71          grid.name=model;          grid.name=model;
72          grid.mygrid=mygrid;          grid.mygrid=mygrid;
73          grid.mytri=mytri;          grid.mytri=mytri;
# Line 77  switch model Line 76  switch model
76          % read T/S Atlas          % read T/S Atlas
77          fldT=mygrid.mskC; fldT(:)=0; fldS=fldT;          fldT=mygrid.mskC; fldT(:)=0; fldS=fldT;
78          for tt=1:12;          for tt=1:12;
79              fldT(:,:,:,tt)=v4_read_bin(fileT,tt).*mygrid.mskC;              fldT(:,:,:,tt)=read_bin(fileT,tt).*mygrid.mskC;
80              fldS(:,:,:,tt)=v4_read_bin(fileS,tt).*mygrid.mskC;              fldS(:,:,:,tt)=read_bin(fileS,tt).*mygrid.mskC;
81          end;          end;
82          atlas.name=model;          atlas.name=model;
83          atlas.T={convert2array(fldT)};  atlas.S={convert2array(fldS)};          atlas.T={convert2array(fldT)};  atlas.S={convert2array(fldS)};
# Line 91  switch model Line 90  switch model
90          dirGrid=[dir_in 'grid/'];          dirGrid=[dir_in 'grid/'];
91    
92          %load grid :          %load grid :
93          global fileFormat; fileFormat='straight';          grid_load(dirGrid,1,'straight'); gcmfaces_bindata;
         grid_load(dirGrid,1); gcmfaces_bindata;  
94          mygrid=rmfield(mygrid,{'XG','YG','RAC','RAZ','DXC','DYC','DXG','DYG'});          mygrid=rmfield(mygrid,{'XG','YG','RAC','RAZ','DXC','DYC','DXG','DYG'});
95          mygrid=rmfield(mygrid,{'hFacC','hFacW','hFacS','Depth','AngleCS','AngleSN'});          mygrid=rmfield(mygrid,{'hFacC','hFacW','hFacS','Depth','AngleCS','AngleSN'});
96          mygrid=rmfield(mygrid,{'hFacCsurf','mskW','mskS','DRC','DRF','RF'});          mygrid=rmfield(mygrid,{'hFacCsurf','mskW','mskS','DRC','DRF','RF'});
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22