params.f

计算线性趋势 回归系数 主要用于气象站点值的线性趋势计算

      data paramlist(213)%grib2dsc /0/
      data paramlist(213)%grib2cat /2/
      data paramlist(213)%grib2num /0/   
      data paramlist(213)%abbrev   /'WDIR    '/
      data paramlist(214)%g1tblver /0/
      data paramlist(214)%grib1val /32/
      data paramlist(214)%grib2dsc /0/
      data paramlist(214)%grib2cat /2/
      data paramlist(214)%grib2num /1/   
      data paramlist(214)%abbrev   /'WIND    '/
      data paramlist(215)%g1tblver /0/
      data paramlist(215)%grib1val /33/
      data paramlist(215)%grib2dsc /0/
      data paramlist(215)%grib2cat /2/
      data paramlist(215)%grib2num /2/   
      data paramlist(215)%abbrev   /'U GRD   '/
      data paramlist(216)%g1tblver /0/
      data paramlist(216)%grib1val /34/
      data paramlist(216)%grib2dsc /0/
      data paramlist(216)%grib2cat /2/
      data paramlist(216)%grib2num /3/   
      data paramlist(216)%abbrev   /'V GRD   '/
      data paramlist(217)%g1tblver /0/
      data paramlist(217)%grib1val /100/
      data paramlist(217)%grib2dsc /10/
      data paramlist(217)%grib2cat /0/
      data paramlist(217)%grib2num /3/   
      data paramlist(217)%abbrev   /'HTSGW   '/
      data paramlist(218)%g1tblver /0/
      data paramlist(218)%grib1val /101/
      data paramlist(218)%grib2dsc /10/
      data paramlist(218)%grib2cat /0/
      data paramlist(218)%grib2num /4/   
      data paramlist(218)%abbrev   /'WVDIR   '/
      data paramlist(219)%g1tblver /0/
      data paramlist(219)%grib1val /103/
      data paramlist(219)%grib2dsc /10/
      data paramlist(219)%grib2cat /0/
      data paramlist(219)%grib2num /6/   
      data paramlist(219)%abbrev   /'WVPER   '/
      data paramlist(220)%g1tblver /0/
      data paramlist(220)%grib1val /107/
      data paramlist(220)%grib2dsc /10/
      data paramlist(220)%grib2cat /0/
      data paramlist(220)%grib2num /10/  
      data paramlist(220)%abbrev   /'DIRPW   '/
      data paramlist(221)%g1tblver /0/
      data paramlist(221)%grib1val /108/
      data paramlist(221)%grib2dsc /10/
      data paramlist(221)%grib2cat /0/
      data paramlist(221)%grib2num /11/  
      data paramlist(221)%abbrev   /'PERPW   '/
      data paramlist(222)%g1tblver /0/
      data paramlist(222)%grib1val /109/
      data paramlist(222)%grib2dsc /10/
      data paramlist(222)%grib2cat /0/
      data paramlist(222)%grib2num /12/  
      data paramlist(222)%abbrev   /'DIRSW   '/
      data paramlist(223)%g1tblver /0/
      data paramlist(223)%grib1val /110/
      data paramlist(223)%grib2dsc /10/
      data paramlist(223)%grib2cat /0/
      data paramlist(223)%grib2num /13/    
      data paramlist(223)%abbrev   /'PERSW   '/
! Added 11/17/2005 - for wave models
      data paramlist(224)%g1tblver /129/
      data paramlist(224)%grib1val /156/
      data paramlist(224)%grib2dsc /0/
      data paramlist(224)%grib2cat /13/
      data paramlist(224)%grib2num /192/ 
      data paramlist(224)%abbrev   /'PMTC    '/
      data paramlist(225)%g1tblver /129/
      data paramlist(225)%grib1val /157/
      data paramlist(225)%grib2dsc /0/
      data paramlist(225)%grib2cat /13/
      data paramlist(225)%grib2num /193/ 
      data paramlist(225)%abbrev   /'PMTF    '/
      data paramlist(226)%g1tblver /3/
      data paramlist(226)%grib1val /11/
      data paramlist(226)%grib2dsc /0/
      data paramlist(226)%grib2cat /0/
      data paramlist(226)%grib2num /0/   
      data paramlist(226)%abbrev   /'TMP     '/
      data paramlist(227)%g1tblver /2/
      data paramlist(227)%grib1val /129/
      data paramlist(227)%grib2dsc /0/
      data paramlist(227)%grib2cat /3/
      data paramlist(227)%grib2num /198/ 
      data paramlist(227)%abbrev   /'MSLMA   '/
      data paramlist(228)%g1tblver /129/
      data paramlist(228)%grib1val /163/
      data paramlist(228)%grib2dsc /0/
      data paramlist(228)%grib2cat /13/
      data paramlist(228)%grib2num /194/ 
      data paramlist(228)%abbrev   /'LPMTF   '/
      data paramlist(229)%g1tblver /129/
      data paramlist(229)%grib1val /164/
      data paramlist(229)%grib2dsc /0/
      data paramlist(229)%grib2cat /13/
      data paramlist(229)%grib2num /195/ 
      data paramlist(229)%abbrev   /'LIPMF   '/
