modd_snow_metamo.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
! ajoutEB
! correction de l'erreur interversion de XVTANG2 et XVTANG3
!######################
      MODULE modd_snow_metamo
!     ######################
!
!!****  *MODD_SNOW_METAMO* - declaration of parameters related
!!                          to snow metamorphism!!
!!    PURPOSE
!!    -------
!       The purpose of this declarative module is to specify  the 
!     parameters related to the metamorphism parameterization of snow.
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!      
!!
!!    AUTHOR
!!    ------
!!      V. Vionnet   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original       02/2008                
!-------------------------------------------------------------------------------
!
!*       0.   DECLARATIONS
!             ------------
!
IMPLICIT NONE
!
!-------------------------------------------------------------------------------
!
! minimum snow layer thickness for thermal calculations.
! Used to prevent numerical problems as snow becomes vanishingly thin.
REAL, PARAMETER                 :: xsnowdzmin = 0.0001
!
! Optical diameter properties
REAL, PARAMETER                 :: xdiaet = 1.e-4
REAL, PARAMETER                 :: xdiagf = 3.e-4
REAL, PARAMETER                 :: xdiafp = 4.e-4
!
! Compaction/Settling Coefficients from Crocus v2.4
!
REAL, PARAMETER        :: xvvisc1 = 7.62237e6   ! pre-exponential viscosity factor (UNIT : N m-2 s)
REAL, PARAMETER        :: xvvisc3 = 0.023       ! density adjustement in the exponential correction for viscosity (UNIT : m3 kg-1)
REAL, PARAMETER        :: xvvisc4 = .1          ! temperature adjustement in the exponential correction for viscosity (UNIT : K-1)
REAL, PARAMETER        :: xvvisc5 = 1.          ! factor for viscosity adjustement to grain type - to be checked
REAL, PARAMETER        :: xvvisc6 = 60.         ! factor for viscosity adjustement to grain type - to be checked 
!                                                             (especially this one ; inconsistency with Crocus v2.4)
REAL, PARAMETER        :: xvvisc7 = 10.         ! factor for viscosity adjustement to grain type - to be checked
REAL, PARAMETER        :: xvro11 = 250.  ! normalization term for density dependence of the viscosity calculation (UNIT : kg m-3)
! 
! Maximum value for TPSNOW%GRAN2
REAL, PARAMETER                 :: xvgran1 = 99.
REAL, PARAMETER                 :: xgran = 99.
!
INTEGER, PARAMETER            :: nvhis1 = 1                
INTEGER, PARAMETER            :: nvhis2 = 2                
INTEGER, PARAMETER            :: nvhis3 = 3                
INTEGER, PARAMETER            :: nvhis4 = 4                
INTEGER, PARAMETER            :: nvhis5 = 5 
!
! Properties of fresh snow
REAL, PARAMETER                 :: xnden1 = 17.12
REAL, PARAMETER                 :: xnden2 = 128.
REAL, PARAMETER                 :: xnden3 = -20.
REAL, PARAMETER                 :: xnsph1 = 7.87
REAL, PARAMETER                 :: xnsph2 = 38.
REAL, PARAMETER                 :: xnsph3 = 50.
REAL, PARAMETER                 :: xnsph4 = 90.
!
REAL, PARAMETER                 :: xuepsi = 1.e-8
REAL, PARAMETER                 :: xepsi = 1.e-8
REAL, PARAMETER                 :: xupourc = 100.
!
! Parameters for Marbouty's function
!
REAL, PARAMETER                 :: xvtang1 = 40.              
REAL, PARAMETER                 :: xvtang2 = 6.              
REAL, PARAMETER                 :: xvtang3 = 22.              
REAL, PARAMETER                 :: xvtang4 = .7               
REAL, PARAMETER                 :: xvtang5 = .3               
REAL, PARAMETER                 :: xvtang6 = 6.               
REAL, PARAMETER                 :: xvtang7 = 1.               
REAL, PARAMETER                 :: xvtang8 = .8               
REAL, PARAMETER                 :: xvtang9 = 16.              
REAL, PARAMETER                 :: xvtanga = .2               
REAL, PARAMETER                 :: xvtangb = .2               
REAL, PARAMETER                 :: xvtangc = 18.             
REAL, PARAMETER                 :: xvrang1 = 400.               
REAL, PARAMETER                 :: xvrang2 = 150.              
REAL, PARAMETER                 :: xvgang1 = 70.               
REAL, PARAMETER                 :: xvgang2 = 25.              
REAL, PARAMETER                 :: xvgang3 = 40.               
REAL, PARAMETER                 :: xvgang4 = 50.               
REAL, PARAMETER                 :: xvgang5 = .1               
REAL, PARAMETER                 :: xvgang6 = 15.              
REAL, PARAMETER                 :: xvgang7 = .1               
REAL, PARAMETER                 :: xvgang8 = .55              
REAL, PARAMETER                 :: xvgang9 = .65              
REAL, PARAMETER                 :: xvganga = .2               
REAL, PARAMETER                 :: xvgangb = .85              
REAL, PARAMETER                 :: xvgangc = .15      
!
! Parameters for snow metamorphism
!
REAL, PARAMETER                 :: xvdent1 = 2314.81481
REAL, PARAMETER                 :: xvdent2 = 7.2338e-7 
REAL, PARAMETER                 :: xvgran6 = 51.              
REAL, PARAMETER                 :: xvvap1 = -6000.            
REAL, PARAMETER                 :: xvvap2 = .4                
REAL, PARAMETER                 :: xvdiam1 = 4.e-4               
REAL, PARAMETER                 :: xvdiam2 = 5.e-4              
REAL, PARAMETER                 :: xvdiam3 = 3.e-4               
REAL, PARAMETER                 :: xvdiam4 = 2.e-4
REAL, PARAMETER                 :: xvdiam5 = 1.e-4
REAL, PARAMETER                 :: xvdiam6 = 1.e-4               
REAL, PARAMETER                 :: xvsphe1 = 1.               
REAL, PARAMETER                 :: xvsphe2 = 11574.074             
REAL, PARAMETER                 :: xvsphe3 = .5               
REAL, PARAMETER                 :: xvsphe4 = .1               
REAL, PARAMETER                 :: xvtail1 = 1.28e-17          
REAL, PARAMETER                 :: xvtail2 = 4.22e-19         
REAL, PARAMETER                 :: xvgrat1 = 5.              
REAL, PARAMETER                 :: xvgrat2 = 15.             
REAL, PARAMETER                 :: xvfi = 1.0417e-9                 
REAL, PARAMETER                 :: xvtelv1 = 0.005 
!
INTEGER,PARAMETER               :: nvdent1 = 3
!
INTEGER :: nvardims !number of dimensions of netcdf input variable
INTEGER :: nlendim1,nlendim2,nlendim3
INTEGER :: nid_var ! Netcdf IDs for  variable
!
INTEGER :: nid_file
REAL, DIMENSION(:,:,:), POINTER :: xdrdt0,xtau,xkappa   ! field read
!
END MODULE modd_snow_metamo