/[MITgcm]/MITgcm/verification/global_ocean.90x40x15/diags_matlab/mit_plotzonalvelocity.m
ViewVC logotype

Annotation of /MITgcm/verification/global_ocean.90x40x15/diags_matlab/mit_plotzonalvelocity.m

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

Revision 1.1 - (hide annotations) (download)
Tue Oct 22 13:30:40 2002 UTC (21 years, 11 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint57t_post, checkpoint58l_post, checkpoint48f_post, checkpoint46k_post, checkpoint51k_post, checkpoint53f_post, checkpoint47j_post, checkpoint54a_pre, checkpoint55c_post, checkpoint53b_pre, checkpoint48d_pre, checkpoint51l_post, checkpoint51j_post, branch-exfmods-tag, checkpoint47e_post, checkpoint57m_post, checkpoint47i_post, checkpoint52l_pre, checkpoint48i_post, checkpoint52e_pre, checkpoint58e_post, checkpoint57v_post, checkpoint57g_pre, checkpoint52j_post, checkpoint47f_post, checkpoint48d_post, checkpoint51o_pre, checkpoint57f_post, checkpoint50e_post, checkpoint52e_post, checkpoint50c_post, checkpoint57s_post, checkpoint51n_pre, checkpoint47d_post, checkpoint57j_post, checkpoint58b_post, checkpoint58m_post, checkpoint47a_post, checkpoint57b_post, checkpoint46l_pre, checkpoint53c_post, checkpoint53d_post, checkpoint57f_pre, checkpoint48a_post, checkpoint55d_pre, checkpoint46n_post, checkpoint51f_pre, checkpoint57g_post, checkpoint48e_post, checkpoint57a_post, checkpoint48h_post, checkpoint55j_post, checkpoint56b_post, checkpoint50c_pre, checkpoint57h_pre, checkpoint57y_post, branchpoint-genmake2, checkpoint46l_post, checkpoint58g_post, checkpoint57x_post, checkpoint52j_pre, checkpoint54a_post, branch-netcdf, checkpoint50d_pre, checkpoint55h_post, checkpoint51r_post, checkpoint52b_pre, checkpoint52n_post, checkpoint54b_post, checkpoint58h_post, checkpoint51i_post, checkpoint57e_post, checkpoint54d_post, checkpoint47h_post, checkpoint48c_post, checkpoint56c_post, checkpoint54e_post, checkpoint58j_post, checkpoint55b_post, checkpoint51e_post, checkpoint51b_post, checkpoint57h_post, checkpoint51l_pre, checkpoint52m_post, checkpoint51c_post, checkpoint57y_pre, checkpoint55, checkpoint53a_post, checkpoint55a_post, checkpoint57c_pre, checkpoint48, checkpoint49, checkpoint47b_post, checkpoint53b_post, checkpoint57o_post, checkpoint55g_post, checkpoint57r_post, checkpoint51o_post, checkpoint48g_post, checkpoint57k_post, checkpoint57d_post, checkpoint55f_post, checkpoint57i_post, checkpoint51q_post, checkpoint52l_post, checkpoint52k_post, checkpoint58, checkpoint57a_pre, checkpoint54, checkpoint57, checkpoint56, checkpoint51, checkpoint50, checkpoint53, checkpoint52, checkpoint50d_post, checkpoint52d_post, checkpoint46m_post, checkpoint51b_pre, checkpoint52a_post, checkpoint57h_done, checkpoint47g_post, checkpoint58f_post, checkpoint52b_post, checkpoint53g_post, checkpoint52f_post, checkpoint57n_post, checkpoint52c_post, checkpoint58d_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint51h_pre, checkpoint58a_post, checkpoint50g_post, checkpoint58i_post, checkpoint57q_post, checkpoint51g_post, ecco_c52_e35, checkpoint57z_post, checkpoint54f_post, checkpoint51f_post, checkpoint48b_post, checkpoint50b_post, eckpoint57e_pre, checkpoint58c_post, checkpoint58k_post, checkpoint57c_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint52a_pre, checkpoint47d_pre, checkpoint51d_post, checkpoint48c_pre, checkpoint51m_post, checkpoint51t_post, checkpoint53d_pre, checkpoint47, checkpoint55e_post, checkpoint54c_post, checkpoint50h_post, checkpoint52i_post, checkpoint51a_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51p_post, checkpoint51n_post, checkpoint55i_post, checkpoint51i_pre, checkpoint57l_post, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, checkpoint50b_pre, checkpoint56a_post, checkpoint51s_post, checkpoint55d_post
Branch point for: netcdf-sm0, branch-genmake2, branch-nonh, tg2-branch, checkpoint51n_branch, release1, branch-exfmods-curt 
 fixed the verification/global_ocean.90x40x15 experiment:
 - new bathymetry (the world according to A., JMC, and M.)
 - new initial fields and forcing fields (*.bin files)
 - new POLY3.COEFFS (for the next release one should switch to a full
   equation of state: JMD95P or MDJWF)
 - fixed several errors and redundancies in the data file
 - experiment uses looped cells
 - added matlab directory with diagnostic scripts for plotting of output
1 mlosch 1.1 % m-file: mit_plotzonalvelocity.m
2    
3     % select timestep
4     k=kmax;
5    
6     up = u(:,:,:,k).*grd.pacific_hfacw;
7     upzm = mit_zonalmean(up,grd.pacific_hfacw,grd.dxg);
8     ua = u(:,:,:,k).*grd.atlantic_hfacw;
9     uazm = mit_zonalmean(ua,grd.atlantic_hfacw,grd.dxg);
10    
11     caxup = [min(upzm(:)) max(upzm(:))];
12     uplev = -1:.1:1;
13     if max(abs(caxup)) < .1
14     uplev = .5*uplev;
15     end
16     if max(abs(caxup)) < .2
17     uplev = .2*uplev;
18     end
19    
20     caxua = [min(uazm(:)) max(uazm(:))];
21     ualev = -1:.1:1;
22     if max(abs(caxua)) < .1
23     ualev = .5*ualev;
24     end
25     if max(abs(caxua)) < .2
26     ualev = .2*ualev;
27     end
28    
29     ixpw = 41;
30     ixpc = 53;
31     ixpe = 65;
32    
33     ixaw = 83;
34     ixac = 85;
35     ixae = 2;
36    
37     zaxis = -grd.zc;
38     zaxis = -grd.zg;
39     yaxis = grd.latg;
40     %zaxis = -[1:grd.nz];
41    
42     figure('PaperPosition',[0.25 0.621429 8 9.75714])
43     sh = subplot(4,2,1);
44     pcol(yaxis,zaxis,upzm')
45     title('Pacific Ocean: zonal average')
46     sh(3) = subplot(4,2,3);
47     pcol(yaxis,zaxis,sq(up(ixpw,:,:))')
48     title(['section at ' num2str(grd.long(ixpw)) '^{\circ}E'])
49     sh(5) = subplot(4,2,5);
50     pcol(yaxis,zaxis,sq(up(ixpc,:,:))')
51     title(['section at ' num2str(grd.long(ixpc)) '^{\circ}E'])
52     sh(7) = subplot(4,2,7);
53     pcol(yaxis,zaxis,sq(up(ixpe,:,:))')
54     title(['section at ' num2str(grd.long(ixpe)) '^{\circ}E'])
55     sh(2) = subplot(4,2,2);
56     pcol(yaxis,zaxis,uazm',ualev)
57     title('Atlantic Ocean: zonal average')
58     sh(4) = subplot(4,2,4);
59     pcol(yaxis,zaxis,sq(ua(ixaw,:,:))')
60     title(['section at ' num2str(grd.long(ixaw)) '^{\circ}E'])
61     sh(6) = subplot(4,2,6);
62     pcol(yaxis,zaxis,sq(ua(ixac,:,:))')
63     title(['section at ' num2str(grd.long(ixac)) '^{\circ}E'])
64     sh(8) = subplot(4,2,8);
65     pcol(yaxis,zaxis,sq(ua(ixae,:,:))')
66     title(['section at ' num2str(grd.long(ixae)) '^{\circ}E'])
67    
68     set(sh,'xlim',[-1 1]*30,'ylim',[-400 0])
69     set(sh,'layer','top')
70     set(sh(1:2:end),'clim',[uplev(1) uplev(end)])
71     set(sh(2:2:end),'clim',[ualev(1) ualev(end)])
72     set(gcf,'currentAxes',sh(end-1));colorbar('h')
73     set(gcf,'currentAxes',sh(end));colorbar('h')
74    
75     suptitle(['experiment ' grd.dname ...
76     ', timestep = ' num2str(timesteps(k)) ...
77     ', ' tuname ' = ' num2str(tim(k)) ', zonal velocity [m/s]'])
  ViewVC Help
Powered by ViewVC 1.1.22