ipxetas.F90

Go to the documentation of this file.

 
 
 
 
! @author Iredell @date 96-04-10
 SUBROUTINE ipxetas(IDIR, IGDTNUMI, IGDTLEN, IGDTMPLI, NPTS_INPUT,  &
                    BITMAP_INPUT, DATA_INPUT, IGDTNUMO, IGDTMPLO, &
                    NPTS_OUTPUT, BITMAP_OUTPUT, DATA_OUTPUT, IRET)
 USE ipolates_mod
 IMPLICIT NONE
!
 INTEGER,         INTENT(IN   )    :: IDIR
 INTEGER,         INTENT(IN   )    :: IGDTNUMI, IGDTLEN
 INTEGER,         INTENT(IN   )    :: IGDTMPLI(IGDTLEN)
 INTEGER,         INTENT(IN   )    :: NPTS_INPUT, NPTS_OUTPUT
 INTEGER,         INTENT(  OUT)    :: IGDTNUMO
 INTEGER,         INTENT(  OUT)    :: IGDTMPLO(IGDTLEN)
 INTEGER,         INTENT(  OUT)    :: IRET
 
 LOGICAL(KIND=1), INTENT(IN   )    :: BITMAP_INPUT(NPTS_INPUT)
 LOGICAL(KIND=1), INTENT(  OUT)    :: BITMAP_OUTPUT(NPTS_OUTPUT)
 
 REAL,            INTENT(IN   )    :: DATA_INPUT(NPTS_INPUT)
 REAL,            INTENT(  OUT)    :: DATA_OUTPUT(NPTS_OUTPUT)
 
 INTEGER                           :: SCAN_MODE, ISCALE, IP, IPOPT(20)
 INTEGER                           :: IBI, IBO, J, KM, NO
 
 REAL                              :: DLONS
 REAL, ALLOCATABLE                 :: OUTPUT_RLAT(:), OUTPUT_RLON(:)
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
 iret = 0
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
! ROUTINE ONLY WORKS FOR ROTATED LAT/LON GRIDS.
 IF (igdtnumi/=1) THEN
   iret=1
   RETURN
 ENDIF
!
 scan_mode=igdtmpli(19)
 IF((scan_mode==68.OR.scan_mode==72).AND.(idir<-2.OR.idir>-1))THEN
   igdtnumo=igdtnumi
   igdtmplo=igdtmpli
   igdtmplo(19)=64
   igdtmplo(8)=igdtmplo(8)*2-1
   IF((igdtmplo(8)*igdtmplo(9))/=npts_output)THEN
     iret=3
     RETURN
   ENDIF
   iscale=igdtmplo(10)*igdtmplo(11)
   IF(iscale==0) iscale=10**6
   dlons=float(igdtmplo(17))/float(iscale)
   dlons=dlons*0.5
   igdtmplo(17)=nint(dlons*float(iscale))
 ELSEIF(scan_mode==64.AND.idir==-1)THEN  ! FULL TO H-GRID
   igdtnumo=igdtnumi
   igdtmplo=igdtmpli
   igdtmplo(19)=68
   igdtmplo(8)=(igdtmplo(8)+1)/2
   IF((igdtmplo(8)*igdtmplo(9))/=npts_output)THEN
     iret=3
     RETURN
   ENDIF
   iscale=igdtmplo(10)*igdtmplo(11)
   IF(iscale==0) iscale=10**6
   dlons=float(igdtmplo(17))/float(iscale)
   dlons=dlons*2.0
   igdtmplo(17)=nint(dlons*float(iscale))
 ELSEIF(scan_mode==64.AND.idir==-2)THEN  ! FULL TO V-GRID
   igdtnumo=igdtnumi
   igdtmplo=igdtmpli
   igdtmplo(19)=72
   igdtmplo(8)=(igdtmplo(8)+1)/2
   IF((igdtmplo(8)*igdtmplo(9))/=npts_output)THEN
     iret=3
     RETURN
   ENDIF
   iscale=igdtmplo(10)*igdtmplo(11)
   IF(iscale==0) iscale=10**6
   dlons=float(igdtmplo(17))/float(iscale)
   dlons=dlons*2.0
   igdtmplo(17)=nint(dlons*float(iscale))
 ELSE
   iret=2
   RETURN
 ENDIF
 
 km=1
 ip=0
 ipopt=0
 ibi=1
 ibo=0
 
 ALLOCATE(output_rlat(npts_output))
 ALLOCATE(output_rlon(npts_output))
 
 CALL ipolates(ip, ipopt, igdtnumi, igdtmpli, igdtlen, &
               igdtnumo, igdtmplo, igdtlen, &
               npts_input, npts_output, km, ibi, bitmap_input, data_input, &
               no, output_rlat, output_rlon, ibo, bitmap_output, data_output, iret)
 
 DEALLOCATE(output_rlat, output_rlon)
 
 IF(iret /= 0)THEN
   print*,'- PROBLEM IN IPOLATES: ', iret
   RETURN
 ENDIF
 
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
! REPLACE ANY UNDEFINED POINTS ALONG THE LEFT AND RIGHT EDGES.
 DO j=1, igdtmplo(9)
   bitmap_output(j*igdtmplo(8))=bitmap_output(j*igdtmplo(8)-1)
   data_output(j*igdtmplo(8))=data_output(j*igdtmplo(8)-1)
   bitmap_output((j-1)*igdtmplo(8)+1)=bitmap_output((j-1)*igdtmplo(8)+2)
   data_output((j-1)*igdtmplo(8)+1)=data_output((j-1)*igdtmplo(8)+2)
 ENDDO
 
 RETURN
! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
 END SUBROUTINE ipxetas