Move the call to clcdrag up from coefkz to clmain (folllowing LMDZ). As both clcdrag and coefkz need zgeop, also move the computation of zgeop from coefkz to clmain.