create density file

[Gabriel Clauzet <gclauzet@xxxxxx>]

Dear Ferret users,

I calculated the potential densities using the function RHO_UN. I want to 
save all  densities referred to the surface (r0), 1000db (r1), 2000db (r2), 
3000db (r3) and 4000 db (r4) in one file. Does any one could help me ?

Thank in advance
Gabriel

___________________________________________

use temp.nc
use sali.nc

  let D0 = RHO_UN(SALI[D=2],TEMP[D=1],0)-1000
  let D1 = RHO_UN(SALI[D=2],TEMP[D=1],1000)-1000
  let D2 = RHO_UN(SALI[D=2],TEMP[D=1],2000)-1000
  let D3 = RHO_UN(SALI[D=2],TEMP[D=1],3000)-1000
  let D4 = RHO_UN(SALI[D=2],TEMP[D=1],4000)-1000

  let  r1 = D0[Z=1:1000,L=@AVE,X=-25]
  let  r2 = D1[Z=1001:2000,L=@AVE,X=-25]
  let  r3 = D2[Z=2001:3000,L=@AVE,X=-25]
  let  r4 = D3[Z=3001:4000,L=@AVE,X=-25]
  let  r5 = D4[Z=4001:6000,L=@AVE,X=-25]

stat r1
Total # of data points: 2394 (1*171*14*1)
 # flagged as bad  data: 339
 Minimum value: 22.296
 Maximum value: 28.556
 Mean    value: 26.672 (unweighted average)
 Standard deviation: 1.2791

stat r2
Total # of data points: 855 (1*171*5*1)
 # flagged as bad  data: 156
 Minimum value: 32.221
 Maximum value: 33.366
 Mean    value: 32.64 (unweighted average)
 Standard deviation: 0.19643

How can I create one file with r1 in the first 1000 meters, r2 between 1000 
and 2000 meters, r3, r4 and r5 ???