diff --git a/org/qmckl_mo.org b/org/qmckl_mo.org
index 66c7b33..b923f24 100644
--- a/org/qmckl_mo.org
+++ b/org/qmckl_mo.org
@@ -768,17 +768,17 @@ assert (rc == QMCKL_SUCCESS);
 rc = qmckl_set_mo_basis_type(context, typ);
 assert (rc == QMCKL_SUCCESS);
 
-const int64_t mo_num = chbrclf_ao_num;
+const int64_t mo_num = chbrclf_mo_num;
 rc = qmckl_set_mo_basis_mo_num(context, mo_num);
 assert (rc == QMCKL_SUCCESS);
 
-double  mo_coefficient[chbrclf_ao_num][mo_num];
-const int64_t min_num = (((chbrclf_ao_num) < (mo_num)) ? (chbrclf_ao_num) : (mo_num));
-for (int i=0; i<min_num; ++i) {
-  mo_coefficient[i][i] = 1.0;
-}
+const double  * mo_coefficient          =  &(chbrclf_mo_coef[0]);
+//const int64_t min_num = (((chbrclf_ao_num) < (mo_num)) ? (chbrclf_ao_num) : (mo_num));
+//for (int i=0; i<min_num; ++i) {
+//  mo_coefficient[i][i] = 1.0;
+//}
 
-rc = qmckl_set_mo_basis_coefficient(context, &(mo_coefficient[0][0]));
+rc = qmckl_set_mo_basis_coefficient(context, mo_coefficient);
 assert (rc == QMCKL_SUCCESS);
 
 assert(qmckl_mo_basis_provided(context));
@@ -787,6 +787,13 @@ double mo_vgl[5][walk_num][elec_num][chbrclf_ao_num];
 rc = qmckl_get_mo_basis_vgl(context, &(mo_vgl[0][0][0][0]));
 assert (rc == QMCKL_SUCCESS);
 
+//double ovl = 0.0;
+//for (int i=0; i<ao_num; ++i) {
+//for (int i=0; i<mo_num; ++i) {
+//  mo_coefficient[i][i] = 1.0;
+//}
+//}
+
 printf("\n");
 printf(" mo_vgl mo_vgl[0][0][26][219] %25.15e\n", mo_vgl[0][0][2][3]);
 printf(" mo_vgl mo_vgl[1][0][26][219] %25.15e\n", mo_vgl[1][0][2][3]);