bin2bin.f File Reference

Go to the source code of this file.

Functions/Subroutines
program	main
subroutine	get1hyb (rew)
subroutine	mergehyb1 (h1)
subroutine	openhyb (icon)

Function/Subroutine Documentation

get1hyb()

subroutine get1hyb

(

logical

rew

)

Definition at line 119 of file bin2bin.f.

References parameter().

      implicit none
      logical rew
      character*128 input0
      integer i, klena

      integer ndepth
      parameter(ndepth= 50)
      integer fn1
      real*8  asdep(ndepth),  munit(ndepth)
      real*8  esize0(ndepth),
     *     age0(ndepth),  cogdep0(ndepth),
     *     seloss0(ndepth),
     *     ng0(ndepth), ne0(ndepth), nmu0(ndepth),
     *     cog0
      integer evno0

      common /zbin2ascii/
     *    asdep, esize0, age0,
     *    cogdep0,  seloss0, munit,
     *    ng0, ne0, nmu0, cog0,
     *    fn1, evno0

      if(rew) rewind fn1

      input0 = "x" 
!////////////
!      write(0,*) ' while'
!////////
      do while (input0(1:10) .ne. " ")
          input0=" "
          read( fn1 ,'(a)') input0
!////////////
!          write(0,*) ' input0=',input0
!//////////////
          if(input0(1:10) .ne. " ")  then
            read(input0(1:klena(input0)), *)
     *        evno0, i, asdep(i), esize0(i), age0(i),
     *        cogdep0(i), seloss0(i),
     *         munit(i), ng0(i), ne0(i), nmu0(i), cog0
!/////////
!          write(0,*) ' input0 read'
!//////////////
          endif
      enddo

Here is the call graph for this function:

main()

program main

(

)

Definition at line 1 of file bin2bin.f.

References copenfw2(), false, kwhistd(), kwhistp(), kwhistr(), kwhists(), kwhistso(), and true.

Here is the call graph for this function:

mergehyb1()

subroutine mergehyb1

(

type(histogram1)

h1

)

Definition at line 166 of file bin2bin.f.

References get1hyb(), and parameter().

      use modhistogram1
      use modhistogram2
      use modhistogram3
      implicit none

      type(histogram1) h1
      type(histogram2) h2
      type(histogram3) h3

      integer ndepth
      integer nc
      parameter(ndepth= 50)
      integer fn1
      real*8  asdep(ndepth),  munit(ndepth)
      real*8  esize0(ndepth),
     *     age0(ndepth),  cogdep0(ndepth),
     *     seloss0(ndepth),
     *     ng0(ndepth), ne0(ndepth), nmu0(ndepth),
     *     cog0
      integer evno0

      common /zbin2ascii/
     *    asdep, esize0, age0,
     *    cogdep0,  seloss0, munit,
     *    ng0, ne0, nmu0, cog0,
     *    fn1, evno0

      integer klena
      integer j

      do while (h1%c%eventno .ne. evno0) 
         call get1hyb( h1%c%eventno .lt. evno0) 
      enddo
!
!       available variables
!     *      idx, ASdep(idx), Esize0(idx), age0(idx), 
!     *      SEloss0(idx), munit(idx), 
!     *      Ng0(idx), Ne0(idx), Nmu0(idx), 
!     *      ASdep(idx)/cog0, cog0
!
!        this part must be consistent with
!        FleshHist/interface.f output  for evid
!//////////////
!      write(0,*) ' id', h1%c%id
!/////////// 

      read(h1%c%id, '(i3)') j
!//////////////
!      write(0,*) ' j=',j
!///////////
      write(h1%c%id, 
     *   '(i3, i5,  f5.2, f5.2,
     *   i5,  i5)')
     *   j, int( asdep(j) ),
     *   age0(j), asdep(j)/cog0,
     *   int(munit(j)), int(cog0)
!//////////////
!      write(0,*) ' j=',j
!///////////
      return
!     *******************
      entry  mergehyb2(h2)
!     *******************
            

      do while (h2%c%eventno .ne. evno0) 
         call get1hyb( h2%c%eventno .lt. evno0) 
      enddo
      read(h2%c%id, '(i3)') j
      write(h2%c%id, 
     *   '(i3, i5,  f5.2, f5.2,
     *   i5,  i5)')
     *   j, int( asdep(j) ),
     *   age0(j), asdep(j)/cog0,
     *   int(munit(j)), int(cog0)

      return
!     *****************      
      entry mergehyb3(h3)
!     ****************
      do while (h3%c%eventno  .ne. evno0) 
         call get1hyb( h3%c%eventno .lt. evno0) 
      enddo

      read(h3%c%id, '(i3)') j
      write(h3%c%id, 
     *   '(i3, i5,  f5.2, f5.2,
     *   i5,  i5)')
     *   j, int( asdep(j) ),
     *   age0(j), asdep(j)/cog0,
     *   int(munit(j)), int(cog0)

Here is the call graph for this function:

openhyb()

subroutine openhyb

(

integer

icon

)

Definition at line 261 of file bin2bin.f.

References copenfw2(), and parameter().

      implicit none
      integer icon  ! output.  1--> hybrid must be read
                    !          0--> hybrid need not be used
      integer leng
      integer ndepth
      parameter(ndepth= 50)
      integer fn1
      real*8  asdep(ndepth),  munit(ndepth)
      real*8  esize0(ndepth),
     *     age0(ndepth),  cogdep0(ndepth),
     *     seloss0(ndepth),
     *     ng0(ndepth), ne0(ndepth), nmu0(ndepth),
     *     cog0
      integer evno0

      common /zbin2ascii/ 
     *    asdep, esize0, age0,
     *    cogdep0,  seloss0, munit,
     *    ng0, ne0, nmu0, cog0,
     *    fn1, evno0

      character*120  hyb0
      integer kgetenv2
      
      fn1= 3
      leng = kgetenv2("HYBFILE0", hyb0)
      call copenfw2(fn1, hyb0, 1, icon)
      if(icon .ne. 1)  then
         write(0,*)
     *    "You haven't given env. var. HYBFILE0"
         write(0,*) 
     *    "or File specified by HYBFILE0"
         if( icon .eq. 0) then
            write(0,*) 'not exists'
         else
            write(0,*) ' cannot be opened '
         endif
         stop 9999
      else
         write(0,*)  hyb0(1:leng), ' opened'
         icon = 1
      endif

      evno0 =0 

Here is the call graph for this function: