forrest@0: ;************************************************************
forrest@0: ; required command line input parameters:
forrest@0: ; ncl 'model_name="10cn" model_grid="T42" dirm="/.../ film="..."' 01.npp.ncl
forrest@0: ;
forrest@0: ; using gsn_table for all
forrest@0: ; output: line plot for each site (4 fields)
forrest@0: ; table for M_score
forrest@0: ;************************************************************
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
forrest@0: load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
forrest@0: ;************************************************************
forrest@0: procedure set_line(lines:string,nline:integer,newlines:string)
forrest@0: begin
forrest@0: ; add line to ascci/html file
forrest@0:
forrest@0: nnewlines = dimsizes(newlines)
forrest@0: if(nline+nnewlines-1.ge.dimsizes(lines))
forrest@0: print("set_line: bad index, not setting anything.")
forrest@0: return
forrest@0: end if
forrest@0: lines(nline:nline+nnewlines-1) = newlines
forrest@0: ; print ("lines = " + lines(nline:nline+nnewlines-1))
forrest@0: nline = nline + nnewlines
forrest@0: return
forrest@0: end
forrest@0: ;*************************************************************
forrest@0: begin
forrest@0:
forrest@0: plot_type = "ps"
forrest@0: plot_type_new = "png"
forrest@0:
forrest@0: ;************************************************
forrest@0: ; read model data
forrest@0: ;************************************************
forrest@0:
forrest@0: ; from command line inputs
forrest@0:
forrest@0: ;--------------------------------------------------
forrest@0: ; edit table.html of current model for movel1_vs_model2
forrest@0:
forrest@0: if (isvar("compare")) then
forrest@0: html_name2 = compare+"/table.html"
forrest@0: html_new2 = html_name2 +".new"
forrest@0:
forrest@0: ; the following only needs to execute once
forrest@0: ; system("sed 1,/model_nameA/s//"+model_name+"/ "+html_name2+" > "+html_new2+";"+ \
forrest@0: ; "mv -f "+html_new2+" "+html_name2+";"+ \
forrest@0: ; "sed 1,/model_nameB/s//"+model_name+"/ "+html_name2+" > "+html_new2+";"+ \
forrest@0: ; "mv -f "+html_new2+" "+html_name2+";"+ \
forrest@0: ; "sed s#"+modeln+"#"+model_name+"# "+html_name2+" > "+html_new2+";"+ \
forrest@0: ; "mv -f "+html_new2+" "+html_name2)
forrest@0: end if
forrest@0:
forrest@0: ;-------------------------------------
forrest@0: ; edit table.html for current model
forrest@0:
forrest@0: html_name = model_name+"/table.html"
forrest@0: html_new = html_name +".new"
forrest@0:
forrest@0: ; the following only needs to execute once
forrest@0: ; system("sed s#model_name#"+model_name+"# "+html_name+" > "+html_new+";"+ \
forrest@0: ; "mv -f "+html_new+" "+html_name)
forrest@0: ;------------------------------------------------
forrest@0: fm = addfile(dirm+film2,"r")
forrest@0:
forrest@0: xm = fm->lon
forrest@0: ym = fm->lat
forrest@0: ;------------------------------------------------
forrest@0: nlat = dimsizes(ym)
forrest@0: nlon = dimsizes(xm)
forrest@0:
forrest@0: ; for 4 fields, 12-monthly
forrest@0: nmon = 12
forrest@0: nfield = 4
forrest@0:
forrest@0: data_mod0 = new ((/nfield,nmon,nlat,nlon/),float)
forrest@0:
forrest@0: ; change to unit of observed (u mol/m2/s)
forrest@0: ; Model_units [=] gC/m2/s
forrest@0: ; 12. = molecular weight of C
forrest@0: ; u mol = 1e-6 mol
forrest@0: factor = 1e6 /12.
forrest@0:
forrest@0: if (ENERGY .eq. "old") then
forrest@0:
forrest@0: data = fm->NEE
forrest@0: data_mod0(0,:,:,:) = data(:,:,:) * factor
forrest@0: delete (data)
forrest@0:
forrest@0: ; data = fm->LATENT
forrest@0: data1 = fm->FCEV
forrest@0: data2 = fm->FCTR
forrest@0: data3 = fm->FGEV
forrest@0: data_mod0(2,:,:,:) = data1(:,:,:)+data2(:,:,:)+data3(:,:,:)
forrest@0: delete (data1)
forrest@0: delete (data2)
forrest@0: delete (data3)
forrest@0:
forrest@0: ; data = fm->SENSIBLE
forrest@0: data = fm->FSH
forrest@0: data_mod0(3,:,:,:) = data(:,:,:)
forrest@0: delete (data)
forrest@0:
forrest@0: ; data = fm->NETRAD
forrest@0: data1 = fm->FSA
forrest@0: data2 = fm->FIRA
forrest@0: data_mod0(1,:,:,:) = data1(:,:,:)-data2(:,:,:)-data_mod0(2,:,:,:)-data_mod0(3,:,:,:)
forrest@0: delete (data1)
forrest@0: delete (data2)
forrest@0:
forrest@0: else
forrest@0:
forrest@0: data = fm->NEE
forrest@0: data_mod0(0,:,:,:) = data(:,:,:) * factor
forrest@0: delete (data)
forrest@0:
forrest@0: data = fm->NETRAD
forrest@0: data_mod0(1,:,:,:) = data(:,:,:)
forrest@0: delete (data)
forrest@0:
forrest@0: data = fm->LATENT
forrest@0: data_mod0(2,:,:,:) = data(:,:,:)
forrest@0: delete (data)
forrest@0:
forrest@0: data = fm->SENSIBLE
forrest@0: data_mod0(3,:,:,:) = data(:,:,:)
forrest@0: delete (data)
forrest@0: end if
forrest@0:
forrest@0: ;************************************************
forrest@0: ; read data: observed
forrest@0: ;************************************************
forrest@0:
forrest@0: station = (/"BOREAS_NSA_OBS" \
forrest@0: ,"CastelPorziano" \
forrest@0: ,"Hyytiala" \
forrest@0: ,"Kaamanen" \
forrest@0: ,"LBA_Tapajos_KM67" \
forrest@0: ,"Lethbridge" \
forrest@0: ,"Tharandt" \
forrest@0: ,"Vielsalm" \
forrest@0: /)
forrest@0:
forrest@0: year_ob = (/"1994-2004" \
forrest@0: ,"1997-2003" \
forrest@0: ,"1996-2003" \
forrest@0: ,"2000-2003" \
forrest@0: ,"2002-2005" \
forrest@0: ,"1999-2004" \
forrest@0: ,"1996-2003" \
forrest@0: ,"1998-2003" \
forrest@0: /)
forrest@0:
forrest@0: field = (/"CO2 Flux" \
forrest@0: ,"Net Radiation" \
forrest@0: ,"Latent Heat" \
forrest@0: ,"Sensible Heat" \
forrest@0: /)
forrest@0:
forrest@0: nstation = dimsizes(station)
forrest@0: nmon = 12
forrest@0: nfield = dimsizes(field)
forrest@0:
forrest@0: data_ob = new ((/nstation, nfield, nmon/),float)
forrest@0: lat_ob = new ((/nstation/),float)
forrest@0: lon_ob = new ((/nstation/),float)
forrest@0:
forrest@0: diri_root = "/fis/cgd/cseg/people/jeff/clamp_data/fluxnet/"
forrest@0:
forrest@0: do n = 0,nstation-1
forrest@0: diri = diri_root + station(n)+"/"
forrest@0: fili = station(n)+"_"+year_ob(n)+"_monthly.nc"
forrest@0: g = addfile (diri+fili,"r")
forrest@0:
forrest@0: lon_ob(n) = g->lon
forrest@0: lat_ob(n) = g->lat
forrest@0:
forrest@0: data = g->CO2_FLUX
forrest@0: data_ob(n,0,:) = dim_avg(data(month|:,year|:))
forrest@0: delete (data)
forrest@0:
forrest@0: data = g->RAD_FLUX
forrest@0: data_ob(n,1,:) = dim_avg(data(month|:,year|:))
forrest@0: delete (data)
forrest@0:
forrest@0: data = g->LH_FLUX
forrest@0: data_ob(n,2,:) = dim_avg(data(month|:,year|:))
forrest@0: delete (data)
forrest@0:
forrest@0: data = g->SH_FLUX
forrest@0: data_ob(n,3,:) = dim_avg(data(month|:,year|:))
forrest@0: delete (data)
forrest@0:
forrest@0: delete (g)
forrest@0: end do
forrest@0:
forrest@0: ;************************************************************
forrest@0: ; interpolate model data into observed station
forrest@0: ; note: model is 0-360E, 90S-90N
forrest@0: ;************************************************************
forrest@0:
forrest@0: ; to be able to handle observation at (-89.98,-24.80)
forrest@0: ym(0) = -90.
forrest@0:
forrest@0: yy = linint2_points_Wrap(xm,ym,data_mod0,True,lon_ob,lat_ob,0)
forrest@0:
forrest@0: delete (data_mod0)
forrest@0: yy!0 = "field"
forrest@0: data_mod = yy(pts|:,field|:,time|:)
forrest@0:
forrest@0: ;************************************************************
forrest@0: ; compute correlation coef and M score
forrest@0: ;************************************************************
forrest@0:
forrest@0: score_max = 5.
forrest@0:
forrest@0: ccr = new ((/nstation, nfield/),float)
forrest@0: M_score = new ((/nstation, nfield/),float)
forrest@0:
forrest@0: do n=0,nstation-1
forrest@0: do m=0,nfield-1
forrest@0: ccr(n,m) = esccr(data_ob(n,m,:),data_mod(n,m,:),0)
forrest@0: bias = sum(abs(data_mod(n,m,:)-data_ob(n,m,:))/(abs(data_mod(n,m,:))+abs(data_ob(n,m,:))))
forrest@0: M_score(n,m) = (1. -(bias/nmon)) * score_max
forrest@0: end do
forrest@0: end do
forrest@0:
forrest@0: M_co2 = avg(M_score(:,0))
forrest@0: M_rad = avg(M_score(:,1))
forrest@0: M_lh = avg(M_score(:,2))
forrest@0: M_sh = avg(M_score(:,3))
forrest@0: M_all = M_co2+ M_rad +M_lh + M_sh
forrest@0:
forrest@0: M_energy_co2 = sprintf("%.2f", M_co2)
forrest@0: M_energy_rad = sprintf("%.2f", M_rad)
forrest@0: M_energy_lh = sprintf("%.2f", M_lh )
forrest@0: M_energy_sh = sprintf("%.2f", M_sh )
forrest@0: M_energy_all = sprintf("%.2f", M_all)
forrest@0:
forrest@0: ;*******************************************************************
forrest@0: ; for station line plot
forrest@0: ;*******************************************************************
forrest@0:
forrest@0: ; for x-axis in xyplot
forrest@0: mon = ispan(1,12,1)
forrest@0: mon@long_name = "month"
forrest@0:
forrest@0: res = True ; plot mods desired
forrest@0: res@xyLineThicknesses = (/2.0,2.0/) ; make 2nd lines thicker
forrest@0: res@xyLineColors = (/"blue","red"/) ; line color (ob,model)
forrest@0: ;-------------------------------------------------------------------------
forrest@0: ; Add a boxed legend using the more simple method
forrest@0:
forrest@0: res@pmLegendDisplayMode = "Always"
forrest@0: ; res@pmLegendWidthF = 0.1
forrest@0: res@pmLegendWidthF = 0.08
forrest@0: res@pmLegendHeightF = 0.06
forrest@0: ; res@pmLegendOrthogonalPosF = -1.17
forrest@0: ; res@pmLegendOrthogonalPosF = -1.00 ;(downward)
forrest@0: res@pmLegendOrthogonalPosF = -0.30 ;(downward)
forrest@0:
forrest@0: ; res@pmLegendParallelPosF = 0.18
forrest@0: res@pmLegendParallelPosF = 0.23 ;(rightward)
forrest@0:
forrest@0: ; res@lgPerimOn = False
forrest@0: res@lgLabelFontHeightF = 0.015
forrest@0: res@xyExplicitLegendLabels = (/"observed",model_name/)
forrest@0: ;-------------------------------------------------------------------
forrest@0: ; for panel plot
forrest@0: res@gsnFrame = False ; Do not draw plot
forrest@0: res@gsnDraw = False ; Do not advance frame
forrest@0:
forrest@0: pres = True ; panel plot mods desired
forrest@0: pres@gsnPanelYWhiteSpacePercent = 5 ; increase white space around
forrest@0: ; indiv. plots in panel
forrest@0: pres@gsnMaximize = True ; fill the page
forrest@0: ;-------------------------------------------------------------------
forrest@0:
forrest@0: plot_data = new((/2,12/),float)
forrest@0: plot_data!0 = "case"
forrest@0: plot_data!1 = "month"
forrest@0:
forrest@0: do n = 0,nstation-1
forrest@0: ;----------------------------
forrest@0: ; for observed
forrest@0:
forrest@0: plot_name = station(n)+"_ob"
forrest@0: title = station(n)+"("+sprintf("%5.2f",lat_ob(n))+","+sprintf("%5.2f",lon_ob(n))+")"
forrest@0: res@tiMainString = title
forrest@0:
forrest@0: wks = gsn_open_wks (plot_type,plot_name)
forrest@0: plot=new(4,graphic) ; create graphic array
forrest@0:
forrest@0: plot_data(0,:) = (/data_ob (n,0,:)/)
forrest@0: plot_data@long_name = field(0)
forrest@0: plot(0)=gsn_csm_xy(wks,mon,plot_data(0,:),res) ; create plot 1
forrest@0:
forrest@0: plot_data(0,:) = (/data_ob (n,1,:)/)
forrest@0: plot_data@long_name = field(1)
forrest@0: plot(1)=gsn_csm_xy(wks,mon,plot_data(0,:),res) ; create plot 2
forrest@0:
forrest@0: plot_data(0,:) = (/data_ob (n,2,:)/)
forrest@0: plot_data@long_name = field(2)
forrest@0: plot(2)=gsn_csm_xy(wks,mon,plot_data(0,:),res) ; create plot 3
forrest@0:
forrest@0: plot_data(0,:) = (/data_ob (n,3,:)/)
forrest@0: plot_data@long_name = field(3)
forrest@0: plot(3)=gsn_csm_xy(wks,mon,plot_data(0,:),res) ; create plot 4
forrest@0:
forrest@0: gsn_panel(wks,plot,(/2,2/),pres) ; create panel plot
forrest@0:
forrest@0: system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
forrest@0: "rm "+plot_name+"."+plot_type)
forrest@0:
forrest@0: clear (wks)
forrest@0: delete (plot)
forrest@0: ;----------------------------
forrest@0: ; for model_vs_ob
forrest@0:
forrest@0: plot_name = station(n)+"_model_vs_ob"
forrest@0: title = station(n)+"("+sprintf("%5.2f",lat_ob(n))+","+sprintf("%5.2f",lon_ob(n))+")"
forrest@0: res@tiMainString = title
forrest@0:
forrest@0: wks = gsn_open_wks (plot_type,plot_name)
forrest@0: plot=new(4,graphic) ; create graphic array
forrest@0:
forrest@0: plot_data(0,:) = (/data_ob (n,0,:)/)
forrest@0: plot_data(1,:) = (/data_mod(n,0,:)/)
forrest@0: plot_data@long_name = field(0)
forrest@0: plot(0)=gsn_csm_xy(wks,mon,plot_data,res) ; create plot 1
forrest@0:
forrest@0: plot_data(0,:) = (/data_ob (n,1,:)/)
forrest@0: plot_data(1,:) = (/data_mod(n,1,:)/)
forrest@0: plot_data@long_name = field(1)
forrest@0: plot(1)=gsn_csm_xy(wks,mon,plot_data,res) ; create plot 2
forrest@0:
forrest@0: plot_data(0,:) = (/data_ob (n,2,:)/)
forrest@0: plot_data(1,:) = (/data_mod(n,2,:)/)
forrest@0: plot_data@long_name = field(2)
forrest@0: plot(2)=gsn_csm_xy(wks,mon,plot_data,res) ; create plot 3
forrest@0:
forrest@0: plot_data(0,:) = (/data_ob (n,3,:)/)
forrest@0: plot_data(1,:) = (/data_mod(n,3,:)/)
forrest@0: plot_data@long_name = field(3)
forrest@0: plot(3)=gsn_csm_xy(wks,mon,plot_data,res) ; create plot 4
forrest@0:
forrest@0: gsn_panel(wks,plot,(/2,2/),pres) ; create panel plot
forrest@0:
forrest@0: system("convert "+plot_name+"."+plot_type+" "+plot_name+"."+plot_type_new+";"+ \
forrest@0: "rm "+plot_name+"."+plot_type)
forrest@0:
forrest@0: clear (wks)
forrest@0: delete (plot)
forrest@0: end do
forrest@0:
forrest@0: ;*******************************************************************
forrest@0: ; html table of site: observed
forrest@0: ;*******************************************************************
forrest@0: output_html = "line_ob.html"
forrest@0:
forrest@0: header = (/"" \
forrest@0: ,"" \
forrest@0: ,"CLAMP metrics" \
forrest@0: ,"" \
forrest@0: ,"Energy at Site: Observation
" \
forrest@0: /)
forrest@0: footer = ""
forrest@0:
forrest@0: table_header = (/ \
forrest@0: "" \
forrest@0: ,"" \
forrest@0: ," | Site Name | " \
forrest@0: ," Latitude | " \
forrest@0: ," Longitude | " \
forrest@0: ," Observed | " \
forrest@0: ,"
" \
forrest@0: /)
forrest@0: table_footer = "
"
forrest@0: row_header = ""
forrest@0: row_footer = "
"
forrest@0:
forrest@0: lines = new(50000,string)
forrest@0: nline = 0
forrest@0:
forrest@0: set_line(lines,nline,header)
forrest@0: set_line(lines,nline,table_header)
forrest@0: ;-----------------------------------------------
forrest@0: ; row of table
forrest@0:
forrest@0: do n = 0,nstation-1
forrest@0: set_line(lines,nline,row_header)
forrest@0:
forrest@0: txt0 = station(n)
forrest@0: txt1 = sprintf("%5.2f", lat_ob(n))
forrest@0: txt2 = sprintf("%5.2f", lon_ob(n))
forrest@0: txt3 = year_ob(n)
forrest@0:
forrest@0: set_line(lines,nline,""+txt0+" | ")
forrest@0: set_line(lines,nline,""+txt1+" | ")
forrest@0: set_line(lines,nline,""+txt2+" | ")
forrest@0: set_line(lines,nline,""+txt3+" | ")
forrest@0:
forrest@0: set_line(lines,nline,row_footer)
forrest@0: end do
forrest@0: ;-----------------------------------------------
forrest@0: set_line(lines,nline,table_footer)
forrest@0: set_line(lines,nline,footer)
forrest@0:
forrest@0: ; Now write to an HTML file.
forrest@0: idx = ind(.not.ismissing(lines))
forrest@0: if(.not.any(ismissing(idx))) then
forrest@0: asciiwrite(output_html,lines(idx))
forrest@0: else
forrest@0: print ("error?")
forrest@0: end if
forrest@0: delete (idx)
forrest@0:
forrest@0: ;*******************************************************************
forrest@0: ; score and line table : model vs observed
forrest@0: ;*******************************************************************
forrest@0: output_html = "score+line_vs_ob.html"
forrest@0:
forrest@0: header = (/"" \
forrest@0: ,"" \
forrest@0: ,"CLAMP metrics" \
forrest@0: ,"" \
forrest@0: ,"Energy at Site: Model "+model_name+"
" \
forrest@0: /)
forrest@0: footer = ""
forrest@0:
forrest@0: delete (table_header)
forrest@0: table_header = (/ \
forrest@0: "" \
forrest@0: ,"" \
forrest@0: ," | Site Name | " \
forrest@0: ," Latitude | " \
forrest@0: ," Longitude | " \
forrest@0: ," Observed | " \
forrest@0: ," CO2 Flux | " \
forrest@0: ," Net Radiation | " \
forrest@0: ," Latent Heat | " \
forrest@0: ," Sensible Heat | " \
forrest@0: ," Average | " \
forrest@0: ,"
" \
forrest@0: /)
forrest@0: table_footer = "
"
forrest@0: row_header = ""
forrest@0: row_footer = "
"
forrest@0:
forrest@0: lines = new(50000,string)
forrest@0: nline = 0
forrest@0:
forrest@0: set_line(lines,nline,header)
forrest@0: set_line(lines,nline,table_header)
forrest@0: ;-----------------------------------------------
forrest@0: ; row of table
forrest@0:
forrest@0: do n = 0,nstation-1
forrest@0: set_line(lines,nline,row_header)
forrest@0:
forrest@0: txt0 = station(n)
forrest@0: txt1 = sprintf("%5.2f", lat_ob(n))
forrest@0: txt2 = sprintf("%5.2f", lon_ob(n))
forrest@0: txt3 = year_ob(n)
forrest@0: txt4 = sprintf("%5.2f", M_score(n,0))
forrest@0: txt5 = sprintf("%5.2f", M_score(n,1))
forrest@0: txt6 = sprintf("%5.2f", M_score(n,2))
forrest@0: txt7 = sprintf("%5.2f", M_score(n,3))
forrest@0: txt8 = sprintf("%5.2f", avg(M_score(n,:)))
forrest@0:
forrest@0: set_line(lines,nline,""+txt0+" | ")
forrest@0: set_line(lines,nline,""+txt1+" | ")
forrest@0: set_line(lines,nline,""+txt2+" | ")
forrest@0: set_line(lines,nline,""+txt3+" | ")
forrest@0: set_line(lines,nline,""+txt4+" | ")
forrest@0: set_line(lines,nline,""+txt5+" | ")
forrest@0: set_line(lines,nline,""+txt6+" | ")
forrest@0: set_line(lines,nline,""+txt7+" | ")
forrest@0: set_line(lines,nline,""+txt8+" | ")
forrest@0:
forrest@0: set_line(lines,nline,row_footer)
forrest@0: end do
forrest@0:
forrest@0: ; last row, summary
forrest@0: set_line(lines,nline,row_header)
forrest@0:
forrest@0: txt0 = "All_"+sprintf("%.0f", nstation)
forrest@0: txt1 = "-"
forrest@0: txt2 = "-"
forrest@0: txt3 = "-"
forrest@0: txt4 = M_energy_co2
forrest@0: txt5 = M_energy_rad
forrest@0: txt6 = M_energy_lh
forrest@0: txt7 = M_energy_sh
forrest@0: txt8 = M_energy_all
forrest@0:
forrest@0: set_line(lines,nline,""+txt0+" | ")
forrest@0: set_line(lines,nline,""+txt1+" | ")
forrest@0: set_line(lines,nline,""+txt2+" | ")
forrest@0: set_line(lines,nline,""+txt3+" | ")
forrest@0: set_line(lines,nline,""+txt4+" | ")
forrest@0: set_line(lines,nline,""+txt5+" | ")
forrest@0: set_line(lines,nline,""+txt6+" | ")
forrest@0: set_line(lines,nline,""+txt7+" | ")
forrest@0: set_line(lines,nline,""+txt8+" | ")
forrest@0:
forrest@0: set_line(lines,nline,row_footer)
forrest@0: ;-----------------------------------------------
forrest@0: set_line(lines,nline,table_footer)
forrest@0: set_line(lines,nline,footer)
forrest@0:
forrest@0: ; Now write to an HTML file.
forrest@0: idx = ind(.not.ismissing(lines))
forrest@0: if(.not.any(ismissing(idx))) then
forrest@0: asciiwrite(output_html,lines(idx))
forrest@0: else
forrest@0: print ("error?")
forrest@0: end if
forrest@0: delete (idx)
forrest@0:
forrest@0: ;**************************************************************************************
forrest@0: ; update score
forrest@0: ;**************************************************************************************
forrest@0:
forrest@0: if (isvar("compare")) then
forrest@0: system("sed 1,/M_energy_co2/s//"+M_energy_co2+"/ "+html_name2+" > "+html_new2+";"+ \
forrest@0: "mv -f "+html_new2+" "+html_name2+";"+ \
forrest@0: "sed 1,/M_energy_rad/s//"+M_energy_rad+"/ "+html_name2+" > "+html_new2+";"+ \
forrest@0: "mv -f "+html_new2+" "+html_name2+";"+ \
forrest@0: "sed 1,/M_energy_lh/s//"+M_energy_lh+"/ "+html_name2+" > "+html_new2+";"+ \
forrest@0: "mv -f "+html_new2+" "+html_name2+";"+ \
forrest@0: "sed 1,/M_energy_sh/s//"+M_energy_sh+"/ "+html_name2+" > "+html_new2+";"+ \
forrest@0: "mv -f "+html_new2+" "+html_name2)
forrest@0: end if
forrest@0:
forrest@0: system("sed s#M_energy_co2#"+M_energy_co2+"# "+html_name+" > "+html_new+";"+ \
forrest@0: "mv -f "+html_new+" "+html_name+";"+ \
forrest@0: "sed s#M_energy_rad#"+M_energy_rad+"# "+html_name+" > "+html_new+";"+ \
forrest@0: "mv -f "+html_new+" "+html_name+";"+ \
forrest@0: "sed s#M_energy_lh#"+M_energy_lh+"# "+html_name+" > "+html_new+";"+ \
forrest@0: "mv -f "+html_new+" "+html_name+";"+ \
forrest@0: "sed s#M_energy_sh#"+M_energy_sh+"# "+html_name+" > "+html_new+";"+ \
forrest@0: "mv -f "+html_new+" "+html_name)
forrest@0:
forrest@0: ;***************************************************************************
forrest@0: ; output plots
forrest@0: ;***************************************************************************
forrest@0: output_dir = model_name+"/energy"
forrest@0:
forrest@0: system("mv *.png *.html " + output_dir)
forrest@0: ;***************************************************************************
forrest@0: end