Ignore:
Timestamp:
11/16/09 15:39:10 (15 years ago)
Author:
pinsard
Message:

add compile_opt idl2, strictarrsubs and subsequent modifications

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/tools/density_binning/density_bin_IDL_gm/msf.pro

    r165 r169  
    1414PRO msf, v, sf $ 
    1515    , MSK=msk 
     16; 
     17  compile_opt idl2, strictarrsubs 
    1618; 
    1719@common 
     
    3436IF n_elements(msk) NE 0 THEN BEGIN  
    3537      zvmask = boundperio( (msk +shift(msk, 0, -1) ) < 1 ) 
    36       zvmask = zvmask(*)#vert 
     38      zvmask = zvmask[*]#vert 
    3739      zv = zv*zvmask 
    3840ENDIF  
     
    4446;  calcul du flux 
    4547; 
    46 FOR i = 0, jpi-1 DO BEGIN for j= 0, jpj-1 do begin ze1v(i,j,*) = replicate(e1v(i,j),jpk) & endfor &endfor 
    47 FOR k = 0, jpk-1 DO BEGIN ze3v(*,*,k)=replicate(e3t(k),jpi*jpj) & endfor 
     48FOR i = 0, jpi-1 DO BEGIN for j= 0, jpj-1 do begin ze1v[i,j,*] = replicate(e1v[i,j],jpk) & endfor &endfor 
     49FOR k = 0, jpk-1 DO BEGIN ze3v[*,*,k]=replicate(e3t[k],jpi*jpj) & endfor 
    4850; 
    4951z= -v*ze1v*ze3v 
     
    5557;  calcul de la msf en integrant depuis le fond 
    5658; 
    57 FOR k = jpk-2, 0, -1 DO begin sf(*, k) = sf(*, k+1)+fm(*, k) & endfor 
     59FOR k = jpk-2, 0, -1 DO begin sf[*, k] = sf[*, k+1]+fm[*, k] & endfor 
    5860; 
    5961;  msfmsk est le masque associe a msf (utilise pour les graphiques) 
Note: See TracChangeset for help on using the changeset viewer.