Changeset 13655 for NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/src/OCE/DOM/phycst.F90
- Timestamp:
- 2020-10-21T16:15:13+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r13648_ASINTER-04_laurent_bulk_ice/src/OCE/DOM/phycst.F90
r12489 r13655 7 7 !! 8.1 ! 1991-11 (G. Madec, M. Imbard) cosmetic changes 8 8 !! NEMO 1.0 ! 2002-08 (G. Madec, C. Ethe) F90, add ice constants 9 !! - ! 2006-08 (G. Madec) style 10 !! 3.2 ! 2006-08 (S. Masson, G. Madec) suppress useless variables + style 11 !! 3.4 ! 2011-11 (C. Harris) minor changes for CICE constants 9 !! - ! 2006-08 (G. Madec) style 10 !! 3.2 ! 2006-08 (S. Masson, G. Madec) suppress useless variables + style 11 !! 3.4 ! 2011-11 (C. Harris) minor changes for CICE constants 12 12 !!---------------------------------------------------------------------- 13 13 … … 26 26 REAL(wp), PUBLIC :: rad = 3.141592653589793_wp / 180._wp !: conversion from degre into radian 27 27 REAL(wp), PUBLIC :: rsmall = 0.5 * EPSILON( 1.e0 ) !: smallest real computer value 28 28 29 29 REAL(wp), PUBLIC :: rday = 24.*60.*60. !: day [s] 30 30 REAL(wp), PUBLIC :: rsiyea !: sideral year [s] … … 36 36 REAL(wp), PUBLIC :: omega !: earth rotation parameter [s-1] 37 37 REAL(wp), PUBLIC :: ra = 6371229._wp !: earth radius [m] 38 REAL(wp), PUBLIC :: grav = 9.80665_wp !: gravity [m/s2] 38 REAL(wp), PUBLIC :: grav = 9.80665_wp !: gravity [m/s2] 39 39 REAL(wp), PUBLIC :: rt0 = 273.15_wp !: freezing point of fresh water [Kelvin] 40 40 … … 43 43 REAL(wp), PUBLIC :: rcp !: ocean specific heat [J/Kelvin] 44 44 REAL(wp), PUBLIC :: r1_rcp !: = 1. / rcp [Kelvin/J] 45 REAL(wp), PUBLIC :: rho0_rcp !: = rho0 * rcp 45 REAL(wp), PUBLIC :: rho0_rcp !: = rho0 * rcp 46 46 REAL(wp), PUBLIC :: r1_rho0_rcp !: = 1. / ( rho0 * rcp ) 47 47 … … 52 52 REAL(wp), PUBLIC :: rLevap = 2.5e+6_wp !: latent heat of evaporation (water) 53 53 REAL(wp), PUBLIC :: vkarmn = 0.4_wp !: von Karman constant 54 REAL(wp), PUBLIC :: stefan = 5.67e-8_wp !: Stefan-Boltzmann constant 54 REAL(wp), PUBLIC :: vkarmn2 = 0.4_wp*0.4_wp !: square of von Karman constant 55 REAL(wp), PUBLIC :: stefan = 5.67e-8_wp !: Stefan-Boltzmann constant 55 56 56 57 REAL(wp), PUBLIC :: rhos = 330._wp !: volumic mass of snow [kg/m3] … … 68 69 !!---------------------------------------------------------------------- 69 70 !! NEMO/OCE 4.0 , NEMO Consortium (2018) 70 !! $Id$ 71 !! $Id$ 71 72 !! Software governed by the CeCILL license (see ./LICENSE) 72 73 !!---------------------------------------------------------------------- 73 74 74 75 CONTAINS 75 76 76 77 SUBROUTINE phy_cst 77 78 !!---------------------------------------------------------------------- … … 86 87 omega = 7.292116e-05 87 88 #else 88 omega = 2._wp * rpi / rsiday 89 omega = 2._wp * rpi / rsiday 89 90 #endif 90 91 … … 125 126 WRITE(numout,*) ' salinity of ice (for pisces) = ', sice , ' psu' 126 127 WRITE(numout,*) ' salinity of sea (for pisces and isf) = ', soce , ' psu' 127 WRITE(numout,*) ' latent heat of evaporation (water) = ', rLevap , ' J/m^3' 128 WRITE(numout,*) ' von Karman constant = ', vkarmn 128 WRITE(numout,*) ' latent heat of evaporation (water) = ', rLevap , ' J/m^3' 129 WRITE(numout,*) ' von Karman constant = ', vkarmn 129 130 WRITE(numout,*) ' Stefan-Boltzmann constant = ', stefan , ' J/s/m^2/K^4' 130 131 WRITE(numout,*)
Note: See TracChangeset
for help on using the changeset viewer.