論文を探す

キャンパス内のネットワーク上でGoogle Scholarなどを使って論文を検索すると大学のアカウントを使ってログインが自動的に行われるため論文を入手することが容易にできます。

でもキャンパス外だと大学の個人情報（ネットリンクIDとパスワード）を入力しないと論文を手に入れることができません。

大学に最初来た時は大学図書館のサイトのホーム画面（http://www.uvic.ca/library/index.php）から検索していたんですがたまに欲しい論文が検索結果に出て来ない場合があります。

そこで図書館の人に訪ねてみると、こちらのデータベース（http://www.uvic.ca/library/find/databases/index.php）から検索するように言われやってみたら全ての論文がでてきました☆自分はGoogle Scholarを使うのでこちらのリンクから大学経由でGoogle Scholarを使うことができます。

UVic学生で論文検索をされる方、要ブックマーク→http://scholar.google.ca.ezproxy.library.uvic.ca/

Plotting for poster.

Figure used for poster

%This is for FABM model output (not GOTM)!

clear all;

filename='resolute.nc';
dep=100;
avgdep=30;
day1=366;
colid='bcgmrk';

figure;

subplot(4,1,1);
n1=squeeze(ncread(filename,'sst'));
n2=squeeze(ncread(filename,'ice_hi'))*100;
xlen=day1:length(n1);
x=linspace(datenum('01-01-1990'),datenum('01-01-1994'),length(xlen));
[ax,h1,h2]=plotyy(x,n1(xlen),x,n2(xlen));
set(h1,'linewidth',2,'color','r');
set(h2,'linewidth',2,'color','b');
set(ax(1),'YLim',[-2 4],'ycolor','r')
set(ax(1),'YTick',-2:1:4)
set(ax(2),'YLim',[0 60],'ycolor','b')
set(ax(2),'YTick',0:10:60)
set(get(ax(1),'Ylabel'),'string','^oC')
set(get(ax(2),'Ylabel'),'string','cm')
axes(ax(2));
title('(a)');
axes(ax(1));
hold on;
datetick('x','mmmyy');
axes(ax(2));
hold on;
datetick('x','mmmyy');
%set(ax(1),'xtick',day1:365:length(xlen),'xlim',[day1 length(xlen)]);set(ax(2),'xtick',day1:365:length(xlen),'xlim',[day1 length(xlen)]);
grid on;
%xlabel('Days since January 1, 1990');
%legend('SST','Ice Thickness');

subplot(4,1,2);
parname={'uvic_eco_ph1' 'uvic_eco_ph2' 'uvic_eco_zo1' 'uvic_eco_zo2'};
n=squeeze(ncread(filename,'uvic_eco_ph1'));
n1=squeeze(mean(n(end-avgdep-1:end,xlen),1))*14/8.75;
n=squeeze(ncread(filename,'uvic_eco_zo1'));
n2=squeeze(mean(n(end-avgdep-1:end,xlen),1))*106/16*12;
[ax,h1,h2]=plotyy(x,n1,x,n2);
set(ax(1),'ycolor','r');
set(ax(2),'ycolor','b');
set(get(ax(1),'Ylabel'),'string','mg-Chl-\alpha/m^3');
set(get(ax(2),'Ylabel'),'string','mg-C/m^3');
set(h1,'linewidth',2,'color','r');
set(h2,'linewidth',2,'color','b');
set(ax(1),'YLim',[0 7],'ycolor','r')
set(ax(1),'YTick',0:1:7)
set(ax(2),'YLim',[0 112],'ycolor','b')
set(ax(2),'YTick',0:16:112)
%xlabel('Days since January 1, 1990');
n=squeeze(ncread(filename,'uvic_eco_ph2'));
n1=squeeze(mean(n(end-avgdep-1:end,xlen),1))*14/8.75;
axes(ax(1));
hold on;
hh=plot(x,n1,'--','color','r');
set(hh,'linewidth',2);
%set(gca,'xtick',1:365:length(n1),'xlim',[1 length(n1)]);
datetick('x','mmmyy');
n=squeeze(ncread(filename,'uvic_eco_zo2'));
n2=squeeze(mean(n(end-avgdep-1:end,xlen),1))*106/16*12;
axes(ax(2));
hold on;
hh=plot(x,n2,'--','color','b');
set(hh,'linewidth',2);
%set(gca,'xtick',1:365:length(n1),'xlim',[1 length(n1)]);
%set(ax(1),'xtick',1:365:length(n1),'xlim',[1 length(n1)]);
%set(ax(2),'xtick',1:365:length(n1),'xlim',[1 length(n1)]);
datetick('x','mmmyy');
grid on;
title('(b)');
%legend('Flagellates','Diatoms','Microzooplankton','Mesozooplankton','orientation','horizontal');

subplot(4,1,3);
n=squeeze(ncread(filename,'uvic_dms_dms'));
[ax, h]=contour(x,-dep:-1,n(1:dep,xlen));
set(h,'linewidth',2);
colorbar('location','west');
ylim([-50 -1]);
title('(c)');
%xlabel('Days since January 1, 1990');
datetick('x','mmmyy');
ylabel('Depth (m)');
%set(gca,'xtick',1:365:length(n1),'xlim',[1 length(n1)]);
grid on;

subplot(4,1,4);
n=squeeze(ncread(filename,'uvic_dms_dmspd'));
[ax, h]=contour(x,-dep:-1,n(1:dep,xlen));
set(h,'linewidth',2);
colorbar('location','west');
ylim([-50 -1]);
title('(d)');
%xlabel('Days since January 1, 1990');
datetick('x','mmmyy');
ylabel('Depth (m)');
%set(gca,'xtick',1:365:length(n1),'xlim',[1 length(n1)]);
grid on;

家賃値上げ

帰宅するとドアのノブのうえに封筒が。光熱費の請求かと思いきやいつもより分厚い。

開けてみると賃貸契約の更新＆九月からの家賃値上げのお知らせだった。

今までは月々$850。これが九月からは$18.70値上がりの$868.70になるとのこと。

こんなん勝手に値上げできるんかと思って調べてみたら州の法律に従っての値上げのようで２０１４年度の上限は2.2%。詳細→http://www.rto.gov.bc.ca/content/news/default.aspx

$850の2.2%は$18.70…

今の家賃での契約は８月までなので９月からの家を探してみるのもいいかもな。

Initial conditions and model behaviour

Initial conditions have effect on the following structure of interannual cycle. For example, the result seems good if I set really high initial zo1 and zo2 values (here, 60 mg-C/m3 for January 1st):

But I am not sure why the initial values drop down so quickly after January 1 until the first bloom. If it remains constant from January 1st, it's more reasonable...

Perhaps setting a background value might solve the problem.

photoinhibition

PI curve with and without photoinhibition