prep_ocean_unif.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
!     #########
SUBROUTINE prep_ocean_unif(KLUOUT,HSURF,PFIELD)
!     #################################################################################
!
!!****  *PREP_OCEAN_UNIF* - prepares oceanic field from prescribed values
!!
!!    PURPOSE
!!    -------
!
!!**  METHOD
!!    ------
!!
!!    REFERENCE
!!    ---------
!!      
!!
!!    AUTHOR
!!    ------
!!     C. Lebeaupin Brossier
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2008
!!------------------------------------------------------------------
!

!
USE modd_csts,       ONLY : xtt
USE modd_prep,       ONLY : cinterp_type
USE modd_prep_seaflux,   ONLY : xsst_unif
USE modd_ocean_csts
USE modd_ocean_grid
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
INTEGER,            INTENT(IN)  :: kluout    ! output listing logical unit
 CHARACTER(LEN=7),   INTENT(IN)  :: hsurf     ! type of field
REAL, POINTER, DIMENSION(:,:,:)   :: pfield    ! field to interpolate horizontally
!
!*      0.2    declarations of local variables
REAL :: zsss_unif=37.6
REAL :: zgradt,zgrads
INTEGER :: jlev
REAL(KIND=JPRB) :: zhook_handle
!
!
!-------------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('PREP_OCEAN_UNIF',0,zhook_handle)
SELECT CASE(hsurf)
!
!*      3.1    oceanic temperature
!
  CASE('TEMP_OC')
    zgradt=1.2/1000.
    ALLOCATE(pfield(1,nockmax-nockmin+1,1))
    pfield(:,1,:) = xsst_unif-xtt
    DO jlev=1,nockmax
      pfield(:,jlev+1,:) = xsst_unif-xtt-zgradt*(xzhoc(jlev)+1)
    ENDDO
!
!*      3.2    oceanic salinity
!
  CASE('SALT_OC')
    zgrads=-0.5/1000.
    ALLOCATE(pfield(1,nockmax-nockmin+1,1))
    pfield(:,1,:) = zsss_unif
    DO jlev=1,nockmax
      pfield(:,jlev+1,:) = zsss_unif+zgrads*(xzhoc(jlev)+1)
    ENDDO
!
!*      3.3    oceanic currents
!
  CASE('UCUR_OC')
    ALLOCATE(pfield(1,nockmax-nockmin+1,1))
    pfield = 0. 
!
  CASE('VCUR_OC')
    ALLOCATE(pfield(1,nockmax-nockmin+1,1))
    pfield = 0. 
!
END SELECT
!
!*      4.     Interpolation method
!              --------------------
!
 cinterp_type='UNIF  '
IF (lhook) CALL dr_hook('PREP_OCEAN_UNIF',1,zhook_handle)
!
!
!-------------------------------------------------------------------------------------
END SUBROUTINE prep_ocean_unif