read_sson.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 read_sso_n (KSIZE_FULL, PSEA, USS, HPROGRAM)
!     ################################
!
!!****  *READ_SSO_n* - routine to read a file for
!!                         physiographic data file of model _n 
!!
!!    PURPOSE
!!    -------
!!       The purpose of this routine is to initialise the 
!!       physiographic data file.
!!
!!
!!**  METHOD
!!    ------
!!      The data are read in the initial surface file :
!!        - 2D physiographic data fields
!!          
!!      It does not read the grid definition. This should have been
!!      read already.
!!
!!    EXTERNAL
!!    --------
!!      
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2003
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
USE modd_sso_n, ONLY : sso_t
!
USE modi_read_surf
!
USE modd_surf_par, ONLY : xundef
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
INTEGER, INTENT(IN) :: KSIZE_FULL
REAL, DIMENSION(:), INTENT(IN) :: PSEA
!
TYPE(sso_t), INTENT(INOUT) :: USS
!
 CHARACTER(LEN=6),  INTENT(IN)  :: HPROGRAM ! calling program
!
!*       0.2   Declarations of local variables
!              -------------------------------
!

!
INTEGER           :: IRESP          ! Error code after redding
! 
 CHARACTER(LEN=12) :: YRECFM         ! Name of the article to be read
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------
!
!*       2.     Orography :
!               ---------
!
!
IF (lhook) CALL dr_hook('READ_SSO_N',0,zhook_handle)
IF(.NOT.ASSOCIATED(uss%XAVG_ZS)) ALLOCATE(uss%XAVG_ZS(ksize_full))
yrecfm='AVG_ZS'
 CALL read_surf(hprogram,yrecfm,uss%XAVG_ZS(:),iresp)
!
IF(.NOT.ASSOCIATED(uss%XSIL_ZS)) ALLOCATE(uss%XSIL_ZS(ksize_full))
yrecfm='SIL_ZS'
 CALL read_surf(hprogram,yrecfm,uss%XSIL_ZS(:),iresp)
!
!
!*       3.     Subgrid Orography :
!               -----------------
!
!
IF(.NOT.ASSOCIATED(uss%XSSO_STDEV)) ALLOCATE(uss%XSSO_STDEV(ksize_full))
yrecfm='SSO_STDEV'
 CALL read_surf(hprogram,yrecfm,uss%XSSO_STDEV(:),iresp)
WHERE (psea(:) == 1.) uss%XSSO_STDEV(:) = xundef
!
IF(.NOT.ASSOCIATED(uss%XMIN_ZS)) ALLOCATE(uss%XMIN_ZS(ksize_full))
yrecfm='MIN_ZS'
 CALL read_surf(hprogram,yrecfm,uss%XMIN_ZS(:),iresp)
!
IF(.NOT.ASSOCIATED(uss%XMAX_ZS)) ALLOCATE(uss%XMAX_ZS(ksize_full))
yrecfm='MAX_ZS'
 CALL read_surf(hprogram,yrecfm,uss%XMAX_ZS(:),iresp)
!
IF(.NOT.ASSOCIATED(uss%XSSO_ANIS)) ALLOCATE(uss%XSSO_ANIS(ksize_full))
yrecfm='SSO_ANIS'
 CALL read_surf(hprogram,yrecfm,uss%XSSO_ANIS(:),iresp)
WHERE (psea(:) == 1.) uss%XSSO_ANIS(:) = xundef
!
IF(.NOT.ASSOCIATED(uss%XSSO_DIR)) ALLOCATE(uss%XSSO_DIR(ksize_full))
yrecfm='SSO_DIR'
 CALL read_surf(hprogram,yrecfm,uss%XSSO_DIR(:),iresp)
WHERE (psea(:) == 1.) uss%XSSO_DIR(:) = xundef
!
IF(.NOT.ASSOCIATED(uss%XSSO_SLOPE)) ALLOCATE(uss%XSSO_SLOPE(ksize_full))
yrecfm='SSO_SLOPE'
 CALL read_surf(hprogram,yrecfm,uss%XSSO_SLOPE(:),iresp)
WHERE (psea(:) == 1.) uss%XSSO_SLOPE(:) = xundef
!
!-------------------------------------------------------------------------------
!
!*       3.     Subgrid Orography roughness:
!               ---------------------------
!
!
IF(.NOT.ASSOCIATED(uss%XHO2IP)) ALLOCATE(uss%XHO2IP(ksize_full))
yrecfm='HO2IP'
 CALL read_surf(hprogram,yrecfm,uss%XHO2IP(:),iresp)
WHERE (psea(:) == 1.) uss%XHO2IP(:) = xundef
!
IF(.NOT.ASSOCIATED(uss%XHO2JP)) ALLOCATE(uss%XHO2JP(ksize_full))
yrecfm='HO2JP'
 CALL read_surf(hprogram,yrecfm,uss%XHO2JP(:),iresp)
WHERE (psea(:) == 1.) uss%XHO2JP(:) = xundef
!
IF(.NOT.ASSOCIATED(uss%XHO2IM)) ALLOCATE(uss%XHO2IM(ksize_full))
yrecfm='HO2IM'
 CALL read_surf(hprogram,yrecfm,uss%XHO2IM(:),iresp)
WHERE (psea(:) == 1.) uss%XHO2IM(:) = xundef
!
IF(.NOT.ASSOCIATED(uss%XHO2JM)) ALLOCATE(uss%XHO2JM(ksize_full))
yrecfm='HO2JM'
 CALL read_surf(hprogram,yrecfm,uss%XHO2JM(:),iresp)
WHERE (psea(:) == 1.) uss%XHO2JM(:) = xundef
!
IF(.NOT.ASSOCIATED(uss%XAOSIP)) ALLOCATE(uss%XAOSIP(ksize_full))
yrecfm='AOSIP'
 CALL read_surf(hprogram,yrecfm,uss%XAOSIP(:),iresp)
WHERE (psea(:) == 1.) uss%XAOSIP(:) = xundef
!
IF(.NOT.ASSOCIATED(uss%XAOSJP)) ALLOCATE(uss%XAOSJP(ksize_full))
yrecfm='AOSJP'
 CALL read_surf(hprogram,yrecfm,uss%XAOSJP(:),iresp)
WHERE (psea(:) == 1.) uss%XAOSJP(:) = xundef
!
IF(.NOT.ASSOCIATED(uss%XAOSIM)) ALLOCATE(uss%XAOSIM(ksize_full))
yrecfm='AOSIM'
 CALL read_surf(hprogram,yrecfm,uss%XAOSIM(:),iresp)
WHERE (psea(:) == 1.) uss%XAOSIM(:) = xundef
!
IF(.NOT.ASSOCIATED(uss%XAOSJM)) ALLOCATE(uss%XAOSJM(ksize_full))
yrecfm='AOSJM'
 CALL read_surf(hprogram,yrecfm,uss%XAOSJM(:),iresp)
WHERE (psea(:) == 1.) uss%XAOSJM(:) = xundef
IF (lhook) CALL dr_hook('READ_SSO_N',1,zhook_handle)
!
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE read_sso_n