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

Contents of /MITgcm_contrib/high_res_cube/input/lookat.m

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

Revision 1.1 - (show annotations) (download)
Wed Dec 10 19:40:13 2003 UTC (21 years, 7 months ago) by dimitri
Branch: MAIN
CVS Tags: hrcube4, hrcube5, hrcube_1, hrcube_2, hrcube_3, HEAD 
o changes to permit hi-res, cubed-sphere, configuration
  - added useSingleCpuIO capability to mdsio_readfield.F and exf_interp.F
  - added "#undef ALLOW_USE_MPI" support to eesupp/src/scatter_2d.F
  - added pkg/exf/exf_set_uv.F for on-the-fly interpolation and rotation
    of surface winds for the cube
  (A verification experiment, based on global_ocean.cs32x15, but with
   pkg/seaice turned on and with on-the-fly interpolation from the NCEP
   Gaussian grid is described in MITgcm_contrib/high_res_cube/README_ice,
   complete with example surface forcing files and matlab scripts to look
   at the output and compare it with that of global_ocean.cs32x15.)
1 figure(1),clf
2 subplot(421),mypcolor(readbin('uwind_192_94_8.bin',[192 94])');colorbar
3 subplot(422),mypcolor(readbin('vwind_192_94_8.bin',[192 94])');colorbar
4 subplot(423),mypcolor(readbin('swdown_192_94_8.bin',[192 94])');colorbar
5 subplot(424),mypcolor(readbin('lwdown_192_94_8.bin',[192 94])');colorbar
6 subplot(425),mypcolor(readbin('precip_192_94_8.bin',[192 94])');colorbar
7 subplot(426),mypcolor(readbin('atemp_192_94_8.bin',[192 94])');colorbar
8 subplot(427),mypcolor(readbin('aqh_192_94_8.bin',[192 94])');colorbar
9
10 figure(2)
11 for fld={'RAZ','RAW','RAS','RA','LONG','LONC','LATG','LATC', ...
12 'DYU','DYG','DYF','DYC','DXV','DXG','DXF','DXC'}
13 tmp=permute(readbin([fld{1} '.bin'],[32,6,32],1,'real*8'),[1 3 2]);
14 cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
15 end
16
17 for fld={'UWIND','VWIND','Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
18 tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);
19 cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
20 end
  ViewVC Help
Powered by ViewVC 1.1.22