%-- 12/08/09  3:06 PM --%
plot([2:5],'*')
%-- 12/08/09  4:21 PM --%
%-- 12/8/09  1:56 PM --%
dir
%-- 12/9/09  9:34 AM --%
dir
edit otps2frc_v3
ls
call_otps2frc
startup
~/matlab/startup
/home/guest/dafu8991/matlab/startup
/home/guest/dafu8991/matlab/startup.m
addpath /home/guest/dafu8991/matlab/
addpath /home/beach/guest/dafu8991/matlab/
startup
call_otps2frc
clear
call_otps2frc
ls
cd OTPS
call_otps2frc
clear
call_otps2frc
clear
call_otps2frc
caclearll_otps2frc
clear
caclearll_otps2frc
call_otps2frc
%-- 12/9/09  2:19 PM --%
pwd
call_otps2frc
nc_dump('tide_interp2.nc')
edit fixtideinterp2
help nc_varput
fixtideinterp2
cd ../
cd INIT
ls
cd CALUS/FORCING/
cd INIT/
ls
fname = 'swf_init.nc'
nc_varput(fname,'hc',3);
nc_varput(fname,'hc',.6096);
%-- 12/9/09  4:58 PM --%
ls
edit create_roms_init
%-- 12/10/09 10:31 AM --%
%-- 4/12/10  7:57 PM --%
%-- 4/13/10 11:57 AM --%
pwd
ls
edit gen_transpathsSWF
clear;gen_transpathsSWF
fname
cd ~/../bevera
cd ~/CALUS/analout
editexit
%-- 4/14/10  8:23 AM --%
ls *.m
edit gen_transpathsSWF
clear;gen_transpathsSWF
edit depths
edit depths2
nc_dump('../FRC2008/swf_grid6.nc')
clear;gen_transpathsSWF
ls
edit make_movie_swf
clear;make_movie_swf
i
edit depths2
%-- 4/15/10 10:52 AM --%
pwd
cd CALUS/analout
ls
clear;make_movie_swf
edit depths2
fname
edit make_movie_swf
tmphis
ii
nc_info(tmphis)
clear;make_movie_swf
ii
t
size(transport_timee)
size(transport_time)
clear;ls *.mat
load sed_trans_paths_wet
figure;pcolor(Ut(100,:,:));shading flat ;colorbar
Ut(100,:,:)
Ut(100,100:125:,100:125)
Ut(100,100:125,100:125)
Ut(200,100:125,100:125)
Ut(210,100:125,100:125)
Vt(210,100:125,100:125)
dir *.m
edit gen_transpathsSWF
clear;gen_transpathsSWF
transport time
transport_time
clear;make_movie_swf
%-- 4/15/10  8:57 PM --%
cd CALUS/analout
ls *.m
edit make_movie_swf
clear;make_movie_swf
hhelp addpath
help addpath
clear;make_movie_swf
which psliceuv
which pslice2
edit psliceuv
clear;make_movie_swf
%-- 4/27/10 11:42 AM --%
ls
edit make_movie_swf
cd ..
mkdir OUTtest
cd OUTstable
ls
cp calus2008_his_003* ../OUTtest
clear;make_movie_swf
jj
nfiles
get(0,'screensize')
nxx
nyy
get(gcf,'screensize')
size(t)
ls
chkt=nc_varget([mdir,MncdfF(3).name],'ocean_time')
size(chekt)
size(chkt)
chkt=nc_varget([mdir,MncdfF(5).name],'ocean_time')
chkt=nc_varget([mdir,MncdfF(1).name],'ocean_time')
size(chkt)
chkt=nc_varget([mdir,
chkt=[mdir,MncdfF(1).name]
[mdir,MncdfF(1).name]
[mdir,MncdfF(2).name]
[mdir,MncdfF(3).name]
[mdir,MncdfF(4).name]
[mdir,MncdfF(5).name]
clear;make_movie_swf
%-- 4/29/10 10:03 AM --%
edit make_movie_swf
clear;make_movie_swf
%-- 5/3/10  1:02 PM --%
dir *.m
clear; make_movie_sef
clear; make_movie_swf
edit make_movie_swf
clear; make_movie_swf
pwd
%-- 5/4/10 11:42 AM --%
cd CALUS/analout
cd ../FORCING/INIT
edit create_roms_init
ls
cd ../../FRC2008
edit create_roms_init
clear;create_roms_init
%-- 5/4/10  8:39 PM --%
ls *.m
edit make_avi_swf
clear;make_avi_swf
%-- 5/5/10  7:05 AM --%
clear;make_avi_swf
mov=close(mov)
cd ../FRC2008
ls
nc_dump('swf_grid6.nc')
checkd=nc_varget('swf_grid6.nc','depthmin')
checkd=nc_varget('swf_grid6.nc','depthmax')
cd ..
ls
ncdump swf_test_init.nc
hcchk=nc_varget('swf_test_init.nc','hc')
cd FRC2008
ls
edit create_roms_init.m
clear;create_roms_init
edit create_roms_init.m
clear;create_roms_init
cd ../analout
clear;make_avi_swf
%-- 5/10/10  9:57 AM --%
ls
edit create_roms_init
ls
clear;create_testwet_roms_init
%-- 5/20/10  9:28 AM --%
dir *.m
edit create _testwet_roms_init
edit create_testwet_roms_init
clear;create_testwet_roms_init
%-- 5/24/10  8:42 AM --%
clear;create_testwet_roms_init
%-- 5/25/10  7:01 AM --%
edit verify_salinity
clear;close all; verify_salinity
print -djpeg salinity
print -djpng salinity
print -dpng salinity
%-- 6/4/10  8:30 PM --%
dir *.m
edit verify_salinity
clear;verify_salinity
figure
subplot(211)
plot(jd,saltsfc,'b.-',jd,saltbot,'m.-',jdFM(infm),mean([salsfc(infm);salbot(infm)]),'r-')
grid
legend('model sfc','model bot','Mean SFWMD')
title('Observed and Predicted Values at Fort Myers')
ylabel('Salinity (psu)')
datetick('x','keeplimits','keepticks')
figure
subplot(211)
plot(jd,saltsfc,'b.-',jd,saltbot,'m.-',jdFM(infm),mean([salsfc(infm) salbot(infm)]'),'r-')
grid
legend('model sfc','model bot','Mean SFWMD')
title('Observed and Predicted Values at Fort Myers')
ylabel('Salinity (psu)')
datetick('x','keeplimits','keepticks')
figure(3);print -djpeg salinity_mean
%-- 6/8/10  8:29 AM --%
ls
edit create_roms_init
create_roms_init
%-- 6/9/10  8:16 AM --%
clear;close all ; verify_salinity
cd ../FRC2008
ls
create_roms_init
edit create_roms_init
create_roms_init
clear;create_roms_init
%-- 6/14/10  9:42 AM --%
ls *.m
clear;verify_salinity
clear;close all;verify_salinity
%-- 6/21/10  8:42 PM --%
dir *.m
clear;gen_transpaths
clear;gen_transpathsSWF
clear;make_movie_swf
%-- 6/22/10 12:03 PM --%
clear;make_movie_swf
t
figure;plot(t);datetick('y')
figure;plot(t,'*');datetick('y')
grid
%-- 6/22/10 12:41 PM --%
clear;make_movie_swf
%-- 7/5/10 10:55 AM --%
ls
verify_salinity
ls
cd ../FRC2008
ls
clear;create_roms_init
%-- 7/12/10  9:46 AM --%
ls
create_roms_init
%-- 7/15/10  9:44 AM --%
dir *.m
create_roms_init
%-- 7/15/10  6:06 PM --%
dir *.m
create_roms_init
%-- 7/21/10  2:40 PM --%
wind_time=nc_varget('swf_wind2008_NW.nc');
wind_time=nc_varget('swf_wind2008_NW.nc','wind_time');
Uwind=nc_varget('swf_wind2008_NW.nc','U_wind(:,200,200)');
Uwind=nc_varget('swf_wind2008_NW.nc','U_wind');
Uwind=nc_varget('swf_wind2008_NW.nc','Uwind');
min(wind_time)
max(wind_time)
figure;plot(wind_time,Uwind(:,200,200),'*-')
Vwind=nc_varget('swf_wind2008_NW.nc','vwind');
Vwind=nc_varget('swf_wind2008_NW.nc','Vwind');
save windforcdat
clear
ls *wind*
wind_time=nc_varget('swf_wind2008.nc','wind_time');
Uwind=nc_varget('swf_wind2008.nc','Uwind');
Vwind=nc_varget('swf_wind2008.nc','Vwind');
Uwind=Uwind(:,200,200);
Vwind=Vwind(:,200,200);
save windforcdat
%-- 7/22/10  3:46 PM --%
wind_time=nc_varget('swf_wind2008_NW.nc','wind_time');
Uwind=nc_varget('swf_wind2008_NW.nc','Uwind');
Uwind=Uwind(:,200,200);
Vwind=nc_varget('swf_wind2008_NW.nc','Vwind');
Vwind=Vwind(:,200,200);
save windNWforcdat
%-- 8/3/10 12:31 PM --%
ls
fn='swf_wind2008_SW.nc';
wind_time=nc_varget(fn,'wind_time');
Uwind=nc_varget(fn,'Uwind');
Vwind=nc_varget(fn,'Vwind');
save windSWforcdat
%-- 8/3/10  2:11 PM --%
fn='swf_wind2008_SW.nc';
wind_time=nc_varget(fn,'wind_time');
Uwind=nc_varget(fn,'Uwind');
Uwind=Uwind(:,200,200);
Vwind=nc_varget(fn,'Vwind');
Vwind=Vwind(:,200,200);
save windSWforcdat
%-- 8/5/10  6:45 PM --%
ls *.m
create_roms_init_wet
%-- 8/7/10 10:10 PM --%
gen_transpathsSWF
clear
ls *.m
make_movie_swf.m
make_movie_swf
clear;close all
make_movie_swf
%-- 8/9/10  9:34 PM --%
dir
create_roms_init_wet
%-- 8/11/10  8:05 PM --%
dir
create_roms_init_dry
dir
%-- 8/17/10 11:03 AM --%
dir
create_roms_init_dry
%-- 8/22/10  3:34 PM --%
ls
chkout
figure(3);print -djpeg w26.jpg
f
time=nc_varget(f,'ocean_time');
close all
clear;dchkout
clear;chkout
w_
%-- 8/24/10  2:35 PM --%
chkout
clear
chkout
edit chkout
close all
clear;chkout
u=nc_varget(f,'u');
v=nc_varget(f,'v');
size(u)
eidit
edit
for t=235:240,
u_=reshape(u(t,z,:,:),200,249);
% if t==1,
%     figure;
% else
%     figure(gcf)
% end
figure(2)
pcolor(u_);
shading flat;
caxis([-2 2])
colorbar;
title(['u, z=' num2str(z) ', t=' ...
num2str(t)])
pause(0.5)
end
for t=235:240,
v_=reshape(v(t,z,:,:),199,250);
% if t==1,
%     figure;
% else
%     figure(gcf)
% end
figure(3)
pcolor(v_);
shading flat;
caxis([-2 2])
colorbar;
title(['v, z=' num2str(z) ', t=' ...
num2str(t)])
pause(0.5)
endexit
%-- 9/9/10 10:50 AM --%
dir *.m
create_roms_initwetlate
create_roms_init_wetlate
%-- 9/9/10  6:46 PM --%
dir
create_roms_init_dry8
%-- 9/14/10  2:01 PM --%
dir
create_roms_init_wet8
%-- 9/23/10  8:13 PM --%
dir *make
ls make*
make_movie_swf8
close all
clear
make_movie_swf8
%-- 9/24/10 11:34 AM --%
make_movie8test
clclear
clear
make_movie8test
clear
closae all
close all
make_movie8test
vi make_movie8tes
k
k:100
nfiles
dir *.mclose all
close all
dir *.m
clear
gen_transpathsSWF
clear
mak_movie8test
make_movie8test
clear;close all;make_movie8test
%-- 9/25/10 12:41 AM --%
make_movie8test
clear
close all
make_movie8test
%-- 9/25/10  8:02 PM --%
make_movie8test
clear;close all
make_movie8test
clear;close all
make_movie8test
clear; close all
make_movie8test
clear;make_movie8test
clear;close all
make_movie8test
clear;close all
make_movie_s0_35
%-- 10/3/10  4:12 PM --%
chkout
close all
clear
chkout
cd ../analout
chkout
clear;close all;chkout
cd ../analout
clear;close all;chkout
cd ../analout
clear;close all;chkout
cd ../analout
%-- 10/3/10  5:22 PM --%
clear;close all;chkout
cd ../analout
edit chkout
f='tmpnc04.nc';
u=nc_varget(f,'u');
v=nc_varget(f,'v');
[T,R,C]=size(u);
cd ../OUTDRY8r
f='tmpnc04.nc';
u=nc_varget(f,'u');
v=nc_varget(f,'v');
[T,R,C]=size(u);
for t=350:T,
u_=reshape(u(t,z,:,:),200,249);
if t==350,
figure;
else
figure(gcf)
end
pcolor(u_);
shading flat;
caxis([-2 2])
colorbar;
title(['u, z=' num2str(z) ', t=' ...
num2str(t)])
pause(0.5)
end
for t=360:T,
v_=reshape(v(t,z,:,:),199,250);
if t==360,
figure;
else
figure(gcf)
end
pcolor(v_);
shading flat;
caxis([-2 2])
colorbar;
title(['v, z=' num2str(z) ', t=' ...
num2str(t)])
pause(0.5)
end
figure(1)
%-- 10/5/10  7:36 PM --%
make_meteorology_all2008
edit make_meteorology_all2008
clear;make_meteorology_all2008
%-- 10/5/10  8:11 PM --%
dir *.m
create_roms_init_dryFeb
%-- 10/6/10  8:33 PM --%
dir
time=nc_varget('calus2008_rst.nc','ocean_time')
format long
time
%-- 10/9/10 10:05 AM --%
dir swf*
create_init_dry6Feb
%-- 10/13/10 10:45 PM --%
gen_transpathsSWF
clear
make_movie_dry10_35
%-- 10/19/10  8:34 PM --%
fn='swf_wind2008_SW.nc';
wind_time=nc_varget(fn,'wind_time');
Uwind=nc_varget(fn,'Uwind');
Uwind=Uwind(:,200,200);
Vwind=nc_varget(fn,'Vwind');
Vwind=Vwind(:,200,200);
save windSWforcdat
fn='swf_wind2008_NW.nc';
wind_time=nc_varget(fn,'wind_time');
Uwind=nc_varget(fn,'Uwind');
Uwind=Uwind(:,200,200);
Vwind=nc_varget(fn,'Vwind');
Vwind=Vwind(:,200,200);
save windNWforcdat
fn='swf_wind2008.nc';
wind_time=nc_varget(fn,'wind_time');
Uwind=nc_varget(fn,'Uwind');
Uwind=Uwind(:,200,200);
Vwind=nc_varget(fn,'Vwind');
Vwind=Vwind(:,200,200);
save windforcdat
%-- 10/20/10  8:48 AM --%
make_movie_dry10_35
%-- 10/20/10  3:52 PM --%
make_movie_dry10_35
%-- 10/21/10  9:24 PM --%
plotfloatlatlondry8
ls plot*
pwd
cd analout
plotfloatlatlondry8
ls
pwd
cd ../analout
plotfloatlatlondry8
clear;dfigs;plotfloatlatlon_dt_dry8ani
clear;plotfloatlatlon_dt_dry8ani
clear;load oldwindforcdatNW
load windNWforcdat
me;
jd_wind=datenum(1968,5,23)wind_time;
figure;plot(jd_windclear;plotfloatlatlon_dt_dry8ani
plotfloatlatlon_dt_dryani
plotfloatlatlon_dt_dry8ani
clear
plotfloatlatlon_dt_dry8ani
clear;close all;plotfloatlatlon8dry
plotfloatlatlondry8
clear ;lotfloatlatlon_dt_dry8ani
plotfloatlatlon_dt_dry8ani
plotfloatlatlon_dt_dry8aniclear;close all;plotfloatlatlon_dt_dry8ani
clear
plotfloatlatlon_dt_dry8ani
ttata
clear
cd ../FRC2008
load windSWforcdat
plot(Uwind
plot(Uwind,Vwind,'*')
clear;makmatwind
cd ../analout
clear;makmatwind
clear;makematwind
clear
load windSWforcdat
jd=datenum(1968,5,24) + windtime;
jd=datenum(1968,5,24) + wind_time;
mangle=nc_varget('swf_grid8.nc','angle',[0 0],[-1 -1]);
w=(Uwind+sqrt(-1).*Vwind).*exp(sqrt(-1).*mangle(200,200));
plot(w,'*')
%-- 11/11/10  8:46 AM --%
dir *.m
makematwind
clear;makematwind
%-- 11/14/10  8:14 AM --%
dir make*10*
make_movie_dry10_35
clear;close all
make_movie_dry10_35
:wqclear;close all makemake_movie_dry20_35
clear ; close all
make_movie_dry20_30
make_movie_dry35
make_movie_dry20_35
%-- 11/16/10 11:34 AM --%
make_movie_dry20_35
%-- 12/22/10  5:47 PM --%
addpath ~/AARONS_TOOLS/PLOTTING/
addpath ~/USGS-toolbox/ecomvis-sea-mat
addpath ~/USGS_toolbox/ecomvis-sea-mat
help psliceuv
addpath ~/USGS_toolbox/ecomviz-sea-mat
help psliceuv
dir make*
edit make_sed_movie
%-- 12/22/10  7:14 PM --%
dir
cd CALUS/analout
dir make*
make_sed_movie
clear;make_sed_movie
close all
make_sed
close all
makemake_sed_movie
make_sed_movie
nfiles
close all
make_sed_movie
l
clear
cd ../FRC2008
create_roms_initAug
%-- 12/26/10 10:04 AM --%
dir *.nc
more on
nc_dump('swf_river2008.nc')
dir
ear
cre
ceate_roms_initAugSedmore off
pwed
pwd
dir
create_init_initAugSed
create_init_
create_roms_initAugSed
cd
n
more off
more on
dir *river*
nc_dump('swf_river2008.nc')
cle
clear
AddRiverSed
clear;AddRiverSed
ls
d
more off
pwd
cd ../analout
dir mke*
%-- 12/26/10  7:46 PM --%
make_sed_movie
clear;make_sed_movie
/home/beach/guest/dafu8991/CALUS/analout
lshelp nc_var_get
clc
help  nc_varget
clc;help squeeze
help depths
clc
dir *.nc
more on;nc_dump('../OUTAUG/calus2008_his001.nc')
more on;nc_dump('../OUTAUG/calus2008_his1.nc')
ls ../OUTAUG/*1*.nc
more on;nc_dump('../OUTAUG/calus2008_his0001.nc')
more on;nc_dump('../OUTAUG/calus2008_his_0001.nc')
more off
clc
mZ
size(mZ)
clc
size(mZ)
help permute
jj
help pclegen2
clear;clc
movieSmallScale
make_sed_movieSmallScale
clc;more off
make_sed_movieSmallScale
clc
clear;make_sed_movieSmallScale
ischclc
clc
size(rin)
size(rn)
min(rn)
ans-max(rn)
clear;make_sed_movieSmallScale
%-- 12/28/10  6:37 PM --%
dir gen*
gen_transpathsSWF
clc
dir *.mat
dir make*
clear;close all;make_sed_movieSmallScale
clc
clear;close all;make_sed_movieSmallScale
length(t)
ii
type tind
size(transport_time)
size(t)
clc
t
clc
gregorian(t(ii))
type tind
plot(transport_time,[1:length(transport_time)],'*');datetick
length(transport_time)
41*64
clc
41*26
gregorian(t(ii))
clear;close all;gen_transpathsSWF
cleaqr
clc
clear
make_sed_movieSmallScale
clc
clear;make_sed_movieSmallScale
clc
clear;make_sed_movieSmallScale
clc
close all
clear;make_sed_movieSmallScale
%-- 12/29/10  9:54 AM --%
clear;make_sed_movieSmallScale
size(mZ)
clc
mZ(:,200,200)
mZ(:,3,3)
clc
clear;make_sed_movieSmallScale_2
size sfcconc
size(sfcconc)
sisize mrlon)
clc
size(mrlon)
clc
clear;make_sed_movieSmallScale_2
