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

     2.1 --- a/all/02.lai.ncl	Mon Jan 26 22:08:20 2009 -0500
     2.2 +++ b/all/02.lai.ncl	Thu Mar 26 14:02:21 2009 -0400
     2.3 @@ -94,8 +94,8 @@
     2.4  ;---------------------------------
     2.5  ; read lai data: observed
     2.6  
     2.7 -  ;ob_name = "MODIS MOD 15A2 2000-2005"
     2.8 -  ob_name = "MODIS MOD 15A2 2000-2004"
     2.9 +  ;ob_name = "MODIS MOD15A2 (2000-2005)"
    2.10 +  ob_name = "MODIS MOD15A2 (2000-2004)"
    2.11  
    2.12    dir_l = diro + "lai/"
    2.13    ;filo = "LAI_2000-2005_MONS_"+model_grid+".nc"
    2.14 @@ -145,7 +145,7 @@
    2.15  
    2.16    delete (wks)
    2.17  
    2.18 -  system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
    2.19 +  system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
    2.20           "rm "+plot_name+"."+plot_type)
    2.21   
    2.22  ;************************************************
    2.23 @@ -167,7 +167,7 @@
    2.24  
    2.25    delete (wks)
    2.26  
    2.27 -  system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
    2.28 +  system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
    2.29           "rm "+plot_name+"."+plot_type)
    2.30   
    2.31  ;*******************************************************************
    2.32 @@ -258,6 +258,7 @@
    2.33    if (m .eq. 0) then
    2.34       data_ob  = dim_avg_Wrap(laiob (lat|:,lon|:,time|:))
    2.35       data_mod = dim_avg_Wrap(laimod(lat|:,lon|:,time|:))
    2.36 +     score_max = 4.0
    2.37    end if
    2.38  
    2.39  ; (B) Max
    2.40 @@ -283,6 +284,8 @@
    2.41          data_mod(j,i) = max(laimod(:,j,i))
    2.42       end do
    2.43       end do
    2.44 +
    2.45 +     score_max = 5.0
    2.46            
    2.47    end if
    2.48  
    2.49 @@ -309,6 +312,8 @@
    2.50          data_mod(j,i) = maxind(laimod(:,j,i)) + 1
    2.51       end do
    2.52       end do
    2.53 +
    2.54 +     score_max = 6.0
    2.55          
    2.56    end if
    2.57  
    2.58 @@ -373,7 +378,9 @@
    2.59  ; compute correlation coef and M score 
    2.60  ;=====================================
    2.61  
    2.62 -  score_max = 5.0
    2.63 +  ; FMH: Max scores are depend on which component is being analyzed, so this
    2.64 +  ; is now set in each if statement above.
    2.65 +  ;score_max = 5.0
    2.66  
    2.67    u       = yvalues(0,:)
    2.68    v       = yvalues(1,:)
    2.69 @@ -489,15 +496,17 @@
    2.70    delete (wks)
    2.71    delete (plot)
    2.72  
    2.73 -  system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
    2.74 +  system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
    2.75           "rm "+plot_name+"."+plot_type)
    2.76 +  ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
    2.77 +  
    2.78  
    2.79  ;============================
    2.80  ; global contour : model
    2.81  ;============================
    2.82  
    2.83    plot_name = "global_"+component(m)+"_model"  
    2.84 -  title     = "Model " + model_name
    2.85 +  title     = "Model " + model_name + " (2000-2004)"
    2.86    resg@tiMainString  = title
    2.87  
    2.88    wks = gsn_open_wks (plot_type,plot_name)
    2.89 @@ -508,8 +517,9 @@
    2.90    delete (wks)
    2.91    delete (plot)
    2.92  
    2.93 -  system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
    2.94 +  system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
    2.95           "rm "+plot_name+"."+plot_type)
    2.96 +  ;system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
    2.97  
    2.98  ;================================
    2.99  ; global contour: model vs ob
   2.100 @@ -544,7 +554,7 @@
   2.101  
   2.102  ; plot model
   2.103  
   2.104 -  title     = "Model "+ model_name
   2.105 +  title     = "Model "+ model_name + " (2000-2004)"
   2.106    resg@tiMainString  = title
   2.107  
   2.108    plot(1) = gsn_csm_contour_map_ce(wks,data_mod,resg) 
   2.109 @@ -571,7 +581,7 @@
   2.110  
   2.111    zz = data_mod
   2.112    zz = data_mod - data_ob
   2.113 -  title = "Model_"+model_name+" - Observed"
   2.114 +  title = "Model_"+model_name+" - MODIS MOD15A2"
   2.115    resg@tiMainString    = title
   2.116  
   2.117    plot(2) = gsn_csm_contour_map_ce(wks,zz,resg) 
   2.118 @@ -586,8 +596,9 @@
   2.119    delete (wks)
   2.120    delete (plot)
   2.121  
   2.122 -  system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
   2.123 +  system("convert -trim -density 150 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
   2.124           "rm "+plot_name+"."+plot_type)
   2.125 +  ;system("convert -trim -density 150 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
   2.126  
   2.127    delete (data_ob)
   2.128    delete (data_mod)
   2.129 @@ -833,7 +844,7 @@
   2.130       delete (wks)
   2.131       delete (plot)
   2.132  
   2.133 -     system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
   2.134 +     system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
   2.135              "rm "+plot_name+"."+plot_type)
   2.136    end do
   2.137  
   2.138 @@ -857,7 +868,7 @@
   2.139       delete (wks)
   2.140       delete (plot)
   2.141  
   2.142 -     system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
   2.143 +     system("convert -trim -density 100 "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
   2.144              "rm "+plot_name+"."+plot_type)
   2.145  
   2.146    delete (plot_data)

     3.1 --- a/all/03.co2.ncl	Mon Jan 26 22:08:20 2009 -0500
     3.2 +++ b/all/03.co2.ncl	Thu Mar 26 14:02:21 2009 -0400
     3.3 @@ -186,7 +186,7 @@
     3.4  
     3.5    if (z .eq. 3) then 
     3.6       zone = "90S-EQ" 
     3.7 -     score_max = 5.0
     3.8 +     score_max = 0.0
     3.9       ind_z = ind(lat_ob .lt. 0. )
    3.10    end if
    3.11  

     4.1 --- a/all/04.biomass.ncl	Mon Jan 26 22:08:20 2009 -0500
     4.2 +++ b/all/04.biomass.ncl	Thu Mar 26 14:02:21 2009 -0400
     4.3 @@ -44,7 +44,7 @@
     4.4  ;--------------------------------------------
     4.5  ; read model data
     4.6  
     4.7 - fm   = addfile (dirm+film1,"r")
     4.8 + fm   = addfile (dirm+film11,"r")
     4.9  
    4.10   if (BGC .eq. "cn") then 
    4.11      data1  = fm->LIVESTEMC
    4.12 @@ -83,7 +83,7 @@
    4.13  ;------------------------------------------------
    4.14  ; read ob data
    4.15  
    4.16 - ob_name = "LC15_Amazon_Biomass"
    4.17 + ob_name = "LC15 Amazon Biomass"
    4.18  
    4.19   dir_b = diro + "biomass/"
    4.20   fil_b = "amazon_biomass_"+model_grid+".nc"
    4.21 @@ -99,9 +99,9 @@
    4.22  ; Units for these variables are:
    4.23  ; dataob   : MgC/ha
    4.24  ; datamod0 : gC/m2
    4.25 -; We want to convert these to KgC/m2
    4.26 +; We want to convert these to kgC/m2
    4.27  ; ha = 100m*100m = 10,000 m2
    4.28 -; MgC/ha*1000/10,000 = KgC/m2
    4.29 +; MgC/ha*1000/10,000 = kgC/m2
    4.30  
    4.31    factor_aboveground = 0.5
    4.32    factor_unit_ob     = 0.1
    4.33 @@ -110,8 +110,8 @@
    4.34    dataob   = dataob * factor_aboveground * factor_unit_ob
    4.35    datamod0 = datamod0 * factor_unit_mod 
    4.36  
    4.37 -  dataob@units      = "KgC/m2"
    4.38 -  datamod0@units    = "KgC/m2"
    4.39 +  dataob@units      = "kg C m~S~-2~N~"
    4.40 +  datamod0@units    = "kg C m~S~-2~N~"
    4.41  
    4.42    dataob@long_name      = "Amazon Biomass"
    4.43    datamod0@long_name    = "Amazon Biomass"
    4.44 @@ -139,7 +139,7 @@
    4.45  ;********************************************************
    4.46  ; sum over amazom_mask area:
    4.47  
    4.48 -; Peta g = 1.e15 g = 1.e12 Kg
    4.49 +; Peta g = 1.e15 g = 1.e12 kg
    4.50    factor_unit = 1.e-12
    4.51  
    4.52  ; mask_amazon = where(mask_amazon .ge. 0.5, mask_amazon ,0.)
    4.53 @@ -204,6 +204,7 @@
    4.54  
    4.55    system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
    4.56           "rm "+plot_name+"."+plot_type)
    4.57 +  ;system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
    4.58  
    4.59  ;------------------------------------------------------------------------
    4.60  ; contour ob
    4.61 @@ -225,6 +226,7 @@
    4.62  
    4.63    system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
    4.64           "rm "+plot_name+"."+plot_type)
    4.65 +  ;system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
    4.66  
    4.67  ;------------------------------------------------------------------------
    4.68  ; contour model
    4.69 @@ -234,7 +236,7 @@
    4.70    resg@cnLevelSpacingF      = 2.              ; interval
    4.71  
    4.72    plot_name = "global_model"
    4.73 -  title     = "Model "+ model_name 
    4.74 +  title     = "Model "+ model_name + " (2000)"
    4.75    resg@tiMainString  = title
    4.76  
    4.77    wks = gsn_open_wks (plot_type,plot_name)   ; open workstation
    4.78 @@ -246,6 +248,7 @@
    4.79  
    4.80    system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
    4.81           "rm "+plot_name+"."+plot_type)
    4.82 +  ;system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
    4.83  
    4.84  ;------------------------------------------------------------------------
    4.85  ; contour model vs ob
    4.86 @@ -273,13 +276,17 @@
    4.87  
    4.88    ccrG = esccr(ug,vg,0)
    4.89  
    4.90 -  score_max = 5.0
    4.91 +  score_max = 10.0
    4.92  
    4.93  ; Miomass = (ccrG*ccrG)* score_max
    4.94  ; new eq
    4.95    bias = sum(abs(ug-vg)/(abs(ug)+abs(vg)))
    4.96 +  print ("bias=" + bias + "dimsizes(ug)=" + dimsizes(ug))
    4.97 +  print ("other bias=" + sum(abs(ug-vg)/(ug+vg)))
    4.98    Mbiomass  = (1. - (bias/dimsizes(ug)))*score_max
    4.99 +  print ("Mbiomass=" + Mbiomass)
   4.100    M_biomass = sprintf("%.2f", Mbiomass)
   4.101 +  print ("M_biomass=" + M_biomass)
   4.102  
   4.103    if (isvar("compare")) then
   4.104       system("sed -e '1,/M_biomass/s/M_biomass/"+M_biomass+"/' "+html_name2+" > "+html_new2+";"+ \ 
   4.105 @@ -309,7 +316,7 @@
   4.106  
   4.107  ;(b) model
   4.108  
   4.109 -  title     = "Model "+ model_name
   4.110 +  title     = "Model "+ model_name + " (2000)"
   4.111    resg@tiMainString  = title
   4.112  
   4.113    plot(1) = gsn_csm_contour_map_ce(wks,datamod,resg) 
   4.114 @@ -318,7 +325,7 @@
   4.115  
   4.116    zz = datamod
   4.117    zz = datamod - dataob
   4.118 -  title = "Model_"+model_name+" - Observed"
   4.119 +  title = "Model "+model_name+" (2000) - Observed"
   4.120  
   4.121    resg@cnMinLevelValF  = -10.          ; Min level
   4.122    resg@cnMaxLevelValF  =  10.          ; Max level
   4.123 @@ -340,6 +347,7 @@
   4.124  
   4.125    system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
   4.126           "rm "+plot_name+"."+plot_type)
   4.127 +  ;system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
   4.128  
   4.129    resg@gsnFrame             = True            ; draw plot 
   4.130    resg@gsnDraw              = True            ; advance frame
   4.131 @@ -363,7 +371,7 @@
   4.132    gRes@txFontHeightF = 0.02
   4.133    gRes@txAngleF = 0
   4.134  
   4.135 -  area_avg_text = "(average over mask area = "+sprintf("%.2f", avg_ob)+" Kg C/m2)"
   4.136 +  area_avg_text = "(average over mask area = "+sprintf("%.2f", avg_ob)+" kg C m~S~-2~N~)"
   4.137  
   4.138    gsn_text_ndc(wks,area_avg_text,0.50,0.81,gRes)
   4.139  ;-----------------------------------------
   4.140 @@ -377,6 +385,7 @@
   4.141  
   4.142    system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
   4.143           "rm "+plot_name+"."+plot_type)
   4.144 +  ;system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
   4.145  
   4.146  ;------------------------------------------------------------------------
   4.147  ; contour model: masked
   4.148 @@ -386,7 +395,7 @@
   4.149    resg@cnLevelSpacingF      = 2.              ; interval
   4.150  
   4.151    plot_name = "global_mask_model"
   4.152 -  title     = "Model "+ model_name 
   4.153 +  title     = "Model "+ model_name + " (2000)"
   4.154    resg@tiMainString  = title
   4.155  
   4.156    wks = gsn_open_wks (plot_type,plot_name)   ; open workstation
   4.157 @@ -398,7 +407,7 @@
   4.158    gRes@txFontHeightF = 0.02
   4.159    gRes@txAngleF = 0
   4.160  
   4.161 -  area_avg_text = "(average over mask area = "+sprintf("%.2f", avg_mod)+" Kg C/m2)"
   4.162 +  area_avg_text = "(average over mask area = "+sprintf("%.2f", avg_mod)+" kg C m~S~-2~N~)"
   4.163  
   4.164    gsn_text_ndc(wks,area_avg_text,0.50,0.81,gRes)
   4.165  ;-----------------------------------------
   4.166 @@ -412,6 +421,7 @@
   4.167  
   4.168    system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
   4.169           "rm "+plot_name+"."+plot_type)
   4.170 +  ;system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
   4.171  
   4.172  ;------------------------------------------------------------------------
   4.173  ; contour model vs ob: masked 
   4.174 @@ -434,7 +444,7 @@
   4.175    delete (ug)
   4.176    delete (vg)  
   4.177  
   4.178 -  score_max = 5.
   4.179 +  score_max = 0.
   4.180  
   4.181    uu = ndtooned(zm)
   4.182    vv = ndtooned(zo)
   4.183 @@ -488,7 +498,7 @@
   4.184  
   4.185  ;(b) model
   4.186  
   4.187 -  title     = "Model "+ model_name
   4.188 +  title     = "Model "+ model_name + " (2000)"
   4.189    resg@tiMainString  = title
   4.190  
   4.191    plot(1) = gsn_csm_contour_map_ce(wks,zm,resg) 
   4.192 @@ -497,7 +507,7 @@
   4.193  
   4.194    zz = zo
   4.195    zz = zm - zo
   4.196 -  title = "Model_"+model_name+" - Observed"
   4.197 +  title = "Model "+model_name+" (2000) - Observed"
   4.198  
   4.199    resg@cnMinLevelValF  = -10.          ; Min level
   4.200    resg@cnMaxLevelValF  =  10.          ; Max level
   4.201 @@ -518,6 +528,7 @@
   4.202  
   4.203    system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
   4.204           "rm "+plot_name+"."+plot_type)
   4.205 +  ;system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new)
   4.206  
   4.207  ;***************************************************************************
   4.208  ; add total score and write to file

     5.1 --- a/all/06.fluxnet.ncl	Mon Jan 26 22:08:20 2009 -0500
     5.2 +++ b/all/06.fluxnet.ncl	Thu Mar 26 14:02:21 2009 -0400
     5.3 @@ -530,7 +530,8 @@
     5.4  ;***************************************************************************
     5.5  ; add total score and write to file
     5.6  ;***************************************************************************
     5.7 -  M_total = M_fluxnet_all
     5.8 +  ;M_total = M_fluxnet_all
     5.9 +  M_total = 0.00 ; we don't actually use these scores
    5.10  
    5.11    asciiwrite("M_save.fluxnet", M_total)
    5.12  

     6.1 --- a/all/07.beta.ncl	Mon Jan 26 22:08:20 2009 -0500
     6.2 +++ b/all/07.beta.ncl	Thu Mar 26 14:02:21 2009 -0400
     6.3 @@ -136,7 +136,7 @@
     6.4  ;---------------------------------------------------
     6.5  ;compute beta_4
     6.6  
     6.7 - score_max = 3.
     6.8 + score_max = 10.
     6.9  
    6.10   beta_4 = new((/n_sta/),float)
    6.11  

     7.1 --- a/all/10.fire.ncl	Mon Jan 26 22:08:20 2009 -0500
     7.2 +++ b/all/10.fire.ncl	Thu Mar 26 14:02:21 2009 -0400
     7.3 @@ -74,7 +74,10 @@
     7.4  
     7.5   fm   = addfile (dirm+film7,"r")
     7.6  
     7.7 + ; FMH: This was Jeff's, but I got an error
     7.8   data_mod = fm->COL_FIRE_CLOSS(18:25,:,:,:)
     7.9 + ; FMH: So I'll try this
    7.10 + ;data_mod = fm->COL_FIRE_CLOSS(17:24,:,:,:)
    7.11  
    7.12   delete (fm)
    7.13  
    7.14 @@ -192,7 +195,7 @@
    7.15  ;----------------------------------------------------
    7.16  ; compute M_global
    7.17  
    7.18 -  score_max = 1.
    7.19 +  score_max = 2.5
    7.20  
    7.21    Mscore1 = cc_space * cc_space * score_max
    7.22  
    7.23 @@ -542,7 +545,7 @@
    7.24  ; unit change from TgC/month to PgC/month
    7.25    unit_factor = 1.e-3
    7.26  
    7.27 -  score_max = 1.
    7.28 +  score_max = 2.5
    7.29  
    7.30    tmp_ob    = new((/ntime/),float)
    7.31    tmp_mod   = new((/ntime/),float)
    7.32 @@ -607,7 +610,7 @@
    7.33  ;------------------------------------------
    7.34  ; data for table : global
    7.35  
    7.36 -  score_max = 1.
    7.37 +  score_max = 2.5
    7.38  
    7.39    tmp_ob  = ndtooned(yvalues_g(:,0,:))
    7.40    tmp_mod = ndtooned(yvalues_g(:,1,:))

     8.1 --- a/all/11.ameriflux.ncl	Mon Jan 26 22:08:20 2009 -0500
     8.2 +++ b/all/11.ameriflux.ncl	Thu Mar 26 14:02:21 2009 -0400
     8.3 @@ -557,7 +557,8 @@
     8.4  ; compute annual cycle correlation coef and M score
     8.5  ;************************************************************
     8.6  
     8.7 - score_max = 1.
     8.8 + ;score_max = 1.
     8.9 + score_max = (/6., 0., 9., 9., 6., 0./)
    8.10  
    8.11   ccr     = new ((/nstation, nfield/),float)
    8.12   M_score = new ((/nstation, nfield/),float) 
    8.13 @@ -566,7 +567,7 @@
    8.14   do m=0,nfield-1   
    8.15      ccr(n,m) = esccr(data_ob_ann(m,:,n),data_mod_ann(m,:,n),0)
    8.16      bias = sum(abs(data_mod_ann(m,:,n)-data_ob_ann(m,:,n))/(abs(data_mod_ann(m,:,n))+abs(data_ob_ann(m,:,n))))
    8.17 -    M_score(n,m) = (1. -(bias/nmonth)) * score_max
    8.18 +    M_score(n,m) = (1. -(bias/nmonth)) * score_max(m)
    8.19   end do
    8.20   end do
    8.21  

     9.1 --- a/all/run_1-model.csh	Mon Jan 26 22:08:20 2009 -0500
     9.2 +++ b/all/run_1-model.csh	Thu Mar 26 14:02:21 2009 -0400
     9.3 @@ -58,7 +58,10 @@
     9.4  else
     9.5     set TEMPLATE = template_1-model_noCO2
     9.6  endif
     9.7 -cp -r $DIR_SCRIPTS/$TEMPLATE $MODEL
     9.8 +#cp -r $DIR_SCRIPTS/$TEMPLATE $MODEL
     9.9 +rm -rf $MODEL
    9.10 +tar xvpf $DIR_SCRIPTS/$TEMPLATE.tar
    9.11 +mv $TEMPLATE $MODEL
    9.12  
    9.13  # add quote, to be usesd in INPUT_TEXT
    9.14  set MODELQ = \"$MODEL\"

    10.1 --- a/all/run_2-model.csh	Mon Jan 26 22:08:20 2009 -0500
    10.2 +++ b/all/run_2-model.csh	Thu Mar 26 14:02:21 2009 -0400
    10.3 @@ -70,8 +70,14 @@
    10.4     set TEMPLATE1 = template_1-model_noCO2
    10.5     set TEMPLATE2 = template_2-model_noCO2
    10.6  endif
    10.7 -cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
    10.8 -cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
    10.9 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
   10.10 +#cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
   10.11 +rm -rf $MODEL1
   10.12 +rm -rf $COMPARE
   10.13 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
   10.14 +mv $TEMPLATE1 $MODEL1
   10.15 +tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
   10.16 +mv $TEMPLATE2 $COMPARE
   10.17  
   10.18  # add quote, to be usesd in INPUT_TEXT
   10.19  set MODELQ = \"$MODEL1\"
   10.20 @@ -147,7 +153,11 @@
   10.21  else
   10.22     set TEMPLATE1 = template_1-model_noCO2
   10.23  endif
   10.24 -cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
   10.25 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
   10.26 +rm -rf $MODEL2
   10.27 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
   10.28 +mv $TEMPLATE1 $MODEL2
   10.29 +
   10.30  
   10.31  # add quote, to be usesd in INPUT_TEXT
   10.32  set MODELQ  = \"$MODEL2\"

    11.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
    11.2 +++ b/all/run_2-model_36cn_37cn_cont.csh	Thu Mar 26 14:02:21 2009 -0400
    11.3 @@ -0,0 +1,247 @@
    11.4 +#!/bin/csh -f
    11.5 +#-------------------------------------------------------------------
    11.6 +# note: user modifies ONLY the "user modification" section 
    11.7 +#
    11.8 +#       COMPARE: model1 vs model2
    11.9 +#       MODELn : model name
   11.10 +#       DIR_M  : directory of model data
   11.11 +#       DIR_O  : directory of observed data
   11.12 +#       DIR_S  : directory of model surface data
   11.13 +#       DIR_SCRIPTS  : directory of run scripts
   11.14 +#       FILE1  : time_mean  climatology from CLM diagnostic package
   11.15 +#       FILE2  : 12-monthly climatology from CLM diagnostic package
   11.16 +#       FILE3  : 12-monthly climatology from ATM diagnostic package
   11.17 +#                leave it blank, if no ATM file:
   11.18 +#                set FILE3 =
   11.19 +#       FILE7 : timeseries file generated from .... 
   11.20 +#       FILE8 : timeseries file generated from .... 
   11.21 +#       GRID   : T31, T42, or 1.9
   11.22 +#       BGC    : cn or casa 
   11.23 +#       ENERGY : new or old (fields in model data)
   11.24 +#-------------------------------------------------------------------
   11.25 +
   11.26 +#*******************************************************
   11.27 +# user modification-(1)
   11.28 +
   11.29 +# directory name of model comparison
   11.30 +#et COMPARE = b30.061n_vs_b30.061m
   11.31 +set COMPARE = i01.36cn_vs_i01.37cn_cont
   11.32 +
   11.33 +#*******************************************************
   11.34 +# user modification-(2)
   11.35 +
   11.36 +# model1 
   11.37 +set MODEL1 = i01.36cn
   11.38 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.36cn_out/
   11.39 +set FILE1  = i01.36cn_out_ANN_climo.nc
   11.40 +set FILE2  = i01.36cn_out_MONS_climo.nc
   11.41 +set FILE3  = 
   11.42 +set FILE4  = $FILE1
   11.43 +set FILE5  = $FILE1
   11.44 +set FILE6  = $FILE1
   11.45 +set FILE7  = i01.36cn_Fire_C_2076-2100_monthly.nc
   11.46 +set FILE8  = i01.36cn_ameriflux_2076-2100_monthly.nc
   11.47 +set FILE9  = $FILE1
   11.48 +set FILE10 = $FILE2
   11.49 +set GRID   = 1.9
   11.50 +set BGC    = cn
   11.51 +set ENERGY = old
   11.52 +
   11.53 +# in the "CLAMP metric processing" section:
   11.54 +#  only 00.initial.ncl and 99.final.ncl are required,
   11.55 +#  user can comment out any one or more of the other ncl scripts,
   11.56 +#  e.g. 
   11.57 +#  #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
   11.58 +
   11.59 +# model surface data
   11.60 +set DIR_S  = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
   11.61 +
   11.62 +# observed data
   11.63 +set DIR_O  = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
   11.64 +
   11.65 +# directory for scripts, templates and ncl files
   11.66 +set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
   11.67 +
   11.68 +#********************************************************
   11.69 +
   11.70 +# create model1 and model1_vs_model2 directory by copying templates
   11.71 +if ($FILE3 != "") then
   11.72 +   set TEMPLATE1 = template_1-model
   11.73 +   set TEMPLATE2 = template_2-model
   11.74 +else
   11.75 +   set TEMPLATE1 = template_1-model_noCO2
   11.76 +   set TEMPLATE2 = template_2-model_noCO2
   11.77 +endif
   11.78 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
   11.79 +#cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
   11.80 +rm -rf $MODEL1
   11.81 +rm -rf $COMPARE
   11.82 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
   11.83 +mv $TEMPLATE1 $MODEL1
   11.84 +tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
   11.85 +mv $TEMPLATE2 $COMPARE
   11.86 +
   11.87 +# add quote, to be usesd in INPUT_TEXT
   11.88 +set MODELQ = \"$MODEL1\"
   11.89 +set DIRMQ  = \"$DIR_M\"
   11.90 +set F1  = \"$FILE1\"
   11.91 +set F2  = \"$FILE2\"
   11.92 +set F3  = \"$FILE3\"
   11.93 +set F4  = \"$FILE4\"
   11.94 +set F5  = \"$FILE5\"
   11.95 +set F6  = \"$FILE6\"
   11.96 +set F7  = \"$FILE7\"
   11.97 +set F8  = \"$FILE8\"
   11.98 +set F9  = \"$FILE9\"
   11.99 +set F10 = \"$FILE10\"
  11.100 +set GRIDQ   = \"$GRID\"
  11.101 +set BGCQ    = \"$BGC\"
  11.102 +set ENERGYQ = \"$ENERGY\"
  11.103 +set DIRSQ   = \"$DIR_S\"
  11.104 +set DIROQ   = \"$DIR_O\"
  11.105 +set DIRCQ   = \"$DIR_SCRIPTS\"
  11.106 +
  11.107 +set COMPAREQ = \"$COMPARE\"
  11.108 +set MODELN   = \"model1\"
  11.109 +
  11.110 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  11.111 + 
  11.112 +echo "`date` - Running diagnostics for model ${MODELQ}"
  11.113 +# CLAMP metric processing of model1
  11.114 +echo "`date` - Running 00.initial.ncl"
  11.115 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  11.116 +echo "`date` - Running 01.npp.ncl"
  11.117 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  11.118 +echo "`date` - Running 01.lai.ncl"
  11.119 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  11.120 +
  11.121 +if ($FILE3 != "") then
  11.122 +echo "`date` - Running 03.co2.ncl"
  11.123 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  11.124 +endif
  11.125 +
  11.126 +echo "`date` - Running 04.biomass.ncl"
  11.127 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  11.128 +echo "`date` - Running 06.fluxnet.ncl"
  11.129 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  11.130 +echo "`date` - Running 07.beta.ncl"
  11.131 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  11.132 +echo "`date` - Running 08.turnover.ncl"
  11.133 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  11.134 +
  11.135 +if ($BGC != "casa") then
  11.136 +echo "`date` - Running 09.carbon_sink.ncl"
  11.137 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  11.138 +else
  11.139 +echo "`date` - Running 09x.carbon_sink.ncl"
  11.140 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  11.141 +endif
  11.142 +
  11.143 +if ($BGC != "casa") then
  11.144 +echo "`date` - Running 10.fire.ncl"
  11.145 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  11.146 +endif
  11.147 +
  11.148 +echo "`date` - Running 11.ameriflux.ncl"
  11.149 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  11.150 +echo "`date` - Running 99.final.ncl"
  11.151 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  11.152 +
  11.153 +#*******************************************************
  11.154 +# user modification-(3)
  11.155 +
  11.156 +# model2
  11.157 +set MODEL2 = i01.37cn_cont
  11.158 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases/i01.37cn_cont_out/
  11.159 +set FILE1  = i01.37cn_cont_out_ANN_climo.nc
  11.160 +set FILE2  = i01.37cn_cont_out_MONS_climo.nc
  11.161 +set FILE3  = 
  11.162 +set FILE4  = $FILE1
  11.163 +set FILE5  = $FILE1
  11.164 +set FILE6  = $FILE1
  11.165 +set FILE7  = i01.37cn_cont_Fire_C_2976-3000_monthly.nc
  11.166 +set FILE8  = i01.37cn_cont_ameriflux_2976-3000_monthly.nc
  11.167 +set FILE9  = $FILE1
  11.168 +set FILE10 = $FILE2
  11.169 +set GRID   = 1.9
  11.170 +set BGC    = cn
  11.171 +set ENERGY = old
  11.172 +#*******************************************************
  11.173 +# create model2 directory by copying templates
  11.174 +if ($FILE3 != "") then
  11.175 +   set TEMPLATE1 = template_1-model
  11.176 +else
  11.177 +   set TEMPLATE1 = template_1-model_noCO2
  11.178 +endif
  11.179 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
  11.180 +rm -rf $MODEL2
  11.181 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
  11.182 +mv $TEMPLATE1 $MODEL2
  11.183 +
  11.184 +# add quote, to be usesd in INPUT_TEXT
  11.185 +set MODELQ  = \"$MODEL2\"
  11.186 +set DIRMQ   = \"$DIR_M\"
  11.187 +set F1  = \"$FILE1\"
  11.188 +set F2  = \"$FILE2\"
  11.189 +set F3  = \"$FILE3\"
  11.190 +set F4  = \"$FILE4\"
  11.191 +set F5  = \"$FILE5\"
  11.192 +set F6  = \"$FILE6\"
  11.193 +set F7  = \"$FILE7\"
  11.194 +set F8  = \"$FILE8\"
  11.195 +set F9  = \"$FILE9\"
  11.196 +set F10 = \"$FILE10\"
  11.197 +set GRIDQ   = \"$GRID\"
  11.198 +set BGCQ    = \"$BGC\"
  11.199 +set ENERGYQ = \"$ENERGY\"
  11.200 +
  11.201 +set COMPAREQ = \"$COMPARE\"
  11.202 +set MODELN   = \"model2\"
  11.203 +
  11.204 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  11.205 + 
  11.206 +echo "`date` - Running diagnostics for model ${MODELQ}"
  11.207 +# CLAMP metric processing of model2
  11.208 +echo "`date` - Running 00.inital.ncl"
  11.209 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  11.210 +echo "`date` - Running 01.npp.ncl"
  11.211 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  11.212 +echo "`date` - Running 02.lai.ncl"
  11.213 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  11.214 +
  11.215 +if ($FILE3 != "") then
  11.216 +echo "`date` - Running 03.co2.ncl"
  11.217 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  11.218 +endif
  11.219 +
  11.220 +echo "`date` - Running 04.biomass.ncl"
  11.221 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  11.222 +echo "`date` - Running 06.fluxnet.ncl"
  11.223 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  11.224 +echo "`date` - Running 07.beta.ncl"
  11.225 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  11.226 +echo "`date` - Running 08.turnover.ncl"
  11.227 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  11.228 +
  11.229 +if ($BGC != "casa") then
  11.230 +echo "`date` - Running 09.carbon_sink.ncl"
  11.231 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  11.232 +else
  11.233 +echo "`date` - Running 09x.carbon_sink.ncl"
  11.234 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  11.235 +endif
  11.236 +
  11.237 +if ($BGC != "casa") then
  11.238 +echo "`date` - Running 10.fire.ncl"
  11.239 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  11.240 +endif
  11.241 +
  11.242 +echo "`date` - Running 11.ameriflux.ncl"
  11.243 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  11.244 +echo "`date` - Running 99.final.ncl"
  11.245 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  11.246 +
  11.247 +# create a tar file from the final output
  11.248 +echo "`date` - tarring up final output into all.tar"
  11.249 +tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
  11.250 +

    12.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
    12.2 +++ b/all/run_2-model_36cn_38cn.csh	Thu Mar 26 14:02:21 2009 -0400
    12.3 @@ -0,0 +1,247 @@
    12.4 +#!/bin/csh -f
    12.5 +#-------------------------------------------------------------------
    12.6 +# note: user modifies ONLY the "user modification" section 
    12.7 +#
    12.8 +#       COMPARE: model1 vs model2
    12.9 +#       MODELn : model name
   12.10 +#       DIR_M  : directory of model data
   12.11 +#       DIR_O  : directory of observed data
   12.12 +#       DIR_S  : directory of model surface data
   12.13 +#       DIR_SCRIPTS  : directory of run scripts
   12.14 +#       FILE1  : time_mean  climatology from CLM diagnostic package
   12.15 +#       FILE2  : 12-monthly climatology from CLM diagnostic package
   12.16 +#       FILE3  : 12-monthly climatology from ATM diagnostic package
   12.17 +#                leave it blank, if no ATM file:
   12.18 +#                set FILE3 =
   12.19 +#       FILE7 : timeseries file generated from .... 
   12.20 +#       FILE8 : timeseries file generated from .... 
   12.21 +#       GRID   : T31, T42, or 1.9
   12.22 +#       BGC    : cn or casa 
   12.23 +#       ENERGY : new or old (fields in model data)
   12.24 +#-------------------------------------------------------------------
   12.25 +
   12.26 +#*******************************************************
   12.27 +# user modification-(1)
   12.28 +
   12.29 +# directory name of model comparison
   12.30 +#et COMPARE = b30.061n_vs_b30.061m
   12.31 +set COMPARE = i01.36cn_vs_i01.38cn
   12.32 +
   12.33 +#*******************************************************
   12.34 +# user modification-(2)
   12.35 +
   12.36 +# model1 
   12.37 +set MODEL1 = i01.36cn
   12.38 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.36cn_out/
   12.39 +set FILE1  = i01.36cn_out_ANN_climo.nc
   12.40 +set FILE2  = i01.36cn_out_MONS_climo.nc
   12.41 +set FILE3  = 
   12.42 +set FILE4  = $FILE1
   12.43 +set FILE5  = $FILE1
   12.44 +set FILE6  = $FILE1
   12.45 +set FILE7  = i01.36cn_Fire_C_2076-2100_monthly.nc
   12.46 +set FILE8  = i01.36cn_ameriflux_2076-2100_monthly.nc
   12.47 +set FILE9  = $FILE1
   12.48 +set FILE10 = $FILE2
   12.49 +set GRID   = 1.9
   12.50 +set BGC    = cn
   12.51 +set ENERGY = old
   12.52 +
   12.53 +# in the "CLAMP metric processing" section:
   12.54 +#  only 00.initial.ncl and 99.final.ncl are required,
   12.55 +#  user can comment out any one or more of the other ncl scripts,
   12.56 +#  e.g. 
   12.57 +#  #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
   12.58 +
   12.59 +# model surface data
   12.60 +set DIR_S  = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
   12.61 +
   12.62 +# observed data
   12.63 +set DIR_O  = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
   12.64 +
   12.65 +# directory for scripts, templates and ncl files
   12.66 +set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
   12.67 +
   12.68 +#********************************************************
   12.69 +
   12.70 +# create model1 and model1_vs_model2 directory by copying templates
   12.71 +if ($FILE3 != "") then
   12.72 +   set TEMPLATE1 = template_1-model
   12.73 +   set TEMPLATE2 = template_2-model
   12.74 +else
   12.75 +   set TEMPLATE1 = template_1-model_noCO2
   12.76 +   set TEMPLATE2 = template_2-model_noCO2
   12.77 +endif
   12.78 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
   12.79 +#cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
   12.80 +rm -rf $MODEL1
   12.81 +rm -rf $COMPARE
   12.82 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
   12.83 +mv $TEMPLATE1 $MODEL1
   12.84 +tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
   12.85 +mv $TEMPLATE2 $COMPARE
   12.86 +
   12.87 +# add quote, to be usesd in INPUT_TEXT
   12.88 +set MODELQ = \"$MODEL1\"
   12.89 +set DIRMQ  = \"$DIR_M\"
   12.90 +set F1  = \"$FILE1\"
   12.91 +set F2  = \"$FILE2\"
   12.92 +set F3  = \"$FILE3\"
   12.93 +set F4  = \"$FILE4\"
   12.94 +set F5  = \"$FILE5\"
   12.95 +set F6  = \"$FILE6\"
   12.96 +set F7  = \"$FILE7\"
   12.97 +set F8  = \"$FILE8\"
   12.98 +set F9  = \"$FILE9\"
   12.99 +set F10 = \"$FILE10\"
  12.100 +set GRIDQ   = \"$GRID\"
  12.101 +set BGCQ    = \"$BGC\"
  12.102 +set ENERGYQ = \"$ENERGY\"
  12.103 +set DIRSQ   = \"$DIR_S\"
  12.104 +set DIROQ   = \"$DIR_O\"
  12.105 +set DIRCQ   = \"$DIR_SCRIPTS\"
  12.106 +
  12.107 +set COMPAREQ = \"$COMPARE\"
  12.108 +set MODELN   = \"model1\"
  12.109 +
  12.110 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  12.111 + 
  12.112 +echo "`date` - Running diagnostics for model ${MODELQ}"
  12.113 +# CLAMP metric processing of model1
  12.114 +echo "`date` - Running 00.initial.ncl"
  12.115 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  12.116 +echo "`date` - Running 01.npp.ncl"
  12.117 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  12.118 +echo "`date` - Running 01.lai.ncl"
  12.119 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  12.120 +
  12.121 +if ($FILE3 != "") then
  12.122 +echo "`date` - Running 03.co2.ncl"
  12.123 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  12.124 +endif
  12.125 +
  12.126 +echo "`date` - Running 04.biomass.ncl"
  12.127 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  12.128 +echo "`date` - Running 06.fluxnet.ncl"
  12.129 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  12.130 +echo "`date` - Running 07.beta.ncl"
  12.131 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  12.132 +echo "`date` - Running 08.turnover.ncl"
  12.133 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  12.134 +
  12.135 +if ($BGC != "casa") then
  12.136 +echo "`date` - Running 09.carbon_sink.ncl"
  12.137 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  12.138 +else
  12.139 +echo "`date` - Running 09x.carbon_sink.ncl"
  12.140 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  12.141 +endif
  12.142 +
  12.143 +if ($BGC != "casa") then
  12.144 +echo "`date` - Running 10.fire.ncl"
  12.145 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  12.146 +endif
  12.147 +
  12.148 +echo "`date` - Running 11.ameriflux.ncl"
  12.149 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  12.150 +echo "`date` - Running 99.final.ncl"
  12.151 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  12.152 +
  12.153 +#*******************************************************
  12.154 +# user modification-(3)
  12.155 +
  12.156 +# model2
  12.157 +set MODEL2 = i01.38cn
  12.158 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.38cn_out/
  12.159 +set FILE1  = i01.38cn_out_ANN_climo.nc
  12.160 +set FILE2  = i01.38cn_out_MONS_climo.nc
  12.161 +set FILE3  = 
  12.162 +set FILE4  = $FILE1
  12.163 +set FILE5  = $FILE1
  12.164 +set FILE6  = $FILE1
  12.165 +set FILE7  = i01.38cn_Fire_C_2151-2175_monthly.nc
  12.166 +set FILE8  = i01.38cn_ameriflux_2151-2175_monthly.nc
  12.167 +set FILE9  = $FILE1
  12.168 +set FILE10 = $FILE2
  12.169 +set GRID   = 1.9
  12.170 +set BGC    = cn
  12.171 +set ENERGY = old
  12.172 +#*******************************************************
  12.173 +# create model2 directory by copying templates
  12.174 +if ($FILE3 != "") then
  12.175 +   set TEMPLATE1 = template_1-model
  12.176 +else
  12.177 +   set TEMPLATE1 = template_1-model_noCO2
  12.178 +endif
  12.179 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
  12.180 +rm -rf $MODEL2
  12.181 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
  12.182 +mv $TEMPLATE1 $MODEL2
  12.183 +
  12.184 +# add quote, to be usesd in INPUT_TEXT
  12.185 +set MODELQ  = \"$MODEL2\"
  12.186 +set DIRMQ   = \"$DIR_M\"
  12.187 +set F1  = \"$FILE1\"
  12.188 +set F2  = \"$FILE2\"
  12.189 +set F3  = \"$FILE3\"
  12.190 +set F4  = \"$FILE4\"
  12.191 +set F5  = \"$FILE5\"
  12.192 +set F6  = \"$FILE6\"
  12.193 +set F7  = \"$FILE7\"
  12.194 +set F8  = \"$FILE8\"
  12.195 +set F9  = \"$FILE9\"
  12.196 +set F10 = \"$FILE10\"
  12.197 +set GRIDQ   = \"$GRID\"
  12.198 +set BGCQ    = \"$BGC\"
  12.199 +set ENERGYQ = \"$ENERGY\"
  12.200 +
  12.201 +set COMPAREQ = \"$COMPARE\"
  12.202 +set MODELN   = \"model2\"
  12.203 +
  12.204 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  12.205 + 
  12.206 +echo "`date` - Running diagnostics for model ${MODELQ}"
  12.207 +# CLAMP metric processing of model2
  12.208 +echo "`date` - Running 00.inital.ncl"
  12.209 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  12.210 +echo "`date` - Running 01.npp.ncl"
  12.211 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  12.212 +echo "`date` - Running 02.lai.ncl"
  12.213 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  12.214 +
  12.215 +if ($FILE3 != "") then
  12.216 +echo "`date` - Running 03.co2.ncl"
  12.217 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  12.218 +endif
  12.219 +
  12.220 +echo "`date` - Running 04.biomass.ncl"
  12.221 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  12.222 +echo "`date` - Running 06.fluxnet.ncl"
  12.223 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  12.224 +echo "`date` - Running 07.beta.ncl"
  12.225 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  12.226 +echo "`date` - Running 08.turnover.ncl"
  12.227 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  12.228 +
  12.229 +if ($BGC != "casa") then
  12.230 +echo "`date` - Running 09.carbon_sink.ncl"
  12.231 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  12.232 +else
  12.233 +echo "`date` - Running 09x.carbon_sink.ncl"
  12.234 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  12.235 +endif
  12.236 +
  12.237 +if ($BGC != "casa") then
  12.238 +echo "`date` - Running 10.fire.ncl"
  12.239 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  12.240 +endif
  12.241 +
  12.242 +echo "`date` - Running 11.ameriflux.ncl"
  12.243 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  12.244 +echo "`date` - Running 99.final.ncl"
  12.245 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  12.246 +
  12.247 +# create a tar file from the final output
  12.248 +echo "`date` - tarring up final output into all.tar"
  12.249 +tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
  12.250 +

    13.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
    13.2 +++ b/all/run_2-model_37cn_cont_42cn.csh	Thu Mar 26 14:02:21 2009 -0400
    13.3 @@ -0,0 +1,247 @@
    13.4 +#!/bin/csh -f
    13.5 +#-------------------------------------------------------------------
    13.6 +# note: user modifies ONLY the "user modification" section 
    13.7 +#
    13.8 +#       COMPARE: model1 vs model2
    13.9 +#       MODELn : model name
   13.10 +#       DIR_M  : directory of model data
   13.11 +#       DIR_O  : directory of observed data
   13.12 +#       DIR_S  : directory of model surface data
   13.13 +#       DIR_SCRIPTS  : directory of run scripts
   13.14 +#       FILE1  : time_mean  climatology from CLM diagnostic package
   13.15 +#       FILE2  : 12-monthly climatology from CLM diagnostic package
   13.16 +#       FILE3  : 12-monthly climatology from ATM diagnostic package
   13.17 +#                leave it blank, if no ATM file:
   13.18 +#                set FILE3 =
   13.19 +#       FILE7 : timeseries file generated from .... 
   13.20 +#       FILE8 : timeseries file generated from .... 
   13.21 +#       GRID   : T31, T42, or 1.9
   13.22 +#       BGC    : cn or casa 
   13.23 +#       ENERGY : new or old (fields in model data)
   13.24 +#-------------------------------------------------------------------
   13.25 +
   13.26 +#*******************************************************
   13.27 +# user modification-(1)
   13.28 +
   13.29 +# directory name of model comparison
   13.30 +#et COMPARE = b30.061n_vs_b30.061m
   13.31 +set COMPARE = i01.37cn_cont_vs_i01.42cn
   13.32 +
   13.33 +#*******************************************************
   13.34 +# user modification-(2)
   13.35 +
   13.36 +# model1 
   13.37 +set MODEL1 = i01.37cn_cont
   13.38 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases/i01.37cn_cont_out/
   13.39 +set FILE1  = i01.37cn_cont_out_ANN_climo.nc
   13.40 +set FILE2  = i01.37cn_cont_out_MONS_climo.nc
   13.41 +set FILE3  = 
   13.42 +set FILE4  = $FILE1
   13.43 +set FILE5  = $FILE1
   13.44 +set FILE6  = $FILE1
   13.45 +set FILE7  = i01.37cn_cont_Fire_C_2976-3000_monthly.nc
   13.46 +set FILE8  = i01.37cn_cont_ameriflux_2976-3000_monthly.nc
   13.47 +set FILE9  = $FILE1
   13.48 +set FILE10 = $FILE2
   13.49 +set GRID   = 1.9
   13.50 +set BGC    = cn
   13.51 +set ENERGY = old
   13.52 +
   13.53 +# in the "CLAMP metric processing" section:
   13.54 +#  only 00.initial.ncl and 99.final.ncl are required,
   13.55 +#  user can comment out any one or more of the other ncl scripts,
   13.56 +#  e.g. 
   13.57 +#  #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
   13.58 +
   13.59 +# model surface data
   13.60 +set DIR_S  = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
   13.61 +
   13.62 +# observed data
   13.63 +set DIR_O  = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
   13.64 +
   13.65 +# directory for scripts, templates and ncl files
   13.66 +set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
   13.67 +
   13.68 +#********************************************************
   13.69 +
   13.70 +# create model1 and model1_vs_model2 directory by copying templates
   13.71 +if ($FILE3 != "") then
   13.72 +   set TEMPLATE1 = template_1-model
   13.73 +   set TEMPLATE2 = template_2-model
   13.74 +else
   13.75 +   set TEMPLATE1 = template_1-model_noCO2
   13.76 +   set TEMPLATE2 = template_2-model_noCO2
   13.77 +endif
   13.78 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
   13.79 +#cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
   13.80 +rm -rf $MODEL1
   13.81 +rm -rf $COMPARE
   13.82 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
   13.83 +mv $TEMPLATE1 $MODEL1
   13.84 +tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
   13.85 +mv $TEMPLATE2 $COMPARE
   13.86 +
   13.87 +# add quote, to be usesd in INPUT_TEXT
   13.88 +set MODELQ = \"$MODEL1\"
   13.89 +set DIRMQ  = \"$DIR_M\"
   13.90 +set F1  = \"$FILE1\"
   13.91 +set F2  = \"$FILE2\"
   13.92 +set F3  = \"$FILE3\"
   13.93 +set F4  = \"$FILE4\"
   13.94 +set F5  = \"$FILE5\"
   13.95 +set F6  = \"$FILE6\"
   13.96 +set F7  = \"$FILE7\"
   13.97 +set F8  = \"$FILE8\"
   13.98 +set F9  = \"$FILE9\"
   13.99 +set F10 = \"$FILE10\"
  13.100 +set GRIDQ   = \"$GRID\"
  13.101 +set BGCQ    = \"$BGC\"
  13.102 +set ENERGYQ = \"$ENERGY\"
  13.103 +set DIRSQ   = \"$DIR_S\"
  13.104 +set DIROQ   = \"$DIR_O\"
  13.105 +set DIRCQ   = \"$DIR_SCRIPTS\"
  13.106 +
  13.107 +set COMPAREQ = \"$COMPARE\"
  13.108 +set MODELN   = \"model1\"
  13.109 +
  13.110 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  13.111 + 
  13.112 +echo "`date` - Running diagnostics for model ${MODELQ}"
  13.113 +# CLAMP metric processing of model1
  13.114 +echo "`date` - Running 00.initial.ncl"
  13.115 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  13.116 +echo "`date` - Running 01.npp.ncl"
  13.117 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  13.118 +echo "`date` - Running 01.lai.ncl"
  13.119 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  13.120 +
  13.121 +if ($FILE3 != "") then
  13.122 +echo "`date` - Running 03.co2.ncl"
  13.123 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  13.124 +endif
  13.125 +
  13.126 +echo "`date` - Running 04.biomass.ncl"
  13.127 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  13.128 +echo "`date` - Running 06.fluxnet.ncl"
  13.129 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  13.130 +echo "`date` - Running 07.beta.ncl"
  13.131 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  13.132 +echo "`date` - Running 08.turnover.ncl"
  13.133 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  13.134 +
  13.135 +if ($BGC != "casa") then
  13.136 +echo "`date` - Running 09.carbon_sink.ncl"
  13.137 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  13.138 +else
  13.139 +echo "`date` - Running 09x.carbon_sink.ncl"
  13.140 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  13.141 +endif
  13.142 +
  13.143 +if ($BGC != "casa") then
  13.144 +echo "`date` - Running 10.fire.ncl"
  13.145 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  13.146 +endif
  13.147 +
  13.148 +echo "`date` - Running 11.ameriflux.ncl"
  13.149 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  13.150 +echo "`date` - Running 99.final.ncl"
  13.151 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  13.152 +
  13.153 +#*******************************************************
  13.154 +# user modification-(3)
  13.155 +
  13.156 +# model2
  13.157 +set MODEL2 = i01.42cn
  13.158 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases/i01.42cn_out/
  13.159 +set FILE1  = i01.42cn_out_ANN_climo.nc
  13.160 +set FILE2  = i01.42cn_out_MONS_climo.nc
  13.161 +set FILE3  = 
  13.162 +set FILE4  = $FILE1
  13.163 +set FILE5  = $FILE1
  13.164 +set FILE6  = $FILE1
  13.165 +set FILE7  = i01.42cn_Fire_C_2976-3000_monthly.nc
  13.166 +set FILE8  = i01.42cn_ameriflux_2976-3000_monthly.nc
  13.167 +set FILE9  = $FILE1
  13.168 +set FILE10 = $FILE2
  13.169 +set GRID   = 1.9
  13.170 +set BGC    = cn
  13.171 +set ENERGY = old
  13.172 +#*******************************************************
  13.173 +# create model2 directory by copying templates
  13.174 +if ($FILE3 != "") then
  13.175 +   set TEMPLATE1 = template_1-model
  13.176 +else
  13.177 +   set TEMPLATE1 = template_1-model_noCO2
  13.178 +endif
  13.179 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
  13.180 +rm -rf $MODEL2
  13.181 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
  13.182 +mv $TEMPLATE1 $MODEL2
  13.183 +
  13.184 +# add quote, to be usesd in INPUT_TEXT
  13.185 +set MODELQ  = \"$MODEL2\"
  13.186 +set DIRMQ   = \"$DIR_M\"
  13.187 +set F1  = \"$FILE1\"
  13.188 +set F2  = \"$FILE2\"
  13.189 +set F3  = \"$FILE3\"
  13.190 +set F4  = \"$FILE4\"
  13.191 +set F5  = \"$FILE5\"
  13.192 +set F6  = \"$FILE6\"
  13.193 +set F7  = \"$FILE7\"
  13.194 +set F8  = \"$FILE8\"
  13.195 +set F9  = \"$FILE9\"
  13.196 +set F10 = \"$FILE10\"
  13.197 +set GRIDQ   = \"$GRID\"
  13.198 +set BGCQ    = \"$BGC\"
  13.199 +set ENERGYQ = \"$ENERGY\"
  13.200 +
  13.201 +set COMPAREQ = \"$COMPARE\"
  13.202 +set MODELN   = \"model2\"
  13.203 +
  13.204 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  13.205 + 
  13.206 +echo "`date` - Running diagnostics for model ${MODELQ}"
  13.207 +# CLAMP metric processing of model2
  13.208 +echo "`date` - Running 00.inital.ncl"
  13.209 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  13.210 +echo "`date` - Running 01.npp.ncl"
  13.211 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  13.212 +echo "`date` - Running 02.lai.ncl"
  13.213 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  13.214 +
  13.215 +if ($FILE3 != "") then
  13.216 +echo "`date` - Running 03.co2.ncl"
  13.217 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  13.218 +endif
  13.219 +
  13.220 +echo "`date` - Running 04.biomass.ncl"
  13.221 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  13.222 +echo "`date` - Running 06.fluxnet.ncl"
  13.223 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  13.224 +echo "`date` - Running 07.beta.ncl"
  13.225 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  13.226 +echo "`date` - Running 08.turnover.ncl"
  13.227 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  13.228 +
  13.229 +if ($BGC != "casa") then
  13.230 +echo "`date` - Running 09.carbon_sink.ncl"
  13.231 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  13.232 +else
  13.233 +echo "`date` - Running 09x.carbon_sink.ncl"
  13.234 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  13.235 +endif
  13.236 +
  13.237 +if ($BGC != "casa") then
  13.238 +echo "`date` - Running 10.fire.ncl"
  13.239 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  13.240 +endif
  13.241 +
  13.242 +echo "`date` - Running 11.ameriflux.ncl"
  13.243 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  13.244 +echo "`date` - Running 99.final.ncl"
  13.245 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  13.246 +
  13.247 +# create a tar file from the final output
  13.248 +echo "`date` - tarring up final output into all.tar"
  13.249 +tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
  13.250 +

    14.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
    14.2 +++ b/all/run_2-model_37cn_cont_43cn.csh	Thu Mar 26 14:02:21 2009 -0400
    14.3 @@ -0,0 +1,247 @@
    14.4 +#!/bin/csh -f
    14.5 +#-------------------------------------------------------------------
    14.6 +# note: user modifies ONLY the "user modification" section 
    14.7 +#
    14.8 +#       COMPARE: model1 vs model2
    14.9 +#       MODELn : model name
   14.10 +#       DIR_M  : directory of model data
   14.11 +#       DIR_O  : directory of observed data
   14.12 +#       DIR_S  : directory of model surface data
   14.13 +#       DIR_SCRIPTS  : directory of run scripts
   14.14 +#       FILE1  : time_mean  climatology from CLM diagnostic package
   14.15 +#       FILE2  : 12-monthly climatology from CLM diagnostic package
   14.16 +#       FILE3  : 12-monthly climatology from ATM diagnostic package
   14.17 +#                leave it blank, if no ATM file:
   14.18 +#                set FILE3 =
   14.19 +#       FILE7 : timeseries file generated from .... 
   14.20 +#       FILE8 : timeseries file generated from .... 
   14.21 +#       GRID   : T31, T42, or 1.9
   14.22 +#       BGC    : cn or casa 
   14.23 +#       ENERGY : new or old (fields in model data)
   14.24 +#-------------------------------------------------------------------
   14.25 +
   14.26 +#*******************************************************
   14.27 +# user modification-(1)
   14.28 +
   14.29 +# directory name of model comparison
   14.30 +#et COMPARE = b30.061n_vs_b30.061m
   14.31 +set COMPARE = i01.37cn_cont_vs_i01.43cn
   14.32 +
   14.33 +#*******************************************************
   14.34 +# user modification-(2)
   14.35 +
   14.36 +# model1 
   14.37 +set MODEL1 = i01.37cn_cont
   14.38 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases/i01.37cn_cont_out/
   14.39 +set FILE1  = i01.37cn_cont_out_ANN_climo.nc
   14.40 +set FILE2  = i01.37cn_cont_out_MONS_climo.nc
   14.41 +set FILE3  = 
   14.42 +set FILE4  = $FILE1
   14.43 +set FILE5  = $FILE1
   14.44 +set FILE6  = $FILE1
   14.45 +set FILE7  = i01.37cn_cont_Fire_C_2976-3000_monthly.nc
   14.46 +set FILE8  = i01.37cn_cont_ameriflux_2976-3000_monthly.nc
   14.47 +set FILE9  = $FILE1
   14.48 +set FILE10 = $FILE2
   14.49 +set GRID   = 1.9
   14.50 +set BGC    = cn
   14.51 +set ENERGY = old
   14.52 +
   14.53 +# in the "CLAMP metric processing" section:
   14.54 +#  only 00.initial.ncl and 99.final.ncl are required,
   14.55 +#  user can comment out any one or more of the other ncl scripts,
   14.56 +#  e.g. 
   14.57 +#  #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
   14.58 +
   14.59 +# model surface data
   14.60 +set DIR_S  = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
   14.61 +
   14.62 +# observed data
   14.63 +set DIR_O  = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
   14.64 +
   14.65 +# directory for scripts, templates and ncl files
   14.66 +set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
   14.67 +
   14.68 +#********************************************************
   14.69 +
   14.70 +# create model1 and model1_vs_model2 directory by copying templates
   14.71 +if ($FILE3 != "") then
   14.72 +   set TEMPLATE1 = template_1-model
   14.73 +   set TEMPLATE2 = template_2-model
   14.74 +else
   14.75 +   set TEMPLATE1 = template_1-model_noCO2
   14.76 +   set TEMPLATE2 = template_2-model_noCO2
   14.77 +endif
   14.78 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
   14.79 +#cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
   14.80 +rm -rf $MODEL1
   14.81 +rm -rf $COMPARE
   14.82 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
   14.83 +mv $TEMPLATE1 $MODEL1
   14.84 +tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
   14.85 +mv $TEMPLATE2 $COMPARE
   14.86 +
   14.87 +# add quote, to be usesd in INPUT_TEXT
   14.88 +set MODELQ = \"$MODEL1\"
   14.89 +set DIRMQ  = \"$DIR_M\"
   14.90 +set F1  = \"$FILE1\"
   14.91 +set F2  = \"$FILE2\"
   14.92 +set F3  = \"$FILE3\"
   14.93 +set F4  = \"$FILE4\"
   14.94 +set F5  = \"$FILE5\"
   14.95 +set F6  = \"$FILE6\"
   14.96 +set F7  = \"$FILE7\"
   14.97 +set F8  = \"$FILE8\"
   14.98 +set F9  = \"$FILE9\"
   14.99 +set F10 = \"$FILE10\"
  14.100 +set GRIDQ   = \"$GRID\"
  14.101 +set BGCQ    = \"$BGC\"
  14.102 +set ENERGYQ = \"$ENERGY\"
  14.103 +set DIRSQ   = \"$DIR_S\"
  14.104 +set DIROQ   = \"$DIR_O\"
  14.105 +set DIRCQ   = \"$DIR_SCRIPTS\"
  14.106 +
  14.107 +set COMPAREQ = \"$COMPARE\"
  14.108 +set MODELN   = \"model1\"
  14.109 +
  14.110 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  14.111 + 
  14.112 +echo "`date` - Running diagnostics for model ${MODELQ}"
  14.113 +# CLAMP metric processing of model1
  14.114 +echo "`date` - Running 00.initial.ncl"
  14.115 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  14.116 +echo "`date` - Running 01.npp.ncl"
  14.117 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  14.118 +echo "`date` - Running 01.lai.ncl"
  14.119 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  14.120 +
  14.121 +if ($FILE3 != "") then
  14.122 +echo "`date` - Running 03.co2.ncl"
  14.123 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  14.124 +endif
  14.125 +
  14.126 +echo "`date` - Running 04.biomass.ncl"
  14.127 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  14.128 +echo "`date` - Running 06.fluxnet.ncl"
  14.129 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  14.130 +echo "`date` - Running 07.beta.ncl"
  14.131 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  14.132 +echo "`date` - Running 08.turnover.ncl"
  14.133 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  14.134 +
  14.135 +if ($BGC != "casa") then
  14.136 +echo "`date` - Running 09.carbon_sink.ncl"
  14.137 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  14.138 +else
  14.139 +echo "`date` - Running 09x.carbon_sink.ncl"
  14.140 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  14.141 +endif
  14.142 +
  14.143 +if ($BGC != "casa") then
  14.144 +echo "`date` - Running 10.fire.ncl"
  14.145 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  14.146 +endif
  14.147 +
  14.148 +echo "`date` - Running 11.ameriflux.ncl"
  14.149 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  14.150 +echo "`date` - Running 99.final.ncl"
  14.151 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  14.152 +
  14.153 +#*******************************************************
  14.154 +# user modification-(3)
  14.155 +
  14.156 +# model2
  14.157 +set MODEL2 = i01.43cn
  14.158 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.43cn_out/
  14.159 +set FILE1  = i01.43cn_out_ANN_climo.nc
  14.160 +set FILE2  = i01.43cn_out_MONS_climo.nc
  14.161 +set FILE3  = 
  14.162 +set FILE4  = $FILE1
  14.163 +set FILE5  = $FILE1
  14.164 +set FILE6  = $FILE1
  14.165 +set FILE7  = i01.43cn_Fire_C_2876-2900_monthly.nc
  14.166 +set FILE8  = i01.43cn_ameriflux_2876-2900_monthly.nc
  14.167 +set FILE9  = $FILE1
  14.168 +set FILE10 = $FILE2
  14.169 +set GRID   = 1.9
  14.170 +set BGC    = cn
  14.171 +set ENERGY = old
  14.172 +#*******************************************************
  14.173 +# create model2 directory by copying templates
  14.174 +if ($FILE3 != "") then
  14.175 +   set TEMPLATE1 = template_1-model
  14.176 +else
  14.177 +   set TEMPLATE1 = template_1-model_noCO2
  14.178 +endif
  14.179 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
  14.180 +rm -rf $MODEL2
  14.181 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
  14.182 +mv $TEMPLATE1 $MODEL2
  14.183 +
  14.184 +# add quote, to be usesd in INPUT_TEXT
  14.185 +set MODELQ  = \"$MODEL2\"
  14.186 +set DIRMQ   = \"$DIR_M\"
  14.187 +set F1  = \"$FILE1\"
  14.188 +set F2  = \"$FILE2\"
  14.189 +set F3  = \"$FILE3\"
  14.190 +set F4  = \"$FILE4\"
  14.191 +set F5  = \"$FILE5\"
  14.192 +set F6  = \"$FILE6\"
  14.193 +set F7  = \"$FILE7\"
  14.194 +set F8  = \"$FILE8\"
  14.195 +set F9  = \"$FILE9\"
  14.196 +set F10 = \"$FILE10\"
  14.197 +set GRIDQ   = \"$GRID\"
  14.198 +set BGCQ    = \"$BGC\"
  14.199 +set ENERGYQ = \"$ENERGY\"
  14.200 +
  14.201 +set COMPAREQ = \"$COMPARE\"
  14.202 +set MODELN   = \"model2\"
  14.203 +
  14.204 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  14.205 + 
  14.206 +echo "`date` - Running diagnostics for model ${MODELQ}"
  14.207 +# CLAMP metric processing of model2
  14.208 +echo "`date` - Running 00.inital.ncl"
  14.209 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  14.210 +echo "`date` - Running 01.npp.ncl"
  14.211 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  14.212 +echo "`date` - Running 02.lai.ncl"
  14.213 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  14.214 +
  14.215 +if ($FILE3 != "") then
  14.216 +echo "`date` - Running 03.co2.ncl"
  14.217 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  14.218 +endif
  14.219 +
  14.220 +echo "`date` - Running 04.biomass.ncl"
  14.221 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  14.222 +echo "`date` - Running 06.fluxnet.ncl"
  14.223 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  14.224 +echo "`date` - Running 07.beta.ncl"
  14.225 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  14.226 +echo "`date` - Running 08.turnover.ncl"
  14.227 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  14.228 +
  14.229 +if ($BGC != "casa") then
  14.230 +echo "`date` - Running 09.carbon_sink.ncl"
  14.231 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  14.232 +else
  14.233 +echo "`date` - Running 09x.carbon_sink.ncl"
  14.234 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  14.235 +endif
  14.236 +
  14.237 +if ($BGC != "casa") then
  14.238 +echo "`date` - Running 10.fire.ncl"
  14.239 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  14.240 +endif
  14.241 +
  14.242 +echo "`date` - Running 11.ameriflux.ncl"
  14.243 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  14.244 +echo "`date` - Running 99.final.ncl"
  14.245 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  14.246 +
  14.247 +# create a tar file from the final output
  14.248 +echo "`date` - tarring up final output into all.tar"
  14.249 +tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
  14.250 +

    15.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
    15.2 +++ b/all/run_2-model_38cn_39cn.csh	Thu Mar 26 14:02:21 2009 -0400
    15.3 @@ -0,0 +1,247 @@
    15.4 +#!/bin/csh -f
    15.5 +#-------------------------------------------------------------------
    15.6 +# note: user modifies ONLY the "user modification" section 
    15.7 +#
    15.8 +#       COMPARE: model1 vs model2
    15.9 +#       MODELn : model name
   15.10 +#       DIR_M  : directory of model data
   15.11 +#       DIR_O  : directory of observed data
   15.12 +#       DIR_S  : directory of model surface data
   15.13 +#       DIR_SCRIPTS  : directory of run scripts
   15.14 +#       FILE1  : time_mean  climatology from CLM diagnostic package
   15.15 +#       FILE2  : 12-monthly climatology from CLM diagnostic package
   15.16 +#       FILE3  : 12-monthly climatology from ATM diagnostic package
   15.17 +#                leave it blank, if no ATM file:
   15.18 +#                set FILE3 =
   15.19 +#       FILE7 : timeseries file generated from .... 
   15.20 +#       FILE8 : timeseries file generated from .... 
   15.21 +#       GRID   : T31, T42, or 1.9
   15.22 +#       BGC    : cn or casa 
   15.23 +#       ENERGY : new or old (fields in model data)
   15.24 +#-------------------------------------------------------------------
   15.25 +
   15.26 +#*******************************************************
   15.27 +# user modification-(1)
   15.28 +
   15.29 +# directory name of model comparison
   15.30 +#et COMPARE = b30.061n_vs_b30.061m
   15.31 +set COMPARE = i01.38cn_vs_i01.39cn
   15.32 +
   15.33 +#*******************************************************
   15.34 +# user modification-(2)
   15.35 +
   15.36 +# model1 
   15.37 +set MODEL1 = i01.38cn
   15.38 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.38cn_out/
   15.39 +set FILE1  = i01.38cn_out_ANN_climo.nc
   15.40 +set FILE2  = i01.38cn_out_MONS_climo.nc
   15.41 +set FILE3  = 
   15.42 +set FILE4  = $FILE1
   15.43 +set FILE5  = $FILE1
   15.44 +set FILE6  = $FILE1
   15.45 +set FILE7  = i01.38cn_Fire_C_2151-2175_monthly.nc
   15.46 +set FILE8  = i01.38cn_ameriflux_2151-2175_monthly.nc
   15.47 +set FILE9  = $FILE1
   15.48 +set FILE10 = $FILE2
   15.49 +set GRID   = 1.9
   15.50 +set BGC    = cn
   15.51 +set ENERGY = old
   15.52 +
   15.53 +# in the "CLAMP metric processing" section:
   15.54 +#  only 00.initial.ncl and 99.final.ncl are required,
   15.55 +#  user can comment out any one or more of the other ncl scripts,
   15.56 +#  e.g. 
   15.57 +#  #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
   15.58 +
   15.59 +# model surface data
   15.60 +set DIR_S  = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
   15.61 +
   15.62 +# observed data
   15.63 +set DIR_O  = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
   15.64 +
   15.65 +# directory for scripts, templates and ncl files
   15.66 +set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
   15.67 +
   15.68 +#********************************************************
   15.69 +
   15.70 +# create model1 and model1_vs_model2 directory by copying templates
   15.71 +if ($FILE3 != "") then
   15.72 +   set TEMPLATE1 = template_1-model
   15.73 +   set TEMPLATE2 = template_2-model
   15.74 +else
   15.75 +   set TEMPLATE1 = template_1-model_noCO2
   15.76 +   set TEMPLATE2 = template_2-model_noCO2
   15.77 +endif
   15.78 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
   15.79 +#cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
   15.80 +rm -rf $MODEL1
   15.81 +rm -rf $COMPARE
   15.82 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
   15.83 +mv $TEMPLATE1 $MODEL1
   15.84 +tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
   15.85 +mv $TEMPLATE2 $COMPARE
   15.86 +
   15.87 +# add quote, to be usesd in INPUT_TEXT
   15.88 +set MODELQ = \"$MODEL1\"
   15.89 +set DIRMQ  = \"$DIR_M\"
   15.90 +set F1  = \"$FILE1\"
   15.91 +set F2  = \"$FILE2\"
   15.92 +set F3  = \"$FILE3\"
   15.93 +set F4  = \"$FILE4\"
   15.94 +set F5  = \"$FILE5\"
   15.95 +set F6  = \"$FILE6\"
   15.96 +set F7  = \"$FILE7\"
   15.97 +set F8  = \"$FILE8\"
   15.98 +set F9  = \"$FILE9\"
   15.99 +set F10 = \"$FILE10\"
  15.100 +set GRIDQ   = \"$GRID\"
  15.101 +set BGCQ    = \"$BGC\"
  15.102 +set ENERGYQ = \"$ENERGY\"
  15.103 +set DIRSQ   = \"$DIR_S\"
  15.104 +set DIROQ   = \"$DIR_O\"
  15.105 +set DIRCQ   = \"$DIR_SCRIPTS\"
  15.106 +
  15.107 +set COMPAREQ = \"$COMPARE\"
  15.108 +set MODELN   = \"model1\"
  15.109 +
  15.110 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  15.111 + 
  15.112 +echo "`date` - Running diagnostics for model ${MODELQ}"
  15.113 +# CLAMP metric processing of model1
  15.114 +echo "`date` - Running 00.initial.ncl"
  15.115 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  15.116 +echo "`date` - Running 01.npp.ncl"
  15.117 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  15.118 +echo "`date` - Running 01.lai.ncl"
  15.119 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  15.120 +
  15.121 +if ($FILE3 != "") then
  15.122 +echo "`date` - Running 03.co2.ncl"
  15.123 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  15.124 +endif
  15.125 +
  15.126 +echo "`date` - Running 04.biomass.ncl"
  15.127 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  15.128 +echo "`date` - Running 06.fluxnet.ncl"
  15.129 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  15.130 +echo "`date` - Running 07.beta.ncl"
  15.131 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  15.132 +echo "`date` - Running 08.turnover.ncl"
  15.133 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  15.134 +
  15.135 +if ($BGC != "casa") then
  15.136 +echo "`date` - Running 09.carbon_sink.ncl"
  15.137 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  15.138 +else
  15.139 +echo "`date` - Running 09x.carbon_sink.ncl"
  15.140 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  15.141 +endif
  15.142 +
  15.143 +if ($BGC != "casa") then
  15.144 +echo "`date` - Running 10.fire.ncl"
  15.145 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  15.146 +endif
  15.147 +
  15.148 +echo "`date` - Running 11.ameriflux.ncl"
  15.149 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  15.150 +echo "`date` - Running 99.final.ncl"
  15.151 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  15.152 +
  15.153 +#*******************************************************
  15.154 +# user modification-(3)
  15.155 +
  15.156 +# model2
  15.157 +set MODEL2 = i01.39cn
  15.158 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.39cn_out/
  15.159 +set FILE1  = i01.39cn_out_ANN_climo.nc
  15.160 +set FILE2  = i01.39cn_out_MONS_climo.nc
  15.161 +set FILE3  = 
  15.162 +set FILE4  = $FILE1
  15.163 +set FILE5  = $FILE1
  15.164 +set FILE6  = $FILE1
  15.165 +set FILE7  = i01.39cn_Fire_C_2101-2125_monthly.nc
  15.166 +set FILE8  = i01.39cn_ameriflux_2101-2125_monthly.nc
  15.167 +set FILE9  = $FILE1
  15.168 +set FILE10 = $FILE2
  15.169 +set GRID   = 1.9
  15.170 +set BGC    = cn
  15.171 +set ENERGY = old
  15.172 +#*******************************************************
  15.173 +# create model2 directory by copying templates
  15.174 +if ($FILE3 != "") then
  15.175 +   set TEMPLATE1 = template_1-model
  15.176 +else
  15.177 +   set TEMPLATE1 = template_1-model_noCO2
  15.178 +endif
  15.179 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
  15.180 +rm -rf $MODEL2
  15.181 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
  15.182 +mv $TEMPLATE1 $MODEL2
  15.183 +
  15.184 +# add quote, to be usesd in INPUT_TEXT
  15.185 +set MODELQ  = \"$MODEL2\"
  15.186 +set DIRMQ   = \"$DIR_M\"
  15.187 +set F1  = \"$FILE1\"
  15.188 +set F2  = \"$FILE2\"
  15.189 +set F3  = \"$FILE3\"
  15.190 +set F4  = \"$FILE4\"
  15.191 +set F5  = \"$FILE5\"
  15.192 +set F6  = \"$FILE6\"
  15.193 +set F7  = \"$FILE7\"
  15.194 +set F8  = \"$FILE8\"
  15.195 +set F9  = \"$FILE9\"
  15.196 +set F10 = \"$FILE10\"
  15.197 +set GRIDQ   = \"$GRID\"
  15.198 +set BGCQ    = \"$BGC\"
  15.199 +set ENERGYQ = \"$ENERGY\"
  15.200 +
  15.201 +set COMPAREQ = \"$COMPARE\"
  15.202 +set MODELN   = \"model2\"
  15.203 +
  15.204 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  15.205 + 
  15.206 +echo "`date` - Running diagnostics for model ${MODELQ}"
  15.207 +# CLAMP metric processing of model2
  15.208 +echo "`date` - Running 00.inital.ncl"
  15.209 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  15.210 +echo "`date` - Running 01.npp.ncl"
  15.211 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  15.212 +echo "`date` - Running 02.lai.ncl"
  15.213 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  15.214 +
  15.215 +if ($FILE3 != "") then
  15.216 +echo "`date` - Running 03.co2.ncl"
  15.217 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  15.218 +endif
  15.219 +
  15.220 +echo "`date` - Running 04.biomass.ncl"
  15.221 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  15.222 +echo "`date` - Running 06.fluxnet.ncl"
  15.223 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  15.224 +echo "`date` - Running 07.beta.ncl"
  15.225 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  15.226 +echo "`date` - Running 08.turnover.ncl"
  15.227 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  15.228 +
  15.229 +if ($BGC != "casa") then
  15.230 +echo "`date` - Running 09.carbon_sink.ncl"
  15.231 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  15.232 +else
  15.233 +echo "`date` - Running 09x.carbon_sink.ncl"
  15.234 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  15.235 +endif
  15.236 +
  15.237 +if ($BGC != "casa") then
  15.238 +echo "`date` - Running 10.fire.ncl"
  15.239 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  15.240 +endif
  15.241 +
  15.242 +echo "`date` - Running 11.ameriflux.ncl"
  15.243 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  15.244 +echo "`date` - Running 99.final.ncl"
  15.245 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  15.246 +
  15.247 +# create a tar file from the final output
  15.248 +echo "`date` - tarring up final output into all.tar"
  15.249 +tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
  15.250 +

    16.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
    16.2 +++ b/all/run_2-model_38cn_41cn.csh	Thu Mar 26 14:02:21 2009 -0400
    16.3 @@ -0,0 +1,247 @@
    16.4 +#!/bin/csh -f
    16.5 +#-------------------------------------------------------------------
    16.6 +# note: user modifies ONLY the "user modification" section 
    16.7 +#
    16.8 +#       COMPARE: model1 vs model2
    16.9 +#       MODELn : model name
   16.10 +#       DIR_M  : directory of model data
   16.11 +#       DIR_O  : directory of observed data
   16.12 +#       DIR_S  : directory of model surface data
   16.13 +#       DIR_SCRIPTS  : directory of run scripts
   16.14 +#       FILE1  : time_mean  climatology from CLM diagnostic package
   16.15 +#       FILE2  : 12-monthly climatology from CLM diagnostic package
   16.16 +#       FILE3  : 12-monthly climatology from ATM diagnostic package
   16.17 +#                leave it blank, if no ATM file:
   16.18 +#                set FILE3 =
   16.19 +#       FILE7 : timeseries file generated from .... 
   16.20 +#       FILE8 : timeseries file generated from .... 
   16.21 +#       GRID   : T31, T42, or 1.9
   16.22 +#       BGC    : cn or casa 
   16.23 +#       ENERGY : new or old (fields in model data)
   16.24 +#-------------------------------------------------------------------
   16.25 +
   16.26 +#*******************************************************
   16.27 +# user modification-(1)
   16.28 +
   16.29 +# directory name of model comparison
   16.30 +#et COMPARE = b30.061n_vs_b30.061m
   16.31 +set COMPARE = i01.38cn_vs_i01.41cn
   16.32 +
   16.33 +#*******************************************************
   16.34 +# user modification-(2)
   16.35 +
   16.36 +# model1 
   16.37 +set MODEL1 = i01.38cn
   16.38 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.38cn_out/
   16.39 +set FILE1  = i01.38cn_out_ANN_climo.nc
   16.40 +set FILE2  = i01.38cn_out_MONS_climo.nc
   16.41 +set FILE3  = 
   16.42 +set FILE4  = $FILE1
   16.43 +set FILE5  = $FILE1
   16.44 +set FILE6  = $FILE1
   16.45 +set FILE7  = i01.38cn_Fire_C_2151-2175_monthly.nc
   16.46 +set FILE8  = i01.38cn_ameriflux_2151-2175_monthly.nc
   16.47 +set FILE9  = $FILE1
   16.48 +set FILE10 = $FILE2
   16.49 +set GRID   = 1.9
   16.50 +set BGC    = cn
   16.51 +set ENERGY = old
   16.52 +
   16.53 +# in the "CLAMP metric processing" section:
   16.54 +#  only 00.initial.ncl and 99.final.ncl are required,
   16.55 +#  user can comment out any one or more of the other ncl scripts,
   16.56 +#  e.g. 
   16.57 +#  #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
   16.58 +
   16.59 +# model surface data
   16.60 +set DIR_S  = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
   16.61 +
   16.62 +# observed data
   16.63 +set DIR_O  = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
   16.64 +
   16.65 +# directory for scripts, templates and ncl files
   16.66 +set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
   16.67 +
   16.68 +#********************************************************
   16.69 +
   16.70 +# create model1 and model1_vs_model2 directory by copying templates
   16.71 +if ($FILE3 != "") then
   16.72 +   set TEMPLATE1 = template_1-model
   16.73 +   set TEMPLATE2 = template_2-model
   16.74 +else
   16.75 +   set TEMPLATE1 = template_1-model_noCO2
   16.76 +   set TEMPLATE2 = template_2-model_noCO2
   16.77 +endif
   16.78 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
   16.79 +#cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
   16.80 +rm -rf $MODEL1
   16.81 +rm -rf $COMPARE
   16.82 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
   16.83 +mv $TEMPLATE1 $MODEL1
   16.84 +tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
   16.85 +mv $TEMPLATE2 $COMPARE
   16.86 +
   16.87 +# add quote, to be usesd in INPUT_TEXT
   16.88 +set MODELQ = \"$MODEL1\"
   16.89 +set DIRMQ  = \"$DIR_M\"
   16.90 +set F1  = \"$FILE1\"
   16.91 +set F2  = \"$FILE2\"
   16.92 +set F3  = \"$FILE3\"
   16.93 +set F4  = \"$FILE4\"
   16.94 +set F5  = \"$FILE5\"
   16.95 +set F6  = \"$FILE6\"
   16.96 +set F7  = \"$FILE7\"
   16.97 +set F8  = \"$FILE8\"
   16.98 +set F9  = \"$FILE9\"
   16.99 +set F10 = \"$FILE10\"
  16.100 +set GRIDQ   = \"$GRID\"
  16.101 +set BGCQ    = \"$BGC\"
  16.102 +set ENERGYQ = \"$ENERGY\"
  16.103 +set DIRSQ   = \"$DIR_S\"
  16.104 +set DIROQ   = \"$DIR_O\"
  16.105 +set DIRCQ   = \"$DIR_SCRIPTS\"
  16.106 +
  16.107 +set COMPAREQ = \"$COMPARE\"
  16.108 +set MODELN   = \"model1\"
  16.109 +
  16.110 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  16.111 + 
  16.112 +echo "`date` - Running diagnostics for model ${MODELQ}"
  16.113 +# CLAMP metric processing of model1
  16.114 +echo "`date` - Running 00.initial.ncl"
  16.115 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  16.116 +echo "`date` - Running 01.npp.ncl"
  16.117 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  16.118 +echo "`date` - Running 01.lai.ncl"
  16.119 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  16.120 +
  16.121 +if ($FILE3 != "") then
  16.122 +echo "`date` - Running 03.co2.ncl"
  16.123 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  16.124 +endif
  16.125 +
  16.126 +echo "`date` - Running 04.biomass.ncl"
  16.127 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  16.128 +echo "`date` - Running 06.fluxnet.ncl"
  16.129 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  16.130 +echo "`date` - Running 07.beta.ncl"
  16.131 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  16.132 +echo "`date` - Running 08.turnover.ncl"
  16.133 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  16.134 +
  16.135 +if ($BGC != "casa") then
  16.136 +echo "`date` - Running 09.carbon_sink.ncl"
  16.137 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  16.138 +else
  16.139 +echo "`date` - Running 09x.carbon_sink.ncl"
  16.140 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  16.141 +endif
  16.142 +
  16.143 +if ($BGC != "casa") then
  16.144 +echo "`date` - Running 10.fire.ncl"
  16.145 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  16.146 +endif
  16.147 +
  16.148 +echo "`date` - Running 11.ameriflux.ncl"
  16.149 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  16.150 +echo "`date` - Running 99.final.ncl"
  16.151 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  16.152 +
  16.153 +#*******************************************************
  16.154 +# user modification-(3)
  16.155 +
  16.156 +# model2
  16.157 +set MODEL2 = i01.41cn
  16.158 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.41cn_out/
  16.159 +set FILE1  = i01.41cn_out_ANN_climo.nc
  16.160 +set FILE2  = i01.41cn_out_MONS_climo.nc
  16.161 +set FILE3  = 
  16.162 +set FILE4  = $FILE1
  16.163 +set FILE5  = $FILE1
  16.164 +set FILE6  = $FILE1
  16.165 +set FILE7  = i01.41cn_Fire_C_2076-2100_monthly.nc
  16.166 +set FILE8  = i01.41cn_ameriflux_2076-2100_monthly.nc
  16.167 +set FILE9  = $FILE1
  16.168 +set FILE10 = $FILE2
  16.169 +set GRID   = 1.9
  16.170 +set BGC    = cn
  16.171 +set ENERGY = old
  16.172 +#*******************************************************
  16.173 +# create model2 directory by copying templates
  16.174 +if ($FILE3 != "") then
  16.175 +   set TEMPLATE1 = template_1-model
  16.176 +else
  16.177 +   set TEMPLATE1 = template_1-model_noCO2
  16.178 +endif
  16.179 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
  16.180 +rm -rf $MODEL2
  16.181 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
  16.182 +mv $TEMPLATE1 $MODEL2
  16.183 +
  16.184 +# add quote, to be usesd in INPUT_TEXT
  16.185 +set MODELQ  = \"$MODEL2\"
  16.186 +set DIRMQ   = \"$DIR_M\"
  16.187 +set F1  = \"$FILE1\"
  16.188 +set F2  = \"$FILE2\"
  16.189 +set F3  = \"$FILE3\"
  16.190 +set F4  = \"$FILE4\"
  16.191 +set F5  = \"$FILE5\"
  16.192 +set F6  = \"$FILE6\"
  16.193 +set F7  = \"$FILE7\"
  16.194 +set F8  = \"$FILE8\"
  16.195 +set F9  = \"$FILE9\"
  16.196 +set F10 = \"$FILE10\"
  16.197 +set GRIDQ   = \"$GRID\"
  16.198 +set BGCQ    = \"$BGC\"
  16.199 +set ENERGYQ = \"$ENERGY\"
  16.200 +
  16.201 +set COMPAREQ = \"$COMPARE\"
  16.202 +set MODELN   = \"model2\"
  16.203 +
  16.204 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  16.205 + 
  16.206 +echo "`date` - Running diagnostics for model ${MODELQ}"
  16.207 +# CLAMP metric processing of model2
  16.208 +echo "`date` - Running 00.inital.ncl"
  16.209 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  16.210 +echo "`date` - Running 01.npp.ncl"
  16.211 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  16.212 +echo "`date` - Running 02.lai.ncl"
  16.213 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  16.214 +
  16.215 +if ($FILE3 != "") then
  16.216 +echo "`date` - Running 03.co2.ncl"
  16.217 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  16.218 +endif
  16.219 +
  16.220 +echo "`date` - Running 04.biomass.ncl"
  16.221 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  16.222 +echo "`date` - Running 06.fluxnet.ncl"
  16.223 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  16.224 +echo "`date` - Running 07.beta.ncl"
  16.225 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  16.226 +echo "`date` - Running 08.turnover.ncl"
  16.227 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  16.228 +
  16.229 +if ($BGC != "casa") then
  16.230 +echo "`date` - Running 09.carbon_sink.ncl"
  16.231 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  16.232 +else
  16.233 +echo "`date` - Running 09x.carbon_sink.ncl"
  16.234 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  16.235 +endif
  16.236 +
  16.237 +if ($BGC != "casa") then
  16.238 +echo "`date` - Running 10.fire.ncl"
  16.239 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  16.240 +endif
  16.241 +
  16.242 +echo "`date` - Running 11.ameriflux.ncl"
  16.243 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  16.244 +echo "`date` - Running 99.final.ncl"
  16.245 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  16.246 +
  16.247 +# create a tar file from the final output
  16.248 +echo "`date` - tarring up final output into all.tar"
  16.249 +tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
  16.250 +

    17.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
    17.2 +++ b/all/run_2-model_39cn_40cn.csh	Thu Mar 26 14:02:21 2009 -0400
    17.3 @@ -0,0 +1,247 @@
    17.4 +#!/bin/csh -f
    17.5 +#-------------------------------------------------------------------
    17.6 +# note: user modifies ONLY the "user modification" section 
    17.7 +#
    17.8 +#       COMPARE: model1 vs model2
    17.9 +#       MODELn : model name
   17.10 +#       DIR_M  : directory of model data
   17.11 +#       DIR_O  : directory of observed data
   17.12 +#       DIR_S  : directory of model surface data
   17.13 +#       DIR_SCRIPTS  : directory of run scripts
   17.14 +#       FILE1  : time_mean  climatology from CLM diagnostic package
   17.15 +#       FILE2  : 12-monthly climatology from CLM diagnostic package
   17.16 +#       FILE3  : 12-monthly climatology from ATM diagnostic package
   17.17 +#                leave it blank, if no ATM file:
   17.18 +#                set FILE3 =
   17.19 +#       FILE7 : timeseries file generated from .... 
   17.20 +#       FILE8 : timeseries file generated from .... 
   17.21 +#       GRID   : T31, T42, or 1.9
   17.22 +#       BGC    : cn or casa 
   17.23 +#       ENERGY : new or old (fields in model data)
   17.24 +#-------------------------------------------------------------------
   17.25 +
   17.26 +#*******************************************************
   17.27 +# user modification-(1)
   17.28 +
   17.29 +# directory name of model comparison
   17.30 +#et COMPARE = b30.061n_vs_b30.061m
   17.31 +set COMPARE = i01.39cn_vs_i01.40cn
   17.32 +
   17.33 +#*******************************************************
   17.34 +# user modification-(2)
   17.35 +
   17.36 +# model1 
   17.37 +set MODEL1 = i01.39cn
   17.38 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.39cn_out/
   17.39 +set FILE1  = i01.39cn_out_ANN_climo.nc
   17.40 +set FILE2  = i01.39cn_out_MONS_climo.nc
   17.41 +set FILE3  = 
   17.42 +set FILE4  = $FILE1
   17.43 +set FILE5  = $FILE1
   17.44 +set FILE6  = $FILE1
   17.45 +set FILE7  = i01.39cn_Fire_C_2101-2125_monthly.nc
   17.46 +set FILE8  = i01.39cn_ameriflux_2101-2125_monthly.nc
   17.47 +set FILE9  = $FILE1
   17.48 +set FILE10 = $FILE2
   17.49 +set GRID   = 1.9
   17.50 +set BGC    = cn
   17.51 +set ENERGY = old
   17.52 +
   17.53 +# in the "CLAMP metric processing" section:
   17.54 +#  only 00.initial.ncl and 99.final.ncl are required,
   17.55 +#  user can comment out any one or more of the other ncl scripts,
   17.56 +#  e.g. 
   17.57 +#  #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
   17.58 +
   17.59 +# model surface data
   17.60 +set DIR_S  = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
   17.61 +
   17.62 +# observed data
   17.63 +set DIR_O  = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
   17.64 +
   17.65 +# directory for scripts, templates and ncl files
   17.66 +set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
   17.67 +
   17.68 +#********************************************************
   17.69 +
   17.70 +# create model1 and model1_vs_model2 directory by copying templates
   17.71 +if ($FILE3 != "") then
   17.72 +   set TEMPLATE1 = template_1-model
   17.73 +   set TEMPLATE2 = template_2-model
   17.74 +else
   17.75 +   set TEMPLATE1 = template_1-model_noCO2
   17.76 +   set TEMPLATE2 = template_2-model_noCO2
   17.77 +endif
   17.78 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
   17.79 +#cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
   17.80 +rm -rf $MODEL1
   17.81 +rm -rf $COMPARE
   17.82 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
   17.83 +mv $TEMPLATE1 $MODEL1
   17.84 +tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
   17.85 +mv $TEMPLATE2 $COMPARE
   17.86 +
   17.87 +# add quote, to be usesd in INPUT_TEXT
   17.88 +set MODELQ = \"$MODEL1\"
   17.89 +set DIRMQ  = \"$DIR_M\"
   17.90 +set F1  = \"$FILE1\"
   17.91 +set F2  = \"$FILE2\"
   17.92 +set F3  = \"$FILE3\"
   17.93 +set F4  = \"$FILE4\"
   17.94 +set F5  = \"$FILE5\"
   17.95 +set F6  = \"$FILE6\"
   17.96 +set F7  = \"$FILE7\"
   17.97 +set F8  = \"$FILE8\"
   17.98 +set F9  = \"$FILE9\"
   17.99 +set F10 = \"$FILE10\"
  17.100 +set GRIDQ   = \"$GRID\"
  17.101 +set BGCQ    = \"$BGC\"
  17.102 +set ENERGYQ = \"$ENERGY\"
  17.103 +set DIRSQ   = \"$DIR_S\"
  17.104 +set DIROQ   = \"$DIR_O\"
  17.105 +set DIRCQ   = \"$DIR_SCRIPTS\"
  17.106 +
  17.107 +set COMPAREQ = \"$COMPARE\"
  17.108 +set MODELN   = \"model1\"
  17.109 +
  17.110 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  17.111 + 
  17.112 +echo "`date` - Running diagnostics for model ${MODELQ}"
  17.113 +# CLAMP metric processing of model1
  17.114 +echo "`date` - Running 00.initial.ncl"
  17.115 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  17.116 +echo "`date` - Running 01.npp.ncl"
  17.117 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  17.118 +echo "`date` - Running 01.lai.ncl"
  17.119 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  17.120 +
  17.121 +if ($FILE3 != "") then
  17.122 +echo "`date` - Running 03.co2.ncl"
  17.123 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  17.124 +endif
  17.125 +
  17.126 +echo "`date` - Running 04.biomass.ncl"
  17.127 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  17.128 +echo "`date` - Running 06.fluxnet.ncl"
  17.129 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  17.130 +echo "`date` - Running 07.beta.ncl"
  17.131 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  17.132 +echo "`date` - Running 08.turnover.ncl"
  17.133 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  17.134 +
  17.135 +if ($BGC != "casa") then
  17.136 +echo "`date` - Running 09.carbon_sink.ncl"
  17.137 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  17.138 +else
  17.139 +echo "`date` - Running 09x.carbon_sink.ncl"
  17.140 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  17.141 +endif
  17.142 +
  17.143 +if ($BGC != "casa") then
  17.144 +echo "`date` - Running 10.fire.ncl"
  17.145 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  17.146 +endif
  17.147 +
  17.148 +echo "`date` - Running 11.ameriflux.ncl"
  17.149 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  17.150 +echo "`date` - Running 99.final.ncl"
  17.151 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  17.152 +
  17.153 +#*******************************************************
  17.154 +# user modification-(3)
  17.155 +
  17.156 +# model2
  17.157 +set MODEL2 = i01.40cn
  17.158 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.40cn_out/
  17.159 +set FILE1  = i01.40cn_out_ANN_climo.nc
  17.160 +set FILE2  = i01.40cn_out_MONS_climo.nc
  17.161 +set FILE3  = 
  17.162 +set FILE4  = $FILE1
  17.163 +set FILE5  = $FILE1
  17.164 +set FILE6  = $FILE1
  17.165 +set FILE7  = i01.40cn_Fire_C_2076-2100_monthly.nc
  17.166 +set FILE8  = i01.40cn_ameriflux_2076-2100_monthly.nc
  17.167 +set FILE9  = $FILE1
  17.168 +set FILE10 = $FILE2
  17.169 +set GRID   = 1.9
  17.170 +set BGC    = cn
  17.171 +set ENERGY = old
  17.172 +#*******************************************************
  17.173 +# create model2 directory by copying templates
  17.174 +if ($FILE3 != "") then
  17.175 +   set TEMPLATE1 = template_1-model
  17.176 +else
  17.177 +   set TEMPLATE1 = template_1-model_noCO2
  17.178 +endif
  17.179 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
  17.180 +rm -rf $MODEL2
  17.181 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
  17.182 +mv $TEMPLATE1 $MODEL2
  17.183 +
  17.184 +# add quote, to be usesd in INPUT_TEXT
  17.185 +set MODELQ  = \"$MODEL2\"
  17.186 +set DIRMQ   = \"$DIR_M\"
  17.187 +set F1  = \"$FILE1\"
  17.188 +set F2  = \"$FILE2\"
  17.189 +set F3  = \"$FILE3\"
  17.190 +set F4  = \"$FILE4\"
  17.191 +set F5  = \"$FILE5\"
  17.192 +set F6  = \"$FILE6\"
  17.193 +set F7  = \"$FILE7\"
  17.194 +set F8  = \"$FILE8\"
  17.195 +set F9  = \"$FILE9\"
  17.196 +set F10 = \"$FILE10\"
  17.197 +set GRIDQ   = \"$GRID\"
  17.198 +set BGCQ    = \"$BGC\"
  17.199 +set ENERGYQ = \"$ENERGY\"
  17.200 +
  17.201 +set COMPAREQ = \"$COMPARE\"
  17.202 +set MODELN   = \"model2\"
  17.203 +
  17.204 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  17.205 + 
  17.206 +echo "`date` - Running diagnostics for model ${MODELQ}"
  17.207 +# CLAMP metric processing of model2
  17.208 +echo "`date` - Running 00.inital.ncl"
  17.209 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  17.210 +echo "`date` - Running 01.npp.ncl"
  17.211 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  17.212 +echo "`date` - Running 02.lai.ncl"
  17.213 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  17.214 +
  17.215 +if ($FILE3 != "") then
  17.216 +echo "`date` - Running 03.co2.ncl"
  17.217 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  17.218 +endif
  17.219 +
  17.220 +echo "`date` - Running 04.biomass.ncl"
  17.221 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  17.222 +echo "`date` - Running 06.fluxnet.ncl"
  17.223 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  17.224 +echo "`date` - Running 07.beta.ncl"
  17.225 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  17.226 +echo "`date` - Running 08.turnover.ncl"
  17.227 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  17.228 +
  17.229 +if ($BGC != "casa") then
  17.230 +echo "`date` - Running 09.carbon_sink.ncl"
  17.231 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  17.232 +else
  17.233 +echo "`date` - Running 09x.carbon_sink.ncl"
  17.234 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  17.235 +endif
  17.236 +
  17.237 +if ($BGC != "casa") then
  17.238 +echo "`date` - Running 10.fire.ncl"
  17.239 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  17.240 +endif
  17.241 +
  17.242 +echo "`date` - Running 11.ameriflux.ncl"
  17.243 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  17.244 +echo "`date` - Running 99.final.ncl"
  17.245 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  17.246 +
  17.247 +# create a tar file from the final output
  17.248 +echo "`date` - tarring up final output into all.tar"
  17.249 +tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
  17.250 +

    18.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
    18.2 +++ b/all/run_2-model_42cn_43cn.csh	Thu Mar 26 14:02:21 2009 -0400
    18.3 @@ -0,0 +1,247 @@
    18.4 +#!/bin/csh -f
    18.5 +#-------------------------------------------------------------------
    18.6 +# note: user modifies ONLY the "user modification" section 
    18.7 +#
    18.8 +#       COMPARE: model1 vs model2
    18.9 +#       MODELn : model name
   18.10 +#       DIR_M  : directory of model data
   18.11 +#       DIR_O  : directory of observed data
   18.12 +#       DIR_S  : directory of model surface data
   18.13 +#       DIR_SCRIPTS  : directory of run scripts
   18.14 +#       FILE1  : time_mean  climatology from CLM diagnostic package
   18.15 +#       FILE2  : 12-monthly climatology from CLM diagnostic package
   18.16 +#       FILE3  : 12-monthly climatology from ATM diagnostic package
   18.17 +#                leave it blank, if no ATM file:
   18.18 +#                set FILE3 =
   18.19 +#       FILE7 : timeseries file generated from .... 
   18.20 +#       FILE8 : timeseries file generated from .... 
   18.21 +#       GRID   : T31, T42, or 1.9
   18.22 +#       BGC    : cn or casa 
   18.23 +#       ENERGY : new or old (fields in model data)
   18.24 +#-------------------------------------------------------------------
   18.25 +
   18.26 +#*******************************************************
   18.27 +# user modification-(1)
   18.28 +
   18.29 +# directory name of model comparison
   18.30 +#et COMPARE = b30.061n_vs_b30.061m
   18.31 +set COMPARE = i01.42cn_vs_i01.43cn
   18.32 +
   18.33 +#*******************************************************
   18.34 +# user modification-(2)
   18.35 +
   18.36 +# model1 
   18.37 +set MODEL1 = i01.42cn
   18.38 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases/i01.42cn_out/
   18.39 +set FILE1  = i01.42cn_out_ANN_climo.nc
   18.40 +set FILE2  = i01.42cn_out_MONS_climo.nc
   18.41 +set FILE3  = 
   18.42 +set FILE4  = $FILE1
   18.43 +set FILE5  = $FILE1
   18.44 +set FILE6  = $FILE1
   18.45 +set FILE7  = i01.42cn_Fire_C_2976-3000_monthly.nc
   18.46 +set FILE8  = i01.42cn_ameriflux_2976-3000_monthly.nc
   18.47 +set FILE9  = $FILE1
   18.48 +set FILE10 = $FILE2
   18.49 +set GRID   = 1.9
   18.50 +set BGC    = cn
   18.51 +set ENERGY = old
   18.52 +
   18.53 +# in the "CLAMP metric processing" section:
   18.54 +#  only 00.initial.ncl and 99.final.ncl are required,
   18.55 +#  user can comment out any one or more of the other ncl scripts,
   18.56 +#  e.g. 
   18.57 +#  #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
   18.58 +
   18.59 +# model surface data
   18.60 +set DIR_S  = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
   18.61 +
   18.62 +# observed data
   18.63 +set DIR_O  = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
   18.64 +
   18.65 +# directory for scripts, templates and ncl files
   18.66 +set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
   18.67 +
   18.68 +#********************************************************
   18.69 +
   18.70 +# create model1 and model1_vs_model2 directory by copying templates
   18.71 +if ($FILE3 != "") then
   18.72 +   set TEMPLATE1 = template_1-model
   18.73 +   set TEMPLATE2 = template_2-model
   18.74 +else
   18.75 +   set TEMPLATE1 = template_1-model_noCO2
   18.76 +   set TEMPLATE2 = template_2-model_noCO2
   18.77 +endif
   18.78 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
   18.79 +#cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
   18.80 +rm -rf $MODEL1
   18.81 +rm -rf $COMPARE
   18.82 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
   18.83 +mv $TEMPLATE1 $MODEL1
   18.84 +tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
   18.85 +mv $TEMPLATE2 $COMPARE
   18.86 +
   18.87 +# add quote, to be usesd in INPUT_TEXT
   18.88 +set MODELQ = \"$MODEL1\"
   18.89 +set DIRMQ  = \"$DIR_M\"
   18.90 +set F1  = \"$FILE1\"
   18.91 +set F2  = \"$FILE2\"
   18.92 +set F3  = \"$FILE3\"
   18.93 +set F4  = \"$FILE4\"
   18.94 +set F5  = \"$FILE5\"
   18.95 +set F6  = \"$FILE6\"
   18.96 +set F7  = \"$FILE7\"
   18.97 +set F8  = \"$FILE8\"
   18.98 +set F9  = \"$FILE9\"
   18.99 +set F10 = \"$FILE10\"
  18.100 +set GRIDQ   = \"$GRID\"
  18.101 +set BGCQ    = \"$BGC\"
  18.102 +set ENERGYQ = \"$ENERGY\"
  18.103 +set DIRSQ   = \"$DIR_S\"
  18.104 +set DIROQ   = \"$DIR_O\"
  18.105 +set DIRCQ   = \"$DIR_SCRIPTS\"
  18.106 +
  18.107 +set COMPAREQ = \"$COMPARE\"
  18.108 +set MODELN   = \"model1\"
  18.109 +
  18.110 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  18.111 + 
  18.112 +echo "`date` - Running diagnostics for model ${MODELQ}"
  18.113 +# CLAMP metric processing of model1
  18.114 +echo "`date` - Running 00.initial.ncl"
  18.115 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  18.116 +echo "`date` - Running 01.npp.ncl"
  18.117 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  18.118 +echo "`date` - Running 01.lai.ncl"
  18.119 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  18.120 +
  18.121 +if ($FILE3 != "") then
  18.122 +echo "`date` - Running 03.co2.ncl"
  18.123 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  18.124 +endif
  18.125 +
  18.126 +echo "`date` - Running 04.biomass.ncl"
  18.127 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  18.128 +echo "`date` - Running 06.fluxnet.ncl"
  18.129 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  18.130 +echo "`date` - Running 07.beta.ncl"
  18.131 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  18.132 +echo "`date` - Running 08.turnover.ncl"
  18.133 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  18.134 +
  18.135 +if ($BGC != "casa") then
  18.136 +echo "`date` - Running 09.carbon_sink.ncl"
  18.137 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  18.138 +else
  18.139 +echo "`date` - Running 09x.carbon_sink.ncl"
  18.140 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  18.141 +endif
  18.142 +
  18.143 +if ($BGC != "casa") then
  18.144 +echo "`date` - Running 10.fire.ncl"
  18.145 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  18.146 +endif
  18.147 +
  18.148 +echo "`date` - Running 11.ameriflux.ncl"
  18.149 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  18.150 +echo "`date` - Running 99.final.ncl"
  18.151 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  18.152 +
  18.153 +#*******************************************************
  18.154 +# user modification-(3)
  18.155 +
  18.156 +# model2
  18.157 +set MODEL2 = i01.43cn
  18.158 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/cases1/i01.43cn_out/
  18.159 +set FILE1  = i01.43cn_out_ANN_climo.nc
  18.160 +set FILE2  = i01.43cn_out_MONS_climo.nc
  18.161 +set FILE3  = 
  18.162 +set FILE4  = $FILE1
  18.163 +set FILE5  = $FILE1
  18.164 +set FILE6  = $FILE1
  18.165 +set FILE7  = i01.43cn_Fire_C_2876-2900_monthly.nc
  18.166 +set FILE8  = i01.43cn_ameriflux_2876-2900_monthly.nc
  18.167 +set FILE9  = $FILE1
  18.168 +set FILE10 = $FILE2
  18.169 +set GRID   = 1.9
  18.170 +set BGC    = cn
  18.171 +set ENERGY = old
  18.172 +#*******************************************************
  18.173 +# create model2 directory by copying templates
  18.174 +if ($FILE3 != "") then
  18.175 +   set TEMPLATE1 = template_1-model
  18.176 +else
  18.177 +   set TEMPLATE1 = template_1-model_noCO2
  18.178 +endif
  18.179 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
  18.180 +rm -rf $MODEL2
  18.181 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
  18.182 +mv $TEMPLATE1 $MODEL2
  18.183 +
  18.184 +# add quote, to be usesd in INPUT_TEXT
  18.185 +set MODELQ  = \"$MODEL2\"
  18.186 +set DIRMQ   = \"$DIR_M\"
  18.187 +set F1  = \"$FILE1\"
  18.188 +set F2  = \"$FILE2\"
  18.189 +set F3  = \"$FILE3\"
  18.190 +set F4  = \"$FILE4\"
  18.191 +set F5  = \"$FILE5\"
  18.192 +set F6  = \"$FILE6\"
  18.193 +set F7  = \"$FILE7\"
  18.194 +set F8  = \"$FILE8\"
  18.195 +set F9  = \"$FILE9\"
  18.196 +set F10 = \"$FILE10\"
  18.197 +set GRIDQ   = \"$GRID\"
  18.198 +set BGCQ    = \"$BGC\"
  18.199 +set ENERGYQ = \"$ENERGY\"
  18.200 +
  18.201 +set COMPAREQ = \"$COMPARE\"
  18.202 +set MODELN   = \"model2\"
  18.203 +
  18.204 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  18.205 + 
  18.206 +echo "`date` - Running diagnostics for model ${MODELQ}"
  18.207 +# CLAMP metric processing of model2
  18.208 +echo "`date` - Running 00.inital.ncl"
  18.209 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  18.210 +echo "`date` - Running 01.npp.ncl"
  18.211 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  18.212 +echo "`date` - Running 02.lai.ncl"
  18.213 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  18.214 +
  18.215 +if ($FILE3 != "") then
  18.216 +echo "`date` - Running 03.co2.ncl"
  18.217 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  18.218 +endif
  18.219 +
  18.220 +echo "`date` - Running 04.biomass.ncl"
  18.221 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  18.222 +echo "`date` - Running 06.fluxnet.ncl"
  18.223 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  18.224 +echo "`date` - Running 07.beta.ncl"
  18.225 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  18.226 +echo "`date` - Running 08.turnover.ncl"
  18.227 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  18.228 +
  18.229 +if ($BGC != "casa") then
  18.230 +echo "`date` - Running 09.carbon_sink.ncl"
  18.231 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  18.232 +else
  18.233 +echo "`date` - Running 09x.carbon_sink.ncl"
  18.234 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  18.235 +endif
  18.236 +
  18.237 +if ($BGC != "casa") then
  18.238 +echo "`date` - Running 10.fire.ncl"
  18.239 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  18.240 +endif
  18.241 +
  18.242 +echo "`date` - Running 11.ameriflux.ncl"
  18.243 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  18.244 +echo "`date` - Running 99.final.ncl"
  18.245 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  18.246 +
  18.247 +# create a tar file from the final output
  18.248 +echo "`date` - tarring up final output into all.tar"
  18.249 +tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
  18.250 +

    19.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
    19.2 +++ b/all/run_2-model_FMH.csh	Thu Mar 26 14:02:21 2009 -0400
    19.3 @@ -0,0 +1,254 @@
    19.4 +#!/bin/csh -f
    19.5 +#-------------------------------------------------------------------
    19.6 +# note: user modifies ONLY the "user modification" section 
    19.7 +#
    19.8 +#       COMPARE: model1 vs model2
    19.9 +#       MODELn : model name
   19.10 +#       DIR_M  : directory of model data
   19.11 +#       DIR_O  : directory of observed data
   19.12 +#       DIR_S  : directory of model surface data
   19.13 +#       DIR_SCRIPTS  : directory of run scripts
   19.14 +#       FILE1  : time_mean  climatology from CLM diagnostic package
   19.15 +#       FILE2  : 12-monthly climatology from CLM diagnostic package
   19.16 +#       FILE3  : 12-monthly climatology from ATM diagnostic package
   19.17 +#                leave it blank, if no ATM file:
   19.18 +#                set FILE3 =
   19.19 +#       FILE7 : timeseries file generated from .... 
   19.20 +#       FILE8 : timeseries file generated from .... 
   19.21 +#       GRID   : T31, T42, or 1.9
   19.22 +#       BGC    : cn or casa 
   19.23 +#       ENERGY : new or old (fields in model data)
   19.24 +#-------------------------------------------------------------------
   19.25 +
   19.26 +#*******************************************************
   19.27 +# user modification-(1)
   19.28 +
   19.29 +# directory name of model comparison
   19.30 +#et COMPARE = b30.061n_vs_b30.061m
   19.31 +set COMPARE = CN_vs_CASA 
   19.32 +
   19.33 +#*******************************************************
   19.34 +# user modification-(2)
   19.35 +
   19.36 +# model1 
   19.37 +set MODEL1 = CASA
   19.38 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/clamp_data/model/
   19.39 +#set FILE1  = i01.10casa_1948-2004_ANN_climo.nc
   19.40 +set FILE1  = i01.10casa_1975-2000_ANN_climo.nc
   19.41 +set FILE2  = i01.10casa_1948-2004_MONS_climo.nc
   19.42 +set FILE3  = f02.03casa_1876-1900_MONS_climo_atm.nc
   19.43 +set FILE4  = i01.06casa_1980-2004_ANN_climo.nc
   19.44 +set FILE5  = i01.10casa_1990-2004_ANN_climo.nc
   19.45 +set FILE6  = i01.07casa_1990-2004_ANN_climo.nc
   19.46 +set FILE7  = i01.10casa_Fire_C_1979-2004_monthly.nc
   19.47 +set FILE8  = i01.10casa_ameriflux_1990-2004_monthly.nc
   19.48 +set FILE9  = i01.10casa_2000-2004_ANN_climo.nc
   19.49 +set FILE10 = i01.10casa_2000-2004_MONS_climo.nc
   19.50 +set FILE11 = i01.10casa_2000_ANN_climo.nc
   19.51 +set GRID   = T42
   19.52 +set BGC    = casa
   19.53 +set ENERGY = new
   19.54 +
   19.55 +# in the "CLAMP metric processing" section:
   19.56 +#  only 00.initial.ncl and 99.final.ncl are required,
   19.57 +#  user can comment out any one or more of the other ncl scripts,
   19.58 +#  e.g. 
   19.59 +#  #ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
   19.60 +
   19.61 +# model surface data
   19.62 +set DIR_S  = /lustre/wolf-ddn/scratch/hof/clamp_data/surface_model/
   19.63 +
   19.64 +# observed data
   19.65 +set DIR_O  = /lustre/wolf-ddn/scratch/hof/clamp_data/observed/
   19.66 +
   19.67 +# directory for scripts, templates and ncl files
   19.68 +set DIR_SCRIPTS = /lustre/wolf-ddn/scratch/hof/clamp/all/
   19.69 +
   19.70 +#********************************************************
   19.71 +
   19.72 +# create model1 and model1_vs_model2 directory by copying templates
   19.73 +if ($FILE3 != "") then
   19.74 +   set TEMPLATE1 = template_1-model
   19.75 +   set TEMPLATE2 = template_2-model
   19.76 +else
   19.77 +   set TEMPLATE1 = template_1-model_noCO2
   19.78 +   set TEMPLATE2 = template_2-model_noCO2
   19.79 +endif
   19.80 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL1
   19.81 +#cp -r $DIR_SCRIPTS/$TEMPLATE2 $COMPARE
   19.82 +rm -rf $MODEL1
   19.83 +rm -rf $COMPARE
   19.84 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
   19.85 +mv $TEMPLATE1 $MODEL1
   19.86 +tar xvpf $DIR_SCRIPTS/$TEMPLATE2.tar
   19.87 +mv $TEMPLATE2 $COMPARE
   19.88 +
   19.89 +# add quote, to be usesd in INPUT_TEXT
   19.90 +set MODELQ = \"$MODEL1\"
   19.91 +set DIRMQ  = \"$DIR_M\"
   19.92 +set F1  = \"$FILE1\"
   19.93 +set F2  = \"$FILE2\"
   19.94 +set F3  = \"$FILE3\"
   19.95 +set F4  = \"$FILE4\"
   19.96 +set F5  = \"$FILE5\"
   19.97 +set F6  = \"$FILE6\"
   19.98 +set F7  = \"$FILE7\"
   19.99 +set F8  = \"$FILE8\"
  19.100 +set F9  = \"$FILE9\"
  19.101 +set F10 = \"$FILE10\"
  19.102 +set F11 = \"$FILE11\"
  19.103 +set GRIDQ   = \"$GRID\"
  19.104 +set BGCQ    = \"$BGC\"
  19.105 +set ENERGYQ = \"$ENERGY\"
  19.106 +set DIRSQ   = \"$DIR_S\"
  19.107 +set DIROQ   = \"$DIR_O\"
  19.108 +set DIRCQ   = \"$DIR_SCRIPTS\"
  19.109 +
  19.110 +set COMPAREQ = \"$COMPARE\"
  19.111 +set MODELN   = \"model1\"
  19.112 +
  19.113 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 film11=$F11 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  19.114 + 
  19.115 +echo "`date` - Running diagnostics for model ${MODELQ}"
  19.116 +# CLAMP metric processing of model1
  19.117 +echo "`date` - Running 00.initial.ncl"
  19.118 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  19.119 +echo "`date` - Running 01.npp.ncl"
  19.120 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  19.121 +echo "`date` - Running 01.lai.ncl"
  19.122 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  19.123 +
  19.124 +if ($FILE3 != "") then
  19.125 +echo "`date` - Running 03.co2.ncl"
  19.126 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  19.127 +endif
  19.128 +
  19.129 +echo "`date` - Running 04.biomass.ncl"
  19.130 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  19.131 +echo "`date` - Running 06.fluxnet.ncl"
  19.132 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  19.133 +echo "`date` - Running 07.beta.ncl"
  19.134 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  19.135 +echo "`date` - Running 08.turnover.ncl"
  19.136 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  19.137 +
  19.138 +if ($BGC != "casa") then
  19.139 +echo "`date` - Running 09.carbon_sink.ncl"
  19.140 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  19.141 +else
  19.142 +echo "`date` - Running 09x.carbon_sink.ncl"
  19.143 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  19.144 +endif
  19.145 +
  19.146 +if ($BGC != "casa") then
  19.147 +echo "`date` - Running 10.fire.ncl"
  19.148 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  19.149 +endif
  19.150 +
  19.151 +echo "`date` - Running 11.ameriflux.ncl"
  19.152 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  19.153 +echo "`date` - Running 99.final.ncl"
  19.154 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  19.155 +
  19.156 +#*******************************************************
  19.157 +# user modification-(3)
  19.158 +
  19.159 +# model2
  19.160 +set MODEL2 = CN
  19.161 +set DIR_M  = /lustre/wolf-ddn/scratch/hof/clamp_data/model/
  19.162 +#set FILE1  = i01.10cn_1948-2004_ANN_climo.nc
  19.163 +set FILE1  = i01.10cn_1975-2000_ANN_climo.nc
  19.164 +set FILE2  = i01.10cn_1948-2004_MONS_climo.nc
  19.165 +set FILE3  = f02.03cn_1901-1925_MONS_climo_atm.nc
  19.166 +set FILE4  = i01.06cn_1980-2004_ANN_climo.nc
  19.167 +set FILE5  = i01.10cn_1990-2004_ANN_climo.nc
  19.168 +set FILE6  = i01.07cn_1990-2004_ANN_climo.nc
  19.169 +set FILE7  = i01.10cn_Fire_C_1979-2004_monthly.nc
  19.170 +set FILE8  = i01.10cn_ameriflux_1990-2004_monthly.nc
  19.171 +set FILE9  = i01.10cn_2000-2004_ANN_climo.nc
  19.172 +set FILE10 = i01.10cn_2000-2004_MONS_climo.nc
  19.173 +set FILE11 = i01.10cn_2000_ANN_climo.nc
  19.174 +set GRID   = T42
  19.175 +set BGC    = cn
  19.176 +set ENERGY = new
  19.177 +
  19.178 +#*******************************************************
  19.179 +# create model2 directory by copying templates
  19.180 +if ($FILE3 != "") then
  19.181 +   set TEMPLATE1 = template_1-model
  19.182 +else
  19.183 +   set TEMPLATE1 = template_1-model_noCO2
  19.184 +endif
  19.185 +#cp -r $DIR_SCRIPTS/$TEMPLATE1 $MODEL2
  19.186 +rm -rf $MODEL2
  19.187 +tar xvpf $DIR_SCRIPTS/$TEMPLATE1.tar
  19.188 +mv $TEMPLATE1 $MODEL2
  19.189 +
  19.190 +# add quote, to be usesd in INPUT_TEXT
  19.191 +set MODELQ  = \"$MODEL2\"
  19.192 +set DIRMQ   = \"$DIR_M\"
  19.193 +set F1  = \"$FILE1\"
  19.194 +set F2  = \"$FILE2\"
  19.195 +set F3  = \"$FILE3\"
  19.196 +set F4  = \"$FILE4\"
  19.197 +set F5  = \"$FILE5\"
  19.198 +set F6  = \"$FILE6\"
  19.199 +set F7  = \"$FILE7\"
  19.200 +set F8  = \"$FILE8\"
  19.201 +set F9  = \"$FILE9\"
  19.202 +set F10 = \"$FILE10\"
  19.203 +set F11 = \"$FILE11\"
  19.204 +set GRIDQ   = \"$GRID\"
  19.205 +set BGCQ    = \"$BGC\"
  19.206 +set ENERGYQ = \"$ENERGY\"
  19.207 +
  19.208 +set COMPAREQ = \"$COMPARE\"
  19.209 +set MODELN   = \"model2\"
  19.210 +
  19.211 +set INPUT_TEXT = "model_name=$MODELQ model_grid=$GRIDQ dirm=$DIRMQ film1=$F1 film2=$F2 film3=$F3 film4=$F4 film5=$F5 film6=$F6 film7=$F7 film8=$F8 film9=$F9 film10=$F10 film11=$F11 BGC=$BGCQ ENERGY=$ENERGYQ dirs=$DIRSQ diro=$DIROQ dirscript=$DIRCQ modeln=$MODELN compare=$COMPAREQ"
  19.212 + 
  19.213 +echo "`date` - Running diagnostics for model ${MODELQ}"
  19.214 +# CLAMP metric processing of model2
  19.215 +echo "`date` - Running 00.inital.ncl"
  19.216 +ncl $INPUT_TEXT $DIR_SCRIPTS/00.initial.ncl
  19.217 +echo "`date` - Running 01.npp.ncl"
  19.218 +ncl $INPUT_TEXT $DIR_SCRIPTS/01.npp.ncl
  19.219 +echo "`date` - Running 02.lai.ncl"
  19.220 +ncl $INPUT_TEXT $DIR_SCRIPTS/02.lai.ncl
  19.221 +
  19.222 +if ($FILE3 != "") then
  19.223 +echo "`date` - Running 03.co2.ncl"
  19.224 +ncl $INPUT_TEXT $DIR_SCRIPTS/03.co2.ncl
  19.225 +endif
  19.226 +
  19.227 +echo "`date` - Running 04.biomass.ncl"
  19.228 +ncl $INPUT_TEXT $DIR_SCRIPTS/04.biomass.ncl
  19.229 +echo "`date` - Running 06.fluxnet.ncl"
  19.230 +ncl $INPUT_TEXT $DIR_SCRIPTS/06.fluxnet.ncl
  19.231 +echo "`date` - Running 07.beta.ncl"
  19.232 +ncl $INPUT_TEXT $DIR_SCRIPTS/07.beta.ncl
  19.233 +echo "`date` - Running 08.turnover.ncl"
  19.234 +ncl $INPUT_TEXT $DIR_SCRIPTS/08.turnover.ncl
  19.235 +
  19.236 +if ($BGC != "casa") then
  19.237 +echo "`date` - Running 09.carbon_sink.ncl"
  19.238 +ncl $INPUT_TEXT $DIR_SCRIPTS/09.carbon_sink.ncl
  19.239 +else
  19.240 +echo "`date` - Running 09x.carbon_sink.ncl"
  19.241 +ncl $INPUT_TEXT $DIR_SCRIPTS/09x.carbon_sink.ncl
  19.242 +endif
  19.243 +
  19.244 +if ($BGC != "casa") then
  19.245 +echo "`date` - Running 10.fire.ncl"
  19.246 +ncl $INPUT_TEXT $DIR_SCRIPTS/10.fire.ncl
  19.247 +endif
  19.248 +
  19.249 +echo "`date` - Running 11.ameriflux.ncl"
  19.250 +ncl $INPUT_TEXT $DIR_SCRIPTS/11.ameriflux.ncl
  19.251 +echo "`date` - Running 99.final.ncl"
  19.252 +ncl $INPUT_TEXT $DIR_SCRIPTS/99.final.ncl
  19.253 +
  19.254 +# create a tar file from the final output
  19.255 +echo "`date` - tarring up final output into all.tar"
  19.256 +tar cf - $MODEL1 $MODEL2 $COMPARE > all.tar
  19.257 +

    20.1 --- /dev/null	Thu Jan 01 00:00:00 1970 +0000
    20.2 +++ b/all/run_clamp_diag_LENS.csh	Thu Mar 26 14:02:21 2009 -0400
    20.3 @@ -0,0 +1,31 @@
    20.4 +#!/bin/csh -f
    20.5 +#PBS -A cli017bgc
    20.6 +#PBS -N clamp_diag
    20.7 +##PBS -q batch
    20.8 +#PBS -l nodes=1:ppn=1
    20.9 +#PBS -l walltime=02:00:00
   20.10 +#PBS -j oe
   20.11 +#PBS -S /bin/csh -V
   20.12 +#PBS -m abe
   20.13 +#PBS -M forrest@climatemodeling.org
   20.14 +##PBS -M forrest@climatemodeling.org,8656843670@messaging.sprintpcs.com
   20.15 +##PBS -M forrest@climatemodeling.org,fmhoffman@ipipi.com
   20.16 +
   20.17 +cd /tmp/work/hof/clamp/all/
   20.18 +
   20.19 +if ($?tcsh) then
   20.20 +	set modules_shell="tcsh"
   20.21 +else
   20.22 +	set modules_shell="csh"
   20.23 +endif
   20.24 +
   20.25 +source /usr/Modules/init/${modules_shell}
   20.26 +
   20.27 +unset modules_shell
   20.28 +
   20.29 +module load nco
   20.30 +module load ncl
   20.31 +
   20.32 +#time ./lnd_template4.6_ORNL_27cn29.csh
   20.33 +#time ./lnd_template4.6_ORNL_27cn31.csh >& lnd_template4.6_ORNL_27cn31.log
   20.34 +time ./run_2-model_FMH.csh >& run_2-model_FMH.log

    21.1 Binary file all/template_1-model.tar has changed

    22.1 --- a/all/template_1-model/index.html	Mon Jan 26 22:08:20 2009 -0500
    22.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    22.3 @@ -1,16 +0,0 @@
    22.4 -<html>
    22.5 -<head>
    22.6 -<title> Score Sheet for CLAMP </title>
    22.7 -<meta name="description" content="score CLAMP">
    22.8 -<meta name="keywords"    content="score BGC CLAMP CN CASA">
    22.9 -<meta name="Authors"     content="Jeff Lee">
   22.10 -</head>
   22.11 -
   22.12 -<!----------------------------------------------------------------------------->
   22.13 -<h2> &nbsp; &nbsp; <a href=http://www.climatemodeling.org/c-lamp/metrics4.pdf>CLAMP Score Sheet for Biogeochemical Model Evaluation</h2>
   22.14 -<!----------------------------------------------------------------------------->
   22.15 -
   22.16 -<!--#include virtual="table.html"-->
   22.17 -
   22.18 -<!----------------------------------------------------------------------------->
   22.19 -</html>

    23.1 --- a/all/template_1-model/table.html	Mon Jan 26 22:08:20 2009 -0500
    23.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    23.3 @@ -1,407 +0,0 @@
    23.4 -
    23.5 -<!----------------------------------------------------------------------------->
    23.6 -<table border=1 cellspacing=0 cellpadding=3 width=100%>
    23.7 -<tr>
    23.8 -   <th bgcolor=DDDDDD width=8% rowspan=2> Metric     </th>
    23.9 -   <th bgcolor=DDDDDD rowspan=2> Metric <br> components</th>
   23.10 -   <th bgcolor=DDDDDD rowspan=2> Observations & <br> comparison protocol</th>
   23.11 -   <th bgcolor=DDDDDD rowspan=2> Model<br> model_name  </th>
   23.12 -   <th bgcolor=DDDDDD colspan=2> Score (points) </th>
   23.13 -</tr>
   23.14 -<tr>
   23.15 -   <th width=6% bgcolor=DDDDDD> Full </th>
   23.16 -   <th width=6% bgcolor=DDDDDD> model_name </th>
   23.17 -</tr>
   23.18 - 
   23.19 -<tr>
   23.20 -   <th rowspan=6>NPP</th>
   23.21 -
   23.22 -   <th rowspan=2>EMDI NPP <br> observations</th>
   23.23 - 
   23.24 -   <th>         
   23.25 -   <a href=./npp/table_site81_ob.html>
   23.26 -   Class_A_table</a> 
   23.27 -   </th>
   23.28 - 
   23.29 -   <th>         
   23.30 -   <a href=./npp/table_site81_model_vs_ob.html>
   23.31 -   table</a> <br>
   23.32 -   <a href=./npp/scatter_model_vs_ob_81.png>
   23.33 -   scatter_plot</a>
   23.34 -   </th>
   23.35 - 
   23.36 -   <th>2.5</th>
   23.37 -   <th>M_npp_S81</th>
   23.38 -</tr>
   23.39 - 
   23.40 -<tr>
   23.41 -   <th>
   23.42 -   <a href=./npp/table_site933_ob.html>
   23.43 -   Class_B_table</a>     
   23.44 -   </th>
   23.45 -    
   23.46 -   <th>
   23.47 -   <a href=./npp/table_site933_model_vs_ob.html>
   23.48 -   table</a> <br>       
   23.49 -   <a href=./npp/scatter_model_vs_ob_933.png>
   23.50 -   scatter_plot</a>
   23.51 -   </th>
   23.52 -
   23.53 -   <th>2.5</th> 
   23.54 -   <th>M_npp_S933</th>
   23.55 -</tr>
   23.56 - 
   23.57 -<tr>
   23.58 -   <th rowspan=2>EMDI NPP<br>normalized by PPT</th>
   23.59 - 
   23.60 -   <th>        
   23.61 -   <a href=./npp/histogram_ob_81.png>
   23.62 -   Class_A_histogram</a>
   23.63 -   </th>
   23.64 - 
   23.65 -   <th>
   23.66 -   <a href=./npp/histogram_model_vs_ob_81.png>
   23.67 -   Class_A_histogram</a>
   23.68 -   </th>
   23.69 - 
   23.70 -   <th>2.5</th>
   23.71 -   <th>M_npp_H81</th>
   23.72 -</tr>
   23.73 - 
   23.74 -<tr>
   23.75 -   <th>
   23.76 -   <a href=./npp/histogram_ob_933.png>
   23.77 -   Class_B_histogram</a>
   23.78 -   </th>
   23.79 - 
   23.80 -   <th>
   23.81 -   <a href=./npp/histogram_model_vs_ob_933.png>
   23.82 -   Class_B_histogram</a>
   23.83 -   </th>
   23.84 -
   23.85 -   <th>2.5</th>
   23.86 -   <th>M_npp_H933</th>
   23.87 -</tr>
   23.88 - 
   23.89 -<tr>
   23.90 -   <th>Correlation with MODIS obs.</th>
   23.91 - 
   23.92 -   <th>
   23.93 -   <a href=./npp/global_ob.png>
   23.94 -   global_map</a>
   23.95 -   </th>
   23.96 - 
   23.97 -   <th>
   23.98 -   <a href=./npp/global_model.png>
   23.99 -   model_map</a> <br>
  23.100 -   <a href=./npp/global_model_vs_ob.png>
  23.101 -   model_vs_obs</a>
  23.102 -   </th>
  23.103 - 
  23.104 -   <th>5</th>
  23.105 -   <th>M_npp_G</th>
  23.106 -</tr>
  23.107 - 
  23.108 -<tr>
  23.109 -   <th>Correlation with MODIS-zonal mean</th>
  23.110 - 
  23.111 -   <th>
  23.112 -   <a href=./npp/zonal_ob.png>
  23.113 -   zonal_mean_obs</a>
  23.114 -   </th>
  23.115 - 
  23.116 -   <th>
  23.117 -   <a href=./npp/zonal_model_vs_ob.png>
  23.118 -   zonal_mean <br>model_vs_obs_plot</a>
  23.119 -   </th>
  23.120 - 
  23.121 -   <th>5</th>
  23.122 -   <th>M_npp_Z</th>
  23.123 -</tr>
  23.124 - 
  23.125 -<tr>
  23.126 -   <th rowspan=3>LAI</th>
  23.127 - 
  23.128 -   <th>MODIS Mean</th>
  23.129 - 
  23.130 -   <th>        
  23.131 -   <a href=./lai/global_class_ob.png>
  23.132 -   land_class_obs</a> <br>
  23.133 -   <a href=./lai/global_class_model.png>
  23.134 -   land_class_model</a> <br>       
  23.135 -   <a href=./lai/global_Mean_ob.png>
  23.136 -   global_map</a>
  23.137 -   </th>
  23.138 - 
  23.139 -   <th>
  23.140 -   <a href=./lai/table_model_vs_ob.html>
  23.141 -   model_vs_obs_table</a> <br>
  23.142 -   <a href=./lai/global_Mean_model.png>
  23.143 -   global_map</a> <br>
  23.144 -   <a href=./lai/global_Mean_model_vs_ob.png>
  23.145 -   model_vs_obs</a>
  23.146 -   </th>
  23.147 - 
  23.148 -   <th>5</th>
  23.149 -   <th>M_lai_Mean</th>
  23.150 -</tr>
  23.151 - 
  23.152 -<tr>
  23.153 -   <th>MODIS Maximum</th>
  23.154 - 
  23.155 -   <th>        
  23.156 -   <a href=./lai/global_Max_ob.png>
  23.157 -   global_map</a>
  23.158 -   </th>
  23.159 - 
  23.160 -   <th>       
  23.161 -   <a href=./lai/global_Max_model.png>
  23.162 -   global_map</a> <br>
  23.163 -   <a href=./lai/global_Max_model_vs_ob.png>
  23.164 -   model_vs_obs</a>
  23.165 -   </th>
  23.166 - 
  23.167 -   <th>5</th>
  23.168 -   <th>M_lai_Max</th>
  23.169 -</tr>
  23.170 -<tr>
  23.171 -   <th>MODIS Phase</th>
  23.172 - 
  23.173 -   <th>        
  23.174 -   <a href=./lai/global_Phase_ob.png>
  23.175 -   global_map</a>
  23.176 -   </th>
  23.177 - 
  23.178 -   <th>       
  23.179 -   <a href=./lai/global_Phase_model.png>
  23.180 -   global_map</a> <br>
  23.181 -   <a href=./lai/global_Phase_model_vs_ob.png>
  23.182 -   model_vs_obs</a>
  23.183 -   </th>
  23.184 - 
  23.185 -   <th>5</th>
  23.186 -   <th>M_lai_Phase</th>
  23.187 -</tr>
  23.188 -
  23.189 -<tr>
  23.190 -   <th>CO2 <br> Seasonal <br> Cycle</th>
  23.191 -   <th>Comparison with NOAA observations <br> phase and amplitude</th>
  23.192 -
  23.193 -   <th>--</th>
  23.194 -
  23.195 -   <th>
  23.196 -   <a href=./co2/score+line_vs_ob.html>
  23.197 -   model_vs_obs</a>
  23.198 -   </th>
  23.199 -
  23.200 -   <th >20</th>
  23.201 -   <th >M_co2</th>
  23.202 -</tr>
  23.203 - 
  23.204 -<tr>
  23.205 -   <th rowspan=2>Carbon<br> Stocks</th>
  23.206 -
  23.207 -   <th>Aboveground live biomass <br> in South America</th>
  23.208 - 
  23.209 -   <th>        
  23.210 -   <a href=./biomass/global_ob.png>
  23.211 -   obs_amazon</a>
  23.212 -   </th>
  23.213 - 
  23.214 -   <th>        
  23.215 -   <a href=./biomass/global_model.png>
  23.216 -   model_amazon</a> <br>
  23.217 -   <a href=./biomass/global_model_vs_ob.png>
  23.218 -   model_vs_obs</a>
  23.219 -   </th>
  23.220 -
  23.221 -   <th>5</th> 
  23.222 -   <th>M_biomass</th>
  23.223 -</tr>
  23.224 - 
  23.225 -<tr>
  23.226 -
  23.227 -   <th>Aboveground live <br> biomass <br> within Amazon Basin <br>(sum within Legal Amazon)</th>
  23.228 -
  23.229 -   <th>
  23.230 -   <a href=./biomass/mask_ob.png>        
  23.231 -   mask</a> <br>
  23.232 -   <a href=./biomass/global_mask_ob.png>
  23.233 -   obs_masked</a> <br>
  23.234 -   Sum_biomass_ob (Pg C)
  23.235 -   </th>
  23.236 -
  23.237 -   <th>
  23.238 -   <a href=./biomass/global_mask_model.png>
  23.239 -   model_masked</a> <br>
  23.240 -   <a href=./biomass/global_mask_vs_ob.png>
  23.241 -   model_vs_obs</a> <br>       
  23.242 -   Sum_biomass_mod (Pg C)
  23.243 -   </th>
  23.244 -
  23.245 -   <th>5</th>
  23.246 -   <th>M_biomask</th>
  23.247 -</tr>
  23.248 - 
  23.249 -<tr>
  23.250 -   <th rowspan=4>Energy and<br>C Fluxes<br>from Fluxnet</th>
  23.251 - 
  23.252 -   <th>NEE</th>
  23.253 -
  23.254 -   <th rowspan=4>     
  23.255 -   <a href=./fluxnet/line_ob.html> 
  23.256 -    line_plot</a>
  23.257 -   </th>
  23.258 -
  23.259 -   <th rowspan=4>     
  23.260 -   <a href=./fluxnet/score+line_vs_ob.html> 
  23.261 -    model_vs_obs</a>
  23.262 -   </th>
  23.263 - 
  23.264 -   <th>5</th>
  23.265 -   <th>M_fluxnet_nee</th>
  23.266 -</tr>
  23.267 - 
  23.268 -<tr>
  23.269 -   <th>Net radiation</th>
  23.270 -   <th>5</th>
  23.271 -   <th>M_fluxnet_rad</th>
  23.272 -</tr>
  23.273 - 
  23.274 -<tr>
  23.275 -   <th>Latent heat</th>
  23.276 -   <th>5</th>
  23.277 -   <th>M_fluxnet_lh</th>
  23.278 -</tr>
  23.279 -  
  23.280 -<tr>
  23.281 -   <th>Sensible heat</th>
  23.282 -   <th>5</th>
  23.283 -   <th>M_fluxnet_sh</th>
  23.284 -</tr>
  23.285 -
  23.286 -<tr>
  23.287 -   <th rowspan=6>Energy and<br>C Fluxes<br>from Ameriflux</th>
  23.288 - 
  23.289 -   <th>NEE</th>
  23.290 -
  23.291 -   <th rowspan=6>     
  23.292 -   <a href=./ameriflux/line_ob.html> 
  23.293 -    line_plot</a>
  23.294 -   </th>
  23.295 -
  23.296 -   <th rowspan=6>
  23.297 -  <a href=./ameriflux/score+line_vs_ob.html> 
  23.298 -    model_vs_obs</a> <br>
  23.299 -   <a href=./ameriflux/tseries_vs_ob.html> 
  23.300 -    timeseries_plot</a>
  23.301 -   </th>
  23.302 - 
  23.303 -   <th>1</th>
  23.304 -   <th>M_ameriflux_nee</th>
  23.305 -</tr>
  23.306 - 
  23.307 -<tr>
  23.308 -   <th>Shortwave Incoming</th>
  23.309 -   <th>1</th>
  23.310 -   <th>M_ameriflux_rad</th>
  23.311 -</tr>
  23.312 - 
  23.313 -<tr>
  23.314 -   <th>Latent heat</th>
  23.315 -   <th>1</th>
  23.316 -   <th>M_ameriflux_lh</th>
  23.317 -</tr>
  23.318 -  
  23.319 -<tr>
  23.320 -   <th>Sensible heat</th>
  23.321 -   <th>1</th>
  23.322 -   <th>M_ameriflux_sh</th>
  23.323 -</tr>
  23.324 -
  23.325 -<tr>
  23.326 -   <th>GPP</th>
  23.327 -   <th>1</th>
  23.328 -   <th>M_ameriflux_gpp</th>
  23.329 -</tr>
  23.330 -  
  23.331 -<tr>
  23.332 -   <th>ER</th>
  23.333 -   <th>1</th>
  23.334 -   <th>M_ameriflux_er</th>
  23.335 -</tr>
  23.336 - 
  23.337 -<tr>
  23.338 -   <th rowspan=4>Transient Dynamics</th>
  23.339 -   <th>Beta factor for NPP Stimulation from elevated CO2</th>
  23.340 -   <th>--</th>
  23.341 -   <th>     
  23.342 -   <a href=./beta/table_station.html> 
  23.343 -    FACE_Site_comparison</a> <br>
  23.344 -   <a href=./beta/table_biome.html> 
  23.345 -    biome_table</a>
  23.346 -   </th>
  23.347 -   <th>3</th>
  23.348 -   <th>M_beta</th>
  23.349 -</tr>
  23.350 - 
  23.351 -<tr>
  23.352 -   <th>Turnover times and pool sizes</th>
  23.353 -   <th>--</th>
  23.354 -   <th>     
  23.355 -   <a href=./turnover/table_Leaf.html> 
  23.356 -    Leaf</a> <br>
  23.357 -   <a href=./turnover/table_Wood.html> 
  23.358 -    Wood</a> <br>
  23.359 -   <a href=./turnover/table_Fine_Root.html>
  23.360 -    Fine_Root</a> <br>
  23.361 -   <a href=./turnover/table_Litter.html> 
  23.362 -    Litter</a> <br>
  23.363 -   <a href=./turnover/table_Coarse_Woody_Debris.html> 
  23.364 -    Coarse_Woody_Debris</a> <br>
  23.365 -   <a href=./turnover/table_Soil.html> 
  23.366 -    Soil</a>
  23.367 -   </th>
  23.368 -   <th>3</th>
  23.369 -   <th>--</th>
  23.370 -</tr>
  23.371 - 
  23.372 -<tr>
  23.373 -   <th>Carbon Sinks<br>(1990-2004)</th>
  23.374 -   <th>--</th>
  23.375 -
  23.376 -   <th>     
  23.377 -   <a href=./carbon_sink/table_per_m2.html> 
  23.378 -    biome_mean</a> <br>
  23.379 -   <a href=./carbon_sink/table_per_biome.html> 
  23.380 -    biome_total</a>
  23.381 -   </th>
  23.382 -
  23.383 -   <th>2</th>
  23.384 -   <th>--</th>
  23.385 -</tr>
  23.386 - 
  23.387 -<tr>
  23.388 -   <th>Fire Variability</th>
  23.389 -   <th>--</th>
  23.390 -
  23.391 -   <th>     
  23.392 -   <a href=./fire/global_model_vs_ob.png> 
  23.393 -    global_spatial_comparison</a> <br>
  23.394 -   <a href=./fire/table_fire.html> 
  23.395 -    temporal_dynamics_(1997-2004)</a>
  23.396 -   </th>
  23.397 -
  23.398 -   <th>2</th>
  23.399 -   <th>M_fire</th>
  23.400 -</tr>
  23.401 -
  23.402 -<tr>
  23.403 -   <th bgcolor=DDDDDD> Total </th>
  23.404 -   <th colspan=3>  </th>
  23.405 -   <th bgcolor=DDDDDD>100</th>
  23.406 -   <th bgcolor=DDDDDD>M_total</th>
  23.407 -</tr>
  23.408 - 
  23.409 -</table>
  23.410 -<!----------------------------------------------------------------------------->