diff -r 0c6405ab2ff4 -r 4be95183fbcd all/02.lai.ncl
--- a/all/02.lai.ncl	Mon Jan 26 22:08:20 2009 -0500
+++ b/all/02.lai.ncl	Thu Mar 26 14:02:21 2009 -0400
@@ -94,8 +94,8 @@
 ;---------------------------------
 ; read lai data: observed
 
-  ;ob_name = "MODIS MOD 15A2 2000-2005"
-  ob_name = "MODIS MOD 15A2 2000-2004"
+  ;ob_name = "MODIS MOD15A2 (2000-2005)"
+  ob_name = "MODIS MOD15A2 (2000-2004)"
 
   dir_l = diro + "lai/"
   ;filo = "LAI_2000-2005_MONS_"+model_grid+".nc"
@@ -145,7 +145,7 @@
 
   delete (wks)
 
-  system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
+  system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
          "rm "+plot_name+"."+plot_type)
  
 ;************************************************
@@ -167,7 +167,7 @@
 
   delete (wks)
 
-  system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
+  system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
          "rm "+plot_name+"."+plot_type)
  
 ;*******************************************************************
@@ -258,6 +258,7 @@
   if (m .eq. 0) then
      data_ob  = dim_avg_Wrap(laiob (lat|:,lon|:,time|:))
      data_mod = dim_avg_Wrap(laimod(lat|:,lon|:,time|:))
+     score_max = 4.0
   end if
 
 ; (B) Max
@@ -283,6 +284,8 @@
         data_mod(j,i) = max(laimod(:,j,i))
      end do
      end do
+
+     score_max = 5.0
           
   end if
 
@@ -309,6 +312,8 @@
         data_mod(j,i) = maxind(laimod(:,j,i)) + 1
      end do
      end do
+
+     score_max = 6.0
         
   end if
 
@@ -373,7 +378,9 @@
 ; compute correlation coef and M score 
 ;=====================================
 
-  score_max = 5.0
+  ; FMH: Max scores are depend on which component is being analyzed, so this
+  ; is now set in each if statement above.
+  ;score_max = 5.0
 
   u       = yvalues(0,:)
   v       = yvalues(1,:)
@@ -489,15 +496,17 @@
   delete (wks)
   delete (plot)
 
-  system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
+  system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
          "rm "+plot_name+"."+plot_type)
+  ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
+  
 
 ;============================
 ; global contour : model
 ;============================
 
   plot_name = "global_"+component(m)+"_model"  
-  title     = "Model " + model_name
+  title     = "Model " + model_name + " (2000-2004)"
   resg@tiMainString  = title
 
   wks = gsn_open_wks (plot_type,plot_name)
@@ -508,8 +517,9 @@
   delete (wks)
   delete (plot)
 
-  system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
+  system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
          "rm "+plot_name+"."+plot_type)
+  ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
 
 ;================================
 ; global contour: model vs ob
@@ -544,7 +554,7 @@
 
 ; plot model
 
-  title     = "Model "+ model_name
+  title     = "Model "+ model_name + " (2000-2004)"
   resg@tiMainString  = title
 
   plot(1) = gsn_csm_contour_map_ce(wks,data_mod,resg) 
@@ -571,7 +581,7 @@
 
   zz = data_mod
   zz = data_mod - data_ob
-  title = "Model_"+model_name+" - Observed"
+  title = "Model_"+model_name+" - MODIS MOD15A2"
   resg@tiMainString    = title
 
   plot(2) = gsn_csm_contour_map_ce(wks,zz,resg) 
@@ -586,8 +596,9 @@
   delete (wks)
   delete (plot)
 
-  system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
+  system("convert -trim -density 150 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
          "rm "+plot_name+"."+plot_type)
+  ;system("convert -trim -density 150 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
 
   delete (data_ob)
   delete (data_mod)
@@ -833,7 +844,7 @@
      delete (wks)
      delete (plot)
 
-     system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
+     system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
             "rm "+plot_name+"."+plot_type)
   end do
 
@@ -857,7 +868,7 @@
      delete (wks)
      delete (plot)
 
-     system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
+     system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
             "rm "+plot_name+"."+plot_type)
 
   delete (plot_data)