modd_glt_const_evp.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.
!GLT_LIC The GELATO model is a seaice model used in stand-alone or embedded mode. 
!GLT_LIC  It has been developed by Meteo-France. The holder of GELATO is Meteo-France.
!GLT_LIC  
!GLT_LIC  This software is governed by the CeCILL-C license under French law and biding
!GLT_LIC  by the rules of distribution of free software. See the CeCILL-C_V1-en.txt
!GLT_LIC  (English) and CeCILL-C_V1-fr.txt (French) for details. The CeCILL is a free
!GLT_LIC  software license, explicitly compatible with the GNU GPL
!GLT_LIC  (see http://www.gnu.org/licenses/license-list.en.html#CeCILL)
!GLT_LIC  
!GLT_LIC  The CeCILL-C licence agreement grants users the right to modify and re-use the
!GLT_LIC  software governed by this free software license. The exercising of this right
!GLT_LIC  is conditional upon the obligation to make available to the community the
!GLT_LIC  modifications made to the source code of the software so as to contribute to
!GLT_LIC  its evolution.
!GLT_LIC  
!GLT_LIC  In consideration of access to the source code and the rights to copy, modify
!GLT_LIC  and redistribute granted by the license, users are provided only with a limited
!GLT_LIC  warranty and the software's author, the holder of the economic rights, and the
!GLT_LIC  successive licensors only have limited liability. In this respect, the risks
!GLT_LIC  associated with loading, using, modifying and/or developing or reproducing the
!GLT_LIC  software by the user are brought to the user's attention, given its Free
!GLT_LIC  Software status, which may make it complicated to use, with the result that its
!GLT_LIC  use is reserved for developers and experienced professionals having in-depth
!GLT_LIC  computer knowledge. Users are therefore encouraged to load and test the
!GLT_LIC  suitability of the software as regards their requirements in conditions enabling
!GLT_LIC  the security of their systems and/or data to be ensured and, more generally, to
!GLT_LIC  use and operate it in the same conditions of security. 
!GLT_LIC  
!GLT_LIC  The GELATO sofware is cureently distibuted with the SURFEX software, available at 
!GLT_LIC  http://www.cnrm.meteo.fr/surfex. The fact that you download the software deemed that
!GLT_LIC  you had knowledge of the CeCILL-C license and that you accept its terms.
!GLT_LIC  Attempts to use this software in a way not complying with CeCILL-C license
!GLT_LIC  may lead to prosecution. 
!GLT_LIC 
MODULE modd_glt_const_evp
  USE modd_glt_param 
  IMPLICIT NONE

#if ! defined in_surfex
! *** Constants for EVP rheology.
!
  REAL, PARAMETER :: eyc = 0.36      ! E0
  REAL, PARAMETER :: cst = 20.
  REAL, PARAMETER :: pst = 2.75e4    ! N.m-2  [ kg.m-1.s-2 ]

!.. densities of ice, seawater
  REAL, PARAMETER :: rhoi = 910.       ! Hibler 1979, kg/m^3
  REAL, PARAMETER :: rhowdyn = 1030.   ! Pond and Pickard, S=35 ppt, T=10.C, p=0

! .. EVP parameters
  REAL, PARAMETER :: ecc = 4.              ! e^2
  REAL, PARAMETER :: ecci = 1.             ! 4./e^2
  REAL, PARAMETER :: eccm = 2.-0.5*ecci    ! 2*(1-0.25*ecci) = 2*(1-1/e^2) 
  REAL, PARAMETER :: eccp = 1.+ecci*0.25   ! 1+1/e^2

!.. water and air turning angle (converted to radians).
! If geostrophic current is provided, phiwater=25 degree, 0 if surface current
!  REAL, PARAMETER :: phiwater = 25.*3.141592653789/180. 
  REAL, PARAMETER :: phiwater = 0.
  REAL, PARAMETER :: phiair = 0.

!.. ice/water drag coefficient
  REAL, PARAMETER :: dragw = 0.0055*rhowdyn
#endif

!.. earth rotational velocity, rad.s-1 
  REAL, PARAMETER :: omega = 7.292e-5
 
#if ! defined in_surfex
! *** Other constants.
  REAL, PARAMETER :: floor1 = 1.e-5
  REAL, PARAMETER :: floor2 = 1.e-11
  REAL, PARAMETER :: rcon = 4.e-9

!.. time step for sea ice dynamics
  REAL :: dt

!.. time step for sea ice dynamics subcycling
  REAL :: dte
  
!.. one over dte
  REAL :: dtei

!.. wave damping time scale T
  REAL :: tdamp2

!.. constants for stress equation
  REAL :: Tdtei
  REAL :: c2n
  REAL :: ecoef
  REAL :: d2n
  REAL :: ccoef
  REAL :: dcoef
  REAL :: bcoef
  REAL :: acoef

!.. cosine, sine of water and atmosphere turning angles
  REAL :: xcw,xsw,xca,xsa
#endif
!
END MODULE modd_glt_const_evp