SURFEX  V8_0
Surfex V8_0 release
 All Classes Files Functions Variables
modd_snow_metamo.F90
Go to the documentation of this file.
1 !SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
2 !SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
3 !SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
4 !SFX_LIC for details. version 1.
5 ! ajoutEB
6 ! correction de l'erreur interversion de XVTANG2 et XVTANG3
7 !######################
8  MODULE modd_snow_metamo
9 ! ######################
10 !
11 !!**** *MODD_SNOW_METAMO* - declaration of parameters related
12 !! to snow metamorphism!!
13 !! PURPOSE
14 !! -------
15 ! The purpose of this declarative module is to specify the
16 ! parameters related to the metamorphism parameterization of snow.
17 !
18 !!
19 !!** IMPLICIT ARGUMENTS
20 !! ------------------
21 !! None
22 !!
23 !! REFERENCE
24 !! ---------
25 !!
26 !!
27 !! AUTHOR
28 !! ------
29 !! V. Vionnet *Meteo France*
30 !!
31 !! MODIFICATIONS
32 !! -------------
33 !! Original 02/2008
34 !-------------------------------------------------------------------------------
35 !
36 !* 0. DECLARATIONS
37 ! ------------
38 !
39 IMPLICIT NONE
40 !
41 !-------------------------------------------------------------------------------
42 !
43 ! minimum snow layer thickness for thermal calculations.
44 ! Used to prevent numerical problems as snow becomes vanishingly thin.
45 REAL, PARAMETER :: XSNOWDZMIN = 0.0001
46 !
47 ! Optical diameter properties
48 REAL, PARAMETER :: XDIAET = 1.E-4
49 REAL, PARAMETER :: XDIAGF = 3.E-4
50 REAL, PARAMETER :: XDIAFP = 4.E-4
51 !
52 ! Compaction/Settling Coefficients from Crocus v2.4
53 !
54 REAL, PARAMETER :: XVVISC1 = 7.62237E6 ! pre-exponential viscosity factor (UNIT : N m-2 s)
55 REAL, PARAMETER :: XVVISC3 = 0.023 ! density adjustement in the exponential correction for viscosity (UNIT : m3 kg-1)
56 REAL, PARAMETER :: XVVISC4 = .1 ! temperature adjustement in the exponential correction for viscosity (UNIT : K-1)
57 REAL, PARAMETER :: XVVISC5 = 1. ! factor for viscosity adjustement to grain type - to be checked
58 REAL, PARAMETER :: XVVISC6 = 60. ! factor for viscosity adjustement to grain type - to be checked
59 ! (especially this one ; inconsistency with Crocus v2.4)
60 REAL, PARAMETER :: XVVISC7 = 10. ! factor for viscosity adjustement to grain type - to be checked
61 REAL, PARAMETER :: XVRO11 = 250. ! normalization term for density dependence of the viscosity calculation (UNIT : kg m-3)
62 !
63 ! Maximum value for TPSNOW%GRAN2
64 REAL, PARAMETER :: XVGRAN1 = 99.
65 REAL, PARAMETER :: XGRAN = 99.
66 !
67 INTEGER, PARAMETER :: NVHIS1 = 1
68 INTEGER, PARAMETER :: NVHIS2 = 2
69 INTEGER, PARAMETER :: NVHIS3 = 3
70 INTEGER, PARAMETER :: NVHIS4 = 4
71 INTEGER, PARAMETER :: NVHIS5 = 5
72 !
73 ! Properties of fresh snow
74 REAL, PARAMETER :: XNDEN1 = 17.12
75 REAL, PARAMETER :: XNDEN2 = 128.
76 REAL, PARAMETER :: XNDEN3 = -20.
77 REAL, PARAMETER :: XNSPH1 = 7.87
78 REAL, PARAMETER :: XNSPH2 = 38.
79 REAL, PARAMETER :: XNSPH3 = 50.
80 REAL, PARAMETER :: XNSPH4 = 90.
81 !
82 REAL, PARAMETER :: XUEPSI = 1.E-8
83 REAL, PARAMETER :: XEPSI = 1.E-8
84 REAL, PARAMETER :: XUPOURC = 100.
85 !
86 ! Parameters for Marbouty's function
87 !
88 REAL, PARAMETER :: XVTANG1 = 40.
89 REAL, PARAMETER :: XVTANG2 = 6.
90 REAL, PARAMETER :: XVTANG3 = 22.
91 REAL, PARAMETER :: XVTANG4 = .7
92 REAL, PARAMETER :: XVTANG5 = .3
93 REAL, PARAMETER :: XVTANG6 = 6.
94 REAL, PARAMETER :: XVTANG7 = 1.
95 REAL, PARAMETER :: XVTANG8 = .8
96 REAL, PARAMETER :: XVTANG9 = 16.
97 REAL, PARAMETER :: XVTANGA = .2
98 REAL, PARAMETER :: XVTANGB = .2
99 REAL, PARAMETER :: XVTANGC = 18.
100 REAL, PARAMETER :: XVRANG1 = 400.
101 REAL, PARAMETER :: XVRANG2 = 150.
102 REAL, PARAMETER :: XVGANG1 = 70.
103 REAL, PARAMETER :: XVGANG2 = 25.
104 REAL, PARAMETER :: XVGANG3 = 40.
105 REAL, PARAMETER :: XVGANG4 = 50.
106 REAL, PARAMETER :: XVGANG5 = .1
107 REAL, PARAMETER :: XVGANG6 = 15.
108 REAL, PARAMETER :: XVGANG7 = .1
109 REAL, PARAMETER :: XVGANG8 = .55
110 REAL, PARAMETER :: XVGANG9 = .65
111 REAL, PARAMETER :: XVGANGA = .2
112 REAL, PARAMETER :: XVGANGB = .85
113 REAL, PARAMETER :: XVGANGC = .15
114 !
115 ! Parameters for snow metamorphism
116 !
117 REAL, PARAMETER :: XVDENT1 = 2314.81481
118 REAL, PARAMETER :: XVDENT2 = 7.2338E-7
119 REAL, PARAMETER :: XVGRAN6 = 51.
120 REAL, PARAMETER :: XVVAP1 = -6000.
121 REAL, PARAMETER :: XVVAP2 = .4
122 REAL, PARAMETER :: XVDIAM1 = 4.E-4
123 REAL, PARAMETER :: XVDIAM2 = 5.E-4
124 REAL, PARAMETER :: XVDIAM3 = 3.E-4
125 REAL, PARAMETER :: XVDIAM4 = 2.E-4
126 REAL, PARAMETER :: XVDIAM5 = 1.E-4
127 REAL, PARAMETER :: XVDIAM6 = 1.E-4
128 REAL, PARAMETER :: XVSPHE1 = 1.
129 REAL, PARAMETER :: XVSPHE2 = 11574.074
130 REAL, PARAMETER :: XVSPHE3 = .5
131 REAL, PARAMETER :: XVSPHE4 = .1
132 REAL, PARAMETER :: XVTAIL1 = 1.28E-17
133 REAL, PARAMETER :: XVTAIL2 = 4.22E-19
134 REAL, PARAMETER :: XVGRAT1 = 5.
135 REAL, PARAMETER :: XVGRAT2 = 15.
136 REAL, PARAMETER :: XVFI = 1.0417E-9
137 REAL, PARAMETER :: XVTELV1 = 0.005
138 !
139 INTEGER,PARAMETER :: NVDENT1 = 3
140 !
141 INTEGER :: NVARDIMS !number of dimensions of netcdf input variable
142 INTEGER :: NLENDIM1,NLENDIM2,NLENDIM3
143 INTEGER :: NID_VAR ! Netcdf IDs for variable
144 !
145 INTEGER :: NID_FILE
146 REAL, DIMENSION(:,:,:), POINTER :: XDRDT0,XTAU,XKAPPA ! field read
147 !
148 END MODULE modd_snow_metamo
149 
150 
151 
modd_snow_metamo
Definition: modd_snow_metamo.F90:8
Generated on Fri May 13 2016 09:31:04 for SURFEX by   doxygen 1.8.5