gridtemplates.f

Go to the documentation of this file.

 
 
 
      module gridtemplates
 
      integer,parameter :: maxlen=200
      integer,parameter :: maxtemp=31
 
      type gridtemplate
          integer :: template_num
          integer :: mapgridlen
          integer,dimension(MAXLEN) :: mapgrid
          logical :: needext
      end type gridtemplate
 
      type(gridtemplate),dimension(MAXTEMP) :: templates
 
      data templates(1)%template_num /0/     !  Lat/Lon 
      data templates(1)%mapgridlen /19/
      data templates(1)%needext /.false./
      data (templates(1)%mapgrid(j),j=1,19) 
     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1/
 
      data templates(2)%template_num /1/     !  Rotated Lat/Lon 
      data templates(2)%mapgridlen /22/
      data templates(2)%needext /.false./
      data (templates(2)%mapgrid(j),j=1,22) 
     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,4/
 
      data templates(3)%template_num /2/     !  Stretched Lat/Lon 
      data templates(3)%mapgridlen /22/
      data templates(3)%needext /.false./
      data (templates(3)%mapgrid(j),j=1,22) 
     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,-4/
 
      data templates(4)%template_num /3/     !  Stretched & Rotated Lat/Lon 
      data templates(4)%mapgridlen /25/
      data templates(4)%needext /.false./
      data (templates(4)%mapgrid(j),j=1,25) 
     &       /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,4,-4,4,-4/
 
      data templates(5)%template_num /10/     !  Mercator
      data templates(5)%mapgridlen /19/
      data templates(5)%needext /.false./
      data (templates(5)%mapgrid(j),j=1,19)
     &              /1,1,4,1,4,1,4,4,4,-4,4,1,-4,-4,4,1,4,4,4/
 
      data templates(6)%template_num /20/     !  Polar Stereographic
      data templates(6)%mapgridlen /18/
      data templates(6)%needext /.false./
      data (templates(6)%mapgrid(j),j=1,18) 
     &              /1,1,4,1,4,1,4,4,4,-4,4,1,-4,4,4,4,1,1/
 
      data templates(7)%template_num /30/     !  Lambert Conformal
      data templates(7)%mapgridlen /22/
      data templates(7)%needext /.false./
      data (templates(7)%mapgrid(j),j=1,22) 
     &              /1,1,4,1,4,1,4,4,4,-4,4,1,-4,4,4,4,1,1,-4,-4,-4,4/
 
      data templates(8)%template_num /40/     !  Gaussian Lat/Lon
      data templates(8)%mapgridlen /19/
      data templates(8)%needext /.false./
      data (templates(8)%mapgrid(j),j=1,19) 
     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1/
 
      data templates(9)%template_num /41/     !  Rotated Gaussian Lat/Lon
      data templates(9)%mapgridlen /22/
      data templates(9)%needext /.false./
      data (templates(9)%mapgrid(j),j=1,22) 
     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,4/
 
      data templates(10)%template_num /42/     !  Stretched Gaussian Lat/Lon
      data templates(10)%mapgridlen /22/
      data templates(10)%needext /.false./
      data (templates(10)%mapgrid(j),j=1,22) 
     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,-4/
 
      data templates(11)%template_num /43/     !  Strtchd and Rot'd Gaus Lat/Lon
      data templates(11)%mapgridlen /25/
      data templates(11)%needext /.false./
      data (templates(11)%mapgrid(j),j=1,25) 
     &          /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,-4,4,4,-4,4,-4/
 
      data templates(12)%template_num /50/    !  Spherical Harmonic Coefficients
      data templates(12)%mapgridlen /5/
      data templates(12)%needext /.false./
      data (templates(12)%mapgrid(j),j=1,5) /4,4,4,1,1/
 
      data templates(13)%template_num /51/   !  Rotated Spherical Harmonic Coeff
      data templates(13)%mapgridlen /8/
      data templates(13)%needext /.false./
      data (templates(13)%mapgrid(j),j=1,8) /4,4,4,1,1,-4,4,4/
 
      data templates(14)%template_num /52/   !  Stretch Spherical Harmonic Coeff
      data templates(14)%mapgridlen /8/
      data templates(14)%needext /.false./
      data (templates(14)%mapgrid(j),j=1,8) /4,4,4,1,1,-4,4,-4/
 
      data templates(15)%template_num /53/   !  Strch and Rot Spher Harm Coeffs
      data templates(15)%mapgridlen /11/
      data templates(15)%needext /.false./
      data (templates(15)%mapgrid(j),j=1,11) /4,4,4,1,1,-4,4,4,-4,4,-4/
 
      data templates(16)%template_num /90/     !  Space view Perspective
      data templates(16)%mapgridlen /21/
      data templates(16)%needext /.false./
      data (templates(16)%mapgrid(j),j=1,21) 
     &              /1,1,4,1,4,1,4,4,4,-4,4,1,4,4,4,4,1,4,4,4,4/
 
      data templates(17)%template_num /100/    !  Triangular grid (icosahedron)
      data templates(17)%mapgridlen /11/
      data templates(17)%needext /.false./
      data (templates(17)%mapgrid(j),j=1,11) /1,1,2,1,-4,4,4,1,1,1,4/
 
      data templates(18)%template_num /110/ !  Equatorial Azimuthal equidistant
      data templates(18)%mapgridlen /16/
      data templates(18)%needext /.false./
      data (templates(18)%mapgrid(j),j=1,16) 
     &              /1,1,4,1,4,1,4,4,4,-4,4,1,4,4,1,1/
 
       data templates(19)%template_num /120/     !  Azimuth-range 
       data templates(19)%mapgridlen /7/
       data templates(19)%needext /.true./
       data (templates(19)%mapgrid(j),j=1,7) /4,4,-4,4,4,4,1/
 
       data templates(20)%template_num /1000/     !  Cross Section Grid 
       data templates(20)%mapgridlen /20/
       data templates(20)%needext /.true./
       data (templates(20)%mapgrid(j),j=1,20) 
     &              /1,1,4,1,4,1,4,4,4,4,-4,4,1,4,4,1,2,1,1,2/
 
       data templates(21)%template_num /1100/     !  Hovmoller Diagram Grid 
       data templates(21)%mapgridlen /28/
       data templates(21)%needext /.false./
       data (templates(21)%mapgrid(j),j=1,28) 
     &    /1,1,4,1,4,1,4,4,4,4,-4,4,1,-4,4,1,4,1,-4,1,1,-4,2,1,1,1,1,1/
 
       data templates(22)%template_num /1200/     !  Time Section Grid 
       data templates(22)%mapgridlen /16/
       data templates(22)%needext /.true./
       data (templates(22)%mapgrid(j),j=1,16) 
     &              /4,1,-4,1,1,-4,2,1,1,1,1,1,2,1,1,2/
 
      data templates(23)%template_num /31/     !  Albers Equal Area
      data templates(23)%mapgridlen /22/
      data templates(23)%needext /.false./
      data (templates(23)%mapgrid(j),j=1,22) 
     &              /1,1,4,1,4,1,4,4,4,-4,4,1,-4,4,4,4,1,1,-4,-4,-4,4/
 
      data templates(24)%template_num /204/     !  Curilinear Orthogonal Grids
      data templates(24)%mapgridlen /19/
      data templates(24)%needext /.false./
      data (templates(24)%mapgrid(j),j=1,19)
     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1/
 
      data templates(25)%template_num /32768/   ! Rotate Lat/Lon E-grid
      data templates(25)%mapgridlen /19/
      data templates(25)%needext /.false./
      data (templates(25)%mapgrid(j),j=1,19)
     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1/
 
      data templates(26)%template_num /32769/   ! Rotate Lat/Lon Non-E Stagger grid
      data templates(26)%mapgridlen /21/
      data templates(26)%needext /.false./
      data (templates(26)%mapgrid(j),j=1,21)
     &              /1,1,4,1,4,1,4,4,4,4,4,-4,4,1,-4,4,4,4,1,4,4/
