runomp2.m

用MATLAB进行海洋水团的模拟程序

6.883 6.981 6.979 6.981 6.955 6.825 6.755 6.620 6.450 5.979 5.640 5.389 5.057 ...
4.793 4.548 4.409 4.359 4.377 4.428 4.513 4.594 4.674 4.748 4.831 4.908 4.938];
atlsal(27,:) = atl;
atl = [ ...
6.966 7.035 7.036 7.035 7.006 6.887 6.817 6.691 6.527 6.026 5.647 5.373 5.020 ...
4.754 4.517 4.395 4.355 4.382 4.442 4.530 4.614 4.699 4.774 4.848 4.920 4.942];
atlsal(28,:) = atl;
atl = [ ...
7.029 7.072 7.076 7.073 7.045 6.934 6.864 6.748 6.591 6.065 5.649 5.355 4.981 ...
4.716 4.487 4.382 4.353 4.388 4.456 4.547 4.633 4.723 4.798 4.861 4.927 4.944];
atlsal(29,:) = atl;
atl = [ ...
7.074 7.095 7.101 7.099 7.072 6.969 6.898 6.793 6.639 6.095 5.646 5.335 4.943 ...
4.679 4.461 4.374 4.355 4.395 4.469 4.562 4.650 4.744 4.818 4.871 4.931 4.942];
atlsal(30,:) = atl;
atl = [ ...
7.099 7.106 7.114 7.114 7.090 6.995 6.922 6.823 6.669 6.112 5.636 5.312 4.906 ...
4.643 4.441 4.369 4.360 4.404 4.481 4.574 4.665 4.760 4.833 4.878 4.932 4.940];
atlsal(31,:) = atl;
atl = [ ...
7.105 7.104 7.114 7.117 7.098 7.012 6.937 6.838 6.678 6.111 5.616 5.283 4.871 ...
4.611 4.427 4.369 4.368 4.413 4.491 4.585 4.677 4.771 4.842 4.882 4.931 4.936];
atlsal(32,:) = atl;
atl = [ ...
7.088 7.085 7.097 7.104 7.090 7.020 6.943 6.839 6.666 6.091 5.583 5.246 4.836 ...
4.582 4.421 4.373 4.377 4.424 4.501 4.593 4.685 4.777 4.846 4.885 4.928 4.933];
atlsal(33,:) = atl;
atl = [ ...
7.047 7.044 7.059 7.069 7.063 7.015 6.941 6.825 6.634 6.047 5.535 5.202 4.802 ...
4.560 4.421 4.382 4.389 4.435 4.509 4.599 4.691 4.779 4.848 4.888 4.924 4.930];
atlsal(34,:) = atl;
atl = [ ...
6.979 6.980 6.995 7.008 7.012 6.995 6.932 6.799 6.582 5.981 5.474 5.150 4.768 ...
4.545 4.428 4.394 4.401 4.446 4.517 4.603 4.694 4.778 4.848 4.890 4.921 4.930];
atlsal(35,:) = atl;
atl = [ ...
6.884 6.889 6.903 6.918 6.934 6.955 6.915 6.761 6.514 5.892 5.400 5.093 4.737 ...
4.537 4.441 4.409 4.414 4.457 4.524 4.607 4.696 4.777 4.849 4.894 4.919 4.931];
atlsal(36,:) = atl;
atl = [ ...
6.768 6.777 6.788 6.803 6.831 6.893 6.888 6.711 6.428 5.784 5.316 5.032 4.709 ...
4.536 4.456 4.425 4.427 4.467 4.530 4.613 4.699 4.777 4.853 4.899 4.922 4.934];
atlsal(37,:) = atl;
atl = [ ...
6.635 6.647 6.655 6.668 6.706 6.809 6.849 6.646 6.326 5.665 5.229 4.972 4.686 ...
4.541 4.472 4.441 4.439 4.476 4.537 4.620 4.705 4.782 4.860 4.907 4.927 4.938];
atlsal(38,:) = atl;
atl = [ ...
6.495 6.507 6.513 6.525 6.567 6.703 6.790 6.559 6.205 5.540 5.143 4.918 4.672 ...
4.550 4.486 4.454 4.451 4.484 4.544 4.629 4.715 4.790 4.868 4.914 4.935 4.942];
atlsal(39,:) = atl;
atl = [ ...
6.354 6.364 6.371 6.382 6.425 6.581 6.703 6.444 6.066 5.419 5.066 4.874 4.666 ...
4.561 4.497 4.464 4.461 4.492 4.550 4.639 4.727 4.802 4.878 4.921 4.942 4.947];
atlsal(40,:) = atl;
atl = [ ...
6.218 6.227 6.236 6.249 6.291 6.450 6.584 6.298 5.912 5.308 5.005 4.844 4.668 ...
4.571 4.504 4.470 4.469 4.498 4.557 4.648 4.740 4.814 4.886 4.926 4.949 4.950];
atlsal(41,:) = atl;
atl = [ ...
6.097 6.105 6.119 6.136 6.177 6.320 6.435 6.128 5.754 5.220 4.968 4.832 4.677 ...
4.579 4.506 4.472 4.474 4.504 4.563 4.657 4.751 4.826 4.894 4.929 4.954 4.954];
atlsal(42,:) = atl;
atl = [ ...
5.999 6.007 6.027 6.049 6.092 6.205 6.266 5.948 5.607 5.160 4.957 4.838 4.691 ...
4.585 4.506 4.472 4.477 4.510 4.570 4.665 4.760 4.836 4.899 4.932 4.957 4.957];
atlsal(43,:) = atl;
atl = [ ...
5.929 5.940 5.966 5.994 6.041 6.114 6.099 5.782 5.486 5.130 4.970 4.862 4.708 ...
4.588 4.504 4.472 4.479 4.517 4.578 4.673 4.767 4.843 4.903 4.934 4.959 4.959];
atlsal(44,:) = atl;
atl = [ ...
5.889 5.905 5.934 5.968 6.019 6.053 5.957 5.651 5.401 5.128 4.999 4.894 4.722 ...
4.589 4.503 4.473 4.482 4.525 4.588 4.683 4.774 4.849 4.906 4.935 4.960 4.962];
atlsal(45,:) = atl;
atl = [ ...
5.873 5.893 5.926 5.962 6.018 6.022 5.855 5.566 5.355 5.145 5.035 4.928 4.733 ...
4.589 4.504 4.476 4.487 4.534 4.598 4.695 4.782 4.855 4.909 4.937 4.962 4.963];
atlsal(46,:) = atl;
atl = [ ...
5.871 5.892 5.925 5.964 6.024 6.010 5.796 5.524 5.342 5.171 5.071 4.958 4.739 ...
4.588 4.506 4.481 4.493 4.542 4.609 4.706 4.790 4.859 4.911 4.939 4.964 4.964];
atlsal(47,:) = atl;
atl = [ ...
5.868 5.888 5.920 5.959 6.024 6.004 5.769 5.513 5.348 5.199 5.100 4.980 4.743 ...
4.589 4.511 4.488 4.501 4.549 4.618 4.716 4.797 4.864 4.913 4.941 4.966 4.964];
atlsal(48,:) = atl;
atl = [ ...
5.850 5.870 5.899 5.939 6.012 5.993 5.759 5.517 5.364 5.223 5.123 4.997 4.748 ...
4.593 4.516 4.497 4.509 4.555 4.626 4.723 4.804 4.868 4.915 4.943 4.966 4.963];
atlsal(49,:) = atl;
atl = [ ...
5.750 5.774 5.801 5.853 5.951 5.938 5.747 5.530 5.387 5.246 5.142 5.014 4.767 ...
4.611 4.533 4.514 4.524 4.564 4.635 4.728 4.810 4.872 4.915 4.943 4.964 4.959];
atlsal(50,:) = atl;
atl = [ ...
5.679 5.708 5.737 5.800 5.915 5.903 5.738 5.529 5.387 5.242 5.138 5.015 4.778 ...
4.625 4.544 4.523 4.531 4.568 4.638 4.728 4.811 4.873 4.914 4.942 4.963 4.957];
atlsal(51,:) = atl;
atl = [ ...
5.616 5.650 5.683 5.759 5.886 5.873 5.728 5.520 5.377 5.229 5.123 5.008 4.789 ...
4.642 4.557 4.534 4.539 4.572 4.639 4.726 4.810 4.872 4.912 4.940 4.963 4.955];
atlsal(52,:) = atl;
atl = [ ...
5.579 5.615 5.654 5.741 5.871 5.850 5.715 5.503 5.356 5.205 5.100 4.993 4.797 ...
4.662 4.574 4.546 4.547 4.578 4.642 4.726 4.810 4.872 4.911 4.939 4.963 4.953];
atlsal(53,:) = atl;
atl = [ ...
5.580 5.615 5.659 5.749 5.868 5.831 5.693 5.474 5.325 5.173 5.069 4.974 4.806 ...
4.687 4.596 4.562 4.560 4.588 4.647 4.729 4.813 4.874 4.912 4.939 4.964 4.952];
atlsal(54,:) = atl;
atlsal=30 + atlsal';
figure('Color',[0.2 0.6 0.8],'MenuBar','none', ...
	'Name','OMP Analysis version 2.0', ...
	'NumberTitle','off','Position',[10 20 600 420],'Resize','off');

