v4_3deg/input/gen_init.m

% Generate initial condition files
pin='../../../../MITgcm/verification/tutorial_global_oce_biogeo/input/';

% Uvel, Vvel, Theta, Salt, Eta
fin=[pin 'pickup.0005184000.data'];
tmp=readbin(fin,[128,64,123],1,'real*8');
writebin('Uvel.0005184000',tmp(:,:,1:15));
writebin('Vvel.0005184000',tmp(:,:,16:30));
writebin('Theta.0005184000',tmp(:,:,31:45));
writebin('Salt.0005184000',tmp(:,:,46:60));
writebin('Eta.0005184000',tmp(:,:,121));

% DIC, Alk, PO4, DOP, O2
fin=[pin 'pickup_ptracers.0005184000.data'];
tmp=readbin(fin,[128,64,75],1,'real*8');
writebin('DIC.0005184000',tmp(:,:,1:15));
writebin('Alk.0005184000',tmp(:,:,16:30));
writebin('PO4.0005184000',tmp(:,:,31:45));
writebin('DOP.0005184000',tmp(:,:,46:60));
writebin('O2.0005184000',tmp(:,:,61:75));
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	writebin('Uvel.0005184000',tmp(:,:,1:15));
8	writebin('Vvel.0005184000',tmp(:,:,16:30));
9	writebin('Theta.0005184000',tmp(:,:,31:45));
10	writebin('Salt.0005184000',tmp(:,:,46:60));
11	writebin('Eta.0005184000',tmp(:,:,121));
12
13	% DIC, Alk, PO4, DOP, O2
14	fin=[pin 'pickup_ptracers.0005184000.data'];
15	tmp=readbin(fin,[128,64,75],1,'real*8');
16	writebin('DIC.0005184000',tmp(:,:,1:15));
17	writebin('Alk.0005184000',tmp(:,:,16:30));
18	writebin('PO4.0005184000',tmp(:,:,31:45));
19	writebin('DOP.0005184000',tmp(:,:,46:60));
20	writebin('O2.0005184000',tmp(:,:,61:75));