!
!     GDT 3.4 Added (08/05/2013)
!
      data templates(27)%template_num /4/     !  Variable resolution Latitude/Longitude
      data templates(27)%mapgridlen /13/
      data templates(27)%needext /.true./
      data (templates(27)%mapgrid(j),j=1,13)
     &              /1,1,4,1,4,1,4,4,4,4,4,1,1/
!
!     GDT 3.5 Added (08/05/2013)
!
      data templates(28)%template_num /5/     !  Variable resolution rotate Latitude/Longitude
      data templates(28)%mapgridlen /16/
      data templates(28)%needext /.true./
      data (templates(28)%mapgrid(j),j=1,16)
     &              /1,1,4,1,4,1,4,4,4,4,4,1,1,-4,4,4/
!
!     GDT 3.12 Added (08/05/2013)
!
      data templates(29)%template_num /12/     !  Transverse Mercator
      data templates(29)%mapgridlen /22/
      data templates(29)%needext /.false./
      data (templates(29)%mapgrid(j),j=1,22)
     &              /1,1,4,1,4,1,4,4,4,-4,4,1,-4,4,4,1,4,4,-4,-4,-4,-4/
!
!     GDT 3.101 Added (08/05/2013)
!
      data templates(30)%template_num /101/     !  General unstructured grid
      data templates(30)%mapgridlen /4/
      data templates(30)%needext /.false./
      data (templates(30)%mapgrid(j),j=1,4)
     &              /1,4,1,-4/
