function merra2grb2(arg) * * simple script to convert a specific merra file to grib2 * * assume the file is already open * 'set x 1 540' 'set y 1 361' 'set t 1 last' '!rm test test.inv' 'g2write ps test PRES:surface:anl' 'g2write slp test PRES:mean sea level:anl' * * levels='1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 10 7 5 3 2 1' i=1 while (i <= 21) lev=subwrd(levels,i) 'set lev ' lev say lev 'g2write h test HGT:' lev ' mb:anl' 'g2write t test TMP:' lev ' mb:anl' 'g2write u test UGRD:' lev ' mb:anl' 'g2write v test VGRD:' lev ' mb:anl' 'g2write qv test SPFH:' lev ' mb:anl' 'g2write o3 test O3MR:' lev ' mb:anl' i=i+1 endwhile 'disable fwrite' '! wgrib2 -i