NCEPLIBS-ip 4.0.0
Private Member Functions
spectral_interp_mod::polates4 Interface Reference

Private Member Functions

subroutine polates4_grib1 (IPOPT, KGDSI, KGDSO, MI, MO, KM, IBI, GI, NO, RLAT, RLON, IBO, LO, GO, IRET)
 THIS SUBPROGRAM PERFORMS SPECTRAL INTERPOLATION FROM ANY GRID TO ANY GRID FOR SCALAR FIELDS. More...
 
subroutine polates4_grib2 (IPOPT, IGDTNUMI, IGDTMPLI, IGDTLENI, IGDTNUMO, IGDTMPLO, IGDTLENO, MI, MO, KM, IBI, GI, NO, RLAT, RLON, IBO, LO, GO, IRET)
 

Detailed Description

Definition at line 16 of file spectral_interp_mod.f90.

Member Function/Subroutine Documentation

◆ polates4_grib1()

subroutine spectral_interp_mod::polates4::polates4_grib1 ( integer, dimension(20), intent(in)  IPOPT,
integer, dimension(200), intent(in)  KGDSI,
integer, dimension(200), intent(in)  KGDSO,
integer, intent(in)  MI,
integer, intent(in)  MO,
integer, intent(in)  KM,
integer, dimension(km), intent(in)  IBI,
real, dimension(mi,km), intent(in)  GI,
integer  NO,
real, dimension(mo), intent(inout)  RLAT,
real, dimension(mo), intent(inout)  RLON,
integer, dimension(km), intent(out)  IBO,
logical*1, dimension(mo,km), intent(out)  LO,
real, dimension(mo,km), intent(out)  GO,
integer, intent(out)  IRET 
)
private

THIS SUBPROGRAM PERFORMS SPECTRAL INTERPOLATION FROM ANY GRID TO ANY GRID FOR SCALAR FIELDS.

IT REQUIRES THAT THE INPUT FIELDS BE UNIFORMLY GLOBAL. OPTIONS ALLOW CHOICES BETWEEN TRIANGULAR SHAPE (IPOPT(1)=0) AND RHOMBOIDAL SHAPE (IPOPT(1)=1) WHICH HAS NO DEFAULT; A SECOND OPTION IS THE TRUNCATION (IPOPT(2)) WHICH DEFAULTS TO A SENSIBLE TRUNCATION FOR THE INPUT GRID (IF OPT(2)=-1). NOTE THAT IF THE OUTPUT GRID IS NOT FOUND IN A SPECIAL LIST, THEN THE TRANSFORM BACK TO GRID IS NOT VERY FAST. THIS SPECIAL LIST CONTAINS GLOBAL CYLINDRICAL GRIDS, POLAR STEREOGRAPHIC GRIDS CENTERED AT THE POLE AND MERCATOR GRIDS. ONLY HORIZONTAL INTERPOLATION IS PERFORMED. THE GRIDS ARE DEFINED BY THEIR GRID DESCRIPTION SECTIONS (PASSED IN INTEGER FORM AS DECODED BY SUBPROGRAM W3FI63). THE CURRENT CODE RECOGNIZES THE FOLLOWING PROJECTIONS: (KGDS(1)=000) EQUIDISTANT CYLINDRICAL (KGDS(1)=001) MERCATOR CYLINDRICAL (KGDS(1)=003) LAMBERT CONFORMAL CONICAL (KGDS(1)=004) GAUSSIAN CYLINDRICAL (SPECTRAL NATIVE) (KGDS(1)=005) POLAR STEREOGRAPHIC AZIMUTHAL (KGDS(1)=203) ROTATED EQUIDISTANT CYLINDRICAL (E-STAGGER) (KGDS(1)=205) ROTATED EQUIDISTANT CYLINDRICAL (B-STAGGER) WHERE KGDS COULD BE EITHER INPUT KGDSI OR OUTPUT KGDSO. AS AN ADDED BONUS THE NUMBER OF OUTPUT GRID POINTS AND THEIR LATITUDES AND LONGITUDES ARE ALSO RETURNED. ON THE OTHER HAND, THE OUTPUT CAN BE A SET OF STATION POINTS IF KGDSO(1)<0, IN WHICH CASE THE NUMBER OF POINTS AND THEIR LATITUDES AND LONGITUDES MUST BE INPUT. OUTPUT BITMAPS WILL NOT BE CREATED.

