Locked History Actions

attachment:sm.kumac of Programs Memo

Attachment 'sm.kumac'

Download

 macro main
   
* **************************************
* **************************************
* Modify according to your environment
   
   leveldir = kumac/sm/levels
   
   file1 = 26Al.lpt
   file2 = 26Mg.lpt
   file3 = 26Na.lpt
   file4 = 26Ne.lpt
   file5 = 26F.lpt
   file6 = 26O.lpt
   
   nucl1 = '^26!Al'
   nucl2 = '^26!Mg'
   nucl3 = '^26!Na'
   nucl4 = '^26!Ne'
   nucl5 = '^26!F'
   nucl6 = '^26!O'
   
* ************************************** 
* **************************************
   
   maxj = 9
   maxjp = [maxj] + 2
   
   opt ngri
   opt ndat
   set xtic 0.001
   set ndvy 515
   set yval 999
   set lwid 3
   set xlab 0.8
   
   do i = 1, 100
     if ([file[i]] .eq. '[file[i]]') then
       nuclides = [i] - 1
       breakl
     endif
   enddo
   
   null -5 $sigma([nuclides]*[maxjp]-1) -0.99 10 
   atit ! 'E?x! (MeV)' ! 220
   
   do k=1,[nuclides]
     
     sh rm -f [file[k]].part
     sh awk '''NR>8{print $5,$6,substr($2,4,1) }''' [leveldir]/[file[k]] > [leveldir]/[file[k]].part
     
     vex = $vexist(j[k])
     if ([vex] .eq. 1) then
       vec/del lvl[k]
       vec/del j[k]
       vec/del pi[k]
     endif
     
     vec/read lvl[k],j[k],pi[k] [leveldir]/[file[k]].part
     lvl0_[k] = $sigma(lvl[k](1))
     sigma lvl[k] = lvl[k]-(lvl[k](1))
     sigma j[k] = j[k] + $sigma(([k]-1)*[maxjp])
     xoffset[k]=$sigma(([k]-1)*[maxjp])-1
     vlen = $vlen(j[k],1)
     
     do i=1,[vlen]
       if ($sigma(lvl[k]([i])) .lt. 10) then
         dline [xoffset[k]] $sigma(j[k]([i])) $sigma(lvl[k]([i])) $sigma(lvl[k]([i]))
       endif
     enddo  
     itx $sigma(([k]-1)*[maxjp]+1) -0.7 [nucl[k]]
     itx $sigma(([k]-1)*[maxjp]+1) -0.3 $format([lvl0_[k]],F.2)
   enddo
 return

Attached Files

To refer to attachments on a page, use attachment:filename, as shown below in the list of files. Do NOT use the URL of the [get] link, since this is subject to change and can break easily.
  • [get | view] (2011-11-12 13:19:19, 14.4 KB) [[attachment:a26_sm.png]]
  • [get | view] (2011-11-12 13:20:04, 1.7 KB) [[attachment:sm.kumac]]

You are not allowed to attach a file to this page.