! Added 8/25/2006 - JFB NCAR/MMM
      data paramlist(230)%g1tblver /2/
      data paramlist(230)%grib1val /189/
      data paramlist(230)%grib2dsc /0/
      data paramlist(230)%grib2cat /0/
      data paramlist(230)%grib2num /15/  
      data paramlist(230)%abbrev   /'VPTMP   '/
      data paramlist(231)%g1tblver /2/
      data paramlist(231)%grib1val /178/
      data paramlist(231)%grib2dsc /0/
      data paramlist(231)%grib2cat /1/
      data paramlist(231)%grib2num /23/  
      data paramlist(231)%abbrev   /'ICMR    '/
      data paramlist(232)%g1tblver /2/
      data paramlist(232)%grib1val /179/
      data paramlist(232)%grib2dsc /0/
      data paramlist(232)%grib2cat /1/
      data paramlist(232)%grib2num /32/  
      data paramlist(232)%abbrev   /'GRMR    '/
      data paramlist(233)%g1tblver /2/
      data paramlist(233)%grib1val /198/
      data paramlist(233)%grib2dsc /0/
      data paramlist(233)%grib2cat /1/
      data paramlist(233)%grib2num /207/ 
      data paramlist(233)%abbrev   /'NCIP    '/
      data paramlist(234)%g1tblver /2/
      data paramlist(234)%grib1val /186/
      data paramlist(234)%grib2dsc /0/
      data paramlist(234)%grib2cat /1/
      data paramlist(234)%grib2num /206/ 
      data paramlist(234)%abbrev   /'TIPD    '/
      data paramlist(235)%g1tblver /2/
      data paramlist(235)%grib1val /188/
      data paramlist(235)%grib2dsc /2/
      data paramlist(235)%grib2cat /0/
      data paramlist(235)%grib2num /206/ 
      data paramlist(235)%abbrev   /'RDRIP   '/
      data paramlist(236)%g1tblver /2/
      data paramlist(236)%grib1val /239/
      data paramlist(236)%grib2dsc /0/
      data paramlist(236)%grib2cat /1/
      data paramlist(236)%grib2num /208/ 
      data paramlist(236)%abbrev   /'SNO T   '/
      data paramlist(237)%g1tblver /130/
      data paramlist(237)%grib1val /171/
      data paramlist(237)%grib2dsc /2/
      data paramlist(237)%grib2cat /3/
      data paramlist(237)%grib2num /193/ 
      data paramlist(237)%abbrev   /'RLYRS   '/
      data paramlist(238)%g1tblver /2/
      data paramlist(238)%grib1val /187/
      data paramlist(238)%grib2dsc /0/
      data paramlist(238)%grib2cat /17/
      data paramlist(238)%grib2num /192/ 
      data paramlist(238)%abbrev   /'LTNG    '/
      data paramlist(239)%g1tblver /2/
      data paramlist(239)%grib1val /137/
      data paramlist(239)%grib2dsc /0/
      data paramlist(239)%grib2cat /3/
      data paramlist(239)%grib2num /199/ 
      data paramlist(239)%abbrev   /'TSLSA   '/

      contains


         subroutine param_g1_to_g2(g1val,g1ver,g2disc,g2cat,g2num)
!$$$  SUBPROGRAM DOCUMENTATION BLOCK
!                .      .    .                                       .
! SUBPROGRAM:    param_g1_to_g2 
!   PRGMMR: Gilbert         ORG: W/NP11    DATE: 2001-06-05
!
! ABSTRACT: This subroutine returns the corresponding GRIB2 Discipline
!   Category and Number for a given GRIB1 parameter value and table version.
!
! PROGRAM HISTORY LOG:
! 2000-05-11  Gilbert
!
! USAGE:    CALL param_g1_to_g2(g1val,g1ver,g2disc,g2cat,g2num)
!   INPUT ARGUMENT LIST:
!     g1val    - GRIB1 parameter number for which discipline is requested
!     g1ver    - GRIB1 parameter table version number
!
!   OUTPUT ARGUMENT LIST:      
!     g2disc   - corresponding GRIB2 Discipline number
!     g2cat    - corresponding GRIB2 Category number
!     g2num    - corresponding GRIB2 Parameter number within Category g2cat
!
! ATTRIBUTES:
!   LANGUAGE: Fortran 90
!   MACHINE:  IBM SP
!
!$$$
           integer,intent(in) :: g1val,g1ver
           integer,intent(out) :: g2disc,g2cat,g2num

           g2disc=255
           g2cat=255
           g2num=255
