| 1 | ;pro try |
|---|
| 2 | |
|---|
| 3 | |
|---|
| 4 | restore, filename='getturb.dat' |
|---|
| 5 | |
|---|
| 6 | ;;wt, tke, ztke, h, ht, t, p, pt, stst, localtime, xtke, ytke, wmax, wmin, depressions, psmin, filename='getturb.dat' |
|---|
| 7 | |
|---|
| 8 | |
|---|
| 9 | nz = n_elements(wt(*,0)) |
|---|
| 10 | nt = n_elements(wt(0,*)) |
|---|
| 11 | |
|---|
| 12 | ; |
|---|
| 13 | ; nz: colonnes |
|---|
| 14 | ; nt: lignes |
|---|
| 15 | ; |
|---|
| 16 | |
|---|
| 17 | openw, lun, 'LES_MERIDIANI_vert_levels.txt', /GET_LUN |
|---|
| 18 | printf, lun, h, format='(F12.6)' |
|---|
| 19 | close, lun |
|---|
| 20 | free_lun, lun |
|---|
| 21 | |
|---|
| 22 | openw, lun, 'LES_MERIDIANI_local_times.txt', /GET_LUN |
|---|
| 23 | printf, lun, localtime, format='(F12.6)' |
|---|
| 24 | close, lun |
|---|
| 25 | free_lun, lun |
|---|
| 26 | |
|---|
| 27 | openw, lun, 'LES_MERIDIANI_TKE.txt', /GET_LUN |
|---|
| 28 | printf, lun, tke, format='('+string(nz,'(I0)')+'F12.6)' |
|---|
| 29 | close, lun |
|---|
| 30 | free_lun, lun |
|---|
| 31 | |
|---|
| 32 | openw, lun, 'LES_MERIDIANI_xTKE.txt', /GET_LUN |
|---|
| 33 | printf, lun, xtke, format='('+string(nz,'(I0)')+'F12.6)' |
|---|
| 34 | close, lun |
|---|
| 35 | free_lun, lun |
|---|
| 36 | |
|---|
| 37 | openw, lun, 'LES_MERIDIANI_yTKE.txt', /GET_LUN |
|---|
| 38 | printf, lun, ytke, format='('+string(nz,'(I0)')+'F12.6)' |
|---|
| 39 | close, lun |
|---|
| 40 | free_lun, lun |
|---|
| 41 | |
|---|
| 42 | openw, lun, 'LES_MERIDIANI_zTKE.txt', /GET_LUN |
|---|
| 43 | printf, lun, ztke, format='('+string(nz,'(I0)')+'F12.6)' |
|---|
| 44 | close, lun |
|---|
| 45 | free_lun, lun |
|---|
| 46 | |
|---|
| 47 | openw, lun, 'LES_MERIDIANI_Wmax.txt', /GET_LUN |
|---|
| 48 | printf, lun, wmax, format='('+string(nz,'(I0)')+'F12.6)' |
|---|
| 49 | close, lun |
|---|
| 50 | free_lun, lun |
|---|
| 51 | |
|---|
| 52 | openw, lun, 'LES_MERIDIANI_Wmin.txt', /GET_LUN |
|---|
| 53 | printf, lun, wmin, format='('+string(nz,'(I0)')+'F12.6)' |
|---|
| 54 | close, lun |
|---|
| 55 | free_lun, lun |
|---|
| 56 | |
|---|
| 57 | openw, lun, 'LES_MERIDIANI_Tpot.txt', /GET_LUN |
|---|
| 58 | printf, lun, t, format='('+string(nz,'(I0)')+'F12.6)' |
|---|
| 59 | close, lun |
|---|
| 60 | free_lun, lun |
|---|
| 61 | |
|---|
| 62 | restore, filename='addturb.dat' |
|---|
| 63 | openw, lun, 'LES_MERIDIANI_VELmax.txt', /GET_LUN |
|---|
| 64 | printf, lun, VELMAX, format='('+string(nz,'(I0)')+'F12.6)' |
|---|
| 65 | close, lun |
|---|
| 66 | free_lun, lun |
|---|
| 67 | |
|---|
| 68 | |
|---|
| 69 | |
|---|
| 70 | ;yeye = fltarr(nz,nt) |
|---|
| 71 | ;openr, lun, 'LES_MERIDIANI_Tpot.txt', /GET_LUN |
|---|
| 72 | ;readf, lun, yeye |
|---|
| 73 | ;close, lun |
|---|
| 74 | ;free_lun, lun |
|---|
| 75 | ;contour, transpose(yeye), nlevels=20 |
|---|
| 76 | |
|---|
| 77 | yeye2 = fltarr(nz,nt) |
|---|
| 78 | openr, lun, 'LES_MERIDIANI_TKE.txt', /GET_LUN |
|---|
| 79 | readf, lun, yeye2 |
|---|
| 80 | close, lun |
|---|
| 81 | free_lun, lun |
|---|
| 82 | contour, transpose(yeye2), nlevels=20 |
|---|
| 83 | |
|---|
| 84 | |
|---|
| 85 | ;end |
|---|
