ini_ssowork.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 ini_ssowork(PMESHLENGTH,PDLAT,PDLON)
!     ###############################################
!
!!**** *INI_SSOWORK* initializes and allocate work arrays for SSO reading
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    10/12/97
!!
!
USE modd_pgdwork,  ONLY : xssqo, lssqo, nsso
USE modd_surf_par, ONLY : nundef, xundef
USE modd_pgd_grid, ONLY : nl
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!----------------------------------------------------------------------------
!
!*    0.1    Declaration of arguments
!            ------------------------
!
REAL, OPTIONAL, INTENT(IN) :: pmeshlength ! average mesh length in degrees
REAL, OPTIONAL, INTENT(IN) :: pdlat       ! input file mesh size (in latitude,  degrees)
REAL, OPTIONAL, INTENT(IN) :: pdlon       ! input file mesh size (in longitude, degrees)
REAL(KIND=JPRB) :: zhook_handle
!
!----------------------------------------------------------------------------
!
!
!*    1.     Adapt subgrid mesh to input file resolution
!            -------------------------------------------
!
IF (lhook) CALL dr_hook('INI_SSOWORK',0,zhook_handle)
IF (present(pmeshlength) .AND. present(pdlat) .AND. present(pdlon)) THEN
  IF (pdlat/= xundef .AND. pdlon /= xundef) THEN
    nsso = nint( 2. * pmeshlength / (pdlat + pdlon) )
    nsso = max(nsso,3)
    nsso = min(nsso,10)
  ELSE
    nsso = 10
  END IF
ELSE
  nsso = 10
END IF
!
!----------------------------------------------------------------------------
!
!*    2.     Allocate subgrid arrays
!            -----------------------
!
IF (ALLOCATED(xssqo)) DEALLOCATE(xssqo)
IF (ALLOCATED(lssqo)) DEALLOCATE(lssqo)
!
ALLOCATE(xssqo(nsso,nsso,nl))
ALLOCATE(lssqo(nsso,nsso,nl))
xssqo(:,:,:) = -99999.
lssqo(:,:,:) = .false.
IF (lhook) CALL dr_hook('INI_SSOWORK',1,zhook_handle)
!
!----------------------------------------------------------------------------
!
END SUBROUTINE ini_ssowork