How to deal with external functions in function "daefun"

[GaglioSSM]

Hi everyone, I am using OpenOCL in the framework of atmospheric re-entry trajectory optimization. Firstly, I would like to thank you for this very powerful tool, which allows the modelling of a wide range of problems. In this regard, I would like to ask you an opinion concerning an issue in my problem I am not able to solve.
One of my ODE (in function daefun) has got a dependence on atmospheric density which is a function of altitude, longitude, latitude, and time (all state functions in my problem). I would like to use "atmosnrlmsise00" (a quite complex function already implemented in Matlab). However, I am meeting some troubles as it requires numbers as inputs and is not compatible with the structure of state variables x.r (the altitude), x.phi (the latitude) etc. Could you suggest me a way to overcome this issue? Maybe to convert the variables to their correspondent values.
Many thanks