PROGRAM HISTORY LOG:

  • 96-04-10 IREDELL
  • 2001-06-18 IREDELL IMPROVE DETECTION OF SPECIAL FAST TRANSFORM
  • 2015-01-27 GAYNO REPLACE CALLS TO GDSWIZ WITH NEW MERGED VERSION OF GDSWZD.
Parameters
IPOPT- INTEGER (20) INTERPOLATION OPTIONS IPOPT(1)=0 FOR TRIANGULAR, IPOPT(1)=1 FOR RHOMBOIDAL; IPOPT(2) IS TRUNCATION NUMBER (DEFAULTS TO SENSIBLE IF IPOPT(2)=-1).
KGDSI- INTEGER (200) INPUT GDS PARAMETERS AS DECODED BY W3FI63
KGDSO- INTEGER (200) OUTPUT GDS PARAMETERS (KGDSO(1)<0 IMPLIES RANDOM STATION POINTS)
MI- INTEGER SKIP NUMBER BETWEEN INPUT GRID FIELDS IF KM>1 OR DIMENSION OF INPUT GRID FIELDS IF KM=1
MO- INTEGER SKIP NUMBER BETWEEN OUTPUT GRID FIELDS IF KM>1 OR DIMENSION OF OUTPUT GRID FIELDS IF KM=1
KM- INTEGER NUMBER OF FIELDS TO INTERPOLATE
IBI- INTEGER (KM) INPUT BITMAP FLAGS (MUST BE ALL 0)
LI- LOGICAL*1 (MI,KM) INPUT BITMAPS (IF SOME IBI(K)=1)
GI- REAL (MI,KM) INPUT FIELDS TO INTERPOLATE
[out]NO- INTEGER NUMBER OF OUTPUT POINTS (ONLY IF KGDSO(1)<0)
[out]RLAT- REAL (NO) OUTPUT LATITUDES IN DEGREES (IF KGDSO(1)<0)
[out]RLON- REAL (NO) OUTPUT LONGITUDES IN DEGREES (IF KGDSO(1)<0)
[out]IBO- INTEGER (KM) OUTPUT BITMAP FLAGS
[out]LO- LOGICAL*1 (MO,KM) OUTPUT BITMAPS (ALWAYS OUTPUT)
[out]GO- REAL (MO,KM) OUTPUT FIELDS INTERPOLATED
[out]IRET- INTEGER RETURN CODE 0 SUCCESSFUL INTERPOLATION 2 UNRECOGNIZED INPUT GRID OR NO GRID OVERLAP 3 UNRECOGNIZED OUTPUT GRID 41 INVALID NONGLOBAL INPUT GRID 42 INVALID SPECTRAL METHOD PARAMETERS

SUBPROGRAMS CALLED:

  • GDSWZD GRID DESCRIPTION SECTION WIZARD
  • SPTRUN SPECTRALLY TRUNCATE GRIDDED SCALAR FIELDS
  • SPTRUNS SPECTRALLY INTERPOLATE SCALARS TO POLAR STEREO.
  • SPTRUNM SPECTRALLY INTERPOLATE SCALARS TO MERCATOR
  • SPTRUNG SPECTRALLY INTERPOLATE SCALARS TO STATIONS

Definition at line 536 of file spectral_interp_mod.f90.

◆ polates4_grib2()

subroutine spectral_interp_mod::polates4::polates4_grib2 ( integer, dimension(20), intent(in)  IPOPT,
integer, intent(in)  IGDTNUMI,
integer, dimension(igdtleni), intent(in)  IGDTMPLI,
integer, intent(in)  IGDTLENI,
integer, intent(in)  IGDTNUMO,
integer, dimension(igdtleno), intent(in)  IGDTMPLO,
integer, intent(in)  IGDTLENO,
integer, intent(in)  MI,
integer, intent(in)  MO,
integer, intent(in)  KM,
integer, dimension(km), intent(in)  IBI,
real, dimension(mi,km), intent(in)  GI,
integer, intent(out)  NO,
real, dimension(mo), intent(inout)  RLAT,
real, dimension(mo), intent(inout)  RLON,
integer, dimension(km), intent(out)  IBO,
logical*1, dimension(mo,km), intent(out)  LO,
real, dimension(mo,km), intent(out)  GO,
integer, intent(out)  IRET 
)
private

Definition at line 102 of file spectral_interp_mod.f90.

The documentation for this interface was generated from the following file: