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