PHY________OUT.f90

Go to the documentation of this file.

      subroutine phy________out(MOTIVATION)

!------------------------------------------------------------------------------+
!                                                         Sat 22-Jun-2013  MAR |
!   MAR          PHY________OUT                                                |
!     subroutine PHY________OUT OUTPUTs    MAR PHYsical parameterizations      |
!                                                                              |
!     version 3.p.4.1 created by H. Gallee,               Tue 12-Mar-2013      |
!           Last Modification by H. Gallee,               Sat 22-Jun-2013      |
!                                                                              |
!------------------------------------------------------------------------------+

      use mod_real
      use mod_phy____dat
      use mod_phy____grd
      use mod_phy_dy_kkl
      use mod_phy_cm_kkl
      use mod_phy_at_kkl
      use mod_sisvat_gpt


      IMPLICIT NONE


      character(len=50)                     ::  MOTIVATION
      integer                               ::       i,     j,   ikl            !
      integer                               ::       k,    kz                   !
      real(kind=real8)                      ::  RelHum




! OUTPUT
! ======


! txt file
! --------

      i      = i_x0
      j      = j_y0
      ikl    = ikl0

      write(4,50)       motivation
 50   format(3(/,1x),'  OUTPUT for VERIFICATION: ',a50,               &
     &         /,1x ,'  ***********************',/,1x)

      write(4,55) day_tu,mon_tu,yeartu,hourtu,minutu,sec_tu,it_exp
 55   format(3x,2(i2,'-'),i4,4x,3(i2,'-'),'  Simulation Iteration No ',i6)

      write(4,62)
      write(4,60)
 60   format(7x,'|','  sigma  ',' |','     Z   ',' |','    T    ',' |' &!
     &     ,'   U  ',' |','   V  ',' |',' RH   ',' |','  Qv   ',' |'   &!
     &                                ,'  Qw   ',' |','  Qi   ',' |'   &!
     &                                ,'  TKE  ',' |','  eps  ',' |'   &!
     &                                ,'  Kzh L',' |','  Kzh  ',' |'   &!
     &,/,'       |','         ',' |','   [km]  ',' |','   [K]   ',' |' &!
     &     ,' [m/s]',' |',' [m/s]',' |',' [-]  ',' |',' [g/kg]',' |'   &!
     &                                ,' [g/kg]',' |',' [g/kg]',' |'   &!
     &                                ,'[m2/s2]',' |','[m3/s2]',' |'   &!
     &                                ,' [m2/s]',' |',' [m2/s]',' |'   &!
     &,/,' ------+','---------','-+','---------','-+','---------','-+' &!
     &     ,'------','-+','------','-+','------','-+','-------','-+'   &!
     &                                ,'-------','-+','-------','-+'   &!
     &                                ,'-------','-+','-------','-+'   &!
     &                                ,'-------','-+','-------','-+')   !
      DO k = 1,mzp
      kz=min(k,mzp)
      relhum =      min(un_1,qv__dy(ikl,k)   /max(qvswcm(ikl,k),eps6))
      write(4,61) k,sigma(k),z___dy(ikl,k)*1.e-3 ,ta__dy(ikl,k)        &!
     &                      ,ua__dy(ikl,k)       ,va__dy(ikl,k)        &!
     &                      ,relhum       ,  1.e3*qv__dy(ikl,k)        &!
     &                 ,1.e3*qw__cm(ikl,kz), 1.e3*qi__cm(ikl,kz)       &!
     &                 ,     tke_at(ikl,kz),      eps_at(ikl,kz)       &!
     &                 ,     kzh_at(ikl,kz),      kzh0at(ikl,kz)
 61   format(i6,' |',f9.6,' |',  f9.4,' |' ,  f9.3,' |' ,2(f6.1,' |')  &!
     &              ,f6.2,' |',  f7.3,' |' ,2(f7.3,' |')               &!
     &                        ,2(f7.2,' |'),2(f7.3,' |')              )
      IF (mod(k,15).EQ.0)                                           THEN
      write(4,62)
 62   format(                                                          &!
     &   ' ------+','---------','-+','---------','-+','---------','-+' &!
     &     ,'------','-+','------','-+','------','-+','-------','-+'   &!
     &                                ,'-------','-+','-------','-+'   &!
     &                                ,'-------','-+','-------','-+'   &!
     &                                ,'-------','-+','-------','-+')   !

      write(4,60)
      END IF
      ENDDO

      write(4,62)
      write(4,63) uts_sv_gpt(ikl),1.e3*uqs_sv_gpt(ikl),us__sv_gpt(ikl) &
     &          , hsensv_gpt(ikl),     hlatsv_gpt(ikl)
 63   format(                                                          &!
     &       6x,' |', 9x ,' |', ( 9x ,' |'),' u*T*     |'              &!
     &                                     ,2( 6x ,' |'), ' u*q*  |'   &!
     &                        ,   7x ,' |' ,2( 7x ,' |'),'  u*    |'   &!
     &                        ,   7x ,' |' ,2( 7x ,' |')               &!
     &    ,/,6x,' |', 9x ,' |', ( 9x ,' |'),  f9.6,' |'                &!
     &                                     ,2( 6x ,' |'),f6.3,' |'     &!
     &                        ,   7x ,' |' ,2( 7x ,' |'),f7.2,' |'     &!
     &                        ,   7x ,' |' ,2( 7x ,' |')               &!
     &    ,/,6x,' |', 9x ,' |', ( 9x ,' |'),  f9.3,' |'                &!
     &                                     ,2( 6x ,' |'),f6.0,' |'     &!
     &                        ,   7x ,' |' ,2( 7x ,' |'), 7x ,' |'     &!
     &                        ,   7x ,' |' ,2( 7x ,' |')               &!
     &    ,/,6x,' |', 9x ,' |', ( 9x ,' |'),' [W/m2]   |'              &!
     &                                     ,2( 6x ,' |'), ' [W/m2]|'   &!
     &                        ,   7x ,' |' ,2( 7x ,' |'),'        |'   &!
     &                        ,   7x ,' |' ,2( 7x ,' |')            )

      k =      mzpp
      write(4,62)
      relhum =      min(un_1,qv__dy(ikl,k)   /max(qvswcm(ikl,k),eps6))
      write(4,64) k,sigma(k),z___dy(ikl,k)*1.e-3 ,ta__dy(ikl,k)        &!
     &                      ,zer0                ,zer0                 &!
     &                      ,relhum       ,  1.e3*qv__dy(ikl,k)
 64   format(i6,' |',f9.6,' |',  f9.4,' |' ,  f9.3,' |' ,2(f6.1,' |')  &!
     &              ,f6.2,' |',  f7.3,' |' ,2( 7x ,' |')               &!
     &                        ,2( 7x ,' |'),2( 7x ,' |')              )
      write(4,62)



! cdf file
! --------


      end subroutine phy________out