SURFEX v7.3
General documentation of Surfex
 All Classes Files Functions Variables Typedefs
/home/dasprezs/EXPORT_v7_3/src/SURFEX/get_coordn.F90
Go to the documentation of this file.
00001 !     #########
00002       SUBROUTINE GET_COORD_n(HPROGRAM,KI,PLON,PLAT)
00003 !     #############################################
00004 !
00005 !!****  *GET_COORD_n* - routine to get some surface fields
00006 !!
00007 !!    PURPOSE
00008 !!    -------
00009 !!
00010 !!**  METHOD
00011 !!    ------
00012 !!
00013 !!    EXTERNAL
00014 !!    --------
00015 !!
00016 !!
00017 !!    IMPLICIT ARGUMENTS
00018 !!    ------------------
00019 !!
00020 !!    REFERENCE
00021 !!    ---------
00022 !!
00023 !!
00024 !!    AUTHOR
00025 !!    ------
00026 !!      P. Le Moigne   *Meteo France*   
00027 !!
00028 !!    MODIFICATIONS
00029 !!    -------------
00030 !!      Original    08/2008
00031 !-------------------------------------------------------------------------------
00032 !
00033 !*       0.    DECLARATIONS
00034 !              ------------
00035 !
00036 USE MODI_GET_LUOUT
00037 !
00038 USE MODD_SURF_ATM_GRID_n,     ONLY : XLON, XLAT
00039 !
00040 !
00041 USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
00042 USE PARKIND1  ,ONLY : JPRB
00043 !
00044 USE MODI_ABOR1_SFX
00045 !
00046 IMPLICIT NONE
00047 !
00048 !*       0.1   Declarations of arguments
00049 !              -------------------------
00050 !
00051  CHARACTER(LEN=6),    INTENT(IN)  :: HPROGRAM
00052 INTEGER,             INTENT(IN)  :: KI         ! horizontal dim. of cover
00053 REAL, DIMENSION(KI), INTENT(OUT) :: PLON, PLAT
00054 !
00055 !
00056 !*       0.2   Declarations of local variables
00057 !              -------------------------------
00058 !
00059 INTEGER :: ILUOUT
00060 REAL(KIND=JPRB) :: ZHOOK_HANDLE
00061 !
00062 !-------------------------------------------------------------------------------
00063 IF (LHOOK) CALL DR_HOOK('GET_COORD_N',0,ZHOOK_HANDLE)
00064  CALL GET_LUOUT(HPROGRAM,ILUOUT)
00065 !-------------------------------------------------------------------------------
00066 !
00067 IF ( SIZE(PLON) /= SIZE(XLON) .OR. SIZE(PLAT) /= SIZE(XLAT) ) THEN
00068   WRITE(ILUOUT,*) 'try to get LON/LAT field from atmospheric model, but size is not correct'
00069   WRITE(ILUOUT,*) 'size of field expected by the atmospheric model (PLON) :', SIZE(PLON)
00070   WRITE(ILUOUT,*) 'size of field expected by the atmospheric model (PLAT) :', SIZE(PLAT)
00071   WRITE(ILUOUT,*) 'size of field in the surface                    (XLAT) :', SIZE(XLAT)
00072   WRITE(ILUOUT,*) 'size of field in the surface                    (XLAT) :', SIZE(XLAT)
00073   CALL ABOR1_SFX('GET_COORDN: LON/LAT SIZE NOT CORRECT')
00074 ELSE
00075   PLON = XLON
00076   PLAT = XLAT
00077 END IF
00078 IF (LHOOK) CALL DR_HOOK('GET_COORD_N',1,ZHOOK_HANDLE)
00079 !
00080 !==============================================================================
00081 !
00082 END SUBROUTINE GET_COORD_n