|
Last change
on this file since 798 was
207,
checked in by aslmd, 14 years ago
|
|
MESOSCALE: A GENERAL CLEAN-UP FOLLOWING UPDATING THE USER MANUAL. EVERYTHING ESSENTIAL IS IN MESOSCALE (much lighter than before). EVERYTHING FOR DEVELOPPERS OR EXPERTS IS IN MESOSCALE_DEV.
|
-
Property svn:executable set to
*
|
|
File size:
2.5 KB
|
| Line | |
|---|
| 1 | 'reinit' |
|---|
| 2 | 'open qss.ctl' |
|---|
| 3 | |
|---|
| 4 | ** script needs height level input data |
|---|
| 5 | ** will look for .75 , 1.5. 4 and 8 km data to plot |
|---|
| 6 | |
|---|
| 7 | 'set display color white' |
|---|
| 8 | 'set xlab off' |
|---|
| 9 | 'set ylab off' |
|---|
| 10 | 'run rgbset.gs' |
|---|
| 11 | |
|---|
| 12 | |
|---|
| 13 | say 'Create gif images as well (1=yes ; 0=no)' |
|---|
| 14 | pull ans |
|---|
| 15 | frame = 1 |
|---|
| 16 | |
|---|
| 17 | 'q file' |
|---|
| 18 | rec=sublin(result,5) |
|---|
| 19 | _endtime=subwrd(rec,12) |
|---|
| 20 | _endlev=subwrd(rec,9) |
|---|
| 21 | |
|---|
| 22 | runscript = 1 |
|---|
| 23 | dis_t = 2 |
|---|
| 24 | |
|---|
| 25 | while(runscript) |
|---|
| 26 | dis_z = 1 |
|---|
| 27 | 'set t ' dis_t |
|---|
| 28 | 'q dims' |
|---|
| 29 | rec=sublin(result,5) |
|---|
| 30 | _analysis=subwrd(rec,6) |
|---|
| 31 | say 'Time is ' _analysis |
|---|
| 32 | |
|---|
| 33 | runlevel = 1 |
|---|
| 34 | while(runlevel) |
|---|
| 35 | 'set z ' dis_z |
|---|
| 36 | 'q dims' |
|---|
| 37 | rec=sublin(result,4) |
|---|
| 38 | newz=subwrd(rec,6) |
|---|
| 39 | say 'Level is ' newz |
|---|
| 40 | |
|---|
| 41 | found_level = 0 |
|---|
| 42 | ** Looking for right level to plot |
|---|
| 43 | if (newz=0.75) |
|---|
| 44 | found_level=1 |
|---|
| 45 | endif |
|---|
| 46 | if (newz=1.5) |
|---|
| 47 | found_level=1 |
|---|
| 48 | endif |
|---|
| 49 | if (newz=4.0) |
|---|
| 50 | found_level=1 |
|---|
| 51 | endif |
|---|
| 52 | if (newz=8.0) |
|---|
| 53 | found_level=1 |
|---|
| 54 | endif |
|---|
| 55 | if (found_level) |
|---|
| 56 | say 'Plotting' |
|---|
| 57 | 'c' |
|---|
| 58 | 'set grads off' |
|---|
| 59 | 'set gxout shaded' |
|---|
| 60 | 'd theta' |
|---|
| 61 | 'set gxout contour' |
|---|
| 62 | 'd p' |
|---|
| 63 | 'set strsiz .2' |
|---|
| 64 | 'set string 1 l 6' |
|---|
| 65 | 'draw string 2.7 8.35 pressure (mb,lines), Theta (color)' |
|---|
| 66 | 'set strsiz .15' |
|---|
| 67 | 'set string 1 l 3' |
|---|
| 68 | 'draw string 6.5 8.0 ' _analysis |
|---|
| 69 | 'draw string 2.5 8.0 ' newz ' km' |
|---|
| 70 | 'run cbar.gs' |
|---|
| 71 | if(ans) |
|---|
| 72 | 'printim qss'frame'.gif gif ' |
|---|
| 73 | frame=frame+1 |
|---|
| 74 | endif |
|---|
| 75 | pull dummy |
|---|
| 76 | 'c' |
|---|
| 77 | 'set grads off' |
|---|
| 78 | 'set gxout shaded' |
|---|
| 79 | 'd w' |
|---|
| 80 | 'set strsiz .2' |
|---|
| 81 | 'set string 1 l 6' |
|---|
| 82 | 'draw string 4.7 8.35 W (color)' |
|---|
| 83 | 'set strsiz .15' |
|---|
| 84 | 'set string 1 l 3' |
|---|
| 85 | 'draw string 6.5 8.0 ' _analysis |
|---|
| 86 | 'draw string 2.5 8.0 ' newz ' km' |
|---|
| 87 | 'run cbar.gs' |
|---|
| 88 | if(ans) |
|---|
| 89 | 'printim qss'frame'.gif gif ' |
|---|
| 90 | frame=frame+1 |
|---|
| 91 | endif |
|---|
| 92 | pull dummy |
|---|
| 93 | 'c' |
|---|
| 94 | 'set grads off' |
|---|
| 95 | 'set gxout shaded' |
|---|
| 96 | 'd qvapor*1000.0' |
|---|
| 97 | 'set strsiz .2' |
|---|
| 98 | 'set string 1 l 6' |
|---|
| 99 | 'draw string 4.2 8.35 QVAPOR (color)' |
|---|
| 100 | 'set strsiz .15' |
|---|
| 101 | 'set string 1 l 3' |
|---|
| 102 | 'draw string 6.5 8.0 ' _analysis |
|---|
| 103 | 'draw string 2.5 8.0 ' newz ' km' |
|---|
| 104 | 'run cbar.gs' |
|---|
| 105 | if(ans) |
|---|
| 106 | 'printim qss'frame'.gif gif ' |
|---|
| 107 | frame=frame+1 |
|---|
| 108 | endif |
|---|
| 109 | pull dummy |
|---|
| 110 | 'c' |
|---|
| 111 | 'set grads off' |
|---|
| 112 | 'set gxout shaded' |
|---|
| 113 | 'd qcloud*1000.0' |
|---|
| 114 | 'set strsiz .2' |
|---|
| 115 | 'set string 1 l 6' |
|---|
| 116 | 'draw string 4.2.7 8.35 QCLOUD (color)' |
|---|
| 117 | 'set strsiz .15' |
|---|
| 118 | 'set string 1 l 3' |
|---|
| 119 | 'draw string 6.5 8.0 ' _analysis |
|---|
| 120 | 'draw string 2.5 8.0 ' newz ' km' |
|---|
| 121 | 'run cbar.gs' |
|---|
| 122 | if(ans) |
|---|
| 123 | 'printim qss'frame'.gif gif ' |
|---|
| 124 | frame=frame+1 |
|---|
| 125 | endif |
|---|
| 126 | pull dummy |
|---|
| 127 | 'c' |
|---|
| 128 | 'set grads off' |
|---|
| 129 | 'set gxout shaded' |
|---|
| 130 | 'd qrain*1000.0' |
|---|
| 131 | 'set strsiz .2' |
|---|
| 132 | 'set string 1 l 6' |
|---|
| 133 | 'draw string 4.2.7 8.35 QRAIN (color)' |
|---|
| 134 | 'set strsiz .15' |
|---|
| 135 | 'set string 1 l 3' |
|---|
| 136 | 'draw string 6.5 8.0 ' _analysis |
|---|
| 137 | 'draw string 2.5 8.0 ' newz ' km' |
|---|
| 138 | 'run cbar.gs' |
|---|
| 139 | pull dummy |
|---|
| 140 | if(ans) |
|---|
| 141 | 'printim qss'frame'.gif gif ' |
|---|
| 142 | frame=frame+1 |
|---|
| 143 | endif |
|---|
| 144 | endif |
|---|
| 145 | |
|---|
| 146 | if ( dis_z=_endlev ) |
|---|
| 147 | runlevel=0 |
|---|
| 148 | endif |
|---|
| 149 | dis_z = dis_z + 1 |
|---|
| 150 | endwhile |
|---|
| 151 | |
|---|
| 152 | |
|---|
| 153 | dis_t = dis_t + 2 |
|---|
| 154 | if ( dis_t>_endtime ) |
|---|
| 155 | runscript=0 |
|---|
| 156 | endif |
|---|
| 157 | |
|---|
| 158 | endwhile |
|---|
| 159 | |
|---|
Note: See
TracBrowser
for help on using the repository browser.
