Changeset 13056 for utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/phycst.F90

Timestamp:

2020-06-07T18:26:09+02:00 (4 years ago)

Author:

rblod

Message:

ticket #2129 : cleaning domcfg

File:

• utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/phycst.F90 (modified) (2 diffs)

utils/tools_dev_r12970_AGRIF_CMEMS/DOMAINcfg/src/phycst.F90

@@ -37 +37 @@
    REAL(wp), PUBLIC ::   ra       = 6371229._wp      !: earth radius                       [m]
    REAL(wp), PUBLIC ::   grav     = 9.80665_wp       !: gravity                            [m/s2]   
-   REAL(wp), PUBLIC ::   rt0      = 273.15_wp        !: freezing point of fresh water [Kelvin]
 
-   REAL(wp), PUBLIC ::   rau0                        !: volumic mass of reference     [kg/m3]
-   REAL(wp), PUBLIC ::   r1_rau0                     !: = 1. / rau0                   [m3/kg]
-   REAL(wp), PUBLIC ::   rcp                         !: ocean specific heat           [J/Kelvin]
-   REAL(wp), PUBLIC ::   r1_rcp                      !: = 1. / rcp                    [Kelvin/J]
-   REAL(wp), PUBLIC ::   rau0_rcp                    !: = rau0 * rcp 
-   REAL(wp), PUBLIC ::   r1_rau0_rcp                 !: = 1. / ( rau0 * rcp )
-
-   REAL(wp), PUBLIC ::   emic     =    0.97_wp       !: emissivity of snow or ice (not used?)
-
-   REAL(wp), PUBLIC ::   sice     =    6.0_wp        !: salinity of ice (for pisces)          [psu]
-   REAL(wp), PUBLIC ::   soce     =   34.7_wp        !: salinity of sea (for pisces and isf)  [psu]
-   REAL(wp), PUBLIC ::   rLevap   =    2.5e+6_wp     !: latent heat of evaporation (water)
-   REAL(wp), PUBLIC ::   vkarmn   =    0.4_wp        !: von Karman constant
-   REAL(wp), PUBLIC ::   stefan   =    5.67e-8_wp    !: Stefan-Boltzmann constant 
-
-   REAL(wp), PUBLIC ::   rhos     =  330._wp         !: volumic mass of snow                                  [kg/m3]
-   REAL(wp), PUBLIC ::   rhoi     =  917._wp         !: volumic mass of sea ice                               [kg/m3]
-   REAL(wp), PUBLIC ::   rhow     = 1000._wp         !: volumic mass of freshwater in melt ponds              [kg/m3]
-   REAL(wp), PUBLIC ::   rcnd_i   =    2.034396_wp   !: thermal conductivity of fresh ice                     [W/m/K]
-   REAL(wp), PUBLIC ::   rcpi     = 2067.0_wp        !: specific heat of fresh ice                            [J/kg/K]
-   REAL(wp), PUBLIC ::   rLsub    =    2.834e+6_wp   !: pure ice latent heat of sublimation                   [J/kg]
-   REAL(wp), PUBLIC ::   rLfus    =    0.334e+6_wp   !: latent heat of fusion of fresh ice                    [J/kg]
-   REAL(wp), PUBLIC ::   rTmlt    =    0.054_wp      !: decrease of seawater meltpoint with salinity
-
-   REAL(wp), PUBLIC ::   r1_rhoi                     !: 1 / rhoi
-   REAL(wp), PUBLIC ::   r1_rhos                     !: 1 / rhos
-   REAL(wp), PUBLIC ::   r1_rcpi                     !: 1 / rcpi
    !!----------------------------------------------------------------------
    !! NEMO/OCE 4.0 , NEMO Consortium (2018)
@@ -89 +61 @@
 #endif
 
-      r1_rhoi = 1._wp / rhoi
-      r1_rhos = 1._wp / rhos
-      r1_rcpi = 1._wp / rcpi
-
-      IF(lwp) THEN
-         WRITE(numout,*)
-         WRITE(numout,*) 'phy_cst : initialization of ocean parameters and constants'
-         WRITE(numout,*) '~~~~~~~'
-         WRITE(numout,*) '      mathematical constant                 rpi = ', rpi
-         WRITE(numout,*) '      day                                rday   = ', rday,   ' s'
-         WRITE(numout,*) '      sideral year                       rsiyea = ', rsiyea, ' s'
-         WRITE(numout,*) '      sideral day                        rsiday = ', rsiday, ' s'
-         WRITE(numout,*) '      omega                              omega  = ', omega,  ' s^-1'
-         WRITE(numout,*)
-         WRITE(numout,*) '      nb of months per year               raamo = ', raamo, ' months'
-         WRITE(numout,*) '      nb of hours per day                 rjjhh = ', rjjhh, ' hours'
-         WRITE(numout,*) '      nb of minutes per hour              rhhmm = ', rhhmm, ' mn'
-         WRITE(numout,*) '      nb of seconds per minute            rmmss = ', rmmss, ' s'
-         WRITE(numout,*)
-         WRITE(numout,*) '      earth radius                         ra   = ', ra   , ' m'
-         WRITE(numout,*) '      gravity                              grav = ', grav , ' m/s^2'
-         WRITE(numout,*)
-         WRITE(numout,*) '      freezing point of water              rt0  = ', rt0  , ' K'
-         WRITE(numout,*)
-         WRITE(numout,*) '   reference density and heat capacity now defined in eosbn2.f90'
-         WRITE(numout,*)
-         WRITE(numout,*) '      thermal conductivity of pure ice          = ', rcnd_i  , ' J/s/m/K'
-         WRITE(numout,*) '      thermal conductivity of snow is defined in a namelist '
-         WRITE(numout,*) '      fresh ice specific heat                   = ', rcpi    , ' J/kg/K'
-         WRITE(numout,*) '      latent heat of fusion of fresh ice / snow = ', rLfus   , ' J/kg'
-         WRITE(numout,*) '      latent heat of subl.  of fresh ice / snow = ', rLsub   , ' J/kg'
-         WRITE(numout,*) '      density of sea ice                        = ', rhoi    , ' kg/m^3'
-         WRITE(numout,*) '      density of snow                           = ', rhos    , ' kg/m^3'
-         WRITE(numout,*) '      density of freshwater (in melt ponds)     = ', rhow    , ' kg/m^3'
-         WRITE(numout,*) '      salinity of ice (for pisces)              = ', sice    , ' psu'
-         WRITE(numout,*) '      salinity of sea (for pisces and isf)      = ', soce    , ' psu'
-         WRITE(numout,*) '      latent heat of evaporation (water)        = ', rLevap  , ' J/m^3' 
-         WRITE(numout,*) '      von Karman constant                       = ', vkarmn 
-         WRITE(numout,*) '      Stefan-Boltzmann constant                 = ', stefan  , ' J/s/m^2/K^4'
-         WRITE(numout,*)
-         WRITE(numout,*) '      conversion: degre ==> radian          rad = ', rad
-         WRITE(numout,*)
-         WRITE(numout,*) '      smallest real computer value       rsmall = ', rsmall
-      ENDIF
-
    END SUBROUTINE phy_cst