| 1 | ;;--------------------------------------- |
|---|
| 2 | ; TRACE |
|---|
| 3 | ; |
|---|
| 4 | ; **** DO NOT MODIFY **** |
|---|
| 5 | ; |
|---|
| 6 | ; set parameters in plot.def then type: |
|---|
| 7 | ; |
|---|
| 8 | ; @trace.idl |
|---|
| 9 | ; |
|---|
| 10 | ;;--------------------------------------- |
|---|
| 11 | |
|---|
| 12 | field1='' & field2='' & winds='' |
|---|
| 13 | backup_data='no' & already_data='no' |
|---|
| 14 | datafolder='./' & plotfolder='./' |
|---|
| 15 | topo=0 |
|---|
| 16 | extract='no' |
|---|
| 17 | |
|---|
| 18 | SPAWN, '\rm user.idl' |
|---|
| 19 | SPAWN, 'grep -v grep plot.def > user.idl' |
|---|
| 20 | |
|---|
| 21 | ;; to avoid dumb and useless warnings from IDL |
|---|
| 22 | space=0. & altitude=0. & lon=0. & lat=0. |
|---|
| 23 | what_I_plot=0. & title_user='' & missing_value=0. |
|---|
| 24 | |
|---|
| 25 | @user.idl |
|---|
| 26 | if (nam5 eq " interp_method = 0 ") then coord='model_level' |
|---|
| 27 | if ((nam5 eq " interp_method = 1 ") and (tabnam(0) lt tabnam(1))) then coord='height' |
|---|
| 28 | |
|---|
| 29 | .compile call_arwpost |
|---|
| 30 | if (extract eq 'yes') then call_arwpost, nam1, nam2, nam3, nam4, nam5, tabnam |
|---|
| 31 | |
|---|
| 32 | denom=plot+'_'+coord+'_'+field1 |
|---|
| 33 | if ((field2 ne '') and (topo eq 0)) then denom=denom+'_'+field2 |
|---|
| 34 | if (topo eq 1) then denom=denom+'_HGT' |
|---|
| 35 | if (winds(0) ne '') then denom=denom+'_'+winds(0)+winds(1) |
|---|
| 36 | |
|---|
| 37 | if (already_data eq 'yes') then SPAWN, 'ln -sf '+datafolder+'/'+coord+'.input.ctl input.ctl' |
|---|
| 38 | if (already_data eq 'yes') then SPAWN, 'ln -sf '+datafolder+'/'+coord+'.input.dat input.dat' |
|---|
| 39 | |
|---|
| 40 | SPAWN, '\rm param_plot.idl' |
|---|
| 41 | SPAWN, 'grep grepALL plot.def > param_plot.idl' |
|---|
| 42 | if ((plot eq 'meridional') or (plot eq 'zonal')) then SPAWN, 'grep grepSEC plot.def >> param_plot.idl' |
|---|
| 43 | if (plot eq 'map') then SPAWN, 'grep grepMAP plot.def >> param_plot.idl' |
|---|
| 44 | .compile map_latlon |
|---|
| 45 | .compile section |
|---|
| 46 | |
|---|
| 47 | !p.charthick = 2.0 |
|---|
| 48 | !p.thick = 3.0 |
|---|
| 49 | !x.thick = 2.0 |
|---|
| 50 | !y.thick = 2.0 |
|---|
| 51 | |
|---|
| 52 | out_wrf, $ |
|---|
| 53 | plot=plot, $ |
|---|
| 54 | field1=field1, $ |
|---|
| 55 | field2=field2, $ |
|---|
| 56 | when=1, $ |
|---|
| 57 | winds=winds, $ |
|---|
| 58 | level=level,$ |
|---|
| 59 | save_data=yeah,$ ;; data is loaded only once |
|---|
| 60 | topo=topo,$ |
|---|
| 61 | ; range=range,$ |
|---|
| 62 | colors=colors, $ |
|---|
| 63 | nlon=nlon, $ |
|---|
| 64 | nlat=nlat, $ |
|---|
| 65 | save_ps=denom+string(1001,'(I0)') |
|---|
| 66 | ; |
|---|
| 67 | ; display the first one as a sample |
|---|
| 68 | ; |
|---|
| 69 | SPAWN, 'display ./'+denom+'1001.ps &' |
|---|
| 70 | |
|---|
| 71 | if (num ge 2) then for i=2,num do $ |
|---|
| 72 | out_wrf, $ |
|---|
| 73 | plot=plot, $ |
|---|
| 74 | field1=field1, $ |
|---|
| 75 | field2=field2, $ |
|---|
| 76 | when=i, $ |
|---|
| 77 | winds=winds, $ |
|---|
| 78 | level=level,$ |
|---|
| 79 | save_data=yeah,$ ;; data is loaded only once |
|---|
| 80 | topo=topo,$ |
|---|
| 81 | range=range,$ |
|---|
| 82 | colors=colors, $ |
|---|
| 83 | nlon=nlon, $ |
|---|
| 84 | nlat=nlat, $ |
|---|
| 85 | save_ps=denom+string(1000+i,'(I0)') |
|---|
| 86 | |
|---|
| 87 | SPAWN, 'mv -f ./'+denom+'????.ps '+plotfolder+'/' |
|---|
| 88 | if (backup_data eq 'yes') then SPAWN, 'cp -f namelist.ARWpost '+datafolder+'/'+coord+'.namelist.ARWpost' |
|---|
| 89 | if (backup_data eq 'yes') then SPAWN, 'cp -f input.ctl '+datafolder+'/'+coord+'.input.ctl' |
|---|
| 90 | if (backup_data eq 'yes') then SPAWN, 'cp -f input.dat '+datafolder+'/'+coord+'.input.dat' |
|---|
| 91 | |
|---|
| 92 | SPAWN, 'cp -f plot.def '+plotfolder+'/'+denom+'.def' |
|---|
| 93 | |
|---|
