macro shift_watch [IDtemp]
   
   global/create ID
   global/create scale
   global/create histMax
   global/create dir

   
   histMax = 12
   dir = ./ps/shift
   
   if [IDtemp] .eq. '' then
     message
     message please input run NO !
     message
     exitm
   endif
   
   ID = [IDtemp]
         
   
*   lnz
*   scale = 'linear'
*   exec main 
*   
*   wait
   
   lgz
   scale = 'log'
   exec main 
   
   mes  'print out ?'
   wait
   
*   scale = 'linear'
*   do i=1,[histMax]
*     sh pwd
*     sh ps2pdf [dir][ID]_[i]_[scale].ps [dir][ID]_[i]_[scale].pdf
*     sh lpr -Pbripslw2 [dir][ID]_[i]_[scale].pdf
** sh lpr [dir][ID]_[i]_[scale].pdf
*   enddo
   
   scale = 'log'
   do i=1,[histMax]
     sh pwd
     sh ps2pdf [dir][ID]_[i]_[scale].ps [dir][ID]_[i]_[scale].pdf
     sh lpr -Pbripslw2 [dir][ID]_[i]_[scale].pdf
* sh lpr [dir][ID]_[i]_[scale].pdf
   enddo

 return  
 
 
 macro main 
   
   import *
   

   
   size 18 24
   
   
   opt date
   opt sta
   zone 2 3
   
   
   set ygti 1.0
   set gsize 0.6
     
      
*  1
   num = 1

   pson [dir][ID]_[num]_[scale].ps   
   title run[ID]
   
   
   ht 101 colz
   hn colz 
   hn colz
   
   do  i = 1,3
     hn colz
     do j = 1,6
       exec dlx [j]+0.5
     enddo
   enddo
   
   do i = 4,6
     sel [i]*10
     itx 1.8 2500 F3
     itx 2.65 2500 F5-1
     itx 3.65 2500 F5-2
     itx 4.8 2500 F7
     itx 5.8 2500 F8
     itx 6.7 2500 F11
   enddo
   
   sel 1
   itx 2 23 shift[num]
   itx 3.5 23 [scale]
   
  
   psoff

   num = [num] + 1
   
   wait
   
*  2   
   pson [dir][ID]_[num]_[scale].ps   
   title run[ID]

   hn colz
   do j = 1,6
     exec dlx [j]+0.5
   enddo
   
   
   itx 1.8 2500 F3
   itx 2.65 2500 F5-1
   itx 3.65 2500 F5-2
   itx 4.8 2500 F7
   itx 5.8 2500 F8
   itx 6.7 2500 F11
   
   
   hn colz 
   
   opt nsta
   zone 1 3 2 s
   hn colz
   exec dlx 2.5   
   exec dlx 4.5   
   exec dlx 8.5   
   exec dlx 9.5   
   exec dlx 13.5   
   exec dlx 14.5   
   exec dlx 18.5   
   exec dlx 22.5   
   exec dlx 26.5   
   exec dlx 30.5   

   hn colz
   exec dlx 2.5   
   exec dlx 4.5   
   exec dlx 8.5   
   exec dlx 9.5   
   exec dlx 13.5   
   exec dlx 14.5   
   exec dlx 18.5   
   exec dlx 22.5   
   exec dlx 26.5   
   exec dlx 30.5   

   do i = 3,4
     sel [i]*10
     itx 1.2 4500 F1
     itx 3 4500 F2
     itx 6.1 4500 F3
     itx 8.6 4500 F4
     itx 11.2 4500 F5
     itx 13.7 4500 F6
     itx 16.2 4500 F7
     itx 20.2 4500 F8
     itx 24.2 4500 F9
     itx 27.7 4500 F10
     itx 32.2 4500 F11
   enddo
   
   sel 1
   itx 2 23 shift[num]
   itx 3.5 23 [scale]

   
   psoff
   wait

   num = [num] + 1
   
*  3
   pson [dir][ID]_[num]_[scale].ps   
   title run[ID]
   
   zone 1 3
   do i = 1,3
     hn colz
     exec dlx 2.5   
     exec dlx 4.5   
     exec dlx 8.5   
     exec dlx 9.5   
     exec dlx 13.5   
     exec dlx 14.5   
     exec dlx 18.5   
     exec dlx 22.5   
     exec dlx 26.5   
     exec dlx 30.5   
   enddo
   
   do i = 1,3
     sel [i]*10
     itx 1.2 4500 F1
     itx 3 4500 F2
     itx 6.1 4500 F3
     itx 8.6 4500 F4
     itx 11.2 4500 F5
     itx 13.7 4500 F6
     itx 16.2 4500 F7
     itx 20.2 4500 F8
     itx 24.2 4500 F9
     itx 27.7 4500 F10
     itx 32.2 4500 F11
   enddo
   
   sel 1
   itx 2 23 shift[num]
   itx 3.5 23 [scale]

   
   psoff
   
   opt sta
   
   num = [num] + 1

   
***
   zone 2 3
   do j = [num],[histMax]
     wait
     pson [dir][ID]_[j]_[scale].ps
     title run[ID]
     do i=1,6
       hn colz
     enddo
     
     sel 1
     itx 2 23 shift[j]
     itx 3.5 23 [scale]
          
     psoff  
   enddo
   
   
   
   return
   
   
 
macro dly

  xmin = $GRAFINFO('WNXMIN')
  xmax = $GRAFINFO('WNXMAX')

  if $OPTION('LOGX') = 1 then
    xmin = $sigma(10**([xmin]))
    xmax = $sigma(10**([xmax]))
  endif

  do i = 1,[#]
    dline [xmin] [xmax] [[i]] [[i]]
  enddo

 return
 
 
macro dlx

  ymin = $GRAFINFO('WNYMIN')
  ymax = $GRAFINFO('WNYMAX')

  if $OPTION('LOGY') = 1 then
    ymin = $sigma(10**([ymin]))
    ymax = $sigma(10**([ymax]))
  endif

  do i = 1,[#]
    dline [[i]] [[i]] [ymin] [ymax]
  enddo

return