1.1 --- a/all/02.lai.ncl Mon Jan 26 22:08:20 2009 -0500
1.2 +++ b/all/02.lai.ncl Thu Mar 26 14:02:21 2009 -0400
1.3 @@ -94,8 +94,8 @@
1.4 ;---------------------------------
1.5 ; read lai data: observed
1.6
1.7 - ;ob_name = "MODIS MOD 15A2 2000-2005"
1.8 - ob_name = "MODIS MOD 15A2 2000-2004"
1.9 + ;ob_name = "MODIS MOD15A2 (2000-2005)"
1.10 + ob_name = "MODIS MOD15A2 (2000-2004)"
1.11
1.12 dir_l = diro + "lai/"
1.13 ;filo = "LAI_2000-2005_MONS_"+model_grid+".nc"
1.14 @@ -145,7 +145,7 @@
1.15
1.16 delete (wks)
1.17
1.18 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.19 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.20 "rm "+plot_name+"."+plot_type)
1.21
1.22 ;************************************************
1.23 @@ -167,7 +167,7 @@
1.24
1.25 delete (wks)
1.26
1.27 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.28 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.29 "rm "+plot_name+"."+plot_type)
1.30
1.31 ;*******************************************************************
1.32 @@ -258,6 +258,7 @@
1.33 if (m .eq. 0) then
1.34 data_ob = dim_avg_Wrap(laiob (lat|:,lon|:,time|:))
1.35 data_mod = dim_avg_Wrap(laimod(lat|:,lon|:,time|:))
1.36 + score_max = 4.0
1.37 end if
1.38
1.39 ; (B) Max
1.40 @@ -283,6 +284,8 @@
1.41 data_mod(j,i) = max(laimod(:,j,i))
1.42 end do
1.43 end do
1.44 +
1.45 + score_max = 5.0
1.46
1.47 end if
1.48
1.49 @@ -309,6 +312,8 @@
1.50 data_mod(j,i) = maxind(laimod(:,j,i)) + 1
1.51 end do
1.52 end do
1.53 +
1.54 + score_max = 6.0
1.55
1.56 end if
1.57
1.58 @@ -373,7 +378,9 @@
1.59 ; compute correlation coef and M score
1.60 ;=====================================
1.61
1.62 - score_max = 5.0
1.63 + ; FMH: Max scores are depend on which component is being analyzed, so this
1.64 + ; is now set in each if statement above.
1.65 + ;score_max = 5.0
1.66
1.67 u = yvalues(0,:)
1.68 v = yvalues(1,:)
1.69 @@ -489,15 +496,17 @@
1.70 delete (wks)
1.71 delete (plot)
1.72
1.73 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.74 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.75 "rm "+plot_name+"."+plot_type)
1.76 + ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.77 +
1.78
1.79 ;============================
1.80 ; global contour : model
1.81 ;============================
1.82
1.83 plot_name = "global_"+component(m)+"_model"
1.84 - title = "Model " + model_name
1.85 + title = "Model " + model_name + " (2000-2004)"
1.86 resg@tiMainString = title
1.87
1.88 wks = gsn_open_wks (plot_type,plot_name)
1.89 @@ -508,8 +517,9 @@
1.90 delete (wks)
1.91 delete (plot)
1.92
1.93 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.94 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.95 "rm "+plot_name+"."+plot_type)
1.96 + ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.97
1.98 ;================================
1.99 ; global contour: model vs ob
1.100 @@ -544,7 +554,7 @@
1.101
1.102 ; plot model
1.103
1.104 - title = "Model "+ model_name
1.105 + title = "Model "+ model_name + " (2000-2004)"
1.106 resg@tiMainString = title
1.107
1.108 plot(1) = gsn_csm_contour_map_ce(wks,data_mod,resg)
1.109 @@ -571,7 +581,7 @@
1.110
1.111 zz = data_mod
1.112 zz = data_mod - data_ob
1.113 - title = "Model_"+model_name+" - Observed"
1.114 + title = "Model_"+model_name+" - MODIS MOD15A2"
1.115 resg@tiMainString = title
1.116
1.117 plot(2) = gsn_csm_contour_map_ce(wks,zz,resg)
1.118 @@ -586,8 +596,9 @@
1.119 delete (wks)
1.120 delete (plot)
1.121
1.122 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.123 + system("convert -trim -density 150 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.124 "rm "+plot_name+"."+plot_type)
1.125 + ;system("convert -trim -density 150 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.126
1.127 delete (data_ob)
1.128 delete (data_mod)
1.129 @@ -833,7 +844,7 @@
1.130 delete (wks)
1.131 delete (plot)
1.132
1.133 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.134 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.135 "rm "+plot_name+"."+plot_type)
1.136 end do
1.137
1.138 @@ -857,7 +868,7 @@
1.139 delete (wks)
1.140 delete (plot)
1.141
1.142 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.143 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.144 "rm "+plot_name+"."+plot_type)
1.145
1.146 delete (plot_data)