prep_ocean_netcdf.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_netcdf(HPROGRAM,HSURF,HFILE,HFILETYPE,&
                              kluout,hncvarname,pfield)  
!     #################################################################################
!
!!****  *PREP_OCEAN_NETCDF* - prepares oceanic fields from Mercator analysis
!!
!!    PURPOSE
!!    -------
!
!!**  METHOD
!!    ------
!!
!!    REFERENCE
!!    ---------
!!      
!!
!!    AUTHOR
!!    ------
!!     C. Lebeaupin Brossier
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2008
!!------------------------------------------------------------------
!
USE modd_prep,       ONLY : cinterp_type, cingrid_type
USE modd_grid_latlonregul, ONLY : nilength, nindepth, xilonarray
!
USE mode_read_netcdf_mercator
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
 CHARACTER(LEN=6),   INTENT(IN)  :: hprogram  ! program calling surf. schemes
 CHARACTER(LEN=7),   INTENT(IN)  :: hsurf     ! type of field
 CHARACTER(LEN=28),  INTENT(IN)  :: hfile     ! file name
 CHARACTER(LEN=6),   INTENT(IN)  :: hfiletype ! file type
INTEGER,            INTENT(IN)  :: kluout    ! output listing logical unit
 CHARACTER(LEN=28),  INTENT(IN), OPTIONAL   :: hncvarname!var to read 
REAL, POINTER, DIMENSION(:,:,:)   :: pfield    ! field to interpolate horizontally
!
!*      0.2    declarations of local variables
REAL,DIMENSION(:), ALLOCATABLE :: zlati
REAL,DIMENSION(:), ALLOCATABLE :: zlong
REAL,DIMENSION(:), ALLOCATABLE :: zdepth
REAL,DIMENSION(:,:,:), ALLOCATABLE :: zfield
REAL(KIND=JPRB) :: zhook_handle
!
!
!-------------------------------------------------------------------------------------
!*      1.     Grid type
!              ---------
IF (lhook) CALL dr_hook('PREP_OCEAN_NETCDF',0,zhook_handle)
 cingrid_type='LATLON'
!
!
IF (.NOT. ALLOCATED(xilonarray)) CALL prep_netcdf_grid(hfile,hncvarname)
!
ALLOCATE(zlati(nilength) )
ALLOCATE(zlong(nilength) )
ALLOCATE(zdepth(nindepth))
!
ALLOCATE(zfield(nilength,nindepth,1))
!
!*      2.     Reading of field
!              ----------------
 CALL read_latlondepval_cdf(hfile,hncvarname,zlong,zlati,zdepth,zfield(:,:,1))
ALLOCATE(pfield(1:SIZE(zfield,1),1:SIZE(zfield,2),1:SIZE(zfield,3)))
pfield=zfield
!
!*      3.     Interpolation method
!              --------------------
!
 cinterp_type='HORIBL'
!
!*      4.     Deallocations
!              -------------
!
IF (ALLOCATED(zlong       ))  DEALLOCATE(zlong  )
IF (ALLOCATED(zlati       ))  DEALLOCATE(zlati  )
IF (ALLOCATED(zdepth      ))  DEALLOCATE(zdepth )
IF (ALLOCATED(zfield      ))  DEALLOCATE(zfield )
IF (lhook) CALL dr_hook('PREP_OCEAN_NETCDF',1,zhook_handle)
!
!-------------------------------------------------------------------------------------
END SUBROUTINE prep_ocean_netcdf