|
Last change
on this file since 170 was
85,
checked in by aslmd, 14 years ago
|
|
LMD_MM_MARS et LMD_LES_MARS: ajout des routines IDL pour tracer les sorties --> voir mesoscale/PLOT
|
|
File size:
428 bytes
|
| Line | |
|---|
| 1 | |
|---|
| 2 | ; |
|---|
| 3 | ; field1='U' |
|---|
| 4 | ; winds=['U','V'] |
|---|
| 5 | ; |
|---|
| 6 | |
|---|
| 7 | if ( (field1 eq 'U') $ |
|---|
| 8 | and (winds(0) eq 'U') $ |
|---|
| 9 | and (winds(1) eq 'V') ) then begin |
|---|
| 10 | |
|---|
| 11 | print, 'computing wind velocity' |
|---|
| 12 | |
|---|
| 13 | title='Horizontal wind velocity (!Nm!N.s!U-1!N)' |
|---|
| 14 | |
|---|
| 15 | u=overvector_x |
|---|
| 16 | v=overvector_y |
|---|
| 17 | |
|---|
| 18 | help, u |
|---|
| 19 | help, v |
|---|
| 20 | |
|---|
| 21 | velocity2=u^2+v^2 |
|---|
| 22 | what_I_plot=sqrt(velocity2) |
|---|
| 23 | |
|---|
| 24 | ;overcontour=sqrt(velocity2) |
|---|
| 25 | |
|---|
| 26 | ;; |
|---|
| 27 | ;; comment if you want wind vectors |
|---|
| 28 | ;; |
|---|
| 29 | ;overvector_x=0 |
|---|
| 30 | ;overvector_y=0 |
|---|
| 31 | |
|---|
| 32 | pal=22 |
|---|
| 33 | |
|---|
| 34 | endif |
|---|
Note: See
TracBrowser
for help on using the repository browser.
