1 | #!/bin/ksh |
---|
2 | # Do not forget to exectute module load R and module load cdo on obelix |
---|
3 | # ${1} can take tree value : all= cdo extract + reformating R script + eval_figure R script |
---|
4 | # rscript= reformating R script + eval_figure R script |
---|
5 | # figure= eval_figure R script |
---|
6 | # ${2} wich step do you want : STOI, ORC-1, FIN, HIST |
---|
7 | # S{3} Output temporal frequency : Daily=D, Monthly=M |
---|
8 | # ${5} A name that describe your simulation experiment ex: reference, V1.X, r5800 etc.... |
---|
9 | # ${4} choose between SRF and SBG |
---|
10 | # $[6...X] A name that describe simulation experiment for wich you want to compare with ${2} |
---|
11 | # |
---|
12 | # Here your can put any variables you would like to extract |
---|
13 | |
---|
14 | if [ ${4} -eq 'SBG' ];then |
---|
15 | var='LAI,RDI,AGE_STAND,DIAMETER,IND,HEIGHT,RDI_TARGET_UPPER,RDI_TARGET_LOWER,GPP,NPP,IND_DOM,HEIGHT_DOM,DIAMETER_DOM,BA,CCTRW_001,CCTRW_002,CCTRW_003,MOISTRESS,HET_RESP,PLANT_STATUS,WOOD_VOLUME_PIX_CUT' |
---|
16 | else |
---|
17 | var='GPP,NPP,NEE,HUMREL,EVAP,ASSIMI,GS,GSMEAN,LAI,TVEG,POTEVAP,VC' |
---|
18 | fi |
---|
19 | |
---|
20 | OUTPATHS=/home/scratch01/$USER/IGCM_OUT/OL2/PROD/ensemble\ |
---|
21 | |
---|
22 | if [ X${1} == 'Xall' ]; then |
---|
23 | |
---|
24 | for i in DE-Bay DE-Hai DE-Tha DE-Wet DK-Sor FI-Hyy FI-Sod FR-Fon FR-Hes FR-LBr FR-Pue IT-Col IT-Cpz IT-Non IT-PT1 IT-Ren IT-SRo NL-Loo SE-Fla SE-Nor SE-Sk1 SE-Sk2 AU-Tum AU-Wac CZ-BK1 ES-ES1 IL-Yat IS-Gun IT-Lav IT-Lec IT-Non IT-Ro1 IT-Ro2 PT-Esp PT-Mi1 RU-Fyo RU-Ha1 RU-Ha2 RU-Zot SK-Tat UK-Gri UK-Ham UK-PL3 US-Blo US-Me4 US-SP1 US-SP2 US-SP3 US-SP4 US-Ho1 US-Ho2 US-Wi0 US-Wi2 US-Wi4 US-Wi5 US-Wi9 AU-Tum AU-Wac US-MMS US-MOz US-WBW US-Bar US-Ha1 US-Oho US-UMB US-WCr US-Wi1 US-Wi8 CA-Man CA-NS1 CA-NS2 CA-NS3 CA-NS4 CA-NS5 CA-Qcu CA-Qfo CA-SF1 CA-SF2 RU-Ha3 US-Aud US-FPe US-Goo US-Var AT-Neu CH-Oe1 DE-Gri DE-Meh DK-Lva ES-VDA FR-Lq1 FR-Lq2 HU-Bug HU-Mat IE-Dri IT-LMa IT-Mal IT-MBo NL-Cal NL-Haa NL-Hor UK-EBu UK-Tad US-Bkg DK-Ris FR-Gri NL-Lut US-Bo1 US-Ne1 US-Ne2 US-Ne3 NL-Lan BE-Vie CH-Lae IT-Cp2 CH-Dav DE-Obe BE-Bra |
---|
25 | do |
---|
26 | if [ ${2} -eq 'HIST' ];then |
---|
27 | if [ ${4} -eq 'SBG' ];then |
---|
28 | PATHS=${OUTPATHS}/FLUXNET/${i}FLUXNETHIST/SBG/Output/MO/ |
---|
29 | else |
---|
30 | PATHS=${OUTPATHS}/FLUXNET/${i}FLUXNETHIST/SRF/Output/MO/ |
---|
31 | fi |
---|
32 | else |
---|
33 | if [ ${4} -eq 'SBG' ];then |
---|
34 | PATHS=${OUTPATHS}/FLUXNET/${i}FLUXNET/SPIN/Output/${i}FLUXNET${2}/SBG/Output/MO/ |
---|
35 | else |
---|
36 | PATHS=${OUTPATHS}/FLUXNET/${i}FLUXNET/SPIN/Output/${i}FLUXNET${2}/SRF/Output/MO/ |
---|
37 | fi |
---|
38 | fi |
---|
39 | |
---|
40 | if [ ! -d ${PATHS} ];then |
---|
41 | continue |
---|
42 | fi |
---|
43 | if [ ${4} -eq 'SBG' ];then |
---|
44 | cdo select,name=${var} $PATHS*1M_stomate_history.nc $OUTPATHS/FLUXNET/${i}_ALL.nc |
---|
45 | else |
---|
46 | cdo select,name=${var} $PATHS*1M_sechiba_history.nc $OUTPATHS/FLUXNET/${i}_ALL.nc |
---|
47 | fi |
---|
48 | done |
---|
49 | Rscript --vanilla convert_nc_txt.R ${OUTPATHS}/FLUXNET/ ${OUTPATHS}/ ${@} |
---|
50 | Rscript --vanilla make_eval_figure.R ${OUTPATHS}/FLUXNET/ ${OUTPATHS}/ ${@} |
---|
51 | tar Jcvf $OUTPATHS/All_${5}.tar.xz $OUTPATHS/FLUXNET/*_ALL.nc |
---|
52 | |
---|
53 | fi |
---|
54 | |
---|
55 | |
---|
56 | if [ X${1} == 'Xrscript' ] ; then |
---|
57 | tar Jxvf $OUTPATHS/All_${5}.tar.xz -C $OUTPATHS/ |
---|
58 | mv ${OUTPATHS}/${OUTPATHS}/FLUXNET/*_ALL.nc $OUTPATHS/ |
---|
59 | Rscript --vanilla convert_nc_txt.R $OUTPATHS/ $OUTPATHS/ ${@} |
---|
60 | Rscript --vanilla make_eval_figure.R $OUTPATHS/ $OUTPATHS/ ${@} |
---|
61 | rm -rf *_ALL.nc |
---|
62 | fi |
---|
63 | |
---|
64 | |
---|
65 | if [ X${1} == 'Xfigure' ] ; then |
---|
66 | tar Jxvf $OUTPATHS/All_${5}.tar.xz -C $OUTPATHS/ |
---|
67 | mv ${OUTPATHS}/${OUTPATHS}/FLUXNET/*_ALL.nc $OUTPATHS/ |
---|
68 | Rscript --vanilla make_eval_figure.R $OUTPATHS/ $OUTPATHS/ ${@} |
---|
69 | rm -rf *_ALL.nc |
---|
70 | fi |
---|
71 | |
---|
72 | |
---|
73 | |
---|
74 | |
---|