|
Last change
on this file since 163 was
163,
checked in by pinsard, 15 years ago
|
|
going on modify unformal header to idldoc 2. header syntax
|
|
File size:
547 bytes
|
| Line | |
|---|
| 1 | ;+ |
|---|
| 2 | ; |
|---|
| 3 | ; @param MSF {in} |
|---|
| 4 | ; msf est une tableau 3d [jpl,jpk,nt] ou nt est le nombre de pas de |
|---|
| 5 | ; |
|---|
| 6 | ; @version |
|---|
| 7 | ; $Id$ |
|---|
| 8 | ; |
|---|
| 9 | ;- |
|---|
| 10 | PRO nadw, msf, nadw, lat, prof |
|---|
| 11 | @common |
|---|
| 12 | ; |
|---|
| 13 | temps |
|---|
| 14 | ; |
|---|
| 15 | N = (size(msf))[3] ; recupere nt |
|---|
| 16 | nadw = fltarr(N) |
|---|
| 17 | lat = fltarr(N) |
|---|
| 18 | prof = fltarr(N) |
|---|
| 19 | |
|---|
| 20 | FOR k = 0, N-1 DO BEGIN |
|---|
| 21 | z = msf(*, *, k) |
|---|
| 22 | nadw(k) = max(z) ; prend le max de msf |
|---|
| 23 | ; |
|---|
| 24 | ; |
|---|
| 25 | ; |
|---|
| 26 | indice = (where(z EQ flow(k)))[0] |
|---|
| 27 | lat(k) = gphit(90, indice MOD jpl) ; recupere la latitude |
|---|
| 28 | prof(k) = zdept(indice/jpl) ; et la profondeur |
|---|
| 29 | ENDFOR |
|---|
| 30 | |
|---|
| 31 | end |
|---|
| 32 | |
|---|
Note: See
TracBrowser
for help on using the repository browser.
