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_surfex_mpi, ONLY : nrank, npio
USE modd_pgdwork,  ONLY : nsso, nsso_all, xsso_all
USE modd_surf_par, ONLY : nundef, xundef
USE modd_surfex_mpi, ONLY : nindex
!
!
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)

INTEGER :: IDIMF
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!----------------------------------------------------------------------------
!
!
!*    1.     Adapt subgrid mesh to input file resolution
!            -------------------------------------------
!
IF (lhook) CALL dr_hook('INI_SSOWORK',0,zhook_handle)
!
idimf = SIZE(nindex)
!
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(xsso_all)) DEALLOCATE(xsso_all)
!
ALLOCATE(xsso_all(idimf,nsso,nsso))
xsso_all(:,:,:) = -xundef
!
IF (ALLOCATED(nsso_all)) DEALLOCATE(nsso_all)
!
ALLOCATE(nsso_all(idimf,nsso,nsso))
nsso_all(:,:,:) = 0
!
IF (lhook) CALL dr_hook('INI_SSOWORK',1,zhook_handle)
!
!----------------------------------------------------------------------------
!
END SUBROUTINE ini_ssowork