SURFEX v7.3: /home/dasprezs/EXPORT_v7_3/src/OFFLIN/oi

Go to the documentation of this file.
00001 SUBROUTINE OI_CAVEGI(VGAT1,VGAT2,VGAT3,VGBT1,VGBT2,VGBT3,VGCT1,VGCT2,       &
00002                       VGAH1,VGAH2,VGAH3,VGBH1,VGBH2,VGBH3,VGCH1,VGCH2,       &
00003                       SIGT2MP,SIGHP2,LSGOBS)   
00004 !****------------------------------------------------------------------------
00005 USE MODD_ASSIM, ONLY : SIGHP1, SIGT2MR, SIGH2MR, SIGT2MO, SIGH2MO, REPS3
00006 !
00007 !
00008 USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
00009 USE PARKIND1  ,ONLY : JPRB
00010 !
00011 IMPLICIT NONE
00012 !
00013 INTEGER :: J
00014 ! 
00015 REAL ,INTENT(OUT) :: VGAT1(24),VGAT2(24),VGAT3(24)
00016 REAL ,INTENT(OUT) :: VGBT1(24),VGBT2(24),VGBT3(24)
00017 REAL ,INTENT(OUT) :: VGCT1(24),VGCT2(24)
00018 REAL ,INTENT(OUT) :: VGAH1(24),VGAH2(24),VGAH3(24)
00019 REAL ,INTENT(OUT) :: VGBH1(24),VGBH2(24),VGBH3(24)
00020 REAL ,INTENT(OUT) :: VGCH1(24),VGCH2(24)
00021 REAL ,INTENT(OUT) :: SIGT2MP(24),SIGHP2(24)
00022 !
00023 LOGICAL :: LSGOBS
00024 REAL(KIND=JPRB) :: ZHOOK_HANDLE
00025 !
00026 !**---------------------------------------------------------------------
00027 !**  1. Initialisation des polynomes bruts et des champs de reference.
00028 !**     -------------------------------------------------------------
00029 !
00030 ! lecture des coefficients polynomiaux
00031 !
00032 IF (LHOOK) CALL DR_HOOK('OI_CAVEGI',0,ZHOOK_HANDLE)
00033 DO J=1,24
00034   READ(61,'(6X,8F10.4)') VGAT1(J),VGAT2(J),VGAT3(J),&
00035       VGBT1(J),VGBT2(J),VGBT3(J),&
00036       VGCT1(J),VGCT2(J)    
00037 ENDDO
00038 !
00039 DO J=1,24
00040   READ(61,'(6X,8F10.4)') VGAH1(J),VGAH2(J),VGAH3(J),&
00041       VGBH1(J),VGBH2(J),VGBH3(J),&
00042       VGCH1(J),VGCH2(J)    
00043 ENDDO
00044 !
00045 DO J = 1,24   
00046   SIGT2MP(J)= MAX(.3 , 2.7*(1.0 - ((REAL(J)-15.)/9.)**2))   
00047   SIGHP2(J)= (REAL(J)*2.0/3. - 15.)*1.E-2
00048 ENDDO
00049 !
00050 !**---------------------------------------------------------------------
00051 !**  3. Initialisation des variables internes.
00052 !**     -------------------------------------
00053 !
00054 LSGOBS = SIGT2MO > 0.0 .AND. SIGH2MO > 0.0 .AND. &
00055     (ABS(SIGH2MO-SIGH2MR) > REPS3 .OR. ABS(SIGT2MO-SIGT2MR) > REPS3)    
00056 IF (LHOOK) CALL DR_HOOK('OI_CAVEGI',1,ZHOOK_HANDLE)
00057 !
00058 !**---------------------------------------------------------------------
00059 END SUBROUTINE OI_CAVEGI
00060 
00061