Project

General

Profile

README_ECOCLIMAP.txt

diane tzanos, 01/07/2020 09:59 AM

 
1
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
2
!ECO_LIC Copyright 2003-2014 Meteo-France 
3
!ECO_LIC This is part of the ECOCLIMAP software and data files 
4
!ECO_LIC governed by the CeCILL-C licence version 1
5
!ECO_LIC See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
6
!ECO_LIC for details. version 1.
7
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
8

    
9
-------------------------------------------------------------------------
10

    
11
How to make its own surface parameters maps with ECOCLIMAP?
12

    
13
-------------------------------------------------------------------------
14

    
15
1) Input data
16
   ----------
17
1.1) In the directory DATA_VEG :
18
the file containing the 215 ecosystems at 1km (ECOCLIMAP_I_GLOBAL_V1.5.dir)
19
    (933120000 words expanded)
20
     and ECOCLIMAP_I_GLOBAL.hdr
21
the file containing the 273 ecosystems at 1km (ECOCLIMAP_II_EUROP_V2.2.dir)
22
    (1866240000 words expanded)
23
    and ECOCLIMAP_II_EUROP.hdr
24

    
25
1.2 ) In the directory DATA_SOIL :
26
the file of clay fraction at 10km             (clay_fao.dir)
27
    (  9331200 words expanded)
28
    and clay_fao.hdr
29
the file of sand fraction at 10km             (sand_fao.dir)
30
    (  9331200 words expanded)
31
    and sand_fao.hdr
32

    
33

    
34
2) programs to compute the parameters are in the directory: SRC
35
   -------------------------------------------------------
36

    
37
3) to run the programs:
38
   -------------------
39

    
40
  You need a FORTRAN90 compiler.
41
  Edit the script script_compile to use your OWN compiler, with its own
42
  options. Choose the write_cover_tex file functions of your compiler. 
43
  3 possible are gfortran, f95, pgf90. It needs to be done once.
44

    
45
  The script script_compile will move the executable in the directory RUN
46

    
47
4) computation of the parameters:
48
   -----------------------------
49

    
50
It is highly recommended to use a separate directory to perform the 
51
computations (the code is supposed to run in the RUN directory).
52

    
53
a) Modify the script_compile file in the directory SRC
54

    
55
c) Edit the ecoclimap.f90 program to choose your final resolution
56
   and your domain. All explaination are given in the program itself.
57

    
58
d) run the script script_compile
59

    
60
e) GO TO the RUN directory and run ecoclimap.exe
61
   The execution time may vary from one computer (and compiler) from
62
   another. For example, it take 1h to have the parameters aggregated
63
   at 1 degre of resolution  at a global scale on a 1GHz Linux PC.
64
   The execution time depends linearly of the number of OUTPUT points
65
   (not input points).
66
   So, at 0.5 degre of resolution at a global scale, it would spend 4h
67
   to run.
68

    
69
f) look-up tables are printed in the class_cover_data.tex file.
70

    
71
5) Output parameters
72
   -----------------
73

    
74
   Results are coded in output files directly in the RUN directory.
75
   There is one file per data and per month. Files are binaries files.
76
   The output file names are defined as follows:
77

    
78
    PARAM.NN.TT.NSCALE
79

    
80
  where
81

    
82
    PARAM=name of the parameter:
83

    
84
     parameters defined for each tile:
85
     - - - - - - - - - - - - - - - - -
86

    
87
       "frac_tile": fraction of each tile on the continental part
88
                    (excluding towns). So, these tiles are not defined
89
		    where there is no continent, and their sum is equal
90
		    to 1 on continent, whatever the water or town covers:
91
		    Sum (frac_tileNN) = 1.
92

    
93
       "lai"      : Leaf Area Index
94
       "veg"      : vegetation fraction (green + dead)
95
       "greeness" : greeness fraction
96

    
97
       "alb"      : albedo (vegetation+soil), WITHOUT snow
98
                    alb field is made using alb_veg and alb_soil fields:
99

    
100
       "alb_veg"  : vegetation albedo (without the bare soil part).
101
       "alb_soil" : bare soil albedo (without the vegetation part).
102

    
103
       "emis"     : emissivity (vegetation+soil), WITHOUT snow
104
       "rsmin"    : Minimal stomathal resistance
105
       "z0"       : roughness length
106
       "z0h"      : roughness length for heat
107

    
108
       "d_root"   : root depth
109
       "d_soil"   : total soil depth
110

    
111
     parameters defined once, for any number of tiles:
112
     - - - - - - - - - - - - - - - - - - - - - - - - -
113

    
114
       "clay"     : clay fraction
115
       "sand"     : sand fraction
116

    
117
       "F_bare"   : fraction of totally bare soil over continent
118
       "F_lowv"   : fraction of low vegetation over continent
119
       "F_tree"   : fraction of woody vegetation over continent
120
                     F_bare + F_lowv + F_tree = 1.
121
       "F_urb"    : urbanization fraction
122
       "F_wat"    : water fraction (seas and lakes)
123

    
124
       "f_*"      : fraction of each vegetation type over continent
125
                      Sum (f_*) = 1.
126

    
127
       "glob_alb" : global albedo      (all tiles, town and water averaged)
128
       "glob_emis": global emissivity  (all tiles, town and water averaged) 
129

    
130

    
131
  and where, if necessary:
132

    
133
    NN=    month  (1=january, 12=december), 
134
        or decade (1=first decade of january, 36=last decade of december)
135

    
136
    TT=    value of the tile described in this file
137

    
138
    NSCALE= averaging value (see program ecoclimap.f90)
139

    
140

    
141

    
142
6) Reading of the output files
143
   ---------------------------
144

    
145
  There is an exemple of short program encapsulated in a shell script
146
  (convert_lai.sh) to read the output files in fortran. Here, the program
147
  write 10*LAI in integer in a ASCII file.
148

    
149
-------------------------------------------------------------------------------
150

    
151
All questions or remarks can be addressed to surfex-support@meteo.fr
152

    
153
-------------------------------------------------------------------------------
154
The scientific references for ECOCLIMAP are : 
155
Masson V., Champeaux J.-L., Chauvin F., Meriguet C., and Lacaze R.: 
156
A Global Database of Land Surface Parameters at 1-km Resolution in 
157
Meteorological and Climate Models, J. Clim., 16,1261-1282, 2003.
158

    
159
Faroux, S., Kaptué Tchuenté, A. T., Roujean, J.-L., Masson, V., Martin, E.,
160
 and Le Moigne, P.: ECOCLIMAP-II/Europe: a twofold database of ecosystems 
161
and surface parameters at 1 km resolution based on satellite information 
162
for use in land surface, meteorological and climate models, Geosci. Model Dev.,
163
 6, 563-582, doi:10.5194/gmd-6-563-2013, 2013.