! for testing
!           g2num=g1val
! for testing

           do n=1,MAXPARAM
              if (paramlist(n)%grib1val.eq.g1val .AND.
     &            paramlist(n)%g1tblver.eq.g1ver ) then
                 g2disc=paramlist(n)%grib2dsc
                 g2cat=paramlist(n)%grib2cat
                 g2num=paramlist(n)%grib2num
                 return
              endif
           enddo

           print *,'param_g1_to_g2:GRIB1 param ',g1val,' not found.',
     &             ' for table version ',g1ver
           return
         end subroutine

         character(len=8) function param_get_abbrev(g2disc,g2cat,g2num)
!$$$  SUBPROGRAM DOCUMENTATION BLOCK
!                .      .    .                                       .
! SUBPROGRAM:    param_get_abbrev 
!   PRGMMR: Gilbert         ORG: W/NP11    DATE: 2002-01-04
!
! ABSTRACT: This function returns the parameter abbreviation for
!   a given GRIB2 Discipline, Category and Parameter number.
!
! PROGRAM HISTORY LOG:
! 2001-06-05  Gilbert
!
! USAGE:     abrev=param_get_abbrev(g2disc,g2cat,g2num)
!   INPUT ARGUMENT LIST:
!     g2disc   - GRIB2 discipline number (See Code Table 0.0)
!     g2cat    - corresponding GRIB2 Category number
!     g2num    - corresponding GRIB2 Parameter number within Category g2cat
!
! RETURNS:  ASCII Paramter Abbreviation
!
! REMARKS: None
!
! ATTRIBUTES:
!   LANGUAGE: Fortran 90
!   MACHINE:  IBM SP
!
!$$$
           integer,intent(in) :: g2disc,g2cat,g2num

           param_get_abbrev='UNKNOWN '

           do n=1,MAXPARAM
              if (paramlist(n)%grib2dsc.eq.g2disc.AND.
     &             paramlist(n)%grib2cat.eq.g2cat.AND.
     &             paramlist(n)%grib2num.eq.g2num) then
                 param_get_abbrev=paramlist(n)%abbrev
                 return
              endif
           enddo

!           print *,'param_get_abbrev:GRIB2 param ',g2disc,g2cat,
!     &              g2num,' not found.'
           return
         end function


         subroutine param_g2_to_g1(g2disc,g2cat,g2num,g1val,g1ver)
!$$$  SUBPROGRAM DOCUMENTATION BLOCK
!                .      .    .                                       .
! SUBPROGRAM:    param_g2_to_g1 
!   PRGMMR: Gilbert         ORG: W/NP11    DATE: 2002-01-04
!
! ABSTRACT: This function returns the GRIB 1 parameter number for 
!   a given GRIB2 Discipline, Category and Parameter number.
!
! PROGRAM HISTORY LOG:
! 2001-06-05  Gilbert
!
! USAGE:     call param_g2_to_g1(g2disc,g2cat,g2num,g1val,g1ver)
!   INPUT ARGUMENT LIST:
!     g2disc   - GRIB2 discipline number (See Code Table 0.0)
!     g2cat    - corresponding GRIB2 Category number
!     g2num    - corresponding GRIB2 Parameter number within Category g2cat
!
!   OUTPUT ARGUMENT LIST:      
!     g1val    - GRIB1 parameter number for which discipline is requested
!     g1ver    - GRIB1 parameter table version number
!
! REMARKS: None
!
! ATTRIBUTES:
!   LANGUAGE: Fortran 90
!   MACHINE:  IBM SP
!
!$$$
           integer,intent(in) :: g2disc,g2cat,g2num
           integer,intent(out) :: g1val,g1ver

           g1val=255
           g1ver=255

! for testing
!           if ( g2disc.eq.255.and.g2cat.eq.255 ) then
!             g1val=g2num
!             g1ver=2
!             return
!           endif
! for testing

           do n=1,MAXPARAM
              if (paramlist(n)%grib2dsc.eq.g2disc.AND.
     &             paramlist(n)%grib2cat.eq.g2cat.AND.
     &             paramlist(n)%grib2num.eq.g2num) then
                 g1val=paramlist(n)%grib1val
                 g1ver=paramlist(n)%g1tblver
                 return
              endif
           enddo

           print *,'param_g2_to_g1:GRIB2 param ',g2disc,g2cat,
     &              g2num,' not found.'
           return
         end subroutine


      end module