| 1 |
ts=354300; |
| 2 |
pn='/nobackupp8/dmenemen/llc/llc_1080/MITgcm/run2/'; |
| 3 |
fld1={'Eta', |
| 4 |
'KPPhbl', |
| 5 |
'PhiBot', |
| 6 |
'SIarea', |
| 7 |
'SIheff', |
| 8 |
'SIhsalt', |
| 9 |
'SIhsnow', |
| 10 |
'SIuice', |
| 11 |
'SIvice', |
| 12 |
'Salt', |
| 13 |
'Theta', |
| 14 |
'U', |
| 15 |
'V', |
| 16 |
'W', |
| 17 |
'oceFWflx', |
| 18 |
'oceQnet', |
| 19 |
'oceQsw', |
| 20 |
'oceSflux', |
| 21 |
'oceTAUX', |
| 22 |
'oceTAUY'}; |
| 23 |
fld2={'_ETAN', |
| 24 |
'_KPPhbl', |
| 25 |
'_PHIBOT', |
| 26 |
'_SIarea', |
| 27 |
'_SIheff', |
| 28 |
'_SIhsalt', |
| 29 |
'_SIhsnow', |
| 30 |
'_SIuice', |
| 31 |
'_SIvice', |
| 32 |
'_SALTanom', |
| 33 |
'_THETA', |
| 34 |
'_UVELMASS', |
| 35 |
'_VVELMASS', |
| 36 |
'_WVELMASS', |
| 37 |
'_oceFWflx', |
| 38 |
'_oceQnet', |
| 39 |
'_oceQsw', |
| 40 |
'_oceSflux', |
| 41 |
'_oceTAUX', |
| 42 |
'_oceTAUY'}; |
| 43 |
|
| 44 |
for f=1:20 |
| 45 |
f1=[pn fld1{f} '.' myint2str(ts-20,10) '.data']; |
| 46 |
f2=[pn fld2{f} '.' myint2str(ts,10) '.data']; |
| 47 |
fl1=quikread_llc(f1,1080); |
| 48 |
fl2=quikread_llc(f2,1080); |
| 49 |
if f==10, fl2=fl2+35; end |
| 50 |
clf |
| 51 |
subplot(311) |
| 52 |
quikplot_llc(fl1) |
| 53 |
thincolorbar |
| 54 |
title(fld1{f}) |
| 55 |
subplot(312) |
| 56 |
quikplot_llc(fl2) |
| 57 |
thincolorbar |
| 58 |
title(fld2{f}) |
| 59 |
subplot(313) |
| 60 |
quikplot_llc(fl2-fl1) |
| 61 |
thincolorbar |
| 62 |
pause |
| 63 |
end |
| 64 |
|
| 65 |
k=88; |
| 66 |
for f=10:13 |
| 67 |
f1=[pn fld1{f} '.' myint2str(ts-20,10) '.data']; |
| 68 |
f2=[pn fld2{f} '.' myint2str(ts,10) '.data']; |
| 69 |
fl1=quikread_llc(f1,1080,k); |
| 70 |
fl2=quikread_llc(f2,1080,k); |
| 71 |
if f==10, fl2=fl2+35; end |
| 72 |
clf |
| 73 |
subplot(311) |
| 74 |
quikplot_llc(fl1) |
| 75 |
thincolorbar |
| 76 |
title(fld1{f}) |
| 77 |
subplot(312) |
| 78 |
quikplot_llc(fl2) |
| 79 |
thincolorbar |
| 80 |
title(fld2{f}) |
| 81 |
subplot(313) |
| 82 |
quikplot_llc(fl2-fl1) |
| 83 |
thincolorbar |
| 84 |
pause |
| 85 |
end |
| 86 |
|
| 87 |
|
| 88 |
============= |
| 89 |
|
| 90 |
ts1=430080; |
| 91 |
ts2=431040; |
| 92 |
u1=quikread_llc(['U.' myint2str(ts1,10) '.data'],1080); |
| 93 |
v1=quikread_llc(['V.' myint2str(ts1,10) '.data'],1080); |
| 94 |
e1=quikread_llc(['Eta.' myint2str(ts1,10) '.data'],1080); |
| 95 |
u2=quikread_llc(['U.' myint2str(ts2,10) '.data'],1080); |
| 96 |
v2=quikread_llc(['V.' myint2str(ts2,10) '.data'],1080); |
| 97 |
e2=quikread_llc(['Eta.' myint2str(ts2,10) '.data'],1080); |
| 98 |
|
| 99 |
clf |
| 100 |
orient tall |
| 101 |
subplot(311) |
| 102 |
quikplot_llc(e1) |
| 103 |
caxis([-1 1]*2.7) |
| 104 |
thincolorbar |
| 105 |
title(['Eta on ' ts2dte(ts1,90,2010)]) |
| 106 |
subplot(312) |
| 107 |
quikplot_llc(e2) |
| 108 |
caxis([-1 1]*2.7) |
| 109 |
thincolorbar |
| 110 |
title(['Eta on ' ts2dte(ts2,90,2010)]) |
| 111 |
subplot(313) |
| 112 |
quikplot_llc(e2-e1) |
| 113 |
caxis([-1 1]*.7) |
| 114 |
thincolorbar |
| 115 |
title('Difference') |
| 116 |
eval(['print -dpsc Eta_day' int2str(ts2*90/60/60/24)]) |
| 117 |
|
| 118 |
clf |
| 119 |
orient tall |
| 120 |
subplot(311) |
| 121 |
quikplot_llc(sqrt(u1.^2+v1.^2)) |
| 122 |
caxis([0 1]*.7) |
| 123 |
thincolorbar |
| 124 |
title(['Surface speed on ' ts2dte(ts1,90,2010)]) |
| 125 |
subplot(312) |
| 126 |
quikplot_llc(sqrt(u2.^2+v2.^2)) |
| 127 |
caxis([0 1]*.7) |
| 128 |
thincolorbar |
| 129 |
title(['Surface speed on ' ts2dte(ts2,90,2010)]) |
| 130 |
subplot(313) |
| 131 |
quikplot_llc(sqrt(u2.^2+v2.^2)-sqrt(u1.^2+v1.^2)) |
| 132 |
caxis([-1 1]*.25) |
| 133 |
thincolorbar |
| 134 |
title('Difference') |
| 135 |
eval(['print -dpsc Speed_day' int2str(ts2*90/60/60/24)]) |
Parent Directory
| 
Revision Log
| 
Revision Graph