change flux unit for STILT module

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="")