;+
;
; @param MSF {in} 
; msf est une tableau 3d [jpl,jpk,nt] ou nt est le nombre de pas de
;
; @uses
; <pro>common</pro>
;
; @version
; $Id$
;
;-
PRO nadw, msf, nadw, lat, prof
;
  compile_opt idl2, strictarrsubs
;
@common
;
temps
;
   N = (size(msf))[3]  ; recupere nt
   nadw = fltarr(N)
   lat = fltarr(N)
   prof = fltarr(N)

   FOR k = 0, N-1 DO BEGIN 
      z = msf[*, *, k]
      nadw[k] = max(z) ; prend le max de msf
;
;
;
      indice = (where(z EQ flow(k)))[0]
; recupere la latitude 
      lat[k] = gphit[90, indice MOD jpl] 
; et la profondeur
      prof[k] = zdept[indice/jpl]
   ENDFOR 

end