r_lwfbrook90 function

Interface function to the LWF-Brook90 model