For the final component of this semester's project, you will be regenerating the plots you've made so far, but for 2006.

1. Make directories capstone/

   2006/monthly/scripts

   2006/monthly/images

   2006/hourly/scripts

   2006/hourly/images

2. Copy over your .csh, .exec, and .gs files from old directories in to appropriate new ones.

3. Data is in /home/gwest/capstone/2006/data/

   -make 2006monthly.ctl using

   grib2ctl.pl /home/gwest/capstone/2006/data/narrmon-a_221_20060201_0000_000.grb > 2006monthly.ctl

   -edit 2006monthly.ctl like you did your previous monthly control file

   -remember to change the tdef line to from 'feb' to 'jan'

   -make index file by doing gribmap -0- i 2006monthly.ctl

   -make 2006hourly00.ctl using

   grib2ctl.pl /home/gwest/capstone/2006/data/narrmonhr-a_221_20060201_0000_000.grb > 2006hourly00.ctl

   -do the same for 03, 06, 09, 12, 15 18, 21

   -edit ctl files as you did for your previous hourly control files

   -make the index files

4. In the .exec files, comment out the lines that run the averaging .gs scripts. Change the variable names to what they are in the control files (remove the 'avg' or 'hravg').

5. For the vertical velocity section, to average the 3 levels, you'll want to define a new variable, newvvelprs as follows:

   define newvvelprs = (vvelprs(lev=500)+vvelprs(lev=600)+vvelprs(lev=700))/3

6. For relative humidity, you'll have to create 4 new RH variables, then average them:

   -set level to 500

   -run rh500.gs, which will have in it

   'define rhprs500 = 100*((lev*100)/(((0.62197*(pow(10,((0.7859+0.03477*(tmpprs-273.15))/(1+0.00412*(tmpprs-273.15))+2))))/spfhprs)+(0.378*(pow(10,((0.7859+0.03477*(tmpprs-273.15))/(1+0.00412*(tmpprs-273.15))+2))))))'

   return(rhprs500)

   -then do that for 600, 700, and 800mb

   -define a new variable, newrhpres, as follows

   define newrhprs = (rhprs500+rhprs600+rhprs700+rhprs800)/4

7. In hourlyanom.exec, in the first part where you do the monthly stuff, in the places where you ran your averaging scripts, you'll want to define monthly variables, for example:

   set lev 850

   define montmpc = tmpprs-273.15

   ....where 'montmpc' is monthly temperature in celsius. This will allow you to difference the monthly and hourly variables later on, for example:

   d (tmpprs-273.15) -montmpc

   ...where tmpprs is the variable from the hourly control file you have open

8. To rotate winds you'll want to use: /home/ldm/oper/models/eta/wind_rotate.gs and /home/ldm/oper/models/eta/wind10m_rotate.gs which will return the variables (ut,vt) and (ut10m,vt10m), respectively, that you will use to display.

9. Also you'll want to update directory paths where needed in .csh and .exec files

10. Debugging will be fairly straight forward. With the number of repetitive changes you're making, you're bound to miss a few things. When you run it, if at the end it says it didnt make the gif, scroll up to where you find the first error. It'll probably indicate the directory or variable that it can't find.