am 4/06
  annual tides in tiegcm-hist-lbc
  before: we couldn't use annual tides since I didn't know
    how to adjust the coefficients in subroutine bndry_annual
    In subroutine bndry_annual was also the background geopotential height
    set
    a quick fix: the model was run w/o annual tides and in subroutine bndry_annual
    80./96.37*CC(0,0) was used to get 80 km at the -10 lower boundary
   
  now: we use msis for the annual tides
  -change msis to feed in not an altitude but a pressure level
  (s. file msis_background) the output is the neutral temperature
  and the geopotenial height
  -longitudinal averaging of geopotential height
  -use global mean if tiegcm run w/o annual tides
  
  to do:
  - what's about if the global mean from MSIS varies over the year?
    Is this a problem?
  - test longer run, different season
  - check if we have geostrophic balance at the equator (no coriolis force)
    therefore the horizontal north-south gradient of the geopotential height 
    and neutral temperature should be zero at the equator 
    the gradient is not zero, but very small (7.e-5 compared to e.g. 6.e-4 off
    the equator) This seems to be ok.
   - was chrashing after 12 hrs- why? try a default run dev5.
     Ray's code is woorking,
     with my modsrc problem arises at the north pole with high values at the
     lowest calculated levels, but not at the LB
     (s. bluesky: tiegcm-dev5-ibm.job + tiegcm-dev5.inp)