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