| [85] | 1 | pro script_con |
|---|
| 2 | |
|---|
| 3 | ;;;;;;; |
|---|
| 4 | file_user='/home/aslmd/EXOMARS/LMD_MM_MARS_Meridiani_S114_Ls247_LT0-23.nc' |
|---|
| 5 | var_user = 'TSURF' |
|---|
| 6 | var_user2 = 'HGT' |
|---|
| 7 | nt = 5 |
|---|
| 8 | format = '(I0)' |
|---|
| 9 | pal=33 |
|---|
| 10 | ;;;;;;; |
|---|
| 11 | |
|---|
| 12 | PS_START, file='toto.ps' |
|---|
| 13 | |
|---|
| 14 | |
|---|
| 15 | cdfid = ncdf_open(file_user) |
|---|
| 16 | varid=ncdf_varid(cdfid,var_user) |
|---|
| 17 | ncdf_varget, cdfid, varid, champ |
|---|
| 18 | |
|---|
| 19 | cdfid = ncdf_open(file_user) |
|---|
| 20 | varid=ncdf_varid(cdfid,var_user2) |
|---|
| 21 | ncdf_varget, cdfid, varid, champ2 |
|---|
| 22 | |
|---|
| 23 | cdfid = ncdf_open(file_user) |
|---|
| 24 | varid=ncdf_varid(cdfid,'XLONG') |
|---|
| 25 | ncdf_varget, cdfid, varid, lon |
|---|
| 26 | |
|---|
| 27 | cdfid = ncdf_open(file_user) |
|---|
| 28 | varid=ncdf_varid(cdfid,'XLAT') |
|---|
| 29 | ncdf_varget, cdfid, varid, lat |
|---|
| 30 | |
|---|
| 31 | |
|---|
| 32 | |
|---|
| 33 | nx = 25 |
|---|
| 34 | ny = 25 |
|---|
| 35 | what_I_plot = reform(champ(nx,ny,*)) |
|---|
| 36 | |
|---|
| 37 | plot, what_I_plot |
|---|
| 38 | |
|---|
| 39 | PS_END, /PNG |
|---|
| 40 | |
|---|
| 41 | end |
|---|
| 42 | |
|---|
| 43 | |
|---|
| 44 | |
|---|
| 45 | what_I_plot = reform(champ(*,*,nt)) |
|---|
| 46 | overcontour = reform(champ2(*,*,nt)) |
|---|
| 47 | longi = reform(lon(*,0)) |
|---|
| 48 | lati = reform(lat(0,*)) |
|---|
| 49 | |
|---|
| 50 | map_latlon, $ |
|---|
| 51 | what_I_plot, $ ; 2D field |
|---|
| 52 | longi, $ ; 1D latitude OR 2D |
|---|
| 53 | lati, $ ; 1D longitude OR 2D |
|---|
| 54 | ; minfield=minfield_init, $ ; minimum value of plotted field (=0: calculate) |
|---|
| 55 | ; maxfield=maxfield_init, $ ; maximum value of plotted field (=0: calculate) |
|---|
| 56 | overcontour=overcontour, $ ; another 2D field to overplot with contour lines (=0: no) |
|---|
| 57 | ; overvector_x=overvector_x, $ ; wind vector - x component (=0: no) |
|---|
| 58 | ; overvector_y=overvector_y, $ ; wind vector - y component (=0: no) |
|---|
| 59 | ct=pal, $ ; color table (33-rainbow is default) |
|---|
| 60 | ; colors=colors, $ ; number of colors/levels (32 is default) |
|---|
| 61 | ; title=title_user, $ ; title of the plot ('' is default) |
|---|
| 62 | format=format ; format of colorbar annotations ('(F6.2)' is default) |
|---|
| 63 | |
|---|
| 64 | |
|---|
| 65 | PS_END, /PNG |
|---|
| 66 | |
|---|
| 67 | end |
|---|
