Do the same calculation for richardson number and f_cdrag between coefcdrag and clcdrag (following revision 2011 of LMDZ). We have then to add input variables f_cdrag_ter and f_cdrag_oce in module clesphys (following revision 1168 of LMDZ), with default values those used for AR4.