/[MITgcm]/MITgcm_contrib/llc_hires/llc_4320/lookat.m
ViewVC logotype

Diff of /MITgcm_contrib/llc_hires/llc_4320/lookat.m

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

revision 1.1 by dimitri, Mon Feb 10 15:34:38 2014 UTC revision 1.3 by dimitri, Sat Apr 2 01:27:56 2016 UTC
# Line 1  Line 1 
1  nx=4320;  nx=4320;
2  pn='/nobackupp8/dmenemen/llc/llc_4320/MITgcm/run/';  pn='~dmenemen/llc_4320/MITgcm/run/';
3    
4    %%%%%%%%%%%%
5    nx=4320;
6    pn='~dmenemen/llc_4320/MITgcm/run/';
7    nm='U';
8    ts=468432;
9    fn=[pn nm '.' myint2str(ts,10) '.data'];
10    fld=read_llc_fkij(fn,nx,3);
11    clf
12    quikpcolor(rot90(fld)');
13    caxis([-1 1]/4)
14    thincolorbar
15    
16    %%%%%%%%%%%%%%
17    % SST of San Pedro basin
18    ts=279360;
19    fn=[pn 'Theta.' myint2str(ts,10) '.data'];
20    f=11;
21    ix=3800:4100;
22    iy=3100:3300;
23    T=quikread_llc(fn,nx,1,'real*4',f);
24    clf reset
25    mypcolor(rot90(T(ix,iy),2)')
26    colormap(cmap)
27    caxis([13.7 16.8])
28    thincolorbar
29    set(gca,'xtick',[],'ytick',[])
30    print -dpsc LA
31    
32  %%%%%%%%%%%%  %%%%%%%%%%%%
33  % compare SST off Peru on November 1, 2011  % compare SST off Peru on November 1, 2011
# Line 52  n=1; Line 80  n=1;
80  clf reset  clf reset
81  colormap(cmap)  colormap(cmap)
82  for ts=10368:144:138384;  for ts=10368:144:138384;
83   mydisp(ts)      mydisp(ts)
84   fn=[pn nm '.' myint2str(ts,10) '.data'];      fn=[pn nm '.' myint2str(ts,10) '.data'];
85   if exist(fn)      if exist(fn)
86    fld=[rot90(quikread_llc(fn,nx,1,'real*4',7),1) quikread_llc(fn,nx,1,'real*4',8)];          % fld=[rot90(quikread_llc(fn,nx,1,'real*4',7),1) quikread_llc(fn,nx,1,'real*4',8)];
87    clf          fld=[ rot90(read_llc_fkij(fn,nx,3),1) read_llc_fkij(fn,nx,4,1,1:nx,(2*nx+1):(3*nx))];
88    quikpcolor(rot90(fld(:,1:5000),1)');          clf
89    caxis([0.9 1])          quikpcolor(rot90(fld(:,1:5000),1)');
90    title('llc4320 fractional sea ice concentration')          caxis([0.9 1])
91    thincolorbar          title('llc4320 fractional sea ice concentration')
92    text(3400,760,ts2dte(ts,25,2011,9,10),'color','w')          thincolorbar
93    eval(['print -dtiff -r135 ' pn 'figs/' nm myint2str(n,4)]);          text(3400,760,ts2dte(ts,25,2011,9,10),'color','w')
94    n=n+1;          eval(['print -dtiff -r135 ' pn 'figs/' nm myint2str(n,4)]);
95   end          n=n+1;
96        end
97  end  end
98    
99  %%%%%%%%%%%%  %%%%%%%%%%%%
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22