drawRigCut2.f

Go to the documentation of this file.

#include "BlockData/cblkManager.h"
!
!   This is to make grid data and physical data (prob) for rigidity cut off
!   of new format
!
      implicit none
#include "Zmanagerp.h"
#include "ZrigCut.h"

      character*120 file
!

      file = '../Contrib/Data/CutoffNew/Kamioka'
      write(errorout, *) ' enter file name=',file
      read(*, *) file

      call crigcut0(file)
      call crigcut2grid(rigtbl2, zensize, azmsize, rigsize)
      call crigcut2data(rigtbl2, zensize, azmsize, rigsize)
      write(errorout, *)
     * 'Note:You have to set "variable aspect ratio" in Nextcontour'
      end
!     ***************************
      subroutine crigcut2data(tbl, izen, iphi, irig)
      implicit none
      integer izen, iphi, irig
      real*4 tbl(izen, iphi, irig), z
      integer i1, i2, i3,icon

      call copenfw(21, "/tmp/nxconData",icon)
      write(21, *) iphi, izen, irig
      do i3 =1, irig
         do i1 = izen, 1, -1 
            do i2 = 1, iphi
               z = tbl(i1, i2, i3)
               if(z .gt. 0.99) z =1.
               write(21,*) z
            enddo
         enddo
      enddo
      close(21)
      end
!     ************************
      subroutine crigcut2grid(tbl, izen, iphi, irig)
      implicit none
#include "ZrigCut.h"
      integer izen, iphi, irig
      real*4 tbl(izen, iphi, irig)
      integer i1, i2, i3,icon

      call copenfw(21, "/tmp/nxconGrid",icon)
      write(21, *) iphi, izen, irig
      do i3 =1, irig
         do i1 =  izen, 1, -1
            do i2 = 1, iphi
               write(21,*) sngl(minazm + (i2-1)*dazm)
            enddo
         enddo
      enddo
      do i3 =1, irig
         do i1 = izen, 1, -1
            do i2 = 1, iphi
               write(21,*) sngl(zenmax + (i1-1)*dzen), 
            enddo
         enddo
      enddo
      do i3 =1, irig
         do i1 = izen, 1, -1
            do i2 = 1, iphi
               write(21,*)
     *          sngl(minrig * 10.**(logdrig*(i3-1)))
            enddo
         enddo
      enddo

      close(21)
      end