source: NEMO/trunk/src/OCE/SBC/ocealb.F90 @ 10069

Last change on this file since 10069 was 10069, checked in by nicolasmartin, 2 years ago

Fix mistakes of previous commit on SVN keywords property

  • Property svn:keywords set to Id
File size: 2.2 KB
Line 
1MODULE ocealb
2   !!======================================================================
3   !!                       ***  MODULE  ocealb  ***
4   !! Ocean forcing:  bulk thermohaline forcing of the ocean
5   !!=====================================================================
6   !! History :
7   !!   NEMO     4.0  ! 2017-07  (C. Rousset) Split ocean and ice albedos
8   !!----------------------------------------------------------------------
9   !!   oce_alb    : albedo for ocean (clear and overcast skies)
10   !!----------------------------------------------------------------------
11   USE phycst         ! physical constants
12   USE in_out_manager ! I/O manager
13   USE lib_mpp        ! MPP library
14   USE lib_fortran    ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 
15
16   IMPLICIT NONE
17   PRIVATE
18
19   PUBLIC   oce_alb   ! routine called by sbccpl
20 
21   !!----------------------------------------------------------------------
22   !! NEMO/OCE 4.0 , NEMO Consortium (2018)
23   !! $Id$
24   !! Software governed by the CeCILL license (see ./LICENSE)
25   !!----------------------------------------------------------------------
26CONTAINS
27
28   SUBROUTINE oce_alb( palb_os , palb_cs )
29      !!----------------------------------------------------------------------
30      !!               ***  ROUTINE oce_alb  ***
31      !!
32      !! ** Purpose :   Computation of the albedo of the ocean
33      !!----------------------------------------------------------------------
34      REAL(wp), DIMENSION(:,:), INTENT(out) ::   palb_os   !  albedo of ocean under overcast sky
35      REAL(wp), DIMENSION(:,:), INTENT(out) ::   palb_cs   !  albedo of ocean under clear sky
36      !!
37      REAL(wp) ::   zcoef 
38      REAL(wp) ::   rmue = 0.40    !  cosine of local solar altitude
39      !!----------------------------------------------------------------------
40      !
41      zcoef = 0.05 / ( 1.1 * rmue**1.4 + 0.15 )   ! Parameterization of Briegled and Ramanathan, 1982
42      palb_cs(:,:) = zcoef 
43      palb_os(:,:) = 0.06                         ! Parameterization of Kondratyev, 1969 and Payne, 1972
44      !
45   END SUBROUTINE oce_alb
46
47   !!======================================================================
48END MODULE ocealb
Note: See TracBrowser for help on using the repository browser.