1.1 --- a/all/01.npp.ncl Mon Jan 26 22:08:20 2009 -0500
1.2 +++ b/all/01.npp.ncl Thu Mar 26 14:02:21 2009 -0400
1.3 @@ -193,7 +193,8 @@
1.4 ;(3) global data, interpolated into model grid
1.5
1.6 dirg = diro + "npp/"
1.7 - filg = "Npp_"+model_grid+"_mean.nc"
1.8 + ;filg = "Npp_"+model_grid+"_mean.nc"
1.9 + filg = "npp_"+model_grid+"_mean_2000-2004.nc"
1.10 fglobe = addfile (dirg+filg,"r")
1.11
1.12 nppglobe0 = fglobe->NPP
1.13 @@ -225,11 +226,11 @@
1.14 ;
1.15 ; rain81 : mm/year
1.16 ; rainmod : mm/s
1.17 -; npp81 : g C/m^2/year
1.18 -; nppmod81: g C/m^2/s
1.19 -; nppglobe: g C/m^2/year
1.20 +; npp81 : gC/m^2/year
1.21 +; nppmod81: gC/m^2/s
1.22 +; nppglobe: gC/m^2/year
1.23 ;
1.24 -; We want to convert these to "m/year" and "g C/m^2/year".
1.25 +; We want to convert these to "m/year" and "gC/m^2/year".
1.26
1.27 nsec_per_year = 60*60*24*365
1.28
1.29 @@ -243,27 +244,27 @@
1.30
1.31 nppmod = nppmod * nsec_per_year
1.32
1.33 - npp81@units = "gC/m^2/yr"
1.34 - nppmod81@units = "gC/m^2/yr"
1.35 - npp933@units = "gC/m^2/yr"
1.36 - nppmod933@units = "gC/m^2/yr"
1.37 - nppmod@units = "gC/m^2/yr"
1.38 - nppglobe@units = "gC/m^2/yr"
1.39 - rain81@units = "m/yr"
1.40 - rainmod81@units = "m/yr"
1.41 - rain933@units = "m/yr"
1.42 - rainmod933@units = "m/yr"
1.43 + npp81@units = "gC m~S~-2~N~ y~S~-1~N~"
1.44 + nppmod81@units = "gC m~S~-2~N~ y~S~-1~N~"
1.45 + npp933@units = "gC m~S~-2~N~ y~S~-1~N~"
1.46 + nppmod933@units = "gC m~S~-2~N~ y~S~-1~N~"
1.47 + nppmod@units = "gC m~S~-2~N~ y~S~-1~N~"
1.48 + nppglobe@units = "gC m~S~-2~N~ y~S~-1~N~"
1.49 + rain81@units = "m y~S~-1~N~"
1.50 + rainmod81@units = "m y~S~-1~N~"
1.51 + rain933@units = "m y~S~-1~N~"
1.52 + rainmod933@units = "m y~S~-1~N~"
1.53
1.54 - npp81@long_name = "Obs. NPP (gC/m2/year)"
1.55 - npp933@long_name = "Obs. NPP (gC/m2/year)"
1.56 - nppmod81@long_name = "Model NPP (gC/m2/year)"
1.57 - nppmod933@long_name = "Model NPP (gC/m2/year)"
1.58 - nppmod@long_name = "Model NPP (gC/m2/year)"
1.59 - nppglobe@long_name = "NPP (gC/m2/year)"
1.60 - rain81@long_name = "PREC (m/year)"
1.61 - rain933@long_name = "PREC (m/year)"
1.62 - rainmod81@long_name = "PREC (m/year)"
1.63 - rainmod933@long_name = "PREC (m/year)"
1.64 + npp81@long_name = "Obs. NPP (gC m~S~-2~N~ y~S~-1~N~)"
1.65 + npp933@long_name = "Obs. NPP (gC m~S~-2~N~ y~S~-1~N~)"
1.66 + nppmod81@long_name = "Model NPP (gC m~S~-2~N~ y~S~-1~N~)"
1.67 + nppmod933@long_name = "Model NPP (gC m~S~-2~N~ y~S~-1~N~)"
1.68 + nppmod@long_name = "Model NPP (gC m~S~-2~N~ y~S~-1~N~)"
1.69 + nppglobe@long_name = "NPP (gC m~S~-2~N~ y~S~-1~N~)"
1.70 + rain81@long_name = "PREC (m y~S~-1~N~)"
1.71 + rain933@long_name = "PREC (m y~S~-1~N~)"
1.72 + rainmod81@long_name = "PREC (m y~S~-1~N~)"
1.73 + rainmod933@long_name = "PREC (m y~S~-1~N~)"
1.74
1.75 ; change longitude from 0-360 to -180-180
1.76 x81 = where(x81 .gt. 180., x81 -360., x81 )
1.77 @@ -289,8 +290,8 @@
1.78 ," <th bgcolor=DDDDDD >Site ID</th>" \
1.79 ," <th bgcolor=DDDDDD >Latitude</th>" \
1.80 ," <th bgcolor=DDDDDD >Longitude</th>" \
1.81 - ," <th bgcolor=DDDDDD >NPP(gC/m2/year)</th>" \
1.82 - ," <th bgcolor=DDDDDD >PPT(m/year)</th>" \
1.83 + ," <th bgcolor=DDDDDD >NPP (gC m<small><sup>-2</sup></small> y<small><sup>-1</sup></small>)</th>" \
1.84 + ," <th bgcolor=DDDDDD >PPT (m y<small><sup>-1</sup></small>)</th>" \
1.85 ,"</tr>" \
1.86 /)
1.87 table_footer = "</table>"
1.88 @@ -354,8 +355,8 @@
1.89 ," <th bgcolor=DDDDDD >Site ID</th>" \
1.90 ," <th bgcolor=DDDDDD >Latitude</th>" \
1.91 ," <th bgcolor=DDDDDD >Longitude</th>" \
1.92 - ," <th bgcolor=DDDDDD >NPP(gC/m2/year)</th>" \
1.93 - ," <th bgcolor=DDDDDD >PPT(m/year)</th>" \
1.94 + ," <th bgcolor=DDDDDD >NPP (gC m<small><sup>-2</sup></small> y<small><sup>-1</sup></small>)</th>" \
1.95 + ," <th bgcolor=DDDDDD >PPT (m y<small><sup>-1</sup></small>)</th>" \
1.96 ,"</tr>" \
1.97 /)
1.98 table_footer = "</table>"
1.99 @@ -425,8 +426,8 @@
1.100 ," <th bgcolor=DDDDDD rowspan=2>Site ID</th>" \
1.101 ," <th bgcolor=DDDDDD rowspan=2>Latitude</th>" \
1.102 ," <th bgcolor=DDDDDD rowspan=2>Longitude</th>" \
1.103 - ," <th bgcolor=DDDDDD colspan=2>NPP(gC/m2.year)</th>" \
1.104 - ," <th bgcolor=DDDDDD colspan=2>RAIN(m/year)</th>" \
1.105 + ," <th bgcolor=DDDDDD colspan=2>NPP (gC m<small><sup>-2</sup></small> y<small><sup>-1</sup></small>)</th>" \
1.106 + ," <th bgcolor=DDDDDD colspan=2>RAIN (m y<small><sup>-1</sup></small>)</th>" \
1.107 ,"</tr>" \
1.108 ,"<tr>" \
1.109 ," <th bgcolor=DDDDDD >observed</th>" \
1.110 @@ -506,8 +507,8 @@
1.111 ," <th bgcolor=DDDDDD rowspan=2>Site ID</th>" \
1.112 ," <th bgcolor=DDDDDD rowspan=2>Latitude</th>" \
1.113 ," <th bgcolor=DDDDDD rowspan=2>Longitude</th>" \
1.114 - ," <th bgcolor=DDDDDD colspan=2>NPP(gC/m2.year)</th>" \
1.115 - ," <th bgcolor=DDDDDD colspan=2>RAIN(m/year)</th>" \
1.116 + ," <th bgcolor=DDDDDD colspan=2>NPP (gC m<small><sup>-2</sup></small> y<small><sup>-1</sup></small>)</th>" \
1.117 + ," <th bgcolor=DDDDDD colspan=2>RAIN (m y<small><sup>-1</sup></small>)</th>" \
1.118 ,"</tr>" \
1.119 ,"<tr>" \
1.120 ," <th bgcolor=DDDDDD >observed</th>" \
1.121 @@ -569,7 +570,8 @@
1.122 ;***************************************************************************
1.123
1.124 plot_name = "scatter_model_vs_ob_81"
1.125 - title = model_name +" vs Class A observations (81 sites)"
1.126 + ;title = model_name +" vs Class A observations (81 sites)"
1.127 + title = model_name +" (1975-2000) vs Class A observations (81 sites)"
1.128
1.129 wks = gsn_open_wks (plot_type,plot_name) ; open workstation
1.130 res = True ; plot mods desired
1.131 @@ -587,7 +589,7 @@
1.132 ccr81 = esccr(nppmod81,npp81,0)
1.133 ;print (ccr81)
1.134
1.135 - score_max = 2.5
1.136 + score_max = 1.0
1.137
1.138 bias = sum(abs(nppmod81-npp81)/(abs(nppmod81)+abs(npp81)))
1.139 M81s = (1. - (bias/dimsizes(y81)))*score_max
1.140 @@ -618,15 +620,16 @@
1.141 delete (wks)
1.142 delete (dum)
1.143
1.144 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.145 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.146 "rm "+plot_name+"."+plot_type)
1.147 + ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.148
1.149 ;***************************************************************************
1.150 ;(A)-6 scatter plot, model vs ob 933
1.151 ;***************************************************************************
1.152
1.153 plot_name = "scatter_model_vs_ob_933"
1.154 - title = model_name +" vs Class B Observations (933 sites)"
1.155 + title = model_name +" (1975-2000) vs Class B Observations (933 sites)"
1.156
1.157 wks = gsn_open_wks (plot_type,plot_name) ; open workstation
1.158 res = True ; plot mods desired
1.159 @@ -644,7 +647,7 @@
1.160
1.161 ccr933 = esccr(nppmod933,npp933,0)
1.162
1.163 - score_max = 2.5
1.164 + score_max = 1.0
1.165
1.166 bias = sum(abs(nppmod933-npp933)/(abs(nppmod933)+abs(npp933)))
1.167 M933s = (1. - (bias/dimsizes(y933)))*score_max
1.168 @@ -675,8 +678,9 @@
1.169 delete (wks)
1.170 delete (dum)
1.171
1.172 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.173 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.174 "rm "+plot_name+"."+plot_type)
1.175 + ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.176
1.177 ;**************************************************************************
1.178 ;(B) histogram 81
1.179 @@ -711,7 +715,7 @@
1.180
1.181 ccr81h = esccr(uu,vv,0)
1.182
1.183 - score_max = 2.5
1.184 + score_max = 2.0
1.185
1.186 bias = sum(abs(vv-uu)/(abs(vv)+abs(uu)))
1.187 M81h = (1.- (bias/dimsizes(uu)))*score_max
1.188 @@ -742,8 +746,8 @@
1.189 resh@trYMinF = min(mn_yvalues) - 10.
1.190 resh@trYMaxF = max(mx_yvalues) + 10.
1.191
1.192 - resh@tiYAxisString = "NPP (g C/m2/year)"
1.193 - resh@tiXAxisString = "Precipitation (m/year)"
1.194 + resh@tiYAxisString = "NPP (gC m~S~-2~N~ y~S~-1~N~)"
1.195 + resh@tiXAxisString = "Precipitation (m y~S~-1~N~)"
1.196
1.197 max_bar = new((/2,nx/),graphic)
1.198 min_bar = new((/2,nx/),graphic)
1.199 @@ -765,6 +769,9 @@
1.200
1.201 xy = gsn_csm_xy(wks,xvalues(0,:),yvalues(0,:),resh)
1.202
1.203 + ; Added by Forrest Hoffman to print out values on Wed Feb 4 14:36:00 EST 2009
1.204 + print("Observations: " + xvalues(0,:) + " " + yvalues(0,:))
1.205 +
1.206 ;-------------------------------
1.207 ;Attach the vertical bar and the horizontal cap line
1.208
1.209 @@ -809,15 +816,16 @@
1.210 draw(xy)
1.211 delete (wks)
1.212
1.213 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.214 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.215 "rm "+plot_name+"."+plot_type)
1.216 + ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.217
1.218 ;****************************************************************************
1.219 ;(B)-2 histogram plot, model vs ob 81 site
1.220 ;****************************************************************************
1.221
1.222 plot_name = "histogram_model_vs_ob_81"
1.223 - title = model_name+ " vs Class A Observations (81 sites)"
1.224 + title = model_name+ " (1975-2000) vs Class A Observations (81 sites)"
1.225 resh@tiMainString = title
1.226
1.227 wks = gsn_open_wks (plot_type,plot_name) ; open workstation
1.228 @@ -847,6 +855,10 @@
1.229 gsn_text_ndc(wks,correlation_text,0.56,0.85,tRes)
1.230
1.231 xy = gsn_csm_xy(wks,xvalues,yvalues,resh)
1.232 +
1.233 + ; Added by Forrest Hoffman to print out values on Wed Feb 4 14:36:00 EST 2009
1.234 + print(model_name + ": " + xvalues(1,:) + " " + yvalues(1,:))
1.235 + ;print("All: " + xvalues(:,:) + " " + yvalues(:,:))
1.236 ;-------------------------------
1.237 ;Attach the vertical bar and the horizontal cap line
1.238
1.239 @@ -883,8 +895,9 @@
1.240 draw(xy)
1.241 delete(wks)
1.242
1.243 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.244 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.245 "rm "+plot_name+"."+plot_type)
1.246 + ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.247
1.248 delete (RAIN1_1D)
1.249 delete (RAIN2_1D)
1.250 @@ -935,7 +948,7 @@
1.251
1.252 ccr933h = esccr(uu,vv,0)
1.253
1.254 - score_max = 2.5
1.255 + score_max = 2.0
1.256
1.257 bias = sum(abs(vv-uu)/(abs(vv)+abs(uu)))
1.258 M933h = (1.- (bias/dimsizes(uu)))*score_max
1.259 @@ -967,8 +980,8 @@
1.260 resh@trYMinF = min(mn_yvalues) - 10.
1.261 resh@trYMaxF = max(mx_yvalues) + 10.
1.262
1.263 - resh@tiYAxisString = "NPP (g C/m2/year)"
1.264 - resh@tiXAxisString = "Precipitation (m/year)"
1.265 + resh@tiYAxisString = "NPP (gC m~S~-2~N~ y~S~-1~N~)"
1.266 + resh@tiXAxisString = "Precipitation (m y~S~-1~N~)"
1.267
1.268 max_bar = new((/2,nx/),graphic)
1.269 min_bar = new((/2,nx/),graphic)
1.270 @@ -990,6 +1003,9 @@
1.271
1.272 xy = gsn_csm_xy(wks,xvalues(0,:),yvalues(0,:),resh)
1.273
1.274 + ; Added by Forrest Hoffman to print out values on Wed Feb 4 14:36:00 EST 2009
1.275 + print("Observations: " + xvalues(0,:) + " " + yvalues(0,:))
1.276 +
1.277 ;-------------------------------
1.278 ;Attach the vertical bar and the horizontal cap line
1.279
1.280 @@ -1027,15 +1043,16 @@
1.281 delete (xy)
1.282 delete (wks)
1.283
1.284 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.285 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.286 "rm "+plot_name+"."+plot_type)
1.287 + ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.288
1.289 ;****************************************************************************
1.290 ;(B)-4 histogram plot, model vs ob 933 site
1.291 ;****************************************************************************
1.292
1.293 plot_name = "histogram_model_vs_ob_933"
1.294 - title = model_name+ " vs Class B Observations (933 sites)"
1.295 + title = model_name+ " (1975-2000) vs Class B Observations (933 sites)"
1.296 resh@tiMainString = title
1.297
1.298 wks = gsn_open_wks (plot_type,plot_name) ; open workstation
1.299 @@ -1065,6 +1082,9 @@
1.300 gsn_text_ndc(wks,correlation_text,0.56,0.85,tRes)
1.301
1.302 xy = gsn_csm_xy(wks,xvalues,yvalues,resh)
1.303 +
1.304 + ; Added by Forrest Hoffman to print out values on Wed Feb 4 14:36:00 EST 2009
1.305 + print("Observations: " + xvalues(1,:) + " " + yvalues(1,:))
1.306 ;-------------------------------
1.307 ;Attach the vertical bar and the horizontal cap line
1.308
1.309 @@ -1102,8 +1122,30 @@
1.310 delete(xy)
1.311 delete(wks)
1.312
1.313 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.314 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.315 "rm "+plot_name+"."+plot_type)
1.316 + ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.317 +
1.318 +;***************************************************************************
1.319 +; Read 2000-2004 dataset for MODIS comparison
1.320 +;***************************************************************************
1.321 +
1.322 +;------------------------------------------------------
1.323 +; read model data
1.324 +
1.325 + fm = addfile (dirm+film9,"r")
1.326 +
1.327 + nppmod0 = fm->NPP
1.328 + xm = fm->lon
1.329 + ym = fm->lat
1.330 +
1.331 + delete (fm)
1.332 +
1.333 + nppmod = nppmod0(0,:,:)
1.334 + delete (nppmod0)
1.335 + nppmod = nppmod * nsec_per_year
1.336 + nppmod@units = "gC m~S~-2~N~ y~S~-1~N~"
1.337 + nppmod@long_name = "Model NPP (gC m~S~-2~N~ y~S~-1~N~)"
1.338
1.339 ;***************************************************************************
1.340 ;(C) global contour
1.341 @@ -1132,7 +1174,7 @@
1.342 nppglobe = where(ismissing(nppglobe).and.(ismissing(nppmod).or.(nppmod.lt.delta)),0.,nppglobe)
1.343
1.344 plot_name = "global_ob"
1.345 - title = "Observed MODIS MOD 17"
1.346 + title = "MODIS MOD17A3 (2000-2004)"
1.347 resg@tiMainString = title
1.348
1.349 wks = gsn_open_wks (plot_type,plot_name) ; open workstation
1.350 @@ -1142,14 +1184,15 @@
1.351
1.352 delete (wks)
1.353
1.354 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.355 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.356 "rm "+plot_name+"."+plot_type)
1.357 + ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.358 ;****************************************************************************
1.359 ;(C)-2 global contour plot, model
1.360 ;****************************************************************************
1.361
1.362 plot_name = "global_model"
1.363 - title = "Model "+ model_name
1.364 + title = "Model "+ model_name + " (2000-2004)"
1.365 resg@tiMainString = title
1.366
1.367 wks = gsn_open_wks (plot_type,plot_name) ; open workstation
1.368 @@ -1159,8 +1202,9 @@
1.369
1.370 delete (wks)
1.371
1.372 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.373 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.374 "rm "+plot_name+"."+plot_type)
1.375 + ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.376
1.377 ;****************************************************************************
1.378 ;(C)-3 global contour plot, model vs ob
1.379 @@ -1190,7 +1234,7 @@
1.380
1.381 ccrG = esccr(ug,vg,0)
1.382
1.383 - score_max = 5.0
1.384 + score_max = 2.0
1.385
1.386 MG = (ccrG*ccrG)* score_max
1.387 M_npp_G = sprintf("%.2f", MG)
1.388 @@ -1215,14 +1259,14 @@
1.389 ;--------------------------------------------------------------------
1.390 ;(a) ob
1.391
1.392 - title = "Observed MODIS MOD 17"
1.393 + title = "MODIS MOD17A3 (2000-2004)"
1.394 resg@tiMainString = title
1.395
1.396 plot(0) = gsn_csm_contour_map_ce(wks,nppglobe,resg)
1.397
1.398 ;(b) model
1.399
1.400 - title = "Model "+ model_name
1.401 + title = "Model "+ model_name + " (2000-2004)"
1.402 resg@tiMainString = title
1.403
1.404 plot(1) = gsn_csm_contour_map_ce(wks,nppmod,resg)
1.405 @@ -1231,7 +1275,7 @@
1.406
1.407 zz = nppmod
1.408 zz = nppmod - nppglobe
1.409 - title = "Model_"+model_name+" - Observed"
1.410 + title = "Model "+model_name+" - MODIS MOD17A3"
1.411
1.412 resg@cnMinLevelValF = -500 ; Min level
1.413 resg@cnMaxLevelValF = 500. ; Max level
1.414 @@ -1250,8 +1294,9 @@
1.415 delete (wks)
1.416 delete (plot)
1.417
1.418 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.419 + system("convert -trim -density 150 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.420 "rm "+plot_name+"."+plot_type)
1.421 + ;system("convert -trim -density 150 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.422
1.423 ;***************************************************************************
1.424 ;(D)-1 zonal line plot, ob
1.425 @@ -1261,7 +1306,7 @@
1.426 vv@long_name = nppglobe@long_name
1.427
1.428 plot_name = "zonal_ob"
1.429 - title = "Observed MODIS MOD 17"
1.430 + title = "MODIS MOD17A3 (2000-2004)"
1.431
1.432 resz = True
1.433 resz@tiMainString = title
1.434 @@ -1272,8 +1317,9 @@
1.435
1.436 delete (wks)
1.437
1.438 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.439 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.440 "rm "+plot_name+"."+plot_type)
1.441 + ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.442
1.443 ;****************************************************************************
1.444 ;(D)-2 zonal line plot, model vs ob
1.445 @@ -1288,7 +1334,7 @@
1.446 ;-------------------------------------------
1.447 ; compute correlation coef and M score
1.448
1.449 - score_max = 5.0
1.450 + score_max = 2.0
1.451
1.452 ccrZ = esccr(s(1,:), s(0,:),0)
1.453 MZ = (ccrZ*ccrZ)* score_max
1.454 @@ -1303,7 +1349,7 @@
1.455 "mv -f "+html_new+" "+html_name)
1.456 ;-------------------------------------------
1.457 plot_name = "zonal_model_vs_ob"
1.458 - title = "Zonal Average"
1.459 + title = "Zonal Average (2000-2004)"
1.460 resz@tiMainString = title
1.461
1.462 wks = gsn_open_wks (plot_type,plot_name)
1.463 @@ -1324,16 +1370,16 @@
1.464 resz@pmLegendDisplayMode = "Always" ; turn on legend
1.465 resz@pmLegendSide = "Top" ; Change location of
1.466 ; resz@pmLegendParallelPosF = .45 ; move units right
1.467 - resz@pmLegendParallelPosF = .82 ; move units right
1.468 + resz@pmLegendParallelPosF = .75 ; move units right
1.469 resz@pmLegendOrthogonalPosF = -0.4 ; move units down
1.470
1.471 resz@pmLegendWidthF = 0.10 ; Change width and
1.472 resz@pmLegendHeightF = 0.10 ; height of legend.
1.473 - resz@lgLabelFontHeightF = .02 ; change font height
1.474 + resz@lgLabelFontHeightF = .015 ; change font height
1.475 ; resz@lgTitleOn = True ; turn on legend title
1.476 ; resz@lgTitleString = "Example" ; create legend title
1.477 ; resz@lgTitleFontHeightF = .025 ; font of legend title
1.478 - resz@xyExplicitLegendLabels = (/"Observed",model_name/) ; explicit labels
1.479 + resz@xyExplicitLegendLabels = (/"MODIS MOD17A3",model_name/) ; explicit labels
1.480 ;--------------------------------------------------------------------
1.481 zRes = True
1.482 zRes@txFontHeightF = 0.025
1.483 @@ -1347,8 +1393,9 @@
1.484
1.485 delete (wks)
1.486
1.487 - system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.488 + system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
1.489 "rm "+plot_name+"."+plot_type)
1.490 + ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
1.491
1.492 ;***************************************************************************
1.493 ; add total score and write to file