topprojectondepth
Without loop
Return value
A 2d array which is the projection of the 3d array following depthes indicated by depthin
Parameters
arrayin
depthin
Examples
we build a possible depthes array
IDL> a=gdept[jpk-1]/(1.*jpi*jpj)*findgen(jpi,jpj)
We build an array to project on these depthes. For the test,
we build a 3d array whose each vector following z is the depth.
IDL> arraytest=replicate(1,jpi*jpj)#gdept
IDL> arraytest=reform(arraytest,jpi,jpj,jpk, /over)
We test the projection of the depth array on the depth...
IDL> plt, 1e6*(a-projectondepth(arraytest,a)),/nocontour
->null field at 1e-6 pres
verifcation projecting the temperature of 20°C for example...
Version history
Version
$Id: projectondepth.pro 142 2006-07-21 12:47:49Z navarro $
History
Sebastien Masson (smasson@lodyc.jussieu.fr)
15/6/2000
Known issues
Restrictions
points at !values.f_nan impossible calculation. Land points masked at valmask.
Other attributes
Uses routines
common.pro