Changeset 169 for trunk/tools/density_binning/density_bin_IDL_gm/msf.pro
- Timestamp:
- 11/16/09 15:39:10 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/tools/density_binning/density_bin_IDL_gm/msf.pro
r165 r169 14 14 PRO msf, v, sf $ 15 15 , MSK=msk 16 ; 17 compile_opt idl2, strictarrsubs 16 18 ; 17 19 @common … … 34 36 IF n_elements(msk) NE 0 THEN BEGIN 35 37 zvmask = boundperio( (msk +shift(msk, 0, -1) ) < 1 ) 36 zvmask = zvmask (*)#vert38 zvmask = zvmask[*]#vert 37 39 zv = zv*zvmask 38 40 ENDIF … … 44 46 ; calcul du flux 45 47 ; 46 FOR i = 0, jpi-1 DO BEGIN for j= 0, jpj-1 do begin ze1v (i,j,*) = replicate(e1v(i,j),jpk) & endfor &endfor47 FOR k = 0, jpk-1 DO BEGIN ze3v (*,*,k)=replicate(e3t(k),jpi*jpj) & endfor48 FOR i = 0, jpi-1 DO BEGIN for j= 0, jpj-1 do begin ze1v[i,j,*] = replicate(e1v[i,j],jpk) & endfor &endfor 49 FOR k = 0, jpk-1 DO BEGIN ze3v[*,*,k]=replicate(e3t[k],jpi*jpj) & endfor 48 50 ; 49 51 z= -v*ze1v*ze3v … … 55 57 ; calcul de la msf en integrant depuis le fond 56 58 ; 57 FOR k = jpk-2, 0, -1 DO begin sf (*, k) = sf(*, k+1)+fm(*, k)& endfor59 FOR k = jpk-2, 0, -1 DO begin sf[*, k] = sf[*, k+1]+fm[*, k] & endfor 58 60 ; 59 61 ; msfmsk est le masque associe a msf (utilise pour les graphiques)
Note: See TracChangeset
for help on using the changeset viewer.