[207] | 1 | 'reinit' |
---|
| 2 | 'open real_z.ctl' |
---|
| 3 | |
---|
| 4 | ** need input data on height levels |
---|
| 5 | |
---|
| 6 | 'set mpdset hires' |
---|
| 7 | 'set display color white' |
---|
| 8 | 'run rgbset.gs' |
---|
| 9 | |
---|
| 10 | say 'Create gif images as well (1=yes ; 0=no)' |
---|
| 11 | pull ans |
---|
| 12 | frame = 1 |
---|
| 13 | |
---|
| 14 | 'q file' |
---|
| 15 | rec=sublin(result,5) |
---|
| 16 | _endtime=subwrd(rec,12) |
---|
| 17 | _endlev=subwrd(rec,9) |
---|
| 18 | |
---|
| 19 | runscript = 1 |
---|
| 20 | dis_t = 1 |
---|
| 21 | |
---|
| 22 | while(runscript) |
---|
| 23 | dis_z = 1 |
---|
| 24 | 'set t ' dis_t |
---|
| 25 | 'q dims' |
---|
| 26 | rec=sublin(result,5) |
---|
| 27 | _analysis=subwrd(rec,6) |
---|
| 28 | say 'Time is ' _analysis |
---|
| 29 | |
---|
| 30 | runlevel = 1 |
---|
| 31 | while(runlevel) |
---|
| 32 | 'set z ' dis_z |
---|
| 33 | 'q dims' |
---|
| 34 | rec=sublin(result,4) |
---|
| 35 | newz=subwrd(rec,6) |
---|
| 36 | say 'Level is ' newz |
---|
| 37 | |
---|
| 38 | ** Looking for 2 km data |
---|
| 39 | if (newz=2) |
---|
| 40 | say 'Plotting' |
---|
| 41 | 'c' |
---|
| 42 | 'set grads off' |
---|
| 43 | 'set gxout shaded' |
---|
| 44 | 'set clevs 10 20 30 40 50 60 70 80 90 100' |
---|
| 45 | 'set ccols 0 0 0 0 33 34 35 36 37 38 39' |
---|
| 46 | 'd rh' |
---|
| 47 | 'set gxout contour' |
---|
| 48 | 'set ccolor 4' |
---|
| 49 | 'set cthick 6' |
---|
| 50 | 'set cint 2' |
---|
| 51 | 'd p' |
---|
| 52 | 'set ccolor 2' |
---|
| 53 | 'set cthick 6' |
---|
| 54 | 'set cint 2' |
---|
| 55 | 'd tc' |
---|
| 56 | 'set strsiz .2' |
---|
| 57 | 'set string 1 l 6' |
---|
| 58 | 'draw string 2.7 8.35 Height (m,blue), T(C,red), RH(color)' |
---|
| 59 | 'set strsiz .15' |
---|
| 60 | 'set string 1 l 3' |
---|
| 61 | 'draw string 7.5 8.0 ' _analysis |
---|
| 62 | 'draw string 1.5 8.0 ' newz ' km' |
---|
| 63 | 'run cbar.gs' |
---|
| 64 | if(ans) |
---|
| 65 | 'printim realz'frame'.gif gif ' |
---|
| 66 | frame=frame+1 |
---|
| 67 | endif |
---|
| 68 | pull dummy |
---|
| 69 | endif |
---|
| 70 | |
---|
| 71 | ** Looking for 5 km data |
---|
| 72 | if (newz=5) |
---|
| 73 | say 'Plotting' |
---|
| 74 | 'c' |
---|
| 75 | 'set ccolor 4' |
---|
| 76 | 'set cthick 6' |
---|
| 77 | 'set cint 2' |
---|
| 78 | 'd p' |
---|
| 79 | 'set ccolor 2' |
---|
| 80 | 'set cthick 6' |
---|
| 81 | 'set cint 2' |
---|
| 82 | 'd tc' |
---|
| 83 | 'set strsiz .2' |
---|
| 84 | 'set string 1 l 6' |
---|
| 85 | 'draw string 3.7 8.35 Height (m,blue), T(C,red)' |
---|
| 86 | 'set strsiz .15' |
---|
| 87 | 'set string 1 l 3' |
---|
| 88 | 'draw string 7.5 8.0 ' _analysis |
---|
| 89 | 'draw string 1.5 8.0 ' newz ' km' |
---|
| 90 | if(ans) |
---|
| 91 | 'printim realz'frame'.gif gif ' |
---|
| 92 | frame=frame+1 |
---|
| 93 | endif |
---|
| 94 | pull dummy |
---|
| 95 | endif |
---|
| 96 | |
---|
| 97 | ** Looking for 10 km data |
---|
| 98 | if (newz=10) |
---|
| 99 | say 'Plotting' |
---|
| 100 | 'c' |
---|
| 101 | 'set ccolor 4' |
---|
| 102 | 'set cthick 6' |
---|
| 103 | 'set cint 2' |
---|
| 104 | 'd p' |
---|
| 105 | 'set ccolor 2' |
---|
| 106 | 'set cthick 6' |
---|
| 107 | 'set cint 2' |
---|
| 108 | 'd tc' |
---|
| 109 | 'set strsiz .2' |
---|
| 110 | 'set string 1 l 6' |
---|
| 111 | 'draw string 3.7 8.35 Height (m,blue), T(C,red)' |
---|
| 112 | 'set strsiz .15' |
---|
| 113 | 'set string 1 l 3' |
---|
| 114 | 'draw string 7.5 8.0 ' _analysis |
---|
| 115 | 'draw string 1.5 8.0 ' newz ' km' |
---|
| 116 | if(ans) |
---|
| 117 | 'printim realz'frame'.gif gif ' |
---|
| 118 | frame=frame+1 |
---|
| 119 | endif |
---|
| 120 | pull dummy |
---|
| 121 | endif |
---|
| 122 | |
---|
| 123 | ** Looking for 16 km data |
---|
| 124 | if (newz=16) |
---|
| 125 | say 'Plotting' |
---|
| 126 | 'c' |
---|
| 127 | 'set ccolor 4' |
---|
| 128 | 'set cthick 6' |
---|
| 129 | 'set cint 2' |
---|
| 130 | 'd p' |
---|
| 131 | 'set ccolor 2' |
---|
| 132 | 'set cthick 6' |
---|
| 133 | 'set cint 2' |
---|
| 134 | 'd tc' |
---|
| 135 | 'set strsiz .2' |
---|
| 136 | 'set string 1 l 6' |
---|
| 137 | 'draw string 3.7 8.35 Height (m,blue), T(C,red)' |
---|
| 138 | 'set strsiz .15' |
---|
| 139 | 'set string 1 l 3' |
---|
| 140 | 'draw string 7.5 8.0 ' _analysis |
---|
| 141 | 'draw string 1.5 8.0 ' newz ' km' |
---|
| 142 | if(ans) |
---|
| 143 | 'printim realz'frame'.gif gif ' |
---|
| 144 | frame=frame+1 |
---|
| 145 | endif |
---|
| 146 | pull dummy |
---|
| 147 | endif |
---|
| 148 | |
---|
| 149 | if ( dis_z=_endlev ) |
---|
| 150 | runlevel=0 |
---|
| 151 | endif |
---|
| 152 | dis_z = dis_z + 1 |
---|
| 153 | endwhile |
---|
| 154 | |
---|
| 155 | if ( dis_t=_endtime ) |
---|
| 156 | runscript=0 |
---|
| 157 | endif |
---|
| 158 | dis_t = dis_t + 1 |
---|
| 159 | endwhile |
---|
| 160 | |
---|