!-------------------------------------- LICENCE BEGIN ------------------------------------ !Environment Canada - Atmospheric Science and Technology License/Disclaimer, ! version 3; Last Modified: May 7, 2008. !This is free but copyrighted software; you can use/redistribute/modify it under the terms !of the Environment Canada - Atmospheric Science and Technology License/Disclaimer !version 3 or (at your option) any later version that should be found at: !http://collaboration.cmc.ec.gc.ca/science/rpn.comm/license.html ! !This software is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; !without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. !See the above mentioned License/Disclaimer for more details. !You should have received a copy of the License/Disclaimer along with this software; !if not, you can write to: EC-RPN COMM Group, 2121 TransCanada, suite 500, Dorval (Quebec), !CANADA, H9P 1J3; or send e-mail to service.rpn@ec.gc.ca !-------------------------------------- LICENCE END -------------------------------------- !subroutine hintvec2(psrcuu,psrcvv,kscol,ksrcgid 6 & ,PTRGUU,PTRGVV,ktcol,ktrggid,klev,cddegint) #if defined (DOC) * ***s/r hintvec * * * Author : S. Pellerin ARMA/AES Avril 2000 * Revision: * * JM Belanger CMDA/SMC Aug 2000 * . 32 bits conversion * S. Pellerin, ARMA, August 2002 * . Fix double declarations ** Purpose: vectorial horizontal interpolation * * * *Arguments * input: * psrcuu(kscol,klev) : source U componant * psrcvv(kscol,klev) : source V componant * kscol : number of collumns of source * ksrcgid : grid id of source * ktcol : number of collumns of target * ktrggid : grid id of target * klev : number of horizontal slabs to interpolate * cddegint : degree of horizontal interpolation * ouput: * ptrguu(ktcol,klev) : target U componant * ptrgvv(ktcol,klev) : target V componant * #endif implicit none *implicits integer ksrcgid,ktrggid integer kscol,ktcol,klev real*8 psrcuu(kscol,klev),ptrguu(ktcol,klev),psrcvv(kscol,klev) real*8 ptrgvv(ktcol,klev) character*(*) cddegint c c---------------------------------------------------------- c integer ier,ezsetopt,vezuvint,ezdefset,ezgprm,iset integer clgrtyp,ini,inj,iig1,iig2,iig3,iig4,ji,itrg,jk real*8 zsrc(1),ztrguu(1),ztrgvv(1) pointer (pzsrc,zsrc) pointer (pztrguu,ztrguu) pointer (pztrgvv,ztrgvv) c c Setting degree of horizontal interpolations ier = ezsetopt('INTERP_DEGREE',cddegint) iset = ezdefset(ktrggid,ksrcgid) do jk = 1,klev ier = vezuvint(ptrguu(1,jk),ptrgvv(1,jk),psrcuu(1,jk) & ,psrcvv(1,jk),ktcol,kscol) enddo return end