15 ngc( :) = (/8, 14, 16, 14, 16, 8, 12, 8, 12, 6, 8, 8, 4, 2, 2, 2 /)
 
   18              &8,  22,  38,  52,  68,  76,  88,  96, &
 
   19            &108, 114, 122, 130, 134, 136, 138, 140/)
 
   21            &1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,             &
 
   22            &1,2,3,4,5,6,7,8,9,10,11,12,13,13,14,14,      &
 
   23            &1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,      &
 
   24            &1,2,3,4,5,6,7,8,9,10,11,12,13,14,14,14,      &
 
   25            &1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,      &
 
   26            &1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,             &
 
   27            &1,1,2,2,3,4,5,6,7,8,9,10,11,11,12,12,        &
 
   28            &1,1,2,2,3,3,4,4,5,5,6,6,7,7,8,8,             &
 
   29            &1,2,3,4,5,6,7,8,9,9,10,10,11,11,12,12,       &
 
   30            &1,1,2,2,3,3,4,4,5,5,5,5,6,6,6,6,             &
 
   31            &1,2,3,3,4,4,5,5,6,6,7,7,7,8,8,8,             &
 
   32            &1,2,3,4,5,5,6,6,7,7,7,7,8,8,8,8,             &
 
   33            &1,1,1,2,2,2,3,3,3,3,4,4,4,4,4,4,             &
 
   34            &1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,             &
 
   35            &1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2,             &
 
   36            &1,1,1,1,1,1,1,1,2,2,2,2,2,2,2,2/)           
 
   40            &1,1,1,1,1,1,1,1,1,1,1,1,2,2,                 &
 
   41            &1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,             &
 
   42            &1,1,1,1,1,1,1,1,1,1,1,1,1,3,                 &
 
   43            &1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,             &
 
   45            &2,2,1,1,1,1,1,1,1,1,2,2,                     &
 
   47            &1,1,1,1,1,1,1,1,2,2,2,2,                     &
 
   58             &2,2,2,2,2,2,2,2,2,2,2,2,2,2,                &
 
   59             &3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,            &
 
   60             &4,4,4,4,4,4,4,4,4,4,4,4,4,4,                &
 
   61             &5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,            &
 
   63             &7,7,7,7,7,7,7,7,7,7,7,7,                    &
 
   65             &9,9,9,9,9,9,9,9,9,9,9,9,                    &
 
   67            &11,11,11,11,11,11,11,11,                     &
 
   68            &12,12,12,12,12,12,12,12,                     &
 
   75      &0.1527534276_jprb,0.1491729617_jprb,0.1420961469_jprb,0.1316886544_jprb,&
 
   76      &0.1181945205_jprb,0.1019300893_jprb,0.0832767040_jprb,0.0626720116_jprb,&
 
   77      &0.0424925_jprb   ,0.0046269894_jprb,0.0038279891_jprb,0.0030260086_jprb,&
 
   78      &0.0022199750_jprb,0.0014140010_jprb,0.000533_jprb    ,0.000075_jprb    /)
 
integer(kind=jpim), dimension(jpgpt) ngb
 
integer(kind=jpim), dimension(jpg *jpband) ngm
 
real(kind=jprb), dimension(jpg) wt
 
integer(kind=jpim), dimension(jpgpt) ngn
 
integer(kind=jpim), dimension(jpband) ngc
 
integer(kind=jpim), dimension(jpband) ngs