!
!     GDT 3.140 Added (08/05/2013)
!
      data templates(31)%template_num /140/     !  Lambert Azimuthal Equal Area Projection
      data templates(31)%mapgridlen /17/
      data templates(31)%needext /.false./
      data (templates(31)%mapgrid(j),j=1,17)
     &              /1,1,4,1,4,1,4,4,4,-4,4,4,4,1,4,4,1/
 
      contains
 
 
         integer function getgridindex(number)
 
           integer,intent(in) :: number
 
           getgridindex=-1
 
           do j=1,maxtemp
              if (number.eq.templates(j)%template_num) then
                 getgridindex=j
                 return
              endif
           enddo
 
         end function
 
 
         subroutine getgridtemplate(number,nummap,map,needext,iret)
 
           integer,intent(in) :: number
           integer,intent(out) :: nummap,map(*),iret
           logical,intent(out) :: needext
 
           iret=0
 
           index=getgridindex(number)
 
           if (index.ne.-1) then
              nummap=templates(index)%mapgridlen
              needext=templates(index)%needext
              map(1:nummap)=templates(index)%mapgrid(1:nummap)
           else
             nummap=0
             needext=.false.
             print *,'getgridtemplate: Grid Template ',number,
     &               ' not defined.'
             iret=1
           endif
 
         end subroutine
 
 
         subroutine extgridtemplate(number,list,nummap,map)
 
           integer,intent(in) :: number,list(*)
           integer,intent(out) :: nummap,map(*)
 
           index=getgridindex(number)
           if (index.eq.-1) return
 
           if ( .not. templates(index)%needext ) return
           nummap=templates(index)%mapgridlen
           map(1:nummap)=templates(index)%mapgrid(1:nummap)
 
           if ( number.eq.120 ) then
              n=list(2)
              do i=1,n
                map(nummap+1)=2
                map(nummap+2)=-2
                nummap=nummap+2
              enddo
           elseif ( number.eq.4 ) then
              ni=list(8)
              do i=1,ni
                map(nummap+1)=4
                nummap=nummap+1
              enddo
              nj=list(9)
              do i=1,nj
                map(nummap+1)=-4
                nummap=nummap+1
              enddo
           elseif ( number.eq.5 ) then
              ni=list(8)
              do i=1,ni
                map(nummap+1)=4
                nummap=nummap+1
              enddo
              nj=list(9)
              do i=1,nj
                map(nummap+1)=-4
                nummap=nummap+1
              enddo
           elseif ( number.eq.1000 ) then
              n=list(20)
              do i=1,n
                map(nummap+1)=4
                nummap=nummap+1
              enddo
           elseif ( number.eq.1200 ) then
              n=list(16)
              do i=1,n
                map(nummap+1)=4
                nummap=nummap+1
              enddo
           endif
 
         end subroutine
 
 
         integer function getgdtlen(number)
 
           integer,intent(in) :: number
 
           getgdtlen=0
 
           index=getgridindex(number)
 
           if (index.ne.-1) then
              getgdtlen=templates(index)%mapgridlen
           endif
 
         end function
 
 
      end