<<prev file | next file >>    view single page | view frames    summary: fields | routine    details: routine

ToBeReviewed/CALCULS/

projectondepth.pro

Allows to project a 3d field following a depth array.

topprojectondepth Without loop

result = projectondepth(arrayin, depthin)

Return value

A 2d array which is the projection of the 3d array following depthes indicated by depthin

Parameters

arrayin       

It is a 3d array whose 3rd dimension must be equal to jpk

depthin       

It is a 2d array indicating for each point n, at which depth to project

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
Produced by IDLdoc 2.0.