|
Last change
on this file since 206 was
11,
checked in by aslmd, 14 years ago
|
|
spiga@svn-planeto:ajoute le modele meso-echelle martien
|
-
Property svn:executable set to
*
|
|
File size:
2.0 KB
|
| Line | |
|---|
| 1 | 'reinit' |
|---|
| 2 | 'open real_z.ctl' |
|---|
| 3 | |
|---|
| 4 | ** script needs z level data as input |
|---|
| 5 | ** will create cross sections NS ; EW, in middle of domain |
|---|
| 6 | |
|---|
| 7 | 'set mpdset hires' |
|---|
| 8 | 'set display color white' |
|---|
| 9 | 'run rgbset.gs' |
|---|
| 10 | |
|---|
| 11 | say 'Create gif images as well (1=yes ; 0=no)' |
|---|
| 12 | pull ans |
|---|
| 13 | frame = 1 |
|---|
| 14 | |
|---|
| 15 | 'q file' |
|---|
| 16 | rec=sublin(result,5) |
|---|
| 17 | _endtime=subwrd(rec,12) |
|---|
| 18 | num_x=subwrd(rec,3) |
|---|
| 19 | num_y=subwrd(rec,6) |
|---|
| 20 | num_z=subwrd(rec,9) |
|---|
| 21 | |
|---|
| 22 | dis_x=num_x/2 |
|---|
| 23 | dis_y=num_y/2 |
|---|
| 24 | |
|---|
| 25 | runscript = 1 |
|---|
| 26 | dis_t = 1 |
|---|
| 27 | 'set z 1 ' num_z |
|---|
| 28 | 'set y ' dis_y |
|---|
| 29 | |
|---|
| 30 | ** DO EW cross section |
|---|
| 31 | while(runscript) |
|---|
| 32 | say 'Plotting EW cross section' |
|---|
| 33 | 'set t ' dis_t |
|---|
| 34 | 'q dims' |
|---|
| 35 | rec=sublin(result,5) |
|---|
| 36 | _analysis=subwrd(rec,6) |
|---|
| 37 | say 'Time is ' _analysis |
|---|
| 38 | |
|---|
| 39 | 'c' |
|---|
| 40 | 'set grads off' |
|---|
| 41 | 'set gxout shaded' |
|---|
| 42 | 'set clevs 10 20 30 40 50 60 70 80 90 100' |
|---|
| 43 | 'set ccols 0 0 0 0 33 34 35 36 37 38 39' |
|---|
| 44 | 'd rh' |
|---|
| 45 | 'set gxout contour' |
|---|
| 46 | 'set ccolor 1' |
|---|
| 47 | 'set cthick 2' |
|---|
| 48 | 'set cint 5' |
|---|
| 49 | 'd tc' |
|---|
| 50 | 'set strsiz .2' |
|---|
| 51 | 'set string 1 l 6' |
|---|
| 52 | 'draw string 2.7 8.35 E-W Relative Humidity and T(C)' |
|---|
| 53 | 'set strsiz .15' |
|---|
| 54 | 'set string 1 l 3' |
|---|
| 55 | 'draw string 8.4 8.2 ' _analysis |
|---|
| 56 | 'draw string 1.5 8.2 y=' dis_y |
|---|
| 57 | 'run cbar.gs' |
|---|
| 58 | if(ans) |
|---|
| 59 | 'printim crossEW'frame'.gif gif ' |
|---|
| 60 | frame=frame+1 |
|---|
| 61 | endif |
|---|
| 62 | pull dummy |
|---|
| 63 | |
|---|
| 64 | if ( dis_t=_endtime ) |
|---|
| 65 | runscript=0 |
|---|
| 66 | endif |
|---|
| 67 | dis_t = dis_t + 1 |
|---|
| 68 | endwhile |
|---|
| 69 | |
|---|
| 70 | 'reset' |
|---|
| 71 | |
|---|
| 72 | runscript = 1 |
|---|
| 73 | dis_t = 1 |
|---|
| 74 | 'set z 1 ' num_z |
|---|
| 75 | 'set x ' dis_x |
|---|
| 76 | |
|---|
| 77 | |
|---|
| 78 | ** DO NS cross section |
|---|
| 79 | say 'Plotting NS cross section' |
|---|
| 80 | while(runscript) |
|---|
| 81 | 'set t ' dis_t |
|---|
| 82 | 'q dims' |
|---|
| 83 | rec=sublin(result,5) |
|---|
| 84 | _analysis=subwrd(rec,6) |
|---|
| 85 | say 'Time is ' _analysis |
|---|
| 86 | |
|---|
| 87 | 'c' |
|---|
| 88 | 'set grads off' |
|---|
| 89 | 'set gxout shaded' |
|---|
| 90 | 'set clevs 10 20 30 40 50 60 70 80 90 100' |
|---|
| 91 | 'set ccols 0 0 0 0 33 34 35 36 37 38 39' |
|---|
| 92 | 'd rh' |
|---|
| 93 | 'set gxout contour' |
|---|
| 94 | 'set ccolor 1' |
|---|
| 95 | 'set cthick 2' |
|---|
| 96 | 'set cint 5' |
|---|
| 97 | 'd tc' |
|---|
| 98 | 'set strsiz .2' |
|---|
| 99 | 'set string 1 l 6' |
|---|
| 100 | 'draw string 2.7 8.35 N-S Relative Humidity and T(C)' |
|---|
| 101 | 'set strsiz .15' |
|---|
| 102 | 'set string 1 l 3' |
|---|
| 103 | 'draw string 8.4 8.2 ' _analysis |
|---|
| 104 | 'draw string 1.5 8.2 x=' dis_x |
|---|
| 105 | 'run cbar.gs' |
|---|
| 106 | if(ans) |
|---|
| 107 | 'printim crossNS'frame'.gif gif ' |
|---|
| 108 | frame=frame+1 |
|---|
| 109 | endif |
|---|
| 110 | pull dummy |
|---|
| 111 | |
|---|
| 112 | if ( dis_t=_endtime ) |
|---|
| 113 | runscript=0 |
|---|
| 114 | endif |
|---|
| 115 | dis_t = dis_t + 1 |
|---|
| 116 | endwhile |
|---|
| 117 | |
|---|
Note: See
TracBrowser
for help on using the repository browser.
