[MITgcm-support] RE: Any example of fortran code forgenerating expanded spherical cube ?.

[Dimitris Menemenlis]

>> You said that  "... there are also other  examples with
>> the cube sphere,  etc., that I can send you.". 
>> 
>> Would you send me the examples for me ?. 

Sung-Dae, instructions for setting up and integrating a 32x32x6x15
cube-sphere, similar to the verification/global_ocean.cs32x15, but
using pkg/exf to read-in the surface forcing from a lat/long grid
are appended.  More cube-sphere examples are described in the file
MITgcm_contrib/high_res_cube/README_ice

The lat/long grid for the input files is specified in data files
called data.exf and data.exf_clim.

D.

##########################################################################
##########################################################################
# getting the code from anonymous CVS server

setenv CVSROOT :pserver:cvsanon at mitgcm.org:/u/u0/gcmpack
cvs login ( CVS password: cvsanon )
cvs co -r hrcube_3 MITgcm_contrib/high_res_cube/README_ice
cvs co -r hrcube_3 MITgcm_contrib/high_res_cube/code-mods
cvs co -r hrcube_3 MITgcm_contrib/high_res_cube/matlab-grid-converter
cvs co -r hrcube_3 MITgcm_contrib/high_res_cube/input
cvs co -r hrcube_3 MITgcm_contrib/high_res_cube/results
cvs co -r hrcube_3 MITgcm_code
cvs co -r hrcube_3 MITgcm/verification/global_ocean.cs32x15


##########################################################################
##########################################################################
# 32*32*6*15, 1-cpu cube sphere

cd MITgcm/verification/global_ocean.cs32x15/code
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/* .
\cp ../../../../MITgcm_contrib/high_res_cube/code-mods/s12t_16x32/* .

cd ../build
../../../tools/genmake2 -mods=../code
make depend
make

cd ..
mkdir run
cd run
\cp ../input/* .
\cp ../../../../MITgcm_contrib/high_res_cube/matlab-grid-converter/tile00* .
\cp ../../../../MITgcm_contrib/high_res_cube/input/* .
../build/mitgcmuv >& output.txt

comparison output is in:
../../../../MITgcm_contrib/high_res_cube/results/output.txt

to use matlab to look at the output
matlab

for fld={'Eta','VICE','UICE','HEFF','AREA','Qsw','Qnet','FV','FU','EmPmR'}
    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(['HEFF.0000000020.data'],[32,6,32]),[1 3 2]);
mask=0*heff; mask(find(heff>.001))=1;
uice=mask.*permute(readbin(['UICE.0000000020.data'],[32,6,32]),[1 3 2]);
vice=mask.*permute(readbin(['VICE.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')




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