axes('Position',[0.1 0.53 0.1 0.2],'color',[1 1 1],'Box','on', ...
	'units','normalized','ylim',[24.5 27],'fontname','helvetica', ...
	'fontsize',fts,'xcol','w','ycol','w');
hold on
plot(100*err,dens,'k.')
axis ij
title('mass conservation','color','w')
xlabel('residual (%)','Color','w','fontname','helvetica','fontsize',fts)
ylabel('sigma-t','Color','w','fontname','helvetica','fontsize',fts)
hold off

axes('Position',[0.1 0.23 0.27 0.2],'Color',[0 0 0],'Box','on', ...
	'Units','normalized','xlim',[min(atllat) max(atllat)],'ylim', ...
	[min(dept) max(dept)], ...
	'xColor','w','yColor','w','fontname','helvetica','fontsize',fts)
hold on
pcolor(atllat,dept,atlsal)
axis ij
colormap(hsv)
shading flat
caxis([34 36.9])
xlabel('latitude (deg.)','fontname','helvetica','fontsize',fts)
ylabel('depth (km)','fontname','helvetica','fontsize',fts,'color','w')
text(-35,2.8,'Atlantic Ocean: salinity','fontname','helvetica','fontsize',fts,'color','w')
hold off

axes('Position',[0.27 0.53 0.12 0.2],'Color',[0 0 0],'Box','on', ...
   'TickDir','out','Units','normalized','xlim',[34.2 35.6],'ylim', ...
     [0 30],'xColor','w','yColor','w', ...
	'fontname','helvetica','fontsize',fts)
