Home Explore Help
Sign In
fmasson
/
forcing-perturbation
1
0
Fork 0
Files Issues 0 Pull Requests 0 Wiki
Branch: main
Branches Tags
forcing-pert...
/
interpolation
/
rotateUVorca_sources
/
phycst.f90

phycst.f90 5.3 KB
Permalink History Raw

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889909192939495 
  1. MODULE phycst
    2. 

  2.    !!======================================================================
    3. 

  3.    !!                    ***  MODULE  phycst  ***
    4. 

  4.    !!     Definition of of both ocean and ice parameters used in the code
    5. 

  5.    !!=====================================================================
    6. 

  6.    !! History :   OPA  !  1990-10  (C. Levy - G. Madec)  Original code
    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. 

  9.    !!              -   !  2006-08  (G. Madec)  style 
    10. 

  10.    !!             3.2  !  2006-08  (S. Masson, G. Madec)  suppress useless variables + style 
    11. 

  11.    !!----------------------------------------------------------------------
    12. 

  12. 

  13.    !!----------------------------------------------------------------------
    14. 

  14.    !!   phy_cst  : define and print physical constant and domain parameters
    15. 

  15.    !!----------------------------------------------------------------------
    16. 

  16. 

  17.    use par_kind 
    18. 

  18. 

  19.    IMPLICIT NONE
    20. 

  20.    PRIVATE
    21. 

  21. 

  22.    PUBLIC   phy_cst     ! routine called by inipar.F90
    23. 

  23. 

  24.    REAL(wp), PUBLIC ::   rpi = 3.141592653589793_wp             !: pi
    25. 

  25.    REAL(wp), PUBLIC ::   rad = 3.141592653589793_wp / 180._wp   !: conversion from degre into radian
    26. 

  26.    REAL(wp), PUBLIC ::   rsmall = 0.5 * EPSILON( 1.e0 )         !: smallest real computer value
    27. 

  27.    
    28. 

  28.    REAL(wp), PUBLIC ::   rday = 24.*60.*60.       !: day (s)
    29. 

  29.    REAL(wp), PUBLIC ::   rsiyea                   !: sideral year (s)
    30. 

  30.    REAL(wp), PUBLIC ::   rsiday                   !: sideral day (s)
    31. 

  31.    REAL(wp), PUBLIC ::   raamo =  12._wp          !: number of months in one year
    32. 

  32.    REAL(wp), PUBLIC ::   rjjhh =  24._wp          !: number of hours in one day
    33. 

  33.    REAL(wp), PUBLIC ::   rhhmm =  60._wp          !: number of minutes in one hour
    34. 

  34.    REAL(wp), PUBLIC ::   rmmss =  60._wp          !: number of seconds in one minute
    35. 

  35. !! REAL(wp), PUBLIC ::   omega = 7.292115083046061e-5_wp ,  &  !: change the last digit!
    36. 

  36.    REAL(wp), PUBLIC ::   omega                    !: earth rotation parameter
    37. 

  37.    REAL(wp), PUBLIC ::   ra    = 6371229._wp      !: earth radius (meter)
    38. 

  38.    REAL(wp), PUBLIC ::   grav  = 9.80665_wp       !: gravity (m/s2)
    39. 

  39.    
    40. 

  40.    REAL(wp), PUBLIC ::   rtt      = 273.16_wp     !: triple point of temperature (Kelvin)
    41. 

  41.    REAL(wp), PUBLIC ::   rt0      = 273.15_wp     !: freezing point of water (Kelvin)
    42. 

  42.    REAL(wp), PUBLIC ::   rt0_snow = 273.15_wp     !: melting point of snow  (Kelvin)
    43. 

  43.    REAL(wp), PUBLIC ::   rt0_ice  = 273.05_wp     !: melting point of ice   (Kelvin)
    44. 

  44. 

  45.    REAL(wp), PUBLIC ::   rau0     = 1020._wp      !: reference volumic mass (density)  (kg/m3)
    46. 

  46.    REAL(wp), PUBLIC ::   rau0r                    !: reference specific volume         (m3/kg)
    47. 

  47.    REAL(wp), PUBLIC ::   rcp      =    4.e+3_wp   !: ocean specific heat
    48. 

  48.    REAL(wp), PUBLIC ::   ro0cpr                   !: = 1. / ( rau0 * rcp )
    49. 

  49. 

  50.    REAL(wp), PUBLIC ::   rcdsn   =   0.22_wp      !: conductivity of the snow
    51. 

  51.    REAL(wp), PUBLIC ::   rcdic   =   2.034396_wp  !: conductivity of the ice
    52. 

  52.    REAL(wp), PUBLIC ::   rcpsn   =   6.9069e+5_wp !: density times specific heat for snow
    53. 

  53.    REAL(wp), PUBLIC ::   rcpic   =   1.8837e+6_wp !: volumetric latent heat fusion of sea ice
    54. 

  54.    REAL(wp), PUBLIC ::   xlsn    = 110.121e+6_wp  !: volumetric latent heat fusion of snow
    55. 

  55.    REAL(wp), PUBLIC ::   xlic    = 300.33e+6_wp   !: volumetric latent heat fusion of ice
    56. 

  56.    REAL(wp), PUBLIC ::   xsn     =   2.8e+6       !: latent heat of sublimation of snow
    57. 

  57.    REAL(wp), PUBLIC ::   rhoic   = 900._wp        !: volumic mass of sea ice (kg/m3)
    58. 

  58.    REAL(wp), PUBLIC ::   rhosn   = 330._wp        !: volumic mass of snow (kg/m3)
    59. 

  59.    REAL(wp), PUBLIC ::   emic    =   0.97_wp      !: emissivity of snow or ice
    60. 

  60.    REAL(wp), PUBLIC ::   sice    =   6.0_wp       !: reference salinity of ice (psu)
    61. 

  61.    REAL(wp), PUBLIC ::   soce    =  34.7_wp       !: reference salinity of sea (psu)
    62. 

  62.    REAL(wp), PUBLIC ::   cevap   =   2.5e+6_wp    !: latent heat of evaporation (water)
    63. 

  63.    REAL(wp), PUBLIC ::   srgamma =   0.9_wp       !: correction factor for solar radiation (Oberhuber, 1974)
    64. 

  64.    REAL(wp), PUBLIC ::   vkarmn  =   0.4_wp       !: von Karman constant
    65. 

  65.    REAL(wp), PUBLIC ::   stefan  =   5.67e-8_wp   !: Stefan-Boltzmann constant 
    66. 

  66.    !!----------------------------------------------------------------------
    67. 

  67.    !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009) 
    68. 

  68.    !! $Id: phycst.F90 1932 2010-06-15 10:28:20Z smasson $ 
    69. 

  69.    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    70. 

  70.    !!----------------------------------------------------------------------
    71. 

  71.    
    72. 

  72. CONTAINS
    73. 

  73.    
    74. 

  74.    SUBROUTINE phy_cst
    75. 

  75.       !!----------------------------------------------------------------------
    76. 

  76.       !!                       ***  ROUTINE phy_cst  ***
    77. 

  77.       !!
    78. 

  78.       !! ** Purpose :   Print model parameters and set and print the constants
    79. 

  79.       !!----------------------------------------------------------------------
    80. 

  80.       CHARACTER (len=64) ::   cform = "(A12, 3(A13, I7) )" 
    81. 

  81.       !!----------------------------------------------------------------------
    82. 

  82. 

  83.       !                                   ! Define additional parameters
    84. 

  84.       rsiyea = 365.25 * rday * 2. * rpi / 6.283076
    85. 

  85.       rsiday = rday / ( 1. + rday / rsiyea )
    86. 

  86.       omega  = 2. * rpi / rsiday 
    87. 

  87. 

  88.       rau0r  = 1. /   rau0  
    89. 

  89.       ro0cpr = 1. / ( rau0 * rcp )
    90. 

  90. 

  91. 

  92.    END SUBROUTINE phy_cst
    93. 

  93. 

  94.    !!======================================================================
    95. 

  95. END MODULE phycst
    96.