SURFEX v7.3
General documentation of Surfex
|
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