New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 719 for trunk/NEMO/OPA_SRC/phycst.F90 – NEMO

Ignore:
Timestamp:
2007-10-16T16:59:56+02:00 (17 years ago)
Author:
ctlod
Message:

get back to the nemo_v2_3 version for trunk

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/NEMO/OPA_SRC/phycst.F90

    • Property svn:keywords changed from Id to Author Date Id Revision
    r703 r719  
    44   !!     Definition of of both ocean and ice parameters used in the code 
    55   !!===================================================================== 
    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    !!---------------------------------------------------------------------- 
     6   !! * Modules used 
    167   USE par_oce          ! ocean parameters 
    178   USE in_out_manager   ! I/O manager 
    … …  
    2011   PRIVATE 
    2112 
    22    PUBLIC   phy_cst     ! routine called by inipar.F90 
     13   !! * Routine accessibility 
     14   PUBLIC phy_cst          ! routine called by inipar.F90 
    2315 
    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 
     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 
    2727    
    2828   REAL(wp), PUBLIC ::          & !: 
    … …  
    6161      xlic    = 300.33e+6_wp  ,   &  !: volumetric latent heat fusion of ice 
    6262      xsn     =   2.8e+6      ,   &  !: latent heat of sublimation of snow 
    63       rhoic   = 900._wp       ,   &  !: volumic mass of sea ice (kg/m3) 
    64       rhosn   = 330._wp       ,   &  !: volumic mass of snow (kg/m3) 
     63      rhoic   = 900._wp       ,   &  !: density of sea ice (kg/m3) 
     64      rhosn   = 330._wp       ,   &  !: density of snow (kg/m3) 
    6565      emic    =   0.97_wp     ,   &  !: emissivity of snow or ice 
    6666      sice    =   6.0_wp      ,   &  !: salinity of ice (psu) 
    … …  
    7070      vkarmn  =   0.4_wp      ,   &  !: von Karman constant 
    7171      stefan  =   5.67e-8_wp         !: Stefan-Boltzmann constant  
    72  
    73    !!---------------------------------------------------------------------- 
    74    !!  OPA 9.0 , LOCEAN-IPSL (2006)  
    75    !! $Id$ 
    76    !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 
    77    !!---------------------------------------------------------------------- 
     72      !!---------------------------------------------------------------------- 
     73      !!  OPA 9.0 , LOCEAN-IPSL (2005)  
     74      !! $Header$  
     75      !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
     76      !!---------------------------------------------------------------------- 
    7877    
    7978CONTAINS 
    … …  
    8483      !! 
    8584      !! ** Purpose :   Print model parameters and set and print the constants 
     85      !! 
     86      !! ** Method  :   no 
     87      !! 
     88      !! History : 
     89      !!        !  90-10  (C. Levy - G. Madec)  Original code 
     90      !!        !  91-11  (G. Madec) 
     91      !!        !  91-12  (M. Imbard) 
     92      !!   8.5  !  02-08  (G. Madec, C. Ethe)  F90, add ice constants  
    8693      !!---------------------------------------------------------------------- 
    87       CHARACTER (len=64) ::   cform = "(A12, 3(A13, I7) )"  
     94      !! * Local variables 
     95      CHARACTER (len=64) ::   cform = "(A9, 3(A13, I7) )"  
    8896      !!---------------------------------------------------------------------- 
    8997 
    … …  
    95103      ! ---------------- 
    96104      IF(lwp) THEN 
    97          WRITE(numout,*) '       Domain info' 
     105         WRITE(numout,*) '       parameter file' 
     106         WRITE(numout,*) 
    98107         WRITE(numout,*) '          dimension of model' 
    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 
     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,*) 
    104114         WRITE(numout,*) '          mpp local domain info (mpp)' 
    105115         WRITE(numout,*) '             jpni    : ', jpni, '   jpreci  : ', jpreci 
    106116         WRITE(numout,*) '             jpnj    : ', jpnj, '   jprecj  : ', jprecj 
    107117         WRITE(numout,*) '             jpnij   : ', jpnij 
     118 
     119         WRITE(numout,*) 
    108120         WRITE(numout,*) '          lateral domain boundary condition type : jperio  = ', jperio 
    109          WRITE(numout,*) '          domain island (use in rigid-lid case)  : jpisl   = ', jpisl, '   jpnisl  = ', jpnisl 
     121         WRITE(numout,*) '          domain island (use in rigid-lid case)  : jpisl   = ', jpisl  
     122         WRITE(numout,*) '                                                   jpnisl  = ', jpnisl 
    110123      ENDIF 
    111124 
    … …  
    113126      ! ---------------- 
    114127      IF(lwp) WRITE(numout,*) 
    115       IF(lwp) WRITE(numout,*) '       Constants' 
     128      IF(lwp) WRITE(numout,*) '       constants' 
    116129 
    117130      IF(lwp) WRITE(numout,*) 
Note: See TracChangeset for help on using the changeset viewer.