If using dye with nesting grids, some of the active tracers may not be written out properly, even with the corresponding switches turn on.
The bug comes from miss using Ltracer under case(‘Hout(idTvar)’) in read_phypar.F. If turn on dye, the size of Ltracer is larger than NAT*Ngrids. However, when calling load_l, the codes in function load_l assumes Ltracer has a size of NAT*Ngrids.
My currently solution is Npts=load_l(Nval,Cval,NAT*Ngrids,Hout(idTvar(1):idTvar(NAT),:)).
Similar solution may apply to Aout.
My version is not up to date, but I think it is still there.
Yan