# Modifying a UM restart file for ACCESS-ESM1.5

**URL:** <https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537>\
**Category:** Technical\
**Tags:** help, land, solved\
**Created:** [20 October 2023 06:04 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537 "2023-10-20T06:04:49Z")\
**Posts on this page:** 20\
**Page:** 2

<div class="post-metadata">

**Author:** ![RachelLaw](https://avatars.discourse-cdn.com/v4/letter/r/da6949/32.png) [@RachelLaw](https://forum.access-hive.org.au/u/RachelLaw)\
**Post date:** [15 December 2023 04:08 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/23 "2023-12-15T04:08:58Z")

</div>

Hi @ben , quick look at the file seems good. It looks like I’ve finally got test runs going so I’ll hopefully start to have some answers early next week.  
Thanks for your ongoing help with this.

---

<div class="post-metadata">

**Author:** ![RachelLaw](https://avatars.discourse-cdn.com/v4/letter/r/da6949/32.png) [@RachelLaw](https://forum.access-hive.org.au/u/RachelLaw)\
**Post date:** [18 December 2023 04:14 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/24 "2023-12-18T04:14:22Z")

</div>

Hi @ben , test runs have given me lots of output to play with. Initial look makes me think I probably didn’t provide the best instructions for averaging the carbon pools in the .astart file. I suspect we need to account for tile fraction, which I’m guessing hasn’t been included. Rachel

---

<div class="post-metadata">

**Author:** ![RachelLaw](https://avatars.discourse-cdn.com/v4/letter/r/da6949/32.png) [@RachelLaw](https://forum.access-hive.org.au/u/RachelLaw)\
**Post date:** [18 December 2023 22:33 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/25 "2023-12-18T22:33:49Z")

</div>

Hi @ben , I’ve tested weighting by vegetation fraction on the wood carbon pool (field 853). There is an example netcdf file here: /g/data/p66/rml599/LandUse-tests/cwood.mean5.nc and the code I used to generate it: /g/data/p66/rml599/LandUse-tests/avTreePFTs.f90 . Please could you make me some new astart files with weighted averages - or show me how to make them. Thanks.

---

<div class="post-metadata">

**Author:** ![ben](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ben/32/282_2.png) [@ben](https://forum.access-hive.org.au/u/ben)\
**Post date:** [18 December 2023 22:38 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/26 "2023-12-18T22:38:20Z")

</div>

Sure thing @RachelLaw, I’ll make the new weighted file and post my code to GitHub.

Standby.

---

<div class="post-metadata">

**Author:** ![ben](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ben/32/282_2.png) [@ben](https://forum.access-hive.org.au/u/ben)\
**Post date:** [18 December 2023 23:11 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/27 "2023-12-18T23:11:47Z")

</div>

Hi @RachelLaw , here is the code that generates the files as is.

> <https://github.com/CABLE-LSM/gists/blob/main/20231219-RachelLaw.ipynb>

(Best to go directly to GitHub to see what it does)

I’ll work on adding your changes for weighting…

---

<div class="post-metadata">

**Author:** ![ben](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ben/32/282_2.png) [@ben](https://forum.access-hive.org.au/u/ben)\
**Post date:** [17 January 2024 00:54 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/28 "2024-01-17T00:54:06Z")

</div>

Hi @RachelLaw ,

As per our discussions yesterday, please find a sample of the fields weighted by vegetation fraction ([Met Office Codes Registry : um/stash/\_m01s00i216](https://reference.metoffice.gov.uk/um/stash/_m01s00i216)) at this location. This is the field you used in your example.

/g/data/p66/bjs581/mean5\_VFW\_851-860\_20240117\_004324.astart

If you could take a look and let me know if this satisfies your requirements that would be great.

If you have any issues with the file, please feel free to email me directly.

Have a great day,

Ben

---

<div class="post-metadata">

**Author:** ![RachelLaw](https://avatars.discourse-cdn.com/v4/letter/r/da6949/32.png) [@RachelLaw](https://forum.access-hive.org.au/u/RachelLaw)\
**Post date:** [17 January 2024 05:17 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/29 "2024-01-17T05:17:56Z")

</div>

Hi @ben ,  
Looks like we are getting close. The only thing I noticed was that areas with no trees now look like they are missing - possibly filled with NaNs e.g. Antarctica. Pools should just be zero in this case.  
Thanks, Rachel

---

<div class="post-metadata">

**Author:** ![ben](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ben/32/282_2.png) [@ben](https://forum.access-hive.org.au/u/ben)\
**Post date:** [17 January 2024 06:44 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/30 "2024-01-17T06:44:19Z")

</div>

Hi Rachel,

I’ve added an explicit step to set the expected 0 values and ensure the fill value is being applied correctly.

New file here: /g/data/p66/bjs581/mean5\_VFW\_851-860\_20240117\_064315.astart

Kind Regards,  
Ben

---

<div class="post-metadata">

**Author:** ![RachelLaw](https://avatars.discourse-cdn.com/v4/letter/r/da6949/32.png) [@RachelLaw](https://forum.access-hive.org.au/u/RachelLaw)\
**Post date:** [18 January 2024 01:15 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/31 "2024-01-18T01:15:50Z")

</div>

Hi @ben  
Thanks - that looks good. I’ll go ahead and set up a model run with it. Please can you now also make me the ‘second step’ file where the vegfrac fields are replaced with vegfrac-dom-tree.nc.  
Regards, Rachel

---

<div class="post-metadata">

**Author:** ![RachelLaw](https://avatars.discourse-cdn.com/v4/letter/r/da6949/32.png) [@RachelLaw](https://forum.access-hive.org.au/u/RachelLaw)\
**Post date:** [22 January 2024 04:18 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/32 "2024-01-22T04:18:07Z")

</div>

hi @ben ,  
The model run with the latest .astart file looks to have performed as expected. It will be interesting to now test the dominant tree case when you are able to make the appropriate .astart (or step me through how to do it).  
Regards, Rachel

---

<div class="post-metadata">

**Author:** ![ben](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ben/32/282_2.png) [@ben](https://forum.access-hive.org.au/u/ben)\
**Post date:** [23 January 2024 03:28 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/33 "2024-01-23T03:28:10Z")

</div>

Will do, I’ll set up new restarts for you.

---

<div class="post-metadata">

**Author:** ![ben](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ben/32/282_2.png) [@ben](https://forum.access-hive.org.au/u/ben)\
**Post date:** [23 January 2024 23:50 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/34 "2024-01-23T23:50:04Z")

</div>

@RachelLaw , here is a restart dropping in your dominant case.  
/g/data/p66/bjs581/mean5\_VFWDOM\_851-860\_20240123\_234638.astart

Please let me know how you go or if there are any issues.

Cheers, Ben

---

<div class="post-metadata">

**Author:** ![RachelLaw](https://avatars.discourse-cdn.com/v4/letter/r/da6949/32.png) [@RachelLaw](https://forum.access-hive.org.au/u/RachelLaw)\
**Post date:** [24 January 2024 01:31 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/35 "2024-01-24T01:31:03Z")

</div>

Hi @ben,  
I had a look at the new .astart file with xconv and it isn’t looking as I expected so maybe my latest request wasn’t very clear.  
I was expecting the 851-860 fields to be the same as in the mean5\_VFW\_ file and the only change in fields would be to the vegetation fraction (fields 216 and 835) to the version with one dominant tree type per grid-cell.  
When I checked the file, the 851-860 fields looked slightly different so I’m wondering whether these are now weighted with the dominant tree type. Possibly another interesting sensitivity test but one I hadn’t considered. Fields 216 and 835 still look like the vegetation fraction for multiple tree types per grid-cell.  
Sorry if I’ve confused things.  
Rachel

---

<div class="post-metadata">

**Author:** ![ben](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ben/32/282_2.png) [@ben](https://forum.access-hive.org.au/u/ben)\
**Post date:** [30 January 2024 23:36 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/36 "2024-01-30T23:36:55Z")

</div>

Hi Rachel,

Can you please check the following file:

/g/data/p66/bjs581/mean5\_VFW\_851-860\_20240130\_233553.astart

It should have the two different options applied…

Cheers, Ben

---

<div class="post-metadata">

**Author:** ![RachelLaw](https://avatars.discourse-cdn.com/v4/letter/r/da6949/32.png) [@RachelLaw](https://forum.access-hive.org.au/u/RachelLaw)\
**Post date:** [31 January 2024 03:27 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/37 "2024-01-31T03:27:38Z")

</div>

Hi @ben,  
Thanks for the update. The 851-860 fields are now what I’d expected (i.e. the same as in mean5\_VFW\_851-860\_20240117\_064315.astart) but the vegetation fraction field isn’t. It should be the same as the vegetation fraction field in mean5\_VF\_8551-860\_20231213\_120638.astart.  
Rachel

---

<div class="post-metadata">

**Author:** ![ben](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ben/32/282_2.png) [@ben](https://forum.access-hive.org.au/u/ben)\
**Post date:** [4 February 2024 10:20 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/38 "2024-02-04T10:20:22Z")

</div>

Hi Rachel,

Please try:  
/g/data/p66/bjs581/mean5\_VFW\_851-860\_20240204\_101944.astart

Cheers, Ben

---

<div class="post-metadata">

**Author:** ![RachelLaw](https://avatars.discourse-cdn.com/v4/letter/r/da6949/32.png) [@RachelLaw](https://forum.access-hive.org.au/u/RachelLaw)\
**Post date:** [4 February 2024 22:47 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/39 "2024-02-04T22:47:58Z")

</div>

hi @ben,  
Looks good so I’ve started the model run with it.  
Thanks, Rachel

---

<div class="post-metadata">

**Author:** ![ben](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ben/32/282_2.png) [@ben](https://forum.access-hive.org.au/u/ben)\
**Post date:** [13 February 2024 03:18 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/40 "2024-02-13T03:18:23Z")

</div>

Hi Rachel,

Good to hear. Here is the latest code that generates your files.

> <https://github.com/CABLE-LSM/gists/blob/main/Rachel-20240116.ipynb>

I’m going to mark this ticket as closed, but please feel free to raise another if you need further assistance.

Cheers, Ben

---

<div class="post-metadata">

**Author:** ![ben](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/ben/32/282_2.png) [@ben](https://forum.access-hive.org.au/u/ben)\
**Post date:** [13 February 2024 03:21 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/41 "2024-02-13T03:21:09Z")

</div>



---

<div class="post-metadata">

**Author:** ![clairecarouge](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/clairecarouge/32/380_2.png) [@clairecarouge](https://forum.access-hive.org.au/u/clairecarouge)\
**Post date:** [23 January 2025 01:12 UTC](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537/42 "2025-01-23T01:12:04Z")

</div>

A post was split to a new topic: [Modify a UM restart file to work with new vegetation distribution](https://forum.access-hive.org.au/t/modify-a-um-restart-file-to-work-with-new-vegetation-distribution/4089)

[Previous page](https://forum.access-hive.org.au/t/modifying-a-um-restart-file-for-access-esm1-5/1537.md?page=1)
