problem: iso_massf_decay in H5_surf does not work properly? Solution 1: Use the python decay routines Solution 2: may iso_massf_decay is ok to use? How to check? need solution 1. If you look at other AGB models (attachment) there is a lot of Li production at low-mass end: Factor 10 - 1e3 in those stars: 1Msun, Z=0.01,Z=0.001,0.0001 1.65Msun Z=0.006,0.001,0.0001 2Msun Z=0.006,0.001,0.0001 4Msun Z=0.01, 0.02,