In procedure cdrag test zri >= 0 rather than zri > 0 and invert the order of the if construction (following LMDZ).