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 |
---|