Project

General

Profile

README_ECOCLIMAP.txt

florian suzat, 06/24/2014 03:23 PM

Download (5.92 KB)

 
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

    
18
the file containing the 215 ecosystems at 1km (ecoclimap_v1)
19
    (933120000 words expanded)
20
the file containing the 273 ecosystems at 1km (ecoclimap_v2)
21
    (1866240000 words expanded)
22
the file of clay fraction at 10km             (clay_10km)
23
    (  9331200 words expanded)
24
the file of sand fraction at 10km             (sand_10km)
25
    (  9331200 words expanded)
26

    
27
These files are in the directory: DATA
28
There are gzipped in the .tar.
29

    
30
2) programs to compute the parameters are in the directory: sources
31
   -------------------------------------------------------
32

    
33
3) to run the programs:
34
   -------------------
35

    
36
  You need a FORTRAN90 compiler.
37
  Edit the script script_compile to use your OWN compiler, with its own
38
  options. Choose the write_cover_tex file functions of your compiler. 3 possible are f95,  
39
  pgf90, g95. It needs to be done once.
40

    
41
4) computation of the parameters:
42
   -----------------------------
43

    
44
It is highly recommended to use a separate directory of the source or
45
DATA directories to perform the computations (for example in a directory
46
named RUN).
47

    
48
a) Put all the source files, the script_compile file and the 3 data files
49
   in your directory RUN
50

    
51
b) go to the RUN directory
52

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

    
56
d) run the script script_compile
57

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

    
67
f) look-up tables are printed in the class_cover_data.tex file.
68

    
69
5) Output parameters
70
   -----------------
71

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

    
76
    PARAM.NN.TT.NSCALE
77

    
78
  where
79

    
80
    PARAM=name of the parameter:
81

    
82
     parameters defined for each tile:
83
     - - - - - - - - - - - - - - - - -
84

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

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

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

    
98
       "alb_veg"  : vegetation albedo (without the bare soil part).
99
       "alb_soil" : bare soil albedo (without the vegetation part).
100

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

    
106
       "d_root"   : root depth
107
       "d_soil"   : total soil depth
108

    
109
     parameters defined once, for any number of tiles:
110
     - - - - - - - - - - - - - - - - - - - - - - - - -
111

    
112
       "clay"     : clay fraction
113
       "sand"     : sand fraction
114

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

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

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

    
128

    
129
  and where, if necessary:
130

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

    
134
    TT=    value of the tile described in this file
135

    
136
    NSCALE= averaging value (see program ecoclimap.f90)
137

    
138

    
139

    
140
6) Reading of the output files
141
   ---------------------------
142

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

    
147
-------------------------------------------------------------------------------
148

    
149
All questions or remarks can be adressed to Stéphanie Faroux & Valery Masson
150
CNRM, TOULOUSE (stephanie.faroux@meteo.fr, valery.masson@meteo.fr) or
151
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.