| 1 | pro mawd |
|---|
| 2 | |
|---|
| 3 | ;;;;;;; |
|---|
| 4 | experiment='a035' |
|---|
| 5 | ;experiment='a040' |
|---|
| 6 | ;experiment='a045' |
|---|
| 7 | ;;f_user='/donnees/aslmd/SVN/trunk/mesoscale/TMPDIR/GCMINI/mawd.nc' |
|---|
| 8 | ;;f_user='/tmp7/aslmd/mawd.nc' |
|---|
| 9 | f_user='/d5/aslmd/WATERCYCLE/mawd_'+experiment+'.nc' |
|---|
| 10 | v_user = 'mtot' |
|---|
| 11 | x_user = 'Time' |
|---|
| 12 | y_user = 'latitude' |
|---|
| 13 | transp = 'yes' |
|---|
| 14 | tes = 'no' |
|---|
| 15 | ;;;;;;; |
|---|
| 16 | ;f_user='/donnees/aslmd/TES/TES.SeasonalClimatology.nc' |
|---|
| 17 | ;v_user = 'water' |
|---|
| 18 | ;x_user = 'time' |
|---|
| 19 | ;y_user = 'latitude' |
|---|
| 20 | ;transp = 'yes' |
|---|
| 21 | ;tes = 'yes' |
|---|
| 22 | ;;;;;;; |
|---|
| 23 | |
|---|
| 24 | |
|---|
| 25 | what_I_plot=0. |
|---|
| 26 | overcontour=0 |
|---|
| 27 | @mawd_inc.pro |
|---|
| 28 | print, lev |
|---|
| 29 | SPAWN, '\rm param_plot.idl ; cp mawd_inc.pro param_plot.idl' |
|---|
| 30 | |
|---|
| 31 | PS_START, file='mawd_'+experiment+'.ps' |
|---|
| 32 | ; !P.Charsize = 1.2 |
|---|
| 33 | ; !p.charthick = 2.0 |
|---|
| 34 | ; !p.thick = 2.0 |
|---|
| 35 | ; !x.thick = 2.0 |
|---|
| 36 | ; !y.thick = 2.0 |
|---|
| 37 | |
|---|
| 38 | cdfid = ncdf_open(f_user) |
|---|
| 39 | |
|---|
| 40 | varid=ncdf_varid(cdfid,v_user) |
|---|
| 41 | ncdf_varget, cdfid, varid, champ |
|---|
| 42 | |
|---|
| 43 | varid=ncdf_varid(cdfid,y_user) |
|---|
| 44 | ncdf_varget, cdfid, varid, yy |
|---|
| 45 | |
|---|
| 46 | varid=ncdf_varid(cdfid,x_user) |
|---|
| 47 | ncdf_varget, cdfid, varid, xx |
|---|
| 48 | |
|---|
| 49 | if (tes eq 'no') then begin |
|---|
| 50 | ;;; en precip-microns |
|---|
| 51 | champ = champ * 1.e6 / 917. |
|---|
| 52 | ;;;; entre 0 et 360 |
|---|
| 53 | ;xx = xx MOD 360 |
|---|
| 54 | xx = xx - 360.*3. |
|---|
| 55 | endif else begin |
|---|
| 56 | xx = xx - 360.*2. |
|---|
| 57 | endelse |
|---|
| 58 | |
|---|
| 59 | help, champ |
|---|
| 60 | if (tes eq 'no') then begin |
|---|
| 61 | what_I_plot = champ |
|---|
| 62 | endif else begin |
|---|
| 63 | what_I_plot = float(reform(champ(0,*,*))) |
|---|
| 64 | endelse |
|---|
| 65 | help, what_I_plot |
|---|
| 66 | overcontour = what_I_plot |
|---|
| 67 | |
|---|
| 68 | if (tes eq 'yes') then begin |
|---|
| 69 | w=where(what_I_plot eq -1.) |
|---|
| 70 | what_I_plot[w] = !VALUES.F_NAN |
|---|
| 71 | overcontour[w] = !VALUES.F_NAN |
|---|
| 72 | endif |
|---|
| 73 | |
|---|
| 74 | if (transp eq 'yes') then begin |
|---|
| 75 | what_I_plot = transpose(what_I_plot) |
|---|
| 76 | overcontour = transpose(overcontour) |
|---|
| 77 | endif |
|---|
| 78 | |
|---|
| 79 | what_I_plot = smooth(what_I_plot,[10,1]) |
|---|
| 80 | |
|---|
| 81 | map_latlon, $ |
|---|
| 82 | what_I_plot, $ ; 2D field |
|---|
| 83 | xx, $ ; 1D latitude OR 2D |
|---|
| 84 | yy, $ ; 1D longitude OR 2D |
|---|
| 85 | ; minfield=minfield_init, $ ; minimum value of plotted field (=0: calculate) |
|---|
| 86 | ; maxfield=maxfield_init, $ ; maximum value of plotted field (=0: calculate) |
|---|
| 87 | overcontour=overcontour, $ ; another 2D field to overplot with contour lines (=0: no) |
|---|
| 88 | ; overvector_x=overvector_x, $ ; wind vector - x component (=0: no) |
|---|
| 89 | ; overvector_y=overvector_y, $ ; wind vector - y component (=0: no) |
|---|
| 90 | ct=pal, $ ; color table (33-rainbow is default) |
|---|
| 91 | ; colors=colors, $ ; number of colors/levels (32 is default) |
|---|
| 92 | ; title=title_user, $ ; title of the plot ('' is default) |
|---|
| 93 | format=format ; format of colorbar annotations ('(F6.2)' is default) |
|---|
| 94 | |
|---|
| 95 | |
|---|
| 96 | PS_END, /PNG |
|---|
| 97 | |
|---|
| 98 | end |
|---|
