| 1 | |
|---|
| 2 | |
|---|
| 3 | SPAWN, '\rm param_plot.idl' |
|---|
| 4 | SPAWN, "echo 'intervaly=5.' > param_plot.idl" |
|---|
| 5 | SPAWN, "echo 'intervalx=5.' >> param_plot.idl" |
|---|
| 6 | |
|---|
| 7 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 8 | file='gcm' & header='' & nlines_header=4 & ncol = 2 |
|---|
| 9 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 10 | nlines = FILE_LINES(file)-nlines_header & data=FLTARR(ncol,nlines) |
|---|
| 11 | OPENR, lun, file, /GET_LUN & READF, lun, header & READF, lun, header & READF, lun, header & READF, lun, header |
|---|
| 12 | READF, lun, data & CLOSE, lun & FREE_LUN, lun |
|---|
| 13 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 14 | altitude = reform(data(0,*)) |
|---|
| 15 | temp = reform(data(1,*)) |
|---|
| 16 | data = 0. |
|---|
| 17 | |
|---|
| 18 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 19 | file='meso' & header='' & nlines_header=4 & ncol = 2 |
|---|
| 20 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 21 | nlines = FILE_LINES(file)-nlines_header & data=FLTARR(ncol,nlines) |
|---|
| 22 | OPENR, lun, file, /GET_LUN & READF, lun, header & READF, lun, header & READF, lun, header & READF, lun, header |
|---|
| 23 | READF, lun, data & CLOSE, lun & FREE_LUN, lun |
|---|
| 24 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 25 | altitude2 = reform(data(0,*)) |
|---|
| 26 | temp2 = reform(data(1,*)) |
|---|
| 27 | data = 0. |
|---|
| 28 | |
|---|
| 29 | |
|---|
| 30 | ; |
|---|
| 31 | ; extremes |
|---|
| 32 | ; |
|---|
| 33 | temp = temp - 12. |
|---|
| 34 | temp2 = temp2 - 12. |
|---|
| 35 | print, 'ATTENTION CACA CACA CACA' |
|---|
| 36 | |
|---|
| 37 | what_I_plot = temp ;- temp2 |
|---|
| 38 | column = altitude |
|---|
| 39 | alt = [20.,60.] |
|---|
| 40 | minfield_init = 100. ;-20 ;100. |
|---|
| 41 | maxfield_init = 200. ;20. ;200. |
|---|
| 42 | overplot = temp2 |
|---|
| 43 | overplot_column = altitude2 |
|---|
| 44 | discrete = 0 |
|---|
| 45 | title_user = 'GCM vs. mesoscale simulations' |
|---|
| 46 | title_axis = ['Temperature (K)','Altitude above MOLA (km)'] |
|---|
| 47 | mention = '' |
|---|
| 48 | |
|---|
| 49 | profile, $ |
|---|
| 50 | what_I_plot, $ ; 1D vertical profile |
|---|
| 51 | column, $ ; altitudes |
|---|
| 52 | alt=alt, $ ; altitude range [altmin, altmax] |
|---|
| 53 | minfield=minfield_init, $ ; minimum value of plotted field (=0: calculate) |
|---|
| 54 | maxfield=maxfield_init, $ ; maximum value of plotted field (=0: calculate) |
|---|
| 55 | inprofile=overplot, $ ; another vertical profile to overplot |
|---|
| 56 | incolumn=overplot_column, $ ; altitudes of the other vertical profile (in case /= column) |
|---|
| 57 | discrete=discrete, $ ; show the profile points (= type of points in !psym) |
|---|
| 58 | title_plot=title_user, $ ; title of the plot ('Profile' is default) |
|---|
| 59 | title_axis=title_axis, $ ; title of the [x,y] axis (['Field','Altitude'] is default) |
|---|
| 60 | mention=mention ; add text precision within the plot window (default is nothing or '') |
|---|
| 61 | |
|---|
| 62 | |
|---|
| 63 | |
|---|
| 64 | ; |
|---|
| 65 | ; condensation |
|---|
| 66 | ; |
|---|
| 67 | latcond=5.9e5 |
|---|
| 68 | tcond1mb=136.27 |
|---|
| 69 | r=192. |
|---|
| 70 | bcond=1./tcond1mb |
|---|
| 71 | acond=r/latcond |
|---|
| 72 | press=610.*exp(-altitude/10.) |
|---|
| 73 | overplot=1./(bcond-acond*alog(.0095*press)) |
|---|
| 74 | overplot_column=altitude |
|---|
| 75 | oplot, overplot, overplot_column, psym=5 |
|---|
| 76 | |
|---|
