Changeset 88 for trunk/NEMO/OPA_SRC/phycst.F90
- Timestamp:
- 2004-04-22T15:50:27+02:00 (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/phycst.F90
r15 r88 91 91 !!---------------------------------------------------------------------- 92 92 !! * Local variables 93 CHARACTER (len=64) :: cform = "( 4(A13, I7) )"93 CHARACTER (len=64) :: cform = "(A9, 3(A13, I7) )" 94 94 !!---------------------------------------------------------------------- 95 95 … … 97 97 IF(lwp) WRITE(numout,*) ' phy_cst : initialization of ocean parameters and constants' 98 98 IF(lwp) WRITE(numout,*) ' ~~~~~~~' 99 IF(lwp) WRITE(numout,*)100 101 99 102 100 ! Ocean Parameters 103 101 ! ---------------- 104 102 IF(lwp) THEN 103 WRITE(numout,*) ' parameter file' 105 104 WRITE(numout,*) 106 WRITE(numout,*) ' parameter file' 107 WRITE(numout,*) ' ~~~~~~~~~~~~~~' 105 WRITE(numout,*) ' dimension of model' 106 WRITE(numout,*) ' Local domain Global domain Data domain ' 107 WRITE(numout,cform) ' ',' jpi : ', jpi, ' jpiglo : ', jpiglo, ' jpidta : ', jpidta 108 WRITE(numout,cform) ' ',' jpj : ', jpj, ' jpjglo : ', jpjglo, ' jpjdta : ', jpjdta 109 WRITE(numout,cform) ' ',' jpk : ', jpk, ' jpk : ', jpk , ' jpkdta : ', jpkdta 110 WRITE(numout,*) ' ',' jpij : ', jpij 111 WRITE(numout,*) 112 WRITE(numout,*) ' mpp local domain info (mpp)' 113 WRITE(numout,*) ' jpni : ', jpni, ' jpreci : ', jpreci 114 WRITE(numout,*) ' jpnj : ', jpnj, ' jprecj : ', jprecj 115 WRITE(numout,*) ' jpnij : ', jpnij 108 116 109 117 WRITE(numout,*) 110 WRITE(numout,*) ' dimension of model' 111 WRITE(numout,*) ' local domain Global domain Data domain ' 112 WRITE(numout,cform) ' jpi : ', jpi, ' jpim1 : ', jpim1, ' jpiglo : ', jpiglo, ' jpidta : ', jpidta 113 WRITE(numout,cform) ' jpj : ', jpj, ' jpjm1 : ', jpjm1, ' jpjglo : ', jpjglo, ' jpjdta : ', jpjdta 114 WRITE(numout,cform) ' jpk : ', jpk, ' jpkm1 : ', jpkm1, ' jpk : ', jpk , ' jpkdta : ', jpkdta 115 WRITE(numout,*) ' jpij : ', jpij 116 WRITE(numout,*) 117 WRITE(numout,*) ' mpp local domain info (mpp)' 118 WRITE(numout,*) ' jpni : ', jpni, ' jpreci : ', jpreci 119 WRITE(numout,*) ' jpnj : ', jpnj, ' jprecj : ', jprecj 120 WRITE(numout,*) ' jpnij : ', jpnij 121 122 WRITE(numout,*) 123 WRITE(numout,*) ' lateral domain boundary condition type : jperio = ', jperio 124 WRITE(numout,*) ' domain island (use in rigid-lid case) : jpisl = ', jpisl 125 WRITE(numout,*) ' jpnisl = ', jpnisl 118 WRITE(numout,*) ' lateral domain boundary condition type : jperio = ', jperio 119 WRITE(numout,*) ' domain island (use in rigid-lid case) : jpisl = ', jpisl 120 WRITE(numout,*) ' jpnisl = ', jpnisl 126 121 ENDIF 127 128 122 129 123 ! Define constants 130 124 ! ---------------- 131 125 IF(lwp) WRITE(numout,*) 132 IF(lwp) WRITE(numout,*) ' constants' 133 IF(lwp) WRITE(numout,*) ' ~~~~~~~~~' 126 IF(lwp) WRITE(numout,*) ' constants' 134 127 135 128 IF(lwp) WRITE(numout,*) 136 IF(lwp) WRITE(numout,*) ' mathematical constantrpi = ', rpi129 IF(lwp) WRITE(numout,*) ' mathematical constant rpi = ', rpi 137 130 138 131 rsiyea = 365.25 * rday * 2. * rpi / 6.283076 … … 140 133 omega = 2. * rpi / rsiday 141 134 IF(lwp) WRITE(numout,*) 142 IF(lwp) WRITE(numout,*) ' dayrday = ', rday, ' s'143 IF(lwp) WRITE(numout,*) ' sideral yearrsiyea = ', rsiyea, ' s'144 IF(lwp) WRITE(numout,*) ' sideral dayrsiday = ', rsiday, ' s'145 IF(lwp) WRITE(numout,*) ' omegaomega = ', omega, ' s-1'135 IF(lwp) WRITE(numout,*) ' day rday = ', rday, ' s' 136 IF(lwp) WRITE(numout,*) ' sideral year rsiyea = ', rsiyea, ' s' 137 IF(lwp) WRITE(numout,*) ' sideral day rsiday = ', rsiday, ' s' 138 IF(lwp) WRITE(numout,*) ' omega omega = ', omega, ' s-1' 146 139 147 140 rjjss = rjjhh * rhhmm * rmmss 148 141 IF(lwp) WRITE(numout,*) 149 IF(lwp) WRITE(numout,*) ' nb of months per yearraamo = ', raamo, ' months'150 IF(lwp) WRITE(numout,*) ' nb of hours per dayrjjhh = ', rjjhh, ' hours'151 IF(lwp) WRITE(numout,*) ' nb of minutes per hourrhhmm = ', rhhmm, ' mn'152 IF(lwp) WRITE(numout,*) ' nb of seconds per minutermmss = ', rmmss, ' s'153 IF(lwp) WRITE(numout,*) ' nb of seconds per dayrjjss = ', rjjss, ' s'142 IF(lwp) WRITE(numout,*) ' nb of months per year raamo = ', raamo, ' months' 143 IF(lwp) WRITE(numout,*) ' nb of hours per day rjjhh = ', rjjhh, ' hours' 144 IF(lwp) WRITE(numout,*) ' nb of minutes per hour rhhmm = ', rhhmm, ' mn' 145 IF(lwp) WRITE(numout,*) ' nb of seconds per minute rmmss = ', rmmss, ' s' 146 IF(lwp) WRITE(numout,*) ' nb of seconds per day rjjss = ', rjjss, ' s' 154 147 155 148 IF(lwp) WRITE(numout,*) 156 IF(lwp) WRITE(numout,*) ' earth radiusra = ', ra, ' m'157 IF(lwp) WRITE(numout,*) ' gravity grav = ', grav , ' m/s2'149 IF(lwp) WRITE(numout,*) ' earth radius ra = ', ra, ' m' 150 IF(lwp) WRITE(numout,*) ' gravity grav = ', grav , ' m/s^2' 158 151 159 152 IF(lwp) WRITE(numout,*) 160 IF(lwp) WRITE(numout,*) ' triple point of temperature rtt = ', rtt , ' K'161 IF(lwp) WRITE(numout,*) ' freezing point of water ( C)rt0 = ', rt0 , ' K'162 IF(lwp) WRITE(numout,*) ' melting point of snow rt0_snow = ', rt0_snow, ' K'163 IF(lwp) WRITE(numout,*) ' melting point of ice rt0_ice = ', rt0_ice , ' K'153 IF(lwp) WRITE(numout,*) ' triple point of temperature rtt = ', rtt , ' K' 154 IF(lwp) WRITE(numout,*) ' freezing point of water rt0 = ', rt0 , ' K' 155 IF(lwp) WRITE(numout,*) ' melting point of snow rt0_snow = ', rt0_snow, ' K' 156 IF(lwp) WRITE(numout,*) ' melting point of ice rt0_ice = ', rt0_ice , ' K' 164 157 165 ro0cpr = 1. / ( rau0 * rcp ) !158 ro0cpr = 1. / ( rau0 * rcp ) 166 159 IF(lwp) WRITE(numout,*) 167 IF(lwp) WRITE(numout,*) ' volumic mass of pure water (kg/m3) rauw = ', rauw, ' kg/m3'168 IF(lwp) WRITE(numout,*) ' volumic mass of reference (kg/m3) rau0 = ', rau0, ' kg/m3'169 IF(lwp) WRITE(numout,*) ' ocean specific heatrcp = ', rcp170 IF(lwp) WRITE(numout,*) ' 1. / ( rau0 * rcp ) = ro0cpr = ', ro0cpr160 IF(lwp) WRITE(numout,*) ' volumic mass of pure water rauw = ', rauw, ' kg/m^3' 161 IF(lwp) WRITE(numout,*) ' volumic mass of reference rau0 = ', rau0, ' kg/m^3' 162 IF(lwp) WRITE(numout,*) ' ocean specific heat rcp = ', rcp 163 IF(lwp) WRITE(numout,*) ' 1. / ( rau0 * rcp ) = ro0cpr = ', ro0cpr 171 164 172 165 IF(lwp) THEN 173 166 WRITE(numout,*) 174 WRITE(numout,*) ' thermal conductivity of the snow = ', rcdsn , ' J.s-1.m-1.K-1'175 WRITE(numout,*) ' thermal conductivity of the ice = ', rcdic , ' J.s-1.m-1.K-1'176 WRITE(numout,*) ' density times specific heat for snow = ', rcpsn , ' J.m-3.K-1'177 WRITE(numout,*) ' density times specific heat for ice = ', rcpic , ' J.m-3.K-1'178 WRITE(numout,*) ' volumetric latent heat fusion of sea ice = ', xlic , ' J.m-3'179 WRITE(numout,*) ' volumetric latent heat fusion of snow = ', xlsn , ' J.m-3'180 WRITE(numout,*) ' latent heat of sublimation of snow = ', xsn , ' J.kg-1'181 WRITE(numout,*) ' density of sea ice = ', rhoic , ' kg.m-3'182 WRITE(numout,*) ' density of snow = ', rhosn , ' kg.m-3'183 WRITE(numout,*) ' emissivity of snow or ice = ', emic184 WRITE(numout,*) ' salinity of ice = ', sice , ' psu'185 WRITE(numout,*) ' salinity of sea = ', soce , ' psu'186 WRITE(numout,*) ' latent heat of evaporation (water) = ', cevap , ' J.m-3'187 WRITE(numout,*) ' correction factor for solar radiation = ', srgamma188 WRITE(numout,*) ' von Karman constant = ', vkarmn189 WRITE(numout,*) ' Stefan-Boltzmann constant = ', stefan , ' J.s-1.m-2.K-4'167 WRITE(numout,*) ' thermal conductivity of the snow = ', rcdsn , ' J/s/m/K' 168 WRITE(numout,*) ' thermal conductivity of the ice = ', rcdic , ' J/s/m/K' 169 WRITE(numout,*) ' density times specific heat for snow = ', rcpsn , ' J/m^3/K' 170 WRITE(numout,*) ' density times specific heat for ice = ', rcpic , ' J/m^3/K' 171 WRITE(numout,*) ' volumetric latent heat fusion of sea ice = ', xlic , ' J/m' 172 WRITE(numout,*) ' volumetric latent heat fusion of snow = ', xlsn , ' J/m' 173 WRITE(numout,*) ' latent heat of sublimation of snow = ', xsn , ' J/kg' 174 WRITE(numout,*) ' density of sea ice = ', rhoic , ' kg/m^3' 175 WRITE(numout,*) ' density of snow = ', rhosn , ' kg/m^3' 176 WRITE(numout,*) ' emissivity of snow or ice = ', emic 177 WRITE(numout,*) ' salinity of ice = ', sice , ' psu' 178 WRITE(numout,*) ' salinity of sea = ', soce , ' psu' 179 WRITE(numout,*) ' latent heat of evaporation (water) = ', cevap , ' J/m^3' 180 WRITE(numout,*) ' correction factor for solar radiation = ', srgamma 181 WRITE(numout,*) ' von Karman constant = ', vkarmn 182 WRITE(numout,*) ' Stefan-Boltzmann constant = ', stefan , ' J/s/m^2/K^4' 190 183 191 184 WRITE(numout,*) 192 WRITE(numout,*) ' conversion: degre ==> radianrad = ', rad185 WRITE(numout,*) ' conversion: degre ==> radian rad = ', rad 193 186 194 187 WRITE(numout,*) 195 WRITE(numout,*) ' smallest real computer value= ', rsmall188 WRITE(numout,*) ' smallest real computer value rsmall = ', rsmall 196 189 ENDIF 197 190
Note: See TracChangeset
for help on using the changeset viewer.