Hello, I need help with the code 139 error mentioned in this post. Apologies, but I couldn’t figure out how to reply directly to that post.
Unfortunately, neither payu sweep nor perturbing the model is resolving the issue. I made some changes to vegetation types in this run, such as removing tundra and adding ice.
The run directory is located at: /home/561/hs2549/access-esm1.5/lgm."
I have made the vegetation changes. My target restart file, after running the update_cable_vegfrac.py script (from /g/data/y99/hs2549/access-esm/scripts/), is located at archive/restart010/atmosphere/test2.astart.
This file differs from restart_dump.astart in the same folder because I have run the perturbation script on it after applying the vegetation updates.
Thanks Himadri,
My guess is that the greatly expanded ice PFT type requires some other land surface variables to be changed in combination to make it stable. I guess the script above worked for the 49 ka case and also the Last Interglacial?
Not exactly sure why it worked in those cases and not this one, but it’s definitely a problem in the Miocene case if I just change the PFT (especially for ice I think), without also changing other land surface variables to be ‘in sync’ with it.
I had a conversation with @mlipson about this… he mentioned that ice gets treated very differently to other PFTs in CABLE. Maybe he could elaborate on this a bit more?
In any case, I think the easiest solution is to remap all of the land surface variables from a working analogue… so that ice PFT cells have corresponding land surface restart variables that have already been working. I can help set that up.
Hi again @HIMADRI_SAINI,
Just looking at your old restart file vs your new restart file. In your restart folder: /home/561/hs2549/access-esm1.5/lgm/archive/restart010/atmosphere
There are files: restart_dump.astart.bak restart_dump.astart
I looked at the PFT variable in each of these files, and it appears that the land-sea mask has changed slightly. The new restart file has more NaNs in the PFT variable around the edge of Australia and New Guinea.
I made netcdfs of these in this folder: /g/data/y99/dkh157/access_esm/LGM/veg_pert
The files old_restart.pft.nc and new_restart.pft.nc are the PFTs in your restart dir, and: mask_diff.nc
Is a difference of the two masks. This difference files shows a smattering of non-zeros, so I think in the first instance the run is crashing because there are NaNs in the PFT restart variable.
Update: I have kept the land-sea mask in my new restart consistent with the old one, but I am still encountering the same error. I would be happy to try your earlier suggestion of remapping all the land surface variables.
Hi @spencerwong thanks for checking in. No, I am not running into this issue anymore. In my simulation, large extent of the Laurentide ice sheet was leading to a strong decrease in salinity near Greenland, which led to the model crashing after a few months. This has been resolved by adjusting the runoff to match my modified topography.