hold on
plot(sal,temp,'.r')
title('Indian Ocean','color','w')
ylabel('temperature','fontname','helvetica','fontsize',fts,'color','w')
xlabel('salinity','fontname','helvetica','fontsize',fts,'color','w')
hold off
set(gca,'fontname','helvetica','fontsize',fts,'color','w')
for i = 1:4
tx2 = uicontrol('Units','normalized', ...
	'BackgroundColor',[0.2 0.2 0.2], ...
	'ForegroundColor',[1 1 1], ...
	'Position',[0.0 1-0.035*i 1.0 0.05], ...
	'HorizontalAlignment','center', ...
	'String',ht3(i,:), ...
	'Style','text','fontsize',fts+1);
end

graphmode = uicontrol('Units','normalized', ...
	'BackgroundColor',[1 0.8 0.5], ...
	'Callback',[	'clear atllat atlsal dens dept err sal temp,' ...
	'omp2gui'], ...
	'Position',[0.23 0.05 0.2 0.066], ...
	'String','interactive (GUI)', ...
	'Style','pushbutton');
	
intermode = uicontrol('Units','normalized', ...
	'BackgroundColor',[1 0.8 0.5], ...
	'Callback',['clear atllat atlsal dens dept err sal temp,' ...
	'omp2int'], ...
	'Position',[0.02 0.05 0.2 0.066], ...
	'String','interactive (listing)', ...
	'Style','pushbutton');

