Changeset 161 for trunk/tools/density_binning/density_bin_IDL_gm/msfs.pro
- Timestamp:
- 11/12/09 17:15:09 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/tools/density_binning/density_bin_IDL_gm/msfs.pro
r2 r161 10 10 tmask=replicate(1b,jpi,jpj,jpk) 11 11 domdef,0,jpi-1,0,jpj-1,0,jpk-1,/index 12 13 14 12 if keyword_set(mask) then begin 15 13 print, 'subdomain: use of mask argument' … … 20 18 msf=total(msf*(e1v[*]#replicate(1,jpk)),1,/nan) 21 19 msf=1e-6*msf 22 23 20 mmsk=msf 24 21 indm=where(mmsk NE 0 ) 25 22 mmsk(indm)=1 26 27 23 msf=total(msf,2,/cumulative) 28 24 ;msf=total(msf,2,/cumulative,/nan) … … 36 32 endfor 37 33 endif 38 39 34 ; mask the msf 40 35 msf=msf*mmsk 41 42 36 varname='Meridional Stream Function' 43 37 varunit='Sv' 44 38 vargrid='W' 45 46 47 39 if keyword_set(orca) then begin 48 40 case orca of … … 61 53 ;if n_elements(valmask) eq 0 then valmask=1e20 62 54 ;if terre[0] ne -1 then msf[terre]=valmask 63 64 55 ;------------------------ 65 56 return, msf
Note: See TracChangeset
for help on using the changeset viewer.