| 1 | ; |
|---|
| 2 | ; time series of rms deviation on a sigma surface (Paul William) |
|---|
| 3 | ; |
|---|
| 4 | |
|---|
| 5 | sizefld=size(fld) |
|---|
| 6 | fld=fld(*,*,round((sig_surface-sig_min)/sig_del),*) |
|---|
| 7 | fld = reform(fld, sizefld(1), sizefld(2), sizefld(4)) |
|---|
| 8 | fld2 = fltarr(sizefld(4)) |
|---|
| 9 | ;FOR count = 0, sizefld(4)-1 DO fld2(count) = mean(fld(*,*,count), /nan) |
|---|
| 10 | FOR count = 0, sizefld(4)-1 DO fld2(count) = stddev(fld(*,*,count), /nan) |
|---|
| 11 | fld = fld2 |
|---|
| 12 | mask_z, fld, cmd, boite_plt1d, dimplot, legz |
|---|
| 13 | niveau = 1 |
|---|
| 14 | sin = '' |
|---|
| 15 | take_log = ',take_log=0' |
|---|
| 16 | pltcmd = 'plt1d,{a:fld, d:'''+date_txt+''', n:'''+varlegend+' Stddev '+legz+legendsurf+''', u:'''+field.units+''', g:'''+vargrid+'''},'''+plt1dtyp+''''+mincmd+maxcmd+intcmd+', boite=boite_plt1d'+overc+sin+com_strplt+take_log |
|---|
| 17 | IF cmd.timave EQ '1y' AND cmd.exp NE 'SQ2' THEN BEGIN ; curve-fitting and over-plot |
|---|
| 18 | time_years = FLOAT(INDGEN(sizefld(4))) |
|---|
| 19 | weights = replicate(1.0, sizefld(4)) |
|---|
| 20 | params = [fld(0),-0.1,0.] ; initial guess |
|---|
| 21 | fldfit = CURVEFIT(time_years, fld, weights, params, error, FUNCTION_NAME='gfunct') |
|---|
| 22 | PRINT, 'Best fit: sigma = ', params(0), '*exp(', params(1), 't) + ', params(2) |
|---|
| 23 | PRINT, 'errors: ', error |
|---|
| 24 | PRINT, 'e-folding time (years) = ', -1.0/params(1) |
|---|
| 25 | look = ', linestyle=0,xthick=2,ythick=2,zthick=2' |
|---|
| 26 | com_strplt = ',petit = ['+string(win[0])+','+string(win[1])+','+string(win[2])+']'+filltxt+nocoltxt+', LANDSCAPE = '+string(landscape)+', NOCOLORBAR = '+string(nocolorbar)+', NOERASE = '+string(noerase)+look+labstr+c_anot+colbar+cal_typ+',cont_thick=2'+type_yz+window_number+tit_str+xlogax |
|---|
| 27 | pltcmd = 'plt1d,{a:fld, d:'''+date_txt+''', n:'''+varlegend+' Stddev (tau ='+strcompress(-1.0/params(1))+'yr) '+legz+legendsurf+''', u:'''+field.units+''', g:'''+vargrid+'''},'''+plt1dtyp+''''+mincmd+maxcmd+intcmd+', boite=boite_plt1d'+overc+sin+com_strplt+take_log |
|---|
| 28 | res = execute(pltcmd) |
|---|
| 29 | noerase = 1 |
|---|
| 30 | look = ', linestyle=2,xthick=2,ythick=2,zthick=2' |
|---|
| 31 | com_strplt = ',petit = ['+string(win[0])+','+string(win[1])+','+string(win[2])+']'+filltxt+nocoltxt+', LANDSCAPE = '+string(landscape)+', NOCOLORBAR = '+string(nocolorbar)+', NOERASE = '+string(noerase)+look+labstr+c_anot+colbar+cal_typ+',cont_thick=2'+type_yz+window_number+tit_str+xlogax |
|---|
| 32 | pltcmd = 'plt1d,{a:fldfit, d:'''+date_txt+''', n:'''+varlegend+' Stddev (tau ='+strcompress(-1.0/params(1))+'yr) '+legz+legendsurf+''', u:'''+field.units+''', g:'''+vargrid+'''},'''+plt1dtyp+''''+mincmd+maxcmd+intcmd+', boite=boite_plt1d'+overc+sin+com_strplt+take_log |
|---|
| 33 | ENDIF |
|---|
