write_budget_coupl_rout.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 write_budget_coupl_rout
!     ##########################
!
!!
!!    PURPOSE
!!    -------
!        
!     
!!**  METHOD
!!    ------
!
!!    EXTERNAL
!!    --------
!!
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------ 
!!      
!!    REFERENCE
!!    ---------
!!     
!!    AUTHOR
!!    ------
!!
!!      L. Bouilloud & B. Vincendon     * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!
!!      Original  03/2008 
!!                03/2014 (B. Vincendon) add new control variables
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modd_topd_par, ONLY : nunit
USE modd_topodyn,       ONLY : ccat, nncat, nnb_topd_step
USE modd_budget_coupl_rout
!
USE modi_open_file
USE modi_close_file
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!      none
!*      0.2    declarations of local variables
INTEGER           :: jcat,jstp ! loop control
INTEGER           :: inb_var   ! number of variable to write
 CHARACTER(LEN=28) :: yfile     ! file name
 CHARACTER(LEN=40) :: yform     ! Writing format
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('WRITE_BUDGET_COUPL_ROUT',0,zhook_handle)
!
!*       1.     WRITING BUDGET FILES
!               ------------
inb_var=12
!
DO jcat=1,nncat
  !
  yform='(i6,12f15.1)'
  yfile=trim('bilan_bv_')//trim(ccat(jcat))//trim('.txt')
  !
  CALL open_file('ASCII ',nunit,hfile=yfile,hform='FORMATTED',haction='WRITE')
  !
  WRITE(nunit,*) '     T','         ',yb_var(1),'         ',yb_var(2),'         ',&
                                      yb_var(3),'         ',yb_var(4),'         ',&
                                      yb_var(5),'         ',yb_var(6),'         ',&
                                      yb_var(7),'         ',yb_var(8),'         ',&
                                      yb_var(9),'         ',yb_var(10),'         ',&
                                      yb_var(11),'         ',yb_var(12)
  !
  DO jstp=1,nnb_topd_step
    WRITE(nunit,yform) jstp,xb_var_bv(jstp,jcat,1:inb_var)
  ENDDO
  !
  CALL close_file('ASCII ',nunit)
  ! 
  yfile=trim('bilan_nobv_')//trim(ccat(jcat))//trim('.txt')
  !
  CALL open_file('ASCII ',nunit,hfile=yfile,hform='FORMATTED',haction='WRITE')
  ! 
  WRITE(nunit,*) '     T','         ',yb_var(1),'         ',yb_var(2),'         ',&
                                      yb_var(3),'         ',yb_var(4),'         ',&
                                      yb_var(5),'         ',yb_var(6),'         ',&
                                      yb_var(7),'         ',yb_var(8),'         ',&
                                      yb_var(9),'         ',yb_var(10),'         ',&
                                      yb_var(11),'         ',yb_var(12)

  !
  DO jstp=1,nnb_topd_step
    WRITE(nunit,yform) jstp,xb_var_nobv(jstp,jcat,1:inb_var)
  ENDDO
  !
  CALL close_file('ASCII ',nunit)
  !
  yform='(i6,5f15.1)'
  yfile=trim('bilan_q.txt')
  !
  CALL open_file('ASCII ',nunit,hfile=yfile,hform='FORMATTED',haction='WRITE')
  !
  WRITE(nunit,*) '     T','         ',yb_varq(1),'         ',yb_varq(2),'         ',&
                                      yb_varq(3),'         ',yb_varq(4),'         ',&
                                      yb_varq(5),'         '
  ! 
  DO jstp=1,nnb_topd_step
    WRITE(nunit,yform) jstp,xb_var_q(jstp,jcat,1:5)
  ENDDO
  !
  CALL close_file('ASCII ',nunit)
!  ENDIF 
  !
ENDDO !JCAT

yform='(i6,12f15.1)'
yfile=trim('bilan_tot.txt')
!
 CALL open_file('ASCII ',nunit,hfile=yfile,hform='FORMATTED',haction='WRITE')
!   
WRITE(nunit,*) '     T','         ',yb_var(1),'         ',yb_var(2),'         ',&
                                    yb_var(3),'         ',yb_var(4),'         ',&
                                    yb_var(5),'         ',yb_var(6),'         ',&
                                    yb_var(7),'         ',yb_var(8),'         ',&
                                    yb_var(9),'         ',yb_var(10),'         ',&
                                      yb_var(11),'         ',yb_var(12)

!
DO jstp=1,nnb_topd_step
  WRITE(nunit,yform) jstp,xb_var_tot(jstp,1:inb_var)
ENDDO
!
 CALL close_file('ASCII ',nunit)
!
IF (lhook) CALL dr_hook('WRITE_BUDGET_COUPL_ROUT',1,zhook_handle)
!
END SUBROUTINE write_budget_coupl_rout