Changeset 9935 for NEMO/trunk/src/OCE/DOM/phycst.F90
- Timestamp:
- 2018-07-12T16:12:48+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/DOM/phycst.F90
r9929 r9935 46 46 REAL(wp), PUBLIC :: r1_rau0_rcp !: = 1. / ( rau0 * rcp ) 47 47 48 !clem: not sure these are needed for cice 49 #if defined key_cice 50 REAL(wp), PUBLIC :: rt0_snow = 273.15_wp !: melting point of snow [Kelvin] 51 REAL(wp), PUBLIC :: rt0_ice = 273.05_wp !: melting point of ice [Kelvin] 52 REAL(wp), PUBLIC :: emic = 0.97_wp !: emissivity of snow or ice 53 REAL(wp), PUBLIC :: xlsn !: = lfus*rhosn (volumetric latent heat fusion of snow) [J/m3] 54 #endif 48 REAL(wp), PUBLIC :: emic = 0.97_wp !: emissivity of snow or ice (not used?) 55 49 56 50 REAL(wp), PUBLIC :: sice = 6.0_wp !: salinity of ice (for pisces) [psu] 57 51 REAL(wp), PUBLIC :: soce = 34.7_wp !: salinity of sea (for pisces and isf) [psu] 58 REAL(wp), PUBLIC :: cevap = 2.5e+6_wp !: latent heat of evaporation (water) 59 REAL(wp), PUBLIC :: srgamma = 0.9_wp !: correction factor for solar radiation (Oberhuber, 1974) 52 REAL(wp), PUBLIC :: rLevap = 2.5e+6_wp !: latent heat of evaporation (water) 60 53 REAL(wp), PUBLIC :: vkarmn = 0.4_wp !: von Karman constant 61 54 REAL(wp), PUBLIC :: stefan = 5.67e-8_wp !: Stefan-Boltzmann constant 62 55 63 REAL(wp), PUBLIC :: rhosn = 330._wp !: volumic mass of snow [kg/m3] 64 REAL(wp), PUBLIC :: rhoic = 917._wp !: volumic mass of sea ice [kg/m3] 65 REAL(wp), PUBLIC :: rhofw = 1000._wp !: volumic mass of freshwater in melt ponds [kg/m3] 66 REAL(wp), PUBLIC :: rcdic = 2.034396_wp !: thermal conductivity of fresh ice [W/m/K] 67 #if defined key_cice 68 REAL(wp), PUBLIC :: rcdsn = 0.31_wp !: thermal conductivity of snow [W/m/K] 69 #endif 70 REAL(wp), PUBLIC :: cpic = 2067.0_wp !: specific heat of fresh ice [J/kg/K] 71 REAL(wp), PUBLIC :: lsub = 2.834e+6_wp !: pure ice latent heat of sublimation [J/kg] 72 REAL(wp), PUBLIC :: lfus = 0.334e+6_wp !: latent heat of fusion of fresh ice [J/kg] 73 REAL(wp), PUBLIC :: tmut = 0.054_wp !: decrease of seawater meltpoint with salinity 56 REAL(wp), PUBLIC :: rhos = 330._wp !: volumic mass of snow [kg/m3] 57 REAL(wp), PUBLIC :: rhoi = 917._wp !: volumic mass of sea ice [kg/m3] 58 REAL(wp), PUBLIC :: rhow = 1000._wp !: volumic mass of freshwater in melt ponds [kg/m3] 59 REAL(wp), PUBLIC :: rcnd_i = 2.034396_wp !: thermal conductivity of fresh ice [W/m/K] 60 REAL(wp), PUBLIC :: rcpi = 2067.0_wp !: specific heat of fresh ice [J/kg/K] 61 REAL(wp), PUBLIC :: rLsub = 2.834e+6_wp !: pure ice latent heat of sublimation [J/kg] 62 REAL(wp), PUBLIC :: rLfus = 0.334e+6_wp !: latent heat of fusion of fresh ice [J/kg] 63 REAL(wp), PUBLIC :: rTmlt = 0.054_wp !: decrease of seawater meltpoint with salinity 74 64 75 REAL(wp), PUBLIC :: r1_rhoi c !: 1 / rhoic76 REAL(wp), PUBLIC :: r1_rhos n !: 1 / rhosn77 REAL(wp), PUBLIC :: r1_ cpic !: 1 / cpic65 REAL(wp), PUBLIC :: r1_rhoi !: 1 / rhoi 66 REAL(wp), PUBLIC :: r1_rhos !: 1 / rhos 67 REAL(wp), PUBLIC :: r1_rcpi !: 1 / rcpi 78 68 !!---------------------------------------------------------------------- 79 69 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 99 89 #endif 100 90 101 #if defined key_cice 102 xlsn = lfus * rhosn 103 #endif 104 105 r1_rhoic = 1._wp / rhoic 106 r1_rhosn = 1._wp / rhosn 107 r1_cpic = 1._wp / cpic 91 r1_rhoi = 1._wp / rhoi 92 r1_rhos = 1._wp / rhos 93 r1_rcpi = 1._wp / rcpi 108 94 109 95 IF(lwp) THEN … … 129 115 WRITE(numout,*) ' reference density and heat capacity now defined in eosbn2.f90' 130 116 WRITE(numout,*) 131 #if defined key_cice 132 WRITE(numout,*) ' thermal conductivity of the snow = ', rcdsn , ' J/s/m/K' 133 #endif 134 WRITE(numout,*) ' thermal conductivity of pure ice = ', rcdic , ' J/s/m/K' 135 WRITE(numout,*) ' fresh ice specific heat = ', cpic , ' J/kg/K' 136 WRITE(numout,*) ' latent heat of fusion of fresh ice / snow = ', lfus , ' J/kg' 137 WRITE(numout,*) ' latent heat of subl. of fresh ice / snow = ', lsub , ' J/kg' 138 WRITE(numout,*) ' density of sea ice = ', rhoic , ' kg/m^3' 139 WRITE(numout,*) ' density of snow = ', rhosn , ' kg/m^3' 140 WRITE(numout,*) ' density of freshwater (in melt ponds) = ', rhofw , ' kg/m^3' 117 WRITE(numout,*) ' thermal conductivity of pure ice = ', rcnd_i , ' J/s/m/K' 118 WRITE(numout,*) ' thermal conductivity of snow is defined in a namelist ' 119 WRITE(numout,*) ' fresh ice specific heat = ', rcpi , ' J/kg/K' 120 WRITE(numout,*) ' latent heat of fusion of fresh ice / snow = ', rLfus , ' J/kg' 121 WRITE(numout,*) ' latent heat of subl. of fresh ice / snow = ', rLsub , ' J/kg' 122 WRITE(numout,*) ' density of sea ice = ', rhoi , ' kg/m^3' 123 WRITE(numout,*) ' density of snow = ', rhos , ' kg/m^3' 124 WRITE(numout,*) ' density of freshwater (in melt ponds) = ', rhow , ' kg/m^3' 141 125 WRITE(numout,*) ' salinity of ice (for pisces) = ', sice , ' psu' 142 126 WRITE(numout,*) ' salinity of sea (for pisces and isf) = ', soce , ' psu' 143 WRITE(numout,*) ' latent heat of evaporation (water) = ', cevap , ' J/m^3' 144 WRITE(numout,*) ' correction factor for solar radiation = ', srgamma 127 WRITE(numout,*) ' latent heat of evaporation (water) = ', rLevap , ' J/m^3' 145 128 WRITE(numout,*) ' von Karman constant = ', vkarmn 146 129 WRITE(numout,*) ' Stefan-Boltzmann constant = ', stefan , ' J/s/m^2/K^4'
Note: See TracChangeset
for help on using the changeset viewer.