# UM Nesting suite - LAM\_INCLUSION grids not aligning

**URL:** <https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745>\
**Category:** rAM3/Regional Nesting Suite\
**Created:** [8 May 2023 06:03 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745 "2023-05-08T06:03:18Z")\
**Posts on this page:** 20\
**Page:** 1

<div class="post-metadata">

**Author:** ![sonyafiddes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/sonyafiddes/32/603_2.png) [@sonyafiddes](https://forum.access-hive.org.au/u/sonyafiddes)\
**Post date:** [8 May 2023 06:03 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/1 "2023-05-08T06:03:18Z")

</div>

Hi team,  
Running the UM nesting suite vn13.1, I am getting the following error in the Regn1\_resn\_1\_RA2M\_um\_createbc\_000 task:

? Error code: 4  
? Error from routine: LAM\_INCLUSION  
? Error message: Eastern boundary at 506.1782500 not within source LAM with 146.1647500  
? Error from processor: 0  
? Error number: 2

It appears that 360 degrees is being added to my eastern boundary somewhere, but I’m not sure why, as this routine has worked in the past for me and as far as I can tell has not changed (from the met office log files). I’ve checked the output files so far (inc. the ancillary files, the ERA5 grib files and the ec\_cb000 files) and they all look correct.

If anyone can offer up any advise I’d be most grateful!

Also, I am wondering if this is the most suitable place for help or should I be posting this to the access help desk?

Cheers,  
Sonya

---

<div class="post-metadata">

**Author:** ![ZhangchengPei](https://avatars.discourse-cdn.com/v4/letter/z/8e8cbc/32.png) [@ZhangchengPei](https://forum.access-hive.org.au/u/ZhangchengPei)\
**Post date:** [8 May 2023 06:17 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/2 "2023-05-08T06:17:03Z")

</div>

Hi team,

This also happens to my suite when I try to run it over the Davis Station, Antarctica with a rotated pole.

It failed in the ec\_um\_recon\_000 task:

? Error code: 4  
? Error from routine: LAM\_INCLUSION  
? Error message: Eastern boundary at 442.3114347 not within source LAM with 182.2500000  
? Error from processor: 90  
? Error number: 6

The difference of eastern boundary is not 360, which may be caused by the rotated pole?

Cheers,  
Zhangcheng

---

<div class="post-metadata">

**Author:** ![Paola-CMS](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/paola-cms/32/57_2.png) [@Paola-CMS](https://forum.access-hive.org.au/u/Paola-CMS)\
**Post date:** [8 May 2023 07:03 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/3 "2023-05-08T07:03:55Z")

</div>

Hi Sonya,

just pointing out that I mentioned the helpdesk specifically for the error with era5grib as we (CMS-CLEX) maintain that and it was just by accident that I noticed your question on the forum. So it’s a different case from your current question. 🙂  
Cheers,  
Paola

---

<div class="post-metadata">

**Author:** ![sonyafiddes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/sonyafiddes/32/603_2.png) [@sonyafiddes](https://forum.access-hive.org.au/u/sonyafiddes)\
**Post date:** [8 May 2023 23:07 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/4 "2023-05-08T23:07:58Z")

</div>

Thank Paola - I had gathered that! 🙂

---

<div class="post-metadata">

**Author:** ![Scott](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/scott/32/37_2.png) [@Scott](https://forum.access-hive.org.au/u/Scott)\
**Post date:** [9 May 2023 00:56 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/5 "2023-05-09T00:56:40Z")

</div>

Hi Sonya, what are the boundaries of both your lam domain and of the ERA5 grib data? There is an option in the ERA5grib settings to include all longitudes which may be needed if working at the date line or near the poles

---

<div class="post-metadata">

**Author:** ![Aidan](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/aidan/32/42_2.png) [@Aidan](https://forum.access-hive.org.au/u/Aidan)\
**Post date:** [9 May 2023 01:59 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/6 "2023-05-09T01:59:27Z")

</div>

I don’t know if this helpful, but Monika Krysta had a similar error code

[https://code.metoffice.gov.uk/trac/roses-u/wiki/ticket/182/nci\_3dvar\_hybrid\_lam\_aw264\_sy](https://code.metoffice.gov.uk/trac/roses-u/wiki/ticket/182/nci_3dvar_hybrid_lam_aw264_sy)

---

<div class="post-metadata">

**Author:** ![Scott](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/scott/32/37_2.png) [@Scott](https://forum.access-hive.org.au/u/Scott)\
**Post date:** [9 May 2023 02:06 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/7 "2023-05-09T02:06:57Z")

</div>

Ah yes the other thing to check is that the era5 domain needs to be a bit larger than the nest domain for the boundaries to work - say 10 to 20 grid lengths on each side, so if you had a 0.1 degree grid size on your nest then you’d want the era5 domain to be 1 to 2 degrees larger than the nest domain on each side.

This would show up in the error message as the boundaries only having a small difference though, so not likely to be the case here

---

<div class="post-metadata">

**Author:** ![sonyafiddes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/sonyafiddes/32/603_2.png) [@sonyafiddes](https://forum.access-hive.org.au/u/sonyafiddes)\
**Post date:** [10 May 2023 00:02 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/8 "2023-05-10T00:02:29Z")

</div>

Hi Scott,

I am trying to run near Cape Grim in Tassie (though I did stuff up my lats here - so not quite right), so there shouldn’t be an issue being near the poles or date line (though I note Zhangcheng is running close to the pole with a similar issue).

In the more recent versions, the nesting suite doesn’t ask you to define the boundaries itself, rather it takes them from a file generated from the Ancillary suite: grid.nl :

&grid  
delta\_lambda\_targ=0.0135,  
delta\_phi\_targ=0.0135,  
global=.false.,  
igrid\_targ=6,  
inwsw=1,  
lambda\_origin\_targ=143.215,  
lambda\_pole=180.0,  
phi\_origin\_targ=-46.185,  
phi\_pole=90.0,  
points\_lambda\_targ=220,  
points\_phi\_targ=220,  
rotated=.false.,

I am wondering if I should just change the points lambda/phi target to 200 to reduce the domain the nesting suite is after, in a similar way to what we did previously?

Thanks for taking the time on this.

---

<div class="post-metadata">

**Author:** ![Scott](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/scott/32/37_2.png) [@Scott](https://forum.access-hive.org.au/u/Scott)\
**Post date:** [10 May 2023 00:32 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/9 "2023-05-10T00:32:17Z")

</div>

Looks fine, what are the grid details for the driving era5 data?

---

<div class="post-metadata">

**Author:** ![sonyafiddes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/sonyafiddes/32/603_2.png) [@sonyafiddes](https://forum.access-hive.org.au/u/sonyafiddes)\
**Post date:** [10 May 2023 00:59 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/10 "2023-05-10T00:59:39Z")

</div>

Well thats the thing possibly - I think they are the same, as in at least the GUI there is now nowhere to set them differently that I can see? I’ve looked that the grib files created by the ERAgrib task and they look ok - ie. they are slightly bigger than the domain I have set.

 ![Screenshot 2023-05-10 at 10.57.10 am](https://us1.discourse-cdn.com/flex020/uploads/access1/original/1X/0c1221205abbb1907f0eb491a77cd5b188069613.png)

---

<div class="post-metadata">

**Author:** ![Scott](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/scott/32/37_2.png) [@Scott](https://forum.access-hive.org.au/u/Scott)\
**Post date:** [10 May 2023 01:25 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/11 "2023-05-10T01:25:42Z")

</div>

Ah right. So how this works is it starts off with the GRIB era5 data, it will then convert that GRIB data to UM format by running it through the reconfiguration, then the CREATEBC task will pull out the boundary data from the UM format file.

When the GRIB to UM step is run the reconfiguration needs a set of ancillary files to fill in extra fields, which is set at dm\_ec\_lam\_ancil\_dir:  
 ![image](https://us1.discourse-cdn.com/flex020/uploads/access1/original/1X/8a002ada59bea51dcf72faae3f7ff6a160f55ca0.jpeg)

These ancillaries should be created by the regional ancillary suite, at around 0.1 degree resolution (in case the UM grid doesn’t exactly match up with the ERA5 grid) and covering a somewhat larger area than the target domain as I said above which will be important for the CREATEBC step.

---

<div class="post-metadata">

**Author:** ![sonyafiddes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/sonyafiddes/32/603_2.png) [@sonyafiddes](https://forum.access-hive.org.au/u/sonyafiddes)\
**Post date:** [10 May 2023 01:36 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/12 "2023-05-10T01:36:03Z")

</div>

Yep - I’ve got all that set I think correctly?

 ![Screenshot 2023-05-10 at 11.35.09 am](https://us1.discourse-cdn.com/flex020/uploads/access1/original/1X/ae42dd57bfafc3f7b89f3970776d86ac7c1eaadf.png)

---

<div class="post-metadata">

**Author:** ![sonyafiddes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/sonyafiddes/32/603_2.png) [@sonyafiddes](https://forum.access-hive.org.au/u/sonyafiddes)\
**Post date:** [10 May 2023 01:36 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/13 "2023-05-10T01:36:25Z")

</div>

Oh wait - should it be EC-forecast?

---

<div class="post-metadata">

**Author:** ![Scott](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/scott/32/37_2.png) [@Scott](https://forum.access-hive.org.au/u/Scott)\
**Post date:** [10 May 2023 01:40 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/14 "2023-05-10T01:40:44Z")

</div>

I couldn’t say sorry, will depend on your personal configuration. You can always check the grid namelist in the ancillary directory

---

<div class="post-metadata">

**Author:** ![sonyafiddes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/sonyafiddes/32/603_2.png) [@sonyafiddes](https://forum.access-hive.org.au/u/sonyafiddes)\
**Post date:** [10 May 2023 06:20 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/15 "2023-05-10T06:20:03Z")

</div>

Ok I’ve worked out a pretty hacky way of making it run: I’m running the ancillary suite twice - once with a larger grid (220 points) and then with a smaller grid (200 points). In the nested suite I then point the driving model dm\_ec\_lam\_ancil\_dir to the larger grid and the rg01\_rs01\_ancil\_dir to the smaller grid. I’ll see if that solves Zhangchengs problem as well closer to the pole!

---

<div class="post-metadata">

**Author:** ![Scott](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/scott/32/37_2.png) [@Scott](https://forum.access-hive.org.au/u/Scott)\
**Post date:** [10 May 2023 06:31 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/16 "2023-05-10T06:31:23Z")

</div>

You can also increase the number of domains in the ancillary suite to generate them both at the same time

---

<div class="post-metadata">

**Author:** ![sonyafiddes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/sonyafiddes/32/603_2.png) [@sonyafiddes](https://forum.access-hive.org.au/u/sonyafiddes)\
**Post date:** [10 May 2023 06:45 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/17 "2023-05-10T06:45:52Z")

</div>

Oh that would have been much smarter

---

<div class="post-metadata">

**Author:** ![sonyafiddes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/sonyafiddes/32/603_2.png) [@sonyafiddes](https://forum.access-hive.org.au/u/sonyafiddes)\
**Post date:** [17 May 2023 06:38 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/18 "2023-05-17T06:38:58Z")

</div>

> [@Scott](#):
>
> Hi Sonya, what are the boundaries of both your lam domain and of the ERA5 grib data? There is an option in the ERA5grib settings to include all longitudes which may be needed if working at the date line or near the poles

HI Scott - can you elaborate a little more on this ERA5grib setting, for near the poles?

---

<div class="post-metadata">

**Author:** ![Scott](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/scott/32/37_2.png) [@Scott](https://forum.access-hive.org.au/u/Scott)\
**Post date:** [17 May 2023 06:49 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/19 "2023-05-17T06:49:00Z")

</div>

In the `nci_era5grib` app set `POLAR=True`. With this set rather than selecting just the LAM region to convert to GRIB format the tool will restrict only by latitude.

---

<div class="post-metadata">

**Author:** ![sonyafiddes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/sonyafiddes/32/603_2.png) [@sonyafiddes](https://forum.access-hive.org.au/u/sonyafiddes)\
**Post date:** [18 May 2023 01:28 UTC](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745/20 "2023-05-18T01:28:58Z")

</div>

Thanks for that Scott - a few more questions (for now at least) if you or anyone else has the time.

On the nesting suite wiki ([https://git.nci.org.au/bom/ngm/documentation/-/wikis/Supported-Suites#regional-nesting-suite](https://git.nci.org.au/bom/ngm/documentation/-/wikis/Supported-Suites#regional-nesting-suite)) it explicitly says to ensure the ERA5 is not on a rotated pole. I did try to rotate it anyway (for a domain over Davis) and the resulting grib files are definitely garbage. Can you tell me if there is a work around for this, so that we could use ERA5 on a rotated pole? I’ve tried a few things so far, but not having much luck.

I’ve also tried to run the nesting suite driving model using initial conditions from one of my ACCESS runs, but it runs into problems with I think something to do with the land surface, which may not be that surprising given the different land schemes?

? Error from routine: RCF\_CALC\_TILES  
? Error message: Try two-step recon: Change in grid for n-to-1 not currently supported for STASHcode 231

So I guess I want to know if there is either a way to use ACCESS runs for the initial conditions without running into these problems (I’m not sure how to do a two-step recon a suggeted in error?) or if there is a UM global configuration suite that I can set up to generate the initial conditions?

Lastly, does anyone have advice about running the nested suite freely or not? Is it much like the global model in the sense that there are pros and cons for both and it just depends on what science you’re doing, or is there a good reason for one over the other?

[Next page](https://forum.access-hive.org.au/t/um-nesting-suite-lam-inclusion-grids-not-aligning/745.md?page=2)
