FSCINT                    Routines d'interpolation                       FSCINT

NOM
   RGUVINT - Interpolation vectorielle de champs a partir 
             d'une grille source reguliere. 

USAGE
   call rguvint(SPDOUT, DIROUT, li, lj, xlat, xlon, uuin, vvin, ni, nj, 
		grtyp, ig1, ig2, ig3, ig4, sws)

   integer li, lj     - Dimensions de la grille de destination
   real SPDOUT(li,lj) - Champ de destination contenant la vitesse des vents
   real DIROUT(li,lj) - Champ de destination contenant la direction des vents
   real xlat(li,lj)   - Latitudes des points sur la grille de destination
   real xlon(li,lj)   - Longitudes des points sur la grille de destination

   integer ni, nj     - Dimensions de la grille source
   real uuin(ni, nj)  - Champ contenant la composante U des vents sur la
 			grille source
   real vvin(ni, nj)  - Champ contenant la composante V des vents sur la
                        grille source
   
   character*1 grtyp  - Type de la grille source ('A','B','G','L','N','S')
   integer ig1, ig2,
           ig3, ig4   - Descripteurs de la grille source
   logical sws        - Drapeau indiquant si l'on desire obtenir la vitesse et
                        la direction des vents ou la vitesse seulement 
			si .true., on calcule la vitesse et la direction
			si .false., on calcule la vitesse seulement.
DESCRIPTION
   Cette routine interpole la composante U et V d'une grille de 
   dimension (ni,nj) sur une grille de dimension (li,lj). Les valeurs
   interpolees reviennent sous forme de vitesse et de direction.

   La direction des vents est definie en degres, par rapport a une grille
   latitude-longitude. La definition des angles suit la convention 
   meteorologique, c'est a dire qu'un vent du nord sera a 0 deg, un vent d'est
   a 90, un vent du sud a 180 et un vent d'ouest a 270 deg.

   Cette routine devrait etre utilisee conjointement avec la routine "rgll2gd",
   qui convertit les vents retournes par "rguvint" en composantes u et v 
   definies sur la grille de destination. 

EXEMPLES D'UTILISATION
      1) On desire obtenir la composante U et V des vents sur la grille 
         de destination.
      call rguvint(UUOUT, VVOUT, li, lj, xlat, xlon, uuin, vvin, ni, nj, 
	           gtypsrc, ig1src, ig2src, ig3src, ig4src, .true.)
      call rgll2gd(UUOUT, VVOUT, xlon, li, lj, 
                   gtypdes, ig1des, ig2des, ig3des, ig4des)

   ou  "gtypsrc, ig1src, ig2src, ig3src, ig4src" sont les descripteurs de la 
   grille source, et "gtypdes, ig1des, ig2des, ig3des, ig4des" sont les 
   descripteurs de la grille de destination.

   Apres l'execution de ces deux routines, UUOUT et VVOUT contiennent les
   composantes U et V des vents sur la grille de destination. 

      2) On desire obtenir seulement le module des vents sur la grille 
         de destination.
      call rguvint(SPDOUT, BIDON, li, lj, xlat, xlon, uuin, vvin, ni, nj, 
	           gtypsrc, ig1src, ig2src, ig3src, ig4src, .false.)

   La methode d'interpolation est de type lagrangienne bi-cubique, 
   bi-lineaire ou plus proche voisin. 

   On doit fournir les latitudes et longitudes (en degres) de la grille 
   de destination dans les champs xlat(li,lj) et xlon(li,lj). 

   Le degre d'interpolation est cubique par defaut; il peut etre modifie
   par un appel aux routines "rgoptc" ou "rgopti". On peut aussi controler, 
   par un appel aux memes routines, le comportement a adopter dans les cas
   d'extrapolation.

COMPILATION
   Pour le moment, le module d'interpolation est contenu dans un fichier ".o"
   residant sur $ARMNLIB/lib. Le module sera eventuellement integre dans
   RMNLIB.

   Pour utiliser le module sur le nec, cidsv02, cidsv07, cidsv08, faire:
   f77 prog.f -o prog $ARMNLIB/lib/fscint2.o $ARMNLIB/lib/rmnxlib.a

   Sur le cray, faire:
   ACCESS, DN=FSCLIB2, ID=ARMNYRC, OWN=ARMN000.
   ...
   LD2, LIB=FSCLIB2:AAAALIB, ...
   
AUTEUR
   Yves Chartier, RPN, 421-4736, ychartier@cid.aes.doe.CA

DERNIERE MISE A JOUR DE LA DOCUMENTATION
   19 decembre 1991

VOIR AUSSI
   igscint, iguvint, rgll2gd, rgscint, rgoptc, rgopti, rgoptr, fstarg1