; basically, this is how to decide on the contour bounds...
; the contour procedure takes as arguments two vectors that 

pro GetContourBounds, ModelChosen, LatMin, LatMax, LongMin, LongMax

LatMin  = 0 & LatMax  = 0 & LongMin = 0 & LongMax = 0

;		print*, "  >  1. obs (Hulme)	 2. ECHam4-OPYC3"
;		print*, "  >  3. HadCM2		 4. HadCM3"
;		print*, "  >  5. CGCM1		 6. CSIRO"
;		print*, "  >  7. CCSR-NIES	 8. ECHam3-LSG"
;		print*, "  >  9. GFDL		10. NCAR-CGCM"
;		print*, "  > 11. NCAR-PCM3	12. 0.5"
;		print*, "  > 13. 10' Europe	14. 0.5 Europe"
;		print*, "  > 15. 10' ATEAM      16. ReGrid x 2"

if            (ModelChosen EQ 3) then begin
  LatMin  =  -90.00 - 1.25
  LatMax  =   90.00 + 1.25
  LongMin = -180.00 - 1.875
  LongMax =  180.00 - 1.875
endif else if (ModelChosen EQ 4) then begin 
  LatMin  =  -90.00 - 1.25
  LatMax  =   90.00 + 1.25
  LongMin = -180.00 - 1.875
  LongMax =  180.00 - 1.875
endif else if (ModelChosen EQ 12) then begin 
  LatMin  =  -90.00 - 0.25
  LatMax  =   90.00 + 0.25
  LongMin = -180.00 - 0.25
  LongMax =  180.00 - 0.25
endif else if (ModelChosen EQ 13) then begin 
  LatMin  =   30.00
  LatMax  =   82.17
  LongMin =  -20.00
  LongMax =   60.00
endif else if (ModelChosen EQ 14) then begin 
  LatMin  =   30.00
  LatMax  =   83.00
  LongMin =  -20.00
  LongMax =   60.00
endif else if (ModelChosen EQ 15) then begin 
  LatMin  =   34.00
  LatMax  =   72.00
  LongMin =  -11.00
  LongMax =   32.00
endif else if (ModelChosen EQ 16) then begin 
  LatMin  =  -89.375
  LatMax  =   89.375
  LongMin = -180.000
  LongMax =  180.000
endif else if (ModelChosen EQ 23) then begin 
  LatMin  =  -90.00 - 1.25
  LatMax  =   90.00 + 1.25
  LongMin = -180.00 - 1.25
  LongMax =  180.00 - 1.25
endif else begin
  print, "  > ##### ERROR: GetContourBounds: no model, gonna crash soon #####"
endelse

end