/[MITgcm]/MITgcm_contrib/high_res_cube/input/lookat_uv.m
ViewVC logotype

Diff of /MITgcm_contrib/high_res_cube/input/lookat_uv.m

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

revision 1.1 by dimitri, Wed Dec 10 19:40:13 2003 UTC revision 1.2 by dimitri, Sat Jan 10 00:50:25 2004 UTC
# Line 11  yg=readbin('YG.data',[32,6,32]);yg=permu Line 11  yg=readbin('YG.data',[32,6,32]);yg=permu
11  %tmp=yg; cx=[min(tmp(:)) max(tmp(:))]; plot_cube  %tmp=yg; cx=[min(tmp(:)) max(tmp(:))]; plot_cube
12    
13  % load input velocity on Gaussian NCEP grid  % load input velocity on Gaussian NCEP grid
14  u=readbin('uwind_192_94_8.bin',[192 94]);  u=readbin('uwind_192_94_12.bin',[192 94]);
15  v=readbin('vwind_192_94_8.bin',[192 94]);  v=readbin('vwind_192_94_12.bin',[192 94]);
16  lon=0:1.875:358.125;  lon=0:1.875:358.125;
17  lat=[-88.5420, -86.6532, -84.7532, -82.8508, -80.9474, -79.0435, ...  lat=[-88.5420, -86.6532, -84.7532, -82.8508, -80.9474, -79.0435, ...
18       -77.1393, -75.2351, -73.3307, -71.4262, -69.5217, -67.6171, ...       -77.1393, -75.2351, -73.3307, -71.4262, -69.5217, -67.6171, ...
# Line 81  figure(2), tmp=v3; plot_cube, title('v f Line 81  figure(2), tmp=v3; plot_cube, title('v f
81  un=permute(readbin('../inp_thsice/ncep_uwind_cs.bin',[32,6,32],1,'real*8'),[1 3 2]);  un=permute(readbin('../inp_thsice/ncep_uwind_cs.bin',[32,6,32],1,'real*8'),[1 3 2]);
82  figure(3), tmp=un; cx=[min(tmp(:)) max(tmp(:))]; plot_cube  figure(3), tmp=un; cx=[min(tmp(:)) max(tmp(:))]; plot_cube
83  title('u from ocean.cs32x15')  title('u from ocean.cs32x15')
84    print -djpeg u_cs32
85  vn=permute(readbin('../inp_thsice/ncep_vwind_cs.bin',[32,6,32],1,'real*8'),[1 3 2]);  vn=permute(readbin('../inp_thsice/ncep_vwind_cs.bin',[32,6,32],1,'real*8'),[1 3 2]);
86  figure(4), tmp=vn; cx=[min(tmp(:)) max(tmp(:))]; plot_cube  figure(4), tmp=vn; cx=[min(tmp(:)) max(tmp(:))]; plot_cube
87  title('v from ocean.cs32x15')  title('v from ocean.cs32x15')
88    print -djpeg v_cs32
89    
90  % load model output files  % load model output files
91  um=permute(readbin('UWIND.0000000020.data',[32,6,32]),[1 3 2]);  um=permute(readbin('UWIND.0000000020.data',[32,6,32]),[1 3 2]);
92  figure(5), tmp=um; plot_cube, title('u from model output')  figure(5), tmp=um; plot_cube, title('u from model output')
93    print -djpeg u_exf
94  vm=permute(readbin('VWIND.0000000020.data',[32,6,32]),[1 3 2]);  vm=permute(readbin('VWIND.0000000020.data',[32,6,32]),[1 3 2]);
95  figure(6), tmp=vm; plot_cube, title('v from model output')  figure(6), tmp=vm; plot_cube, title('v from model output')
96    print -djpeg v_exf
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22