automode = uicontrol('Units','normalized', ...
	'BackgroundColor',[1 0.8 0.5], ...
	'Callback',[ ...
	'clear atllat atlsal dens dept err sal temp,' ...
	'set(intermode,''Visible'',''off''), set(automode,''Visible'',''off''),' ...
	'set(infile,''Visible'',''on''), set(graphmode,''Visible'',''off''),' ...
	'set(whichfile,''Visible'',''on'');'], ...
	'Position',[0.44 0.05 0.2 0.066], ...
	'String','automatic (file)', ...
	'Style','pushbutton');

ax = axes('position',[0.455 0.12 1 0.71],'color',[0.85 0.95 0.95], ...
	'xcolor',[0.7 0.83 0.8],'ycolor',[0.7 0.83 0.8],'xtick',[],'ytick',[]);
hold on
for i = 1:9
	text(0.02,1-0.05*i,ht1(i,:),'color','k','fontname',eval(ftm),'fontsize',ftz);
end

for i = 1:6
	text(0.02,0.53-0.05*i,ht2(i,:),'color','k','fontname',eval(ftm),'fontsize',ftz);
end

ex = uicontrol('Units','normalized', ...
	'BackgroundColor',[0.9 0.9 0.2], ...
	'Callback',['clear all, close all'], ...
	'Position',[0.9 0.05 0.066 0.066], ...
	'String','exit', ...
	'Style','pushbutton');
	
hlp = uicontrol('Units','normalized', ...
	'BackgroundColor',[0.9 0.9 0.2], ...
	'Callback','helpomp', ...
	'Position',[0.82 0.05 0.066 0.066], ...
	'String','help', ...
	'Style','pushbutton');
	
cr = uicontrol('Units','normalized', ...
	'BackgroundColor',[1 0.4 0.2], ...
	'Callback',['crf = figure(''Color'',[0.5 0.6 0.8],''MenuBar'',''none'','...
		'''Name'',''OMP Analysis: Credits'',' ...
		'''Position'',[10 20 600 420],''NumberTitle'',''off'',''Resize'',''off'');' ...
		'displtxt(credit);', ...
		'pause(5),close(crf),' ...
		], ...
	'Position',[0.69 0.05 0.116 0.066], ...
	'String','credits', ...
	'Style','pushbutton');

infile = uicontrol('Units','normalized', ...
	'BackgroundColor',[1 1 1], ...
	'Callback',[ ...
	'newfile = get(infile,''String'');,' ...
	'if exist(newfile) == 0,' ...
	'set(infile,''String'','' file not found''),' ...
	'else,' ...
		'set(infile,''Visible'',''off'');,' ...
		'set(whichfile,''String'',''Reading parameters'');,' ...
		'omp2auto;' ...
	'end'], ...
	'Position',[0.23 0.05 0.21 0.05], ...
	'String','incontr2', ...
	'Style','edit');

whichfile = uicontrol('Units','normalized', ...
	'BackgroundColor',[1 0.9 0.9], ...
	'Callback','omp2int', ...
	'Position',[0.02 0.05 0.18 0.05], ...
	'String','Enter control file:', ...
	'Style','text');
	
set(infile,'Visible','off');
set(whichfile,'Visible','off');