function medsetp(number_of_9set)

 !ls c* > nameUT
 
 !gawk '$1~/c/ {print $1}' nameUT > nameUT1
 
 !mkdir all9
 !cp c*/* all9
 nn=number_of_9set*9;
 cd all9
 pre(nn);
 
 !~/cloud/new/polar

pwd

 load tt

 load hwp

 !cp ../nameUT1 .
 
nodallthwp(tt,hwp)