2009. 3. 27 (금)
휴대폰 구입 LG-SV390 와인폰2
2009년 3월 28일 토요일
FORTRAN julian day
subroutine juld(iy,im,id,ih,mn,jday)
integer iy,im,id,ih,mn
real day(12),jday
data day/31.,28.,31.,30.,31.,30.,31.,31.,30.,31.,30.,31./
if ( mod(iy,4).eq.0 .and. (mod(iy,100).ne.0
& .or. mod(iy,400).eq.0) ) then
day(2)=29.
endif
jday=0.
if (im-1.eq.0) then
jday=0.
else !if (im-1.eq.0) then
do i=1,im-1
jday=jday+day(i)
enddo
endif !if (im-1.eq.0) then
jday=jday+float(id) + (float(ih)/24.) + (float(mn)/1440.)
return
end
integer iy,im,id,ih,mn
real day(12),jday
data day/31.,28.,31.,30.,31.,30.,31.,31.,30.,31.,30.,31./
if ( mod(iy,4).eq.0 .and. (mod(iy,100).ne.0
& .or. mod(iy,400).eq.0) ) then
day(2)=29.
endif
jday=0.
if (im-1.eq.0) then
jday=0.
else !if (im-1.eq.0) then
do i=1,im-1
jday=jday+day(i)
enddo
endif !if (im-1.eq.0) then
jday=jday+float(id) + (float(ih)/24.) + (float(mn)/1440.)
return
end
2009년 3월 11일 수요일
MATLAB tip
matlab tip
2007.03.18. 01:21
figure창 배경색 조절
set(gcf,'color',[1 1 1]); % 배경색을 흰색으로
set(gcf,'color','w'); % 배경색을 흰색으로
stick diagram을 그리면 상단과 오른쪽에 박스가 안그려진다. 이것을 그릴려면
set(gca,'box','on');
figure창 위치 및 크기 조절
figure(1);
set( gcf , 'Position' , [1500 504 560 420] );
활용. 동영상 파일을 만들경우 다른 창이 동영상 작업중인 것 위에 올라가면 동영상에 손상을 준다. 그런경우 창을 완전히 다른곳으로 옮겨 버리면 이런 실수를 막을 수 있다.
1500; 좌측하단 모서리의 x좌표(단위; pixel)
504; 좌측하단 모서리의 y좌표(단위; pixel)
560; 좌측하단 기준 x방향의 길이(단위; pixel)
420; 좌측하단 기준 y방향의 길이(단위; pixel)
선 두껍게 그리기
LinHd(1)=line(temp,-1*depth, 'LineWidth' , 2.0);
plot(x,y,'--rs', 'LineWidth',2,...
'MarkerEdgeColor','k',...
'MarkerFaceColor','g',...
'MarkerSize',10)
set( gca , 'linewidth' , 2 ); % 이경우에는 x,y 축의 선이 두꺼워진다.
Temperature(℃) 표시
‘Temperature(\circC)’
여러 개의 파일을 load하는 방법
files={'05.txt','09.txt','13.txt'}
d=load( char(files(i)) );
grid line을 원하는 곳에 표시하는 방법
grid on;
v = [-0.5 -0.33 0 0.5 4];
set(gca, ‘ytick’, v);
이렇게 하면 y축에 지정된 위치만 grid line이 적용된다.
x도 마찬가지 방법으로 하면 된다.
title에 return문자 넣기
sprintf('\n')
title( ['time=',num2str(time*24*60),'min' ,sprintf('\n'), 'internal wave generator'] );
title이 두줄로 써진다.
숫자 자릿수 지정
num2str(42.67, '%10.2f')
시간 단위로 변화는 동영상 만들 때 단위 자릿수 맞출 때 사용
Scale bar 단위
%%%%%%%%%%%%%%%%%%% Color bar에 단위 넣기 %%%%%%%%%%%%%%%%%%%%%%
h = colorbar;
zlabel = get(h,'YTickLabel');
mm = repmat(' m',length(zlabel),1);
zlabel = [zlabel mm];
set(h, 'YTickLabel' ,zlabel);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Contour 그릴 때 color bar에 단위 넣기 위의 예제는 ‘m’ 단위
동영상 빠르게 재생하게 만들기
movie2avi(K,'temp_1d_tide.avi','compression','Cinepak','fps',30);
fps의 값은 기본 15, 이 값을 크게 하면 된다.
y축 제목
set(gca,'Ydir','reverse');
문자처리하기
원하는 부분 자르기
ofile='/mmrl1/shyou/RTSMVRFY/DAOU/48HR_INT_CORR_TS/BUSN.12H';
% 여기서 BUSN.12H 만 따로 뽑기
name=ofile(end-7:end);
문자열 붙이기
a1='aaa';a2='bbb';
a3=[a1,a2];
x축 또는 y축 원하는 위치에 tick 쓰기
set(gca,'xtick',[1 5 8 11 14 17 20 23]);set(gca,'XTickLabel',{'05/Jun.','05/Oct.','06/Jan.','06/Apr.','06/Jul.','06/Oct.'... ,'07/Jan.','07/Apr.'});
y축 값 증가 거꾸로 하기
y축 0을 중심으로 위쪽이 증가하는 방향으로 그림이 그려진다.
하지만 수심과 같은 경우를 표현할 때는 상부에서 하부로 향하면서 증가하는 숫자가 필요할 경우가 있다.
이때
view(0,-90);
print 문 활용
hh=figure(1);
그림 그리고
print(hh, '-djpeg60', [num2str(i),'.jpg']);close(hh);
위에서 '-djpeg60' 은 파일 형식을 나타낸다.
만약 imagemagik의 convert 를 활용할 수 있다면 ps 파일 형식으로 저장할 것을 추천한다.
-dpsc ; PostScript for color printers
2007.03.18. 01:21
figure창 배경색 조절
set(gcf,'color',[1 1 1]); % 배경색을 흰색으로
set(gcf,'color','w'); % 배경색을 흰색으로
stick diagram을 그리면 상단과 오른쪽에 박스가 안그려진다. 이것을 그릴려면
set(gca,'box','on');
figure창 위치 및 크기 조절
figure(1);
set( gcf , 'Position' , [1500 504 560 420] );
활용. 동영상 파일을 만들경우 다른 창이 동영상 작업중인 것 위에 올라가면 동영상에 손상을 준다. 그런경우 창을 완전히 다른곳으로 옮겨 버리면 이런 실수를 막을 수 있다.
1500; 좌측하단 모서리의 x좌표(단위; pixel)
504; 좌측하단 모서리의 y좌표(단위; pixel)
560; 좌측하단 기준 x방향의 길이(단위; pixel)
420; 좌측하단 기준 y방향의 길이(단위; pixel)
선 두껍게 그리기
LinHd(1)=line(temp,-1*depth, 'LineWidth' , 2.0);
plot(x,y,'--rs', 'LineWidth',2,...
'MarkerEdgeColor','k',...
'MarkerFaceColor','g',...
'MarkerSize',10)
set( gca , 'linewidth' , 2 ); % 이경우에는 x,y 축의 선이 두꺼워진다.
Temperature(℃) 표시
‘Temperature(\circC)’
여러 개의 파일을 load하는 방법
files={'05.txt','09.txt','13.txt'}
d=load( char(files(i)) );
grid line을 원하는 곳에 표시하는 방법
grid on;
v = [-0.5 -0.33 0 0.5 4];
set(gca, ‘ytick’, v);
이렇게 하면 y축에 지정된 위치만 grid line이 적용된다.
x도 마찬가지 방법으로 하면 된다.
title에 return문자 넣기
sprintf('\n')
title( ['time=',num2str(time*24*60),'min' ,sprintf('\n'), 'internal wave generator'] );
title이 두줄로 써진다.
숫자 자릿수 지정
num2str(42.67, '%10.2f')
시간 단위로 변화는 동영상 만들 때 단위 자릿수 맞출 때 사용
Scale bar 단위
%%%%%%%%%%%%%%%%%%% Color bar에 단위 넣기 %%%%%%%%%%%%%%%%%%%%%%
h = colorbar;
zlabel = get(h,'YTickLabel');
mm = repmat(' m',length(zlabel),1);
zlabel = [zlabel mm];
set(h, 'YTickLabel' ,zlabel);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Contour 그릴 때 color bar에 단위 넣기 위의 예제는 ‘m’ 단위
동영상 빠르게 재생하게 만들기
movie2avi(K,'temp_1d_tide.avi','compression','Cinepak','fps',30);
fps의 값은 기본 15, 이 값을 크게 하면 된다.
y축 제목
set(gca,'Ydir','reverse');
문자처리하기
원하는 부분 자르기
ofile='/mmrl1/shyou/RTSMVRFY/DAOU/48HR_INT_CORR_TS/BUSN.12H';
% 여기서 BUSN.12H 만 따로 뽑기
name=ofile(end-7:end);
문자열 붙이기
a1='aaa';a2='bbb';
a3=[a1,a2];
x축 또는 y축 원하는 위치에 tick 쓰기
set(gca,'xtick',[1 5 8 11 14 17 20 23]);set(gca,'XTickLabel',{'05/Jun.','05/Oct.','06/Jan.','06/Apr.','06/Jul.','06/Oct.'... ,'07/Jan.','07/Apr.'});
y축 값 증가 거꾸로 하기
y축 0을 중심으로 위쪽이 증가하는 방향으로 그림이 그려진다.
하지만 수심과 같은 경우를 표현할 때는 상부에서 하부로 향하면서 증가하는 숫자가 필요할 경우가 있다.
이때
view(0,-90);
print 문 활용
hh=figure(1);
그림 그리고
print(hh, '-djpeg60', [num2str(i),'.jpg']);close(hh);
위에서 '-djpeg60' 은 파일 형식을 나타낸다.
만약 imagemagik의 convert 를 활용할 수 있다면 ps 파일 형식으로 저장할 것을 추천한다.
-dpsc ; PostScript for color printers
피드 구독하기:
글 (Atom)