implement GCL

You may need to add a initialization call before the loop in wake_model, that is overridden in GLC to reshape the deficit with an extra dimension for the 40 evaluation points