I have a vn12 nested job on ARCHER2 (u-cs690), which has domains at 12, 4, 1.5, and 0.5 km over a region including the South Orkney Islands in the Southern Ocean.
There is an error making the vegfrac / soil ancillaries at 4, 1.5 and 0.5 km. The error output for the 4 km domain is at: /home/n02/n02/anmcr/cylc-run/u-cs690/log/job/20000101T0000Z/south_orkney_islands_4km_ancil_vegfrac/01/job.out. I have also copied some of the error message below. There seems to be unresolved points in the soil dataset, and actually the qrparm.soil ancillary is not being created. However, these ancillaries seem to be created for the 12 km domain.
Are you able to please suggest how to resolve this issue?
Thanks,
Andrew
UNRESOLVED GRID POINTS IN SOIL DATASET
Number of points unresolved is 35
POINT 4743 LAT -60.5267 LONG 314.1646
POINT 4744 LAT -60.5270 LONG 314.2378
POINT 4841 LAT -60.5620 LONG 314.0175
POINT 4842 LAT -60.5624 LONG 314.0908
POINT 4843 LAT -60.5627 LONG 314.1640
POINT 4844 LAT -60.5630 LONG 314.2372
POINT 4845 LAT -60.5633 LONG 314.3105
POINT 4846 LAT -60.5635 LONG 314.3837
POINT 4847 LAT -60.5637 LONG 314.4570
POINT 4848 LAT -60.5638 LONG 314.5302
POINT 4941 LAT -60.5980 LONG 314.0167
POINT 4942 LAT -60.5984 LONG 314.0900
POINT 4944 LAT -60.5990 LONG 314.2367
POINT 4945 LAT -60.5993 LONG 314.3100
POINT 4946 LAT -60.5995 LONG 314.3833
POINT 4947 LAT -60.5997 LONG 314.4567
POINT 4948 LAT -60.5998 LONG 314.5300
POINT 5045 LAT -60.6353 LONG 314.3095
POINT 5046 LAT -60.6355 LONG 314.3829
POINT 5047 LAT -60.6357 LONG 314.4563
POINT 5048 LAT -60.6358 LONG 314.5298
POINT 5049 LAT -60.6359 LONG 314.6032
POINT 5050 LAT -60.6360 LONG 314.6766
POINT 5051 LAT -60.6360 LONG 314.7500
POINT 5150 LAT -60.6720 LONG 314.6765
POINT 5151 LAT -60.6720 LONG 314.7500
POINT 5154 LAT -60.6718 LONG 314.9705
POINT 5251 LAT -60.7080 LONG 314.7500
POINT 5254 LAT -60.7078 LONG 314.9707
POINT 5259 LAT -60.7067 LONG 315.3386
POINT 5261 LAT -60.7060 LONG 315.4858
POINT 5352 LAT -60.7440 LONG 314.8237
POINT 5358 LAT -60.7430 LONG 315.2656
POINT 5359 LAT -60.7427 LONG 315.3393
POINT 5360 LAT -60.7424 LONG 315.4129
Search radius 1
NO DATA FROM WHICH TO SET UNRESOLVED POINTS
***ERROR: No source data available in target domain
Just to update you that I’m trying to use the RAS (Regional Ancillary Suite) on ARCHER2 to compute the ancillaries. My job is u-cs728, which is a copy of Ella Gilbert’s RAS job for a pan-Antarctic domain (her job is u-cs542).
I ran 12 and 4km RAS suites over S Orkney - vegetation fraction worked, but there is no vegetation. The spiral search failed for soils at 12 km and 4km. Are there any soils on S Orkney?
(u-cs690 didn’t have multiple resolutions when I looked)
Thank you for replying and investigating this. Can I copy your RAS suites? I also tried to run these (taking copies from both Ella Gilbert and Helen Burns, but both suite failed.)
I switched off the multiple domains in u-cs690 as my aim was to get the 12 km domain working in the first instance. I don’t think there is any soil over S. Orkney, and any soil likely to be frozen. I think the main consideration would be the surface type, and how that impacts surface roughness.
I simply ran your u-cs728 - as I mentioned, the spiral search failed for soils, (and lai, and topographic index.) That’s maybe not a surprise given the domain.
Thanks for the update. That was also what I got (the failure for some ancils).
I used the ancils that were created in by model suite u-cs690 but its failing with the error ‘NaNs in error term in BiCGstab’. I also switched on the advanced ‘print’ options. Is this likely an issue with the ancillaries?
Thanks,
Andrew
Error code: 1
? Error from routine: EG_BICGSTAB
? Error message: NaNs in error term in BiCGstab after 1 iterations
? This is a common point for the model to fail if it
? has ingested or developed NaNs or infinities
? elsewhere in the code.
? See the following URL for more information:
? Met Office Science Repository Service
? Error from processor: 17
? Error number: 47
As you mentioned, the RAS suite failed to make a soil ancillary (qrparm.soil). Am I able to run my model without including this in the reconfiguration? Is one option to set this to zero? I am aware that there are issues with the soil source file for latitudes higher than 60S.
I copied all the ancillaries made by the RAS suite to /work/n02/n02/anmcr/ancils/south_orkney_islands/12km. I also copied a qrparm.soil ancillary to this location, that I made using the ancillary generation of the u-cs690 suite. I will see if this helps.
I’ve gone back to trying to make the ancillaries. Ella Gilbert has also been working on this for Antarctica, so I took a copy of her recent RAS suite. This is u-cs864. This fails for ‘ancil_soils_hydr’, with an error that states ‘The provided source doesn’t appear to have any valid data.’ See also below? Are you able to please advise?
I looked at app/ancil_soils_hydr.conf , which gives the source of the input data as source=${ANCIL_PREPROC_PATH}/soils_hydrology.nc. ANCIL_PREPROC_PATH = $ROSE_DATA/etc/ants_preproc. But I’m not sure what $ROSE_DATA refers to?
Thanks,
Andrew
The following have been reloaded with a version change:
cce/11.0.4 => cce/12.0.3
/opt/conda/envs/ants/lib/python3.7/site-packages/mule-2020.1.1-py3.7.egg/mule/stashmaster.py:284: UserWarning: Ancillary files do not define the UM version number in the Fixed Length Header. No STASHmaster file loaded: Fields will not have STASH entries attached.
warnings.warn(msg)
Traceback (most recent call last):
File “/work/n02/n02/anmcr/cylc-run/u-cs864/share/RAS_contrib/SoilParameters/ancil_soils.py”, line 789, in
main(args.sources, args.lct_ancillary, args.soils_lookup, args.output)
File “/work/n02/n02/anmcr/cylc-run/u-cs864/share/RAS_contrib/SoilParameters/ancil_soils.py”, line 768, in main
make_consistent_lct(parameters, lct_fractions)
File “/work/n02/n02/anmcr/cylc-run/u-cs864/share/RAS_contrib/SoilParameters/ancil_soils.py”, line 651, in make_consistent_lct
parameters[0], target_mask=target_mask_noice
File “/opt/conda/envs/ants/lib/python3.7/site-packages/ants/analysis/_merge.py”, line 865, in init
self._call_spiral_search(source)
File “/opt/conda/envs/ants/lib/python3.7/site-packages/ants/analysis/_merge.py”, line 920, in _call_spiral_search
raise ValueError(msg)
ValueError: The provided source doesn’t appear to have any valid data.
[FAIL] python_env ${CONTRIB_PATH}/SoilParameters/ancil_soils.py ${source}
[FAIL] --lct-ancillary ${vegfrac} --soils-lookup ${soils_lookup}
[FAIL] -o ${output} --ants-config ${config} <<‘STDIN’
[FAIL]
[FAIL] ‘STDIN’ # return-code=1
2022-12-14T16:23:47Z CRITICAL - failed/EXIT
I should add that my suite now runs fine for the Arctic but not over the Antarctic - it also fails on the soil / veg frac steps , which I imagine is something to do with the Antarctic having no soil or vegetation in the source dataset.