Re: Calculating total carbon of different vegetation types

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Thu Nov 07 2013 - 11:00:03 MST

Lifen,

It doesn't seem there's any way to tell if you are using the correct vegetation values, because your "Vegetation_types.nc" file contains no metadata that explains the values.

Did you look at your "t2" data at all?

I added the line:

   printMinMax(t2,0)

and noticed that "t2" goes from 1 to 10, but you are checking for values 0 to 9 in your do loop. Perhaps this is where your error is:

bc_density = where(t2.eq. i, t1, bc_density@_FillValue)

Change "t2.eq.i" to "t2.eq.(i+1)"

--Mary

On Nov 6, 2013, at 8:12 AM, "Jiang, Lifen" <lfjiang@ou.edu> wrote:

> Hi all,
>
> I have global carbon data on lands (entitled "Carbon") and global vegetation type data (entitled "Vegetation_types "). I was trying to obtain total carbon of different vegetation types using code in the file "Code".
>
> However, the results (as in excel) seemed not realistic because I would expect tropical area should have highest values judging from the map of carbon.
>
> I am worried if I confused the vegetation names but I couldn't find where was wrong.
>
> I appreciate if anybody could help check if it is matched correctly between total carbon and vegetation names.
>
> Thanks,
>
> Lifen
> <Carbon.nc><Vegetation_types.nc><Code.ncl><Carbon of vegetation types .xlsx>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Nov 7 11:00:14 2013

This archive was generated by hypermail 2.1.8 : Mon Nov 11 2013 - 09:45:33 MST