Ignore:
Timestamp:
2007-10-10T10:14:32+02:00 (12 years ago)
Author:
smasson
Message:

code modifications associated with the new routines, see ticket:4

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/phycst.F90

    r699 r703  
    44   !!     Definition of of both ocean and ice parameters used in the code 
    55   !!===================================================================== 
    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   !!---------------------------------------------------------------------- 
    716   USE par_oce          ! ocean parameters 
    817   USE in_out_manager   ! I/O manager 
     
    1120   PRIVATE 
    1221 
    13    !! * Routine accessibility 
    14    PUBLIC phy_cst          ! routine called by inipar.F90 
     22   PUBLIC   phy_cst     ! routine called by inipar.F90 
    1523 
    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 
    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       ,   &  !: density of sea ice (kg/m3) 
    64       rhosn   = 330._wp       ,   &  !: density of snow (kg/m3) 
     63      rhoic   = 900._wp       ,   &  !: volumic mass of sea ice (kg/m3) 
     64      rhosn   = 330._wp       ,   &  !: volumic mass 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       !!  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   !!---------------------------------------------------------------------- 
    7778    
    7879CONTAINS 
     
    8384      !! 
    8485      !! ** 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  
    9386      !!---------------------------------------------------------------------- 
    94       !! * Local variables 
    95       CHARACTER (len=64) ::   cform = "(A9, 3(A13, I7) )"  
     87      CHARACTER (len=64) ::   cform = "(A12, 3(A13, I7) )"  
    9688      !!---------------------------------------------------------------------- 
    9789 
     
    10395      ! ---------------- 
    10496      IF(lwp) THEN 
    105          WRITE(numout,*) '       parameter file' 
    106          WRITE(numout,*) 
     97         WRITE(numout,*) '       Domain info' 
    10798         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 
    114104         WRITE(numout,*) '          mpp local domain info (mpp)' 
    115105         WRITE(numout,*) '             jpni    : ', jpni, '   jpreci  : ', jpreci 
    116106         WRITE(numout,*) '             jpnj    : ', jpnj, '   jprecj  : ', jprecj 
    117107         WRITE(numout,*) '             jpnij   : ', jpnij 
    118  
    119          WRITE(numout,*) 
    120108         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 
    123110      ENDIF 
    124111 
     
    126113      ! ---------------- 
    127114      IF(lwp) WRITE(numout,*) 
    128       IF(lwp) WRITE(numout,*) '       constants' 
     115      IF(lwp) WRITE(numout,*) '       Constants' 
    129116 
    130117      IF(lwp) WRITE(numout,*) 
Note: See TracChangeset for help on using the changeset viewer.