source: LMDZ6/branches/Optimisation_LMDZ/libf/phymar/Mod_SISVAT_TRV.f90 @ 5304

Last change on this file since 5304 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR


Integration of MAR physics

File size: 3.7 KB
Line 
1      module Mod_SISVAT_TRV
2
3
4!--------------------------------------------------------------------------+
5!                                                     Thu 28-Feb-2013  MAR |
6!     module Mod_SISVAT_TRV contains the constants of the                  |
7!                radiative tranfer model through Vegetation                |
8!                                                                          |
9!     version 3.p.4.1 created by H. Gallee,           Wed 30-Jan-2013      |
10!                    modified by H. Gallee,           Thu 28-Feb-2013      |
11!                                                                          |
12!--------------------------------------------------------------------------+
13
14
15! Global   Variables
16! ==================
17
18      use Mod_Real
19      use Mod_SISVAT_dim
20
21
22      IMPLICIT NONE
23
24
25
26! Internal Variables
27! ==================
28
29      integer, SAVE  ::   ivg
30
31
32      real(kind=real8), SAVE    ::  reVisL(0:nvgt)      ! Reflectivity  / Visible / Live Leaves
33      real(kind=real8), SAVE    ::  renIRL(0:nvgt)      ! Reflectivity  / Near IR / Live Leaves
34      real(kind=real8), SAVE    ::  trVisL(0:nvgt)      ! Transmitivity / Visible / Live Leaves
35      real(kind=real8), SAVE    ::  trnIRL(0:nvgt)      ! Transmitivity / Near IR / Live Leaves
36      real(kind=real8), SAVE    ::  reVisD(0:nvgt)      ! Reflectivity  / Visible / Dead Leaves
37      real(kind=real8), SAVE    ::  renIRD(0:nvgt)      ! Reflectivity  / Near IR / Dead Leaves
38      real(kind=real8), SAVE    ::  trVisD(0:nvgt)      ! Transmitivity / Visible / Dead Leaves
39      real(kind=real8), SAVE    ::  trnIRD(0:nvgt)      ! Transmitivity / Near IR / Dead Leaves
40
41      DATA (reVisL(ivg),renIRL(ivg),trVisL(ivg),trnIRL(ivg),           &
42     &      reVisD(ivg),renIRD(ivg),trVisD(ivg),trnIRD(ivg),ivg=0,nvgt)&
43
44!     reVisL renIRL trVisL trnIRL reVisD renIRD trVisD trnIRD    SVAT     CLASSES
45!     ------ ------ ------ ------ ------ ------ ------ --------+ ----------------
46     &/0.11,  0.58,  0.07,  0.25,  0.36,  0.58,  0.22,  0.38, &!  0 NO VEGETATION
47     & 0.11,  0.58,  0.07,  0.25,  0.36,  0.58,  0.22,  0.38, &!  1 CROPS LOW
48     & 0.11,  0.58,  0.07,  0.25,  0.36,  0.58,  0.22,  0.38, &!  2 CROPS MEDIUM
49     & 0.11,  0.58,  0.07,  0.25,  0.36,  0.58,  0.22,  0.38, &!  3 CROPS HIGH
50     & 0.11,  0.58,  0.07,  0.25,  0.36,  0.58,  0.22,  0.38, &!  4 GRASS LOW
51     & 0.11,  0.58,  0.07,  0.25,  0.36,  0.58,  0.22,  0.38, &!  5 GRASS MEDIUM
52     & 0.11,  0.58,  0.07,  0.25,  0.36,  0.58,  0.22,  0.38, &!  6 GRASS HIGH
53     & 0.10,  0.45,  0.05,  0.25,  0.16,  0.39,  0.01,  0.01, &!  7 BROADL LOW
54     & 0.10,  0.45,  0.05,  0.25,  0.16,  0.39,  0.01,  0.01, &!  8 BROADL MEDIUM
55     & 0.10,  0.45,  0.05,  0.25,  0.16,  0.39,  0.01,  0.01, &!  9 BROADL HIGH
56     & 0.07,  0.35,  0.05,  0.10,  0.10,  0.39,  0.01,  0.01, &! 10 NEEDLE LOW
57     & 0.07,  0.35,  0.05,  0.10,  0.10,  0.39,  0.01,  0.01, &! 11 NEEDLE MEDIUM
58     & 0.07,  0.35,  0.05,  0.10,  0.10,  0.39,  0.01,  0.01/  ! 12 NEEDLE HIGH
59
60
61
62      real(kind=real8), SAVE    ::  reVisS = 0.85       ! Reflectivity  / Visible / Canopy Snow
63      real(kind=real8), SAVE    ::  renIRS = 0.85       ! Reflectivity  / Near IR / Canopy Snow
64      real(kind=real8), SAVE    ::  trVisS = 0.00       ! Transmitivity / Visible / Canopy Snow
65      real(kind=real8), SAVE    ::  trnIRS = 0.00       ! Transmitivity / Near IR / Canopy Snow
66!     REMARK: Possible Refinement by taking actual Surface Snow Reflectivities
67!     ^^^^^^
68
69      real(kind=real8), SAVE    ::  snCaMx = 0.5        ! Canopy Snow Thickness for having Snow
70                                                  ! Snow Reflectivity and Transmitivity
71      real(kind=real8), SAVE    ::  CriStR = 25.        ! Critical Radiation Stomatal Resistance
72
73
74
75      end module Mod_SISVAT_TRV
Note: See TracBrowser for help on using the repository browser.