Hello there,
starting from Berkhoff's equation, it is possible to derive a new form from it that takes into account the presence of currents (Liu 1990).Nevertheless this is not a trivial procedure, and as far as I know, there is not yet a fortran subroutine that modifies the ARTEMIS module so as to consider currents (in fact, I think of this as a new release of the ARTEMIS featuring the 'currents' form,rather than a subroutine).
Regarding the wind, you should bear in mind that ARTEMIS is a phase-averaged model. Historically, although attemps have been made to include generation or decay, the results have been disappointing, and therefore not included in this kind of models (those based on Berkhoff's equation). During the last years it seems that there have been some success including wind as forcing to a phase-averaged model, but for Boussinesq-type models only (see Chen 2004).
To sum up: currents can be added, but it seems to be more a revision of the whole ARTEMIS module rather than a single subroutine (if I am wrong, then it would be very thoughtful to share that code!). With regards to the inclusion of winds...I wouldn't bet for that.
Hope that it helps.
Ruben.