[MITgcm-support] cubed sphere: init format, subdomains and tiles
Jean-Michel Campin
jmc at ocean.mit.edu
Mon May 7 10:14:55 EDT 2012
Hi Maura,
Few things:
> I have compiled the 2 experiments with the corresponding W2_EXCH2_SIZE.h and
> w2_e2setup.F files, obtained by driver.m.
You don't need to generate theses files (see top of the file:
MITgcm/utils/exch2/code-mods/README); might need a parameter file "data.exch2"
only if you have blank tiles or non regular (i.e., with faces of different
size) cube-sphere grid.
The way you generate input files seems OK to me (one of the simplest example
of a matlab script that is used to generate input file on regular CS-grid is:
MITgcm/verification/adjustment.cs-32x32x1/input.nlfs/gendata.m).
And if you use MITgcm/utils/matlab/rdmds.m to read in your output files,
in both cases, the loaded array should come with a [60*6 60 nr] shape.
The reason why your output files have not the same shape (after using rdmds.m)
might be related to the old W2_EXCH2_SIZE.h + w2_e2setup.F files you are using.
Cheers,
Jean-Michel
On Mon, May 07, 2012 at 02:41:32PM +0200, Maura BRUNETTI wrote:
> Dear MITgcm users,
>
> I am testing different cubed sphere configurations. I have constructed initial conditions in the
> format (60x6x60) (i.e. bathymetry [60 6 60], Levitus climatology [60 6 60 nr=15], ncep files
> and Levitus relaxation [60 6 60 nt=12]): in the first experiment, I use 6 procs (sNx =sNy =60,
> nSx =1, nPx=6), while in the second experiment I use 24 procs (sNx = sNy = 30, nSx =1, nPx =24).
> Thus, the subdomain dimension (sNx,sNy) and the number of tiles change in the 2 experiments.
> I have compiled the 2 experiments with the corresponding W2_EXCH2_SIZE.h and
> w2_e2setup.F files, obtained by driver.m.
>
> Since I obtain different results with the 2 configurations, I am trying to understand where the error
> is ... Is it correct to use the same initial files in the format (60x6x60) for the 2 experiments?
>
> At the end of the second experiment, I obtain arrays of the format (720x30) (for example,
> the temperature has the following dimension: [720 30 15 4] after 4 output). So they are not in
> the standard format (60*6 60). I have found a possible way to convert this format in the standard
> one, but I am not sure that it is correct... Is there a matlab file which properly reshape the array in
> the standard format (just using t = reshape(t,[360 60 15 4]) is not sufficient...)
>
> Cheers,
>
> Maura
>
>
> --
> Dr. Maura Brunetti
> Institute for Environmental Sciences (ISE)
> University of Geneva -- Switzerland
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list