; basically, this gives you the contour vectors...
; the contour procedure takes as arguments two vectors that we supply here

pro GetContourVec, ModelChosen, Lons, Lats

;		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   2) then begin	
  Corners = [-181.40625,-90.0,178.59375,90.0]
  LonN = 96 & LatN = 73
endif else if (ModelChosen EQ   3) then begin	; LonN = 96 LatN = 73
  Corners = [-181.125,-91.25,178.125,91.25]
  LonN = 96 & LatN = 73
endif else if (ModelChosen EQ   4) then begin 
  Corners = [-181.125,-91.25,178.125,91.25]
  LonN = 96 & LatN = 73
endif else if (ModelChosen EQ   5) then begin 
  Corners = [-181.125,-90.0,178.125,90.0]
  LonN = 96 & LatN = 48
endif else if (ModelChosen EQ   6) then begin 
  Corners = [-182.8125,-90.0,177.1875,90.0]
  LonN = 64 & LatN = 56
endif else if (ModelChosen EQ  11) then begin 
  Corners = [-181.40625,-90.0,178.59375,90.0]
  LonN = 128 & LatN = 64
endif else if (ModelChosen EQ  12) then begin 
  Corners = [-180.25,-90.25,179.75,90.25]
  LonN = 720 & LatN = 360
endif else if (ModelChosen EQ  13) then begin 
  Corners = [-20.0,30.0,60.0,82.16666]
  LonN = 480 & LatN = 313
endif else if (ModelChosen EQ  14) then begin 
  Corners = [-20.0,30.0,60.0,83.0]
  LonN = 160 & LatN = 106
endif else if (ModelChosen EQ  15) then begin 
  Corners = [-11.0,34.0,32.0,72.0]
  LonN = 258 & LatN = 228
endif else if (ModelChosen EQ  16) then begin 
  Corners = [-181.125,-88.75,178.125,88.75]
  LonN = 96 & LatN = 71
endif else if (ModelChosen EQ  17) then begin 
  Corners = [-11.0,2.33333,49.0,61.0]
  LonN = 80 & LatN = 72
endif else if (ModelChosen EQ  18) then begin 
  Corners = [-180.625,-91.25,179.375,91.25]
  LonN = 144 & LatN = 73
endif else if (ModelChosen EQ  20) then begin 
  Corners = [-10.019,47.819,3.5551,61.018]
  LonN = 180 & LatN = 290
endif else if (ModelChosen EQ  21) then begin 
  Corners = [-6.927,49.732,2.881,60.747]
  LonN = 130 & LatN = 242
endif else if (ModelChosen EQ  22) then begin 
  Corners = [-182.5,-92.5,177.5,92.5]
  LonN = 72 & LatN = 36
endif else if (ModelChosen EQ  23) then begin 
  Corners = [-181.125,-91.175,178.875,91.175]
  LonN = 144 & LatN =  72
endif else if (ModelChosen EQ 102) then begin
  Corners = [-181.40625,-90.0,181.40625,90.0]
  LonN = 129 & LatN = 64
endif else if (ModelChosen EQ 103) then begin
  Corners = [-181.125,-91.25,181.125,91.25]
  LonN = 97 & LatN = 73
endif else if (ModelChosen EQ 105) then begin 
  Corners = [-181.125,-90.0,181.125,90.0]
  LonN = 97 & LatN = 48
endif else if (ModelChosen EQ 106) then begin 
  Corners = [-182.8125,-90.0,182.1875,90.0]
  LonN = 65 & LatN = 56
endif else if (ModelChosen EQ 111) then begin 
  Corners = [-181.40625,-90.0,181.40625,90.0]
  LonN = 129 & LatN = 65
endif else if (ModelChosen EQ 112) then begin 
  Corners = [-180.25,-90.25,180.25,90.25]
  LonN = 721 & LatN = 360
endif else if (ModelChosen EQ 115) then begin 
  Corners = [-11.0,34.0,32.0,72.0]
  LonN = 258 & LatN = 228
endif else if (ModelChosen EQ 117) then begin 
  Corners = [-11.0,2.33333,49.0,61.0]
  LonN = 80 & LatN = 72
endif else if (ModelChosen EQ 118) then begin 
  Corners = [-180.625,-91.25,180.625,91.25]
  LonN = 145 & LatN = 73
endif else if (ModelChosen EQ 120) then begin 
  Corners = [-10.019,47.819,3.5551,61.018]
  LonN = 180 & LatN = 290
endif else if (ModelChosen EQ 121) then begin 
  Corners = [-6.927,49.732,2.881,60.747]
  LonN = 130 & LatN = 242
endif else if (ModelChosen EQ 122) then begin 
  Corners = [-182.5,-92.5,182.5,92.5]
  LonN = 73 & LatN = 36
endif else if (ModelChosen EQ 123) then begin 
  Corners = [-181.125,-91.175,181.125,91.175]
  LonN = 145 & LatN =  72
endif else begin
  print, "  > ##### ERROR: GetContourVec: no model, this will crash soon #####"
endelse

LonUnit = (Corners [2] - Corners [0]) / float(LonN)
LatUnit = (Corners [3] - Corners [1]) / float(LatN)

Lons = findgen (LonN) & Lons = Lons + 0.5
Lats = findgen (LatN) & Lats = Lats + 0.5

Lons = (Lons * LonUnit) + Corners [0] 
Lats = (Lats * LatUnit) + Corners [1] 

end