[MITgcm-support] Re: Some questions

[Dimitris Menemenlis]

> Hi, menemenlis
> 
> 1. I have run mitgcm as
> http://mitgcm.org/cgi-bin/viewcvs.cgi/MITgcm_contrib/high_res_cube/README_ice
> instructed, but my result(bw.1346.out) is different from yours(output.txt),
> is there any mistake?

output.txt was out of date.  I have just updated experiment to be compatible
with latest MITgcm code

> 2. I can not find calculated ice data, why?

should be there, they are called VICEtave*, UICEtave*, HEFFtave*, and AREAtave*
try downloading a fresh copy of code and rerunning.  follow instruction in
http://mitgcm.org/cgi-bin/viewcvs.cgi/MITgcm_contrib/high_res_cube/README_ice

> 3. How to plot the calculated results?

very crude graphics using matlab from the README_ice instructions:

to use matlab to look at the output
matlab

for fld={'VICEtave','UICEtave','HEFFtave','AREAtave','QSWtave', ...
         'QNETtave','FVtave','FUtave','EmPmRtave'}
    tmp=permute(readbin([fld{1} '.0000000020.data'],[32,6,32]),[1 3 2]);
    cx=[min(tmp(:)) max(tmp(:))]; clf, plot_cube, title(fld{1}), pause
end

heff=permute(readbin(['HEFFtave.0000000020.data'],[32,6,32]),[1 3 2]);
mask=0*heff; mask(find(heff>.001))=1;
uice=mask.*permute(readbin(['UICEtave.0000000020.data'],[32,6,32]),[1 3 2]);
vice=mask.*permute(readbin(['VICEtave.0000000020.data'],[32,6,32]),[1 3 2]);
clf
subplot(221), mypcolor(uice(:,:,3)'); colorbar, title('uice, tile 3')
subplot(222), mypcolor(uice(:,:,6)'); colorbar, title('uice, tile 6')
subplot(223), mypcolor(vice(:,:,3)'); colorbar, title('vice, tile 3')
subplot(224), mypcolor(vice(:,:,6)'); colorbar, title('vice, tile 6')

> 
> Best regards,
>  				
> 
> 　　　　　　　　liuyu
> 　　　　　　　　ly_elvis at yahoo.com.cn
> 　　　　　　　　　　2006-06-02

-- 
Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Institute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
tel: 818-354-1656;  fax: 818-393-6720