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