Opened 22 months ago
Closed 21 months ago
#2867 closed help (answered)
fluxnet headers
Reported by: | awright | Owned by: | pmcguire |
---|---|---|---|
Component: | JULES | Keywords: | FLUXNET, JULES |
Cc: | Platform: | JASMIN | |
UM Version: |
Description
Hello Patrick,
I have downloaded the fluxnet observation datasets (so I can compare them with my model output) from here:
/group_workspaces/jasmin2/jules/pmcguire/fluxnet/kwilliam/suite_data/vn1.0/
As an example I have the directory 'fluxnet' in there, which has AT_Neu-met.dat file in it with 8 columns. The first row is:
265.05, 0, 218.687, 0.00187, 90957, 0, 1.49, 390.43
I want to find out what the column headers are:
In /home/users/azin/roses/u-bh383/app/jules/opt/rose-app-drive-AT_Neu.conf I find out that
there are 7 variables in AT_Neu-met.dat.
file='$FLUXNET2015_DRIVE_DIR/AT_Neu-met.dat'
var='t','sw_down','lw_down','q','pstar','precip','wind'
I was expecting 8 variables. How can I find out more about the column headers (e.g. their units) and why do I have 7 of them instead of 8?
All the best,
Azin
Change History (7)
comment:1 Changed 22 months ago by pmcguire
comment:2 Changed 22 months ago by pmcguire
- Status changed from new to accepted
comment:3 Changed 22 months ago by awright
Hello Patrick,
I am trying to find the headers for observed datasets located in :
/group_workspaces/jasmin2/jules/pmcguire/fluxnet/kwilliam/suite_data/vn1.0/lba_obs and also in …/fluxnet_obs.
1- In lba_obs directory, I have two types of files: carbon and energy.(e.g. BANday-carbon.dat and BANday-energy.dat). How do I find out about the headers?
2- In fluxnet_obs directory, I have one type of file : energy and carbon (e.g. AT_Neu-energyandcarbon-dailyUTC.dat). The headers are already there: YYYYMMDD_UTC, GPP, Reco, NEE, SH, LE. How can I find out more about them?
I have looked at the .conf, .rc and .py files and have not found the information yet.
Many thanks,
Azin
comment:4 Changed 22 months ago by awright
Hello Patrick,
I think I have found where the headers are:
/group_workspaces/jasmin2/jules/pmcguire/fluxnet/kwilliam/suite_data/vn1.0/lba_obs/process_LBA_rose.ncl :
C_hdr = (/"Year_LBAMIP","DoY_LBAMIP","Hour_LBAMIP","NEE","NEEf","NEE_model", \
"Re_5day_ust_Sco2_LUT","GEP_model","par_fill","VPD","mrs"/)
ncolC = dimsizes(C_hdr)
E_hdr = (/"Year_LBAMIP","DoY_LBAMIP","Hour_LBAMIP","FG","Leraw","LE","Hraw","H","Rn"/)
ncolE = dimsizes(E_hdr)
Where C_hdr refers to headers of [sitename, e.g BAN]day_carbon.dat and E_hdr refers to headers of [sitename, e.g BAN]day_energy.dat
All the best,
Azin
comment:5 Changed 22 months ago by pmcguire
Dear Azin:
Excellent! Nice work!
These variables are also listed in the fluxnet_evaluation.py program under daily_energy_file_headings daily_carbon_file_headings, for both the FLUXNET 2015 observations and differently for the LBA observations.
Some of the variables' units are redefined or recomputed in the fluxnet_evaluation.py program (in order to match the modelled output units from JULES) with commands like cube.units = or cube *=.
I found that some of the LBA variables' units are defined (for example for Re_5day_ust_Sco2_LUT) at:
ftp://saleskalab.eebweb.arizona.edu/pub/BrasilFlux_Data/Version.2.0/Documentation/README_Files.txt
and
ftp://saleskalab.eebweb.arizona.edu/pub/BrasilFlux_Data/Version.2.0/Documentation/Brasil%20flux%20network_2.0.pptx
Does this help?
Patrick
comment:6 Changed 22 months ago by awright
Hello Patrick,
Yes, it definitely helps. Thank you Very much!
Azin
comment:7 Changed 21 months ago by pmcguire
- Resolution set to answered
- Status changed from accepted to closed
Hi Azin
To figure out the units of the driving variables, first note that your file /home/users/azin/roses/u-bh383/app/jules/opt/rose-app-drive-AT_Neu.conf defines the settings for the jules_drive namelist of JULES. So then you can go to the JULES Userguide and search for jules_drive. You will then find this namelist described in http://jules-lsm.github.io/vn4.9/namelists/drive.nml.html .
In there, you can find out about the units.
I suspect that the 8th column of the FLUXNET driving data sets is an unused column. The values of the first 7 columns correspond to the values I would expect for the 7 variables: var='t','sw_down','lw_down','q','pstar','precip','wind'.
Patrick