oi_cavegi.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 oi_cavegi(PVGAT1,PVGAT2,PVGAT3,PVGBT1,PVGBT2,PVGBT3,PVGCT1,PVGCT2,       &
                     pvgah1,pvgah2,pvgah3,pvgbh1,pvgbh2,pvgbh3,pvgch1,pvgch2,       &
                     psigt2mp,psighp2,osgobs)   
!****------------------------------------------------------------------------
USE modd_assim, ONLY : xsighp1, xsigt2mr, xsigh2mr, xsigt2mo, xsigh2mo, xreps3
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
INTEGER :: j
! 
REAL ,INTENT(OUT) :: pvgat1(24),pvgat2(24),pvgat3(24)
REAL ,INTENT(OUT) :: pvgbt1(24),pvgbt2(24),pvgbt3(24)
REAL ,INTENT(OUT) :: pvgct1(24),pvgct2(24)
REAL ,INTENT(OUT) :: pvgah1(24),pvgah2(24),pvgah3(24)
REAL ,INTENT(OUT) :: pvgbh1(24),pvgbh2(24),pvgbh3(24)
REAL ,INTENT(OUT) :: pvgch1(24),pvgch2(24)
REAL ,INTENT(OUT) :: psigt2mp(24),psighp2(24)
!
LOGICAL :: osgobs
REAL(KIND=JPRB) :: zhook_handle
!
!**---------------------------------------------------------------------
!**  1. Initialisation des polynomes bruts et des champs de reference.
!**     -------------------------------------------------------------
!
! lecture des coefficients polynomiaux
!
IF (lhook) CALL dr_hook('OI_CAVEGI',0,zhook_handle)
!
DO j=1,24

  READ(61,'(6X,8F10.4)') pvgat1(j),pvgat2(j),pvgat3(j),&
                         pvgbt1(j),pvgbt2(j),pvgbt3(j),&
                         pvgct1(j),pvgct2(j)    

ENDDO
!
DO j=1,24

  READ(61,'(6X,8F10.4)') pvgah1(j),pvgah2(j),pvgah3(j),&
                         pvgbh1(j),pvgbh2(j),pvgbh3(j),&
                         pvgch1(j),pvgch2(j)    

ENDDO
!
DO j = 1,24 

  psigt2mp(j) = max(.3 , 2.7*(1.0 - ((REAL(j)-15.)/9.)**2))   
  psighp2(j) = (REAL(j)*2.0/3. - 15.)*1.e-2

ENDDO
!
!**---------------------------------------------------------------------
!**  3. Initialisation des variables internes.
!**     -------------------------------------
!
osgobs = xsigt2mo > 0.0 .AND. xsigh2mo > 0.0 .AND. &
         (abs(xsigh2mo-xsigh2mr) > xreps3 .OR. abs(xsigt2mo-xsigt2mr) > xreps3)  

IF (lhook) CALL dr_hook('OI_CAVEGI',1,zhook_handle)
!
!**---------------------------------------------------------------------
END SUBROUTINE oi_cavegi