| 1 |
dimitri |
1.1 |
% Generate initial condition files |
| 2 |
|
|
pin='../../../../MITgcm/verification/tutorial_global_oce_biogeo/input/'; |
| 3 |
|
|
|
| 4 |
|
|
% Uvel, Vvel, Theta, Salt, Eta |
| 5 |
|
|
fin=[pin 'pickup.0005184000.data']; |
| 6 |
|
|
tmp=readbin(fin,[128,64,123],1,'real*8'); |
| 7 |
dimitri |
1.2 |
writebin('../input/Uvel.0005184000',tmp(:,:,1:15)); |
| 8 |
|
|
writebin('../input/Vvel.0005184000',tmp(:,:,16:30)); |
| 9 |
|
|
writebin('../input/Theta.0005184000',tmp(:,:,31:45)); |
| 10 |
|
|
writebin('../input/Salt.0005184000',tmp(:,:,46:60)); |
| 11 |
|
|
writebin('../input/Eta.0005184000',tmp(:,:,121)); |
| 12 |
|
|
|
| 13 |
|
|
% Area, Heff |
| 14 |
|
|
Ice=tmp(:,:,31); |
| 15 |
|
|
Ice(find(Ice>-1.8))=0; |
| 16 |
|
|
Ice(find(Ice))=1; |
| 17 |
|
|
writebin('../input/Area.0005184000',Ice); |
| 18 |
|
|
writebin('../input/Heff.0005184000',Ice); |
| 19 |
dimitri |
1.1 |
|
| 20 |
|
|
% DIC, Alk, PO4, DOP, O2 |
| 21 |
|
|
fin=[pin 'pickup_ptracers.0005184000.data']; |
| 22 |
|
|
tmp=readbin(fin,[128,64,75],1,'real*8'); |
| 23 |
dimitri |
1.2 |
writebin('../input/DIC.0005184000',tmp(:,:,1:15)); |
| 24 |
|
|
writebin('../input/Alk.0005184000',tmp(:,:,16:30)); |
| 25 |
|
|
writebin('../input/PO4.0005184000',tmp(:,:,31:45)); |
| 26 |
|
|
writebin('../input/DOP.0005184000',tmp(:,:,46:60)); |
| 27 |
|
|
writebin('../input/O2.0005184000',tmp(:,:,61:75)); |
| 28 |
dimitri |
1.3 |
|
| 29 |
|
|
% individual ptracers extracted from |
| 30 |
|
|
% /nobackup/dcarrol2/v4_3deg/pickup/v4_3deg_pickup_ptracers_optimized.0000000001.data |
| 31 |
|
|
% Type is real*4, dimensions are [128 64 15 40]. |
| 32 |
|
|
for i=1:40 |
| 33 |
|
|
tmp=readbin('v4_3deg_pickup_ptracers_optimized.0000000001.data',[128 64 15],1,'real*4',i-1); |
| 34 |
|
|
writebin(['../input/ptracers_optimized_' myint2str(i) '.0000000001'],tmp); |
| 35 |
|
|
end |
Parent Directory
| 
Revision Log
| 
Revision Graph