Changeset 703 for trunk/NEMO/OPA_SRC/phycst.F90
- Timestamp:
- 2007-10-10T10:14:32+02:00 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/phycst.F90
r699 r703 4 4 !! Definition of of both ocean and ice parameters used in the code 5 5 !!===================================================================== 6 !! * Modules used 6 !! History : ! 90-10 (C. Levy - G. Madec) Original code 7 !! ! 91-11 (G. Madec) 8 !! ! 91-12 (M. Imbard) 9 !! 8.5 ! 02-08 (G. Madec, C. Ethe) F90, add ice constants 10 !! 9.0 ! 06-08 (G. Madec) style 11 !!---------------------------------------------------------------------- 12 13 !!---------------------------------------------------------------------- 14 !! phy_cst : define and print physical constant and domain parameters 15 !!---------------------------------------------------------------------- 7 16 USE par_oce ! ocean parameters 8 17 USE in_out_manager ! I/O manager … … 11 20 PRIVATE 12 21 13 !! * Routine accessibility 14 PUBLIC phy_cst ! routine called by inipar.F90 22 PUBLIC phy_cst ! routine called by inipar.F90 15 23 16 !! * Shared module variables 17 INTEGER, PUBLIC, DIMENSION(12) :: & !: 18 nbiss = (/ 31, 29, 31, 30, 31, 30, & !: number of days per month 19 & 31, 31, 30, 31, 30, 31 /) , & ! (leap-year) 20 nobis = (/ 31, 28, 31, 30, 31, 30, & !: number of days per month 21 & 31, 31, 30, 31, 30, 31 /) ! (365 days a year) 22 23 REAL(wp), PUBLIC :: & !: 24 rpi = 3.141592653589793_wp , & !: pi 25 rad = 3.141592653589793_wp / 180._wp , & !: conversion from degre into radian 26 rsmall = 0.5 * EPSILON( 1. ) !: smallest real computer value 24 REAL(wp), PUBLIC :: rpi = 3.141592653589793_wp !: pi 25 REAL(wp), PUBLIC :: rad = 3.141592653589793_wp / 180._wp !: conversion from degre into radian 26 REAL(wp), PUBLIC :: rsmall = 0.5 * EPSILON( 1. ) !: smallest real computer value 27 27 28 28 REAL(wp), PUBLIC :: & !: … … 61 61 xlic = 300.33e+6_wp , & !: volumetric latent heat fusion of ice 62 62 xsn = 2.8e+6 , & !: latent heat of sublimation of snow 63 rhoic = 900._wp , & !: densityof sea ice (kg/m3)64 rhosn = 330._wp , & !: densityof snow (kg/m3)63 rhoic = 900._wp , & !: volumic mass of sea ice (kg/m3) 64 rhosn = 330._wp , & !: volumic mass of snow (kg/m3) 65 65 emic = 0.97_wp , & !: emissivity of snow or ice 66 66 sice = 6.0_wp , & !: salinity of ice (psu) … … 70 70 vkarmn = 0.4_wp , & !: von Karman constant 71 71 stefan = 5.67e-8_wp !: Stefan-Boltzmann constant 72 !!---------------------------------------------------------------------- 73 !! OPA 9.0 , LOCEAN-IPSL (2005) 74 !! $Id$ 75 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 76 !!---------------------------------------------------------------------- 72 73 !!---------------------------------------------------------------------- 74 !! OPA 9.0 , LOCEAN-IPSL (2006) 75 !! $Id$ 76 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 77 !!---------------------------------------------------------------------- 77 78 78 79 CONTAINS … … 83 84 !! 84 85 !! ** Purpose : Print model parameters and set and print the constants 85 !!86 !! ** Method : no87 !!88 !! History :89 !! ! 90-10 (C. Levy - G. Madec) Original code90 !! ! 91-11 (G. Madec)91 !! ! 91-12 (M. Imbard)92 !! 8.5 ! 02-08 (G. Madec, C. Ethe) F90, add ice constants93 86 !!---------------------------------------------------------------------- 94 !! * Local variables 95 CHARACTER (len=64) :: cform = "(A9, 3(A13, I7) )" 87 CHARACTER (len=64) :: cform = "(A12, 3(A13, I7) )" 96 88 !!---------------------------------------------------------------------- 97 89 … … 103 95 ! ---------------- 104 96 IF(lwp) THEN 105 WRITE(numout,*) ' parameter file' 106 WRITE(numout,*) 97 WRITE(numout,*) ' Domain info' 107 98 WRITE(numout,*) ' dimension of model' 108 WRITE(numout,*) ' Local domain Global domain Data domain ' 109 WRITE(numout,cform) ' ',' jpi : ', jpi, ' jpiglo : ', jpiglo, ' jpidta : ', jpidta 110 WRITE(numout,cform) ' ',' jpj : ', jpj, ' jpjglo : ', jpjglo, ' jpjdta : ', jpjdta 111 WRITE(numout,cform) ' ',' jpk : ', jpk, ' jpk : ', jpk , ' jpkdta : ', jpkdta 112 WRITE(numout,*) ' ',' jpij : ', jpij 113 WRITE(numout,*) 99 WRITE(numout,*) ' Local domain Global domain Data domain ' 100 WRITE(numout,cform) ' ',' jpi : ', jpi, ' jpiglo : ', jpiglo, ' jpidta : ', jpidta 101 WRITE(numout,cform) ' ',' jpj : ', jpj, ' jpjglo : ', jpjglo, ' jpjdta : ', jpjdta 102 WRITE(numout,cform) ' ',' jpk : ', jpk, ' jpk : ', jpk , ' jpkdta : ', jpkdta 103 WRITE(numout,*) ' ',' jpij : ', jpij 114 104 WRITE(numout,*) ' mpp local domain info (mpp)' 115 105 WRITE(numout,*) ' jpni : ', jpni, ' jpreci : ', jpreci 116 106 WRITE(numout,*) ' jpnj : ', jpnj, ' jprecj : ', jprecj 117 107 WRITE(numout,*) ' jpnij : ', jpnij 118 119 WRITE(numout,*)120 108 WRITE(numout,*) ' lateral domain boundary condition type : jperio = ', jperio 121 WRITE(numout,*) ' domain island (use in rigid-lid case) : jpisl = ', jpisl 122 WRITE(numout,*) ' jpnisl = ', jpnisl 109 WRITE(numout,*) ' domain island (use in rigid-lid case) : jpisl = ', jpisl, ' jpnisl = ', jpnisl 123 110 ENDIF 124 111 … … 126 113 ! ---------------- 127 114 IF(lwp) WRITE(numout,*) 128 IF(lwp) WRITE(numout,*) ' constants'115 IF(lwp) WRITE(numout,*) ' Constants' 129 116 130 117 IF(lwp) WRITE(numout,*)
Note: See TracChangeset
for help on using the changeset viewer.