NAME:projectondepth
PURPOSE: routine permettant de projeter un champ 3d suivant un
tableau de profondeurs.
CATEGORY: sans boucles
CALLING SEQUENCE:res=projectondepth(arrayin, depthin)
INPUTS:
arrayin: un tableau 3d dont la 3eme dimension doit etre egale
a jpk
depthin: un tableau 2d indiquant n chaque point a quel
profondeur projeter
KEYWORD PARAMETERS:none
OUTPUTS:res: un tableau 2d projection du tableau 3d suivant les
profondeurs indiquees par depthin
COMMON BLOCKS:common.pro
SIDE EFFECTS: points a !values.f_nan qd calcul impossible. points
terres masques a Valmask.
RESTRICTIONS:
EXAMPLE:
on contruit un tableau de profondeurs possibles
IDL> a=gdept[jpk-1]/(1.*jpi*jpj)*findgen(jpi,jpj)
on contruit un tableau a projeter sur ces profondeurs. pour le test
on construit un tableau 3d dont chaque vecteur suivant z est la
profondeur.
IDL> arraytest=replicate(1,jpi*jpj)#gdept
IDL> arraytest=reform(arraytest,jpi,jpj,jpk, /over)
on test la projection du tabeau profondeur sur la profondeur...
IDL> plt, 1e6*(a-projectondepth(arraytest,a)),/nocontour
->champ nul a 1e-6 pres
verifcation en projettant la temperature sur la profondeur de la 20
degres par exemple...