[MITgcm-support] Defining different regions over the same area with regional diagnostics

[Martin Losch]

Hi Alessandro,

this can be done with defining multi-layer masks. I assume that you 
followed the example in, e.g., 
verification/global_ocean.cs32x15/input.thsice/data.diagnostics? In that 
example, the mask is a 2D field.
You can now add a 2nd layer as a 3rd dimension and in this second layer 
you can assign mask values 4, 5, ... etc. In data.diagnostics, you need 
to set:
  set_regMask(1:3) = 1,  1,  1,  2,  2, ...
  val_regMask(1:3) = 1., 2., 3., 4., 5., ...

in DIAGNOSTICS_SIZE.h you may have to adjust this line:
       PARAMETER( nRegions = 10 , sizRegMsk = 3 , nStats = 4 )
if you want to use more that 3 layers and more than 10 regions.

And yes, it would be nice if this was documented somewhere. If you would 
like to volunteer (after success, of course and while waiting 
impatiently for your new diagnostics output) to write some instructions, 
I'd be happy to add them to the online documentation.

Martin

On 06/21/2016 11:34 AM, gmail wrote:
> Dear MITgcm users,
>
> I am using the regional diagnostic package on my Arctic configuration of the MITgcm.
> I have presently defined 3 disjointed regions in my configuration.
> This works fine but now I would like to add some other regions for the diagnostics.
> In particular I want to define some subregions over the already defined regions.
> Moreover I would like as an output both the large region and the subregion.
>
> Now, because of the way the regional diagnostics are defined it appears it is not possible to have two regions defined over the same area.
> I have searched this a bit but I couldn’t find any possible solution.
>
> Do you have any suggestions?
>
> All the best,
> Alessandro
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>

-- 
Martin Losch
Alfred Wegener Institute for Polar and Marine Research
Postfach 120161, 27515 Bremerhaven, Germany;
Tel./Fax: ++49(0471)4831-1872/1797