- Timestamp:
- 11/12/09 17:15:09 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/tools/density_binning/binning_neutral_and_co/msfs.pro
r2 r161 21 21 pe3v = read_ncdf('e3v_ps',0,/timestep,iodir=ioMESH,/nostruct,/tout,filename=fich_msh) 22 22 vvf = vv * reform(e1v[*]#replicate(1,jpk),jpi,jpj,jpk) * pe3v *vmask() 23 24 23 ;------------------------- 25 24 ; gibraltar cross land adv: 5 Sv … … 32 31 vvf(167,100, 20 ) = vvf(167,100, 20 )+2.5e+6 *vvmask(167,100, 20 ) 33 32 vvf(167,100, 21 ) = vvf(167,100, 21 )+0.8e+6/14.*vvmask(167,100, 21 ) 34 35 33 ;------------------------- 36 34 if keyword_set(mask) then begin &$ … … 38 36 if (size(mask))[0] eq 2 then vvf=vvf*reform(mask[*]#replicate(1,jpk),jpi,jpj,jpk) &$ 39 37 endif &$ 40 41 38 ;------------------------- 42 39 ; crude summation … … 45 42 x1_bin = fltarr(jpj, N_s) ; output array (transport function of density) 46 43 x1_bin(*,*) = 0.e0 47 48 44 FOR is = N_s-1,0,-1 DO begin &$ 49 45 dmsk= density GT s_s(is) &$ … … 51 47 ENDFOR &$ 52 48 msfs=x1_bin *1.e-6 53 54 49 ;------------------------ 55 50 return, msfs
Note: See TracChangeset
for help on using the changeset viewer.