From 2fb7445d5e4089592d8a73a297c1e9570073e1f4 Mon Sep 17 00:00:00 2001 From: weihe <amvdw95@gmail.com> Date: Sun, 5 Jul 2015 17:37:23 +0000 Subject: [PATCH] change flux unit for STILT module --- da/stilt/stilt.co2.simu.2015.May.r | 27 +++++++++++++++++---------- 1 file changed, 17 insertions(+), 10 deletions(-) diff --git a/da/stilt/stilt.co2.simu.2015.May.r b/da/stilt/stilt.co2.simu.2015.May.r index f2e1b24..d3314b1 100755 --- a/da/stilt/stilt.co2.simu.2015.May.r +++ b/da/stilt/stilt.co2.simu.2015.May.r @@ -604,11 +604,18 @@ for(d in 1:ndays) } # calculate the mean values for these fluxes,mean by row/col -gppmean_hour=rowMeans(gpp, na.rm = TRUE, dims = 2) -resmean_hour=rowMeans(res, na.rm = TRUE, dims = 2) -ocnmean_hour=rowMeans(ocn, na.rm = TRUE, dims = 2) -fosmean_hour=rowMeans(fos, na.rm = TRUE, dims = 2) -firmean_hour=rowMeans(fir, na.rm = TRUE, dims = 2) +gppmean_hour=rowMeans(gpp, na.rm = TRUE, dims = 2)*1e-6 +resmean_hour=rowMeans(res, na.rm = TRUE, dims = 2)*1e-6 +ocnmean_hour=rowMeans(ocn, na.rm = TRUE, dims = 2)*1e-6 +fosmean_hour=rowMeans(fos, na.rm = TRUE, dims = 2)*1e-6 +firmean_hour=rowMeans(fir, na.rm = TRUE, dims = 2)*1e-6 + +gppmean_hour=replace(gppmean_hour,is.nan(gppmean_hour),0) +resmean_hour=replace(resmean_hour,is.nan(resmean_hour),0) +ocnmean_hour=replace(ocnmean_hour,is.nan(ocnmean_hour),0) +fosmean_hour=replace(fosmean_hour,is.nan(fosmean_hour),0) +firmean_hour=replace(firmean_hour,is.nan(firmean_hour),0) + xvals <- -179.5:179.5 yvals <- -89.5:89.5 @@ -617,11 +624,11 @@ xdim <- ncdim_def( 'Lon', 'degree', xvals ) ydim <- ncdim_def( 'Lat', 'degree', yvals ) mv <- 0 # missing value -var_gpp <- ncvar_def( name="flux_gpp_prior_mean", units="umol/m2/sec", dim=list(xdim,ydim), longname="Gross Primary Productivity",missval=mv ) -var_res <- ncvar_def( name="flux_res_prior_mean", units="umol/m2/sec", dim=list(xdim,ydim), longname="Total Ecosystem Respiration",missval=mv ) -var_ocn <- ncvar_def( name="flux_ocean_prior_mean", units="umol/m2/sec", dim=list(xdim,ydim), longname="Ocean CO2 assimilation",missval=mv ) -var_fos <- ncvar_def( name="flux_ff_prior_mean", units="umol/m2/sec", dim=list(xdim,ydim), longname="Fossil fuel CO2 emission",missval=mv ) -var_fir <- ncvar_def( name="flux_fires_prior_mean", units="umol/m2/sec", dim=list(xdim,ydim), longname="Fires CO2 emission",missval=mv ) +var_gpp <- ncvar_def( name="flux_gpp_prior_mean", units="mol/m2/sec", dim=list(xdim,ydim), longname="Gross Primary Productivity",missval=mv ) +var_res <- ncvar_def( name="flux_res_prior_mean", units="mol/m2/sec", dim=list(xdim,ydim), longname="Total Ecosystem Respiration",missval=mv ) +var_ocn <- ncvar_def( name="flux_ocean_prior_mean", units="mol/m2/sec", dim=list(xdim,ydim), longname="Ocean CO2 assimilation",missval=mv ) +var_fos <- ncvar_def( name="flux_ff_prior_mean", units="mol/m2/sec", dim=list(xdim,ydim), longname="Fossil fuel CO2 emission",missval=mv ) +var_fir <- ncvar_def( name="flux_fires_prior_mean", units="mol/m2/sec", dim=list(xdim,ydim), longname="Fires CO2 emission",missval=mv ) output_fname=paste(outdir,"/","flux1x1_",tb2,"00","_",tb3,"00",".nc",sep="") -- GitLab