Use named constant f_ri_cd_min in procedure cdrag (following LMDZ). Remove some intermediary variables.