source: LMDZ5/branches/testing/libf/phymar/PHY________OUT.f90 @ 3525

Last change on this file since 3525 was 2160, checked in by Laurent Fairhead, 10 years ago

Merged trunk changes -r2070:2158 into testing branch. Compilation problems introduced by revision r2155 have been corrected by hand

File size: 6.0 KB
Line 
1      subroutine PHY________OUT(MOTIVATION)
2
3!------------------------------------------------------------------------------+
4!                                                         Sat 22-Jun-2013  MAR |
5!   MAR          PHY________OUT                                                |
6!     subroutine PHY________OUT OUTPUTs    MAR PHYsical parameterizations      |
7!                                                                              |
8!     version 3.p.4.1 created by H. Gallee,               Tue 12-Mar-2013      |
9!           Last Modification by H. Gallee,               Sat 22-Jun-2013      |
10!                                                                              |
11!------------------------------------------------------------------------------+
12
13      use Mod_Real
14      use Mod_PHY____dat
15      use Mod_PHY____grd
16      use Mod_PHY_DY_kkl
17      use Mod_PHY_CM_kkl
18      use Mod_PHY_AT_kkl
19      use Mod_SISVAT_gpt
20
21
22      IMPLICIT NONE
23
24
25      character(len=50)                     ::  MOTIVATION
26      integer                               ::       i,     j,   ikl            !
27      integer                               ::       k,    kz                   !
28      real(kind=real8)                      ::  RelHum
29
30
31
32
33! OUTPUT
34! ======
35
36
37! txt file
38! --------
39
40      i      = i_x0
41      j      = j_y0
42      ikl    = ikl0
43
44      write(4,50)       MOTIVATION
45 50   format(3(/,1x),'  OUTPUT for VERIFICATION: ',a50,               &
46     &         /,1x ,'  ***********************',/,1x)
47
48      write(4,55) Day_TU,Mon_TU,YearTU,HourTU,MinuTU,Sec_TU,it_EXP
49 55   format(3x,2(i2,'-'),i4,4x,3(i2,'-'),'  Simulation Iteration No ',i6)
50
51      write(4,62)
52      write(4,60)
53 60   format(7x,'|','  sigma  ',' |','     Z   ',' |','    T    ',' |' &!
54     &     ,'   U  ',' |','   V  ',' |',' RH   ',' |','  Qv   ',' |'   &!
55     &                                ,'  Qw   ',' |','  Qi   ',' |'   &!
56     &                                ,'  TKE  ',' |','  eps  ',' |'   &!
57     &                                ,'  Kzh L',' |','  Kzh  ',' |'   &!
58     &,/,'       |','         ',' |','   [km]  ',' |','   [K]   ',' |' &!
59     &     ,' [m/s]',' |',' [m/s]',' |',' [-]  ',' |',' [g/kg]',' |'   &!
60     &                                ,' [g/kg]',' |',' [g/kg]',' |'   &!
61     &                                ,'[m2/s2]',' |','[m3/s2]',' |'   &!
62     &                                ,' [m2/s]',' |',' [m2/s]',' |'   &!
63     &,/,' ------+','---------','-+','---------','-+','---------','-+' &!
64     &     ,'------','-+','------','-+','------','-+','-------','-+'   &!
65     &                                ,'-------','-+','-------','-+'   &!
66     &                                ,'-------','-+','-------','-+'   &!
67     &                                ,'-------','-+','-------','-+')   !
68      DO k = 1,mzp
69      kz=min(k,mzp)
70      RelHum =      min(un_1,qv__DY(ikl,k)   /max(qvswCM(ikl,k),eps6))
71      write(4,61) k,sigma(k),Z___DY(ikl,k)*1.e-3 ,Ta__DY(ikl,k)        &!
72     &                      ,Ua__DY(ikl,k)       ,Va__DY(ikl,k)        &!
73     &                      ,RelHum       ,  1.e3*qv__DY(ikl,k)        &!
74     &                 ,1.e3*qw__CM(ikl,kz), 1.e3*qi__CM(ikl,kz)       &!
75     &                 ,     TKE_AT(ikl,kz),      eps_AT(ikl,kz)       &!
76     &                 ,     Kzh_AT(ikl,kz),      Kzh0AT(ikl,kz)
77 61   format(i6,' |',f9.6,' |',  f9.4,' |' ,  f9.3,' |' ,2(f6.1,' |')  &!
78     &              ,f6.2,' |',  f7.3,' |' ,2(f7.3,' |')               &!
79     &                        ,2(f7.2,' |'),2(f7.3,' |')              )
80      IF (mod(k,15).EQ.0)                                           THEN
81      write(4,62)
82 62   format(                                                          &!
83     &   ' ------+','---------','-+','---------','-+','---------','-+' &!
84     &     ,'------','-+','------','-+','------','-+','-------','-+'   &!
85     &                                ,'-------','-+','-------','-+'   &!
86     &                                ,'-------','-+','-------','-+'   &!
87     &                                ,'-------','-+','-------','-+')   !
88
89      write(4,60)
90      END IF
91      ENDDO
92
93      write(4,62)
94      write(4,63) uts_SV_gpt(ikl),1.e3*uqs_SV_gpt(ikl),us__SV_gpt(ikl) &
95     &          , HsenSV_gpt(ikl),     HLatSV_gpt(ikl)
96 63   format(                                                          &!
97     &       6x,' |', 9x ,' |', ( 9x ,' |'),' u*T*     |'              &!
98     &                                     ,2( 6x ,' |'), ' u*q*  |'   &!
99     &                        ,   7x ,' |' ,2( 7x ,' |'),'  u*    |'   &!
100     &                        ,   7x ,' |' ,2( 7x ,' |')               &!
101     &    ,/,6x,' |', 9x ,' |', ( 9x ,' |'),  f9.6,' |'                &!
102     &                                     ,2( 6x ,' |'),f6.3,' |'     &!
103     &                        ,   7x ,' |' ,2( 7x ,' |'),f7.2,' |'     &!
104     &                        ,   7x ,' |' ,2( 7x ,' |')               &!
105     &    ,/,6x,' |', 9x ,' |', ( 9x ,' |'),  f9.3,' |'                &!
106     &                                     ,2( 6x ,' |'),f6.0,' |'     &!
107     &                        ,   7x ,' |' ,2( 7x ,' |'), 7x ,' |'     &!
108     &                        ,   7x ,' |' ,2( 7x ,' |')               &!
109     &    ,/,6x,' |', 9x ,' |', ( 9x ,' |'),' [W/m2]   |'              &!
110     &                                     ,2( 6x ,' |'), ' [W/m2]|'   &!
111     &                        ,   7x ,' |' ,2( 7x ,' |'),'        |'   &!
112     &                        ,   7x ,' |' ,2( 7x ,' |')            )
113
114      k =      mzpp
115      write(4,62)
116      RelHum =      min(un_1,qv__DY(ikl,k)   /max(qvswCM(ikl,k),eps6))
117      write(4,64) k,sigma(k),Z___DY(ikl,k)*1.e-3 ,Ta__DY(ikl,k)        &!
118     &                      ,zer0                ,zer0                 &!
119     &                      ,RelHum       ,  1.e3*qv__DY(ikl,k)
120 64   format(i6,' |',f9.6,' |',  f9.4,' |' ,  f9.3,' |' ,2(f6.1,' |')  &!
121     &              ,f6.2,' |',  f7.3,' |' ,2( 7x ,' |')               &!
122     &                        ,2( 7x ,' |'),2( 7x ,' |')              )
123      write(4,62)
124
125
126
127! cdf file
128! --------
129
130
131      end subroutine PHY________OUT
Note: See TracBrowser for help on using the repository browser.