source: LMDZ5/branches/IPSLCM6.0.11pre/libf/phymar/PHY_MAR___.f90 @ 5434

Last change on this file since 5434 was 2104, checked in by fhourdin, 10 years ago

Corrections de phymar par Gilles Delaygue.

File size: 73.6 KB
Line 
1      subroutine PHY_MAR                                               &
2
3!------------------------------------------------------------------------------+
4!                                                         Mon  1-Jul-2013  MAR |
5!     subroutine PHY_MAR is the MAR PHYsics Driver                             |
6!     interfaces HOST   variables                                              |
7!            and SISVAT variables                                              |
8!                                                                              |
9!     Applied to: MARthusalem               (variables in MAR***.inc files)    |
10!                                                                              |
11!                                                                              |
12! # OPTIONS: #dT  Distinction among  Tendencies of MAR  Physical Parametr.     |
13! # ^^^^^^^^ #cw  Cloud Condensation Nuclei (CCNw) Microphysics Activation     |
14!                                                                              |
15!                                                                              |
16!     version 3.p.4.1 created by H. Gallee,               Tue 12-Mar-2013      |
17!           Last Modification by H. Gallee,               Mon  1-Jul-2013      |
18!                                                                              |
19!------------------------------------------------------------------------------+
20
21     &                  (FlagSV                                        &   ! FLAG  for SISVAT: (T,F) =                    (active OR NOT)
22     &                  ,FlagSV_Veg                                    &   ! FLAG  for SISVAT: (T,F) =       (Variable Vegetation OR NOT)
23     &                  ,FlagSV_SNo                                    &   ! FLAG  for SISVAT: (T,F) =         (Snow Model active OR NOT)
24     &                  ,FlagSV_BSn                                    &   ! FLAG  for SISVAT: (T,F) = (Blowing Snow Model active OR NOT)
25     &                  ,FlagSV_KzT                                    &   ! FLAG  for SISVAT: (T,F) = (pkt Turb.Transfert active OR NOT in SISVAT)
26     &                  ,FlagSV_SWD                                    &   ! FLAG  for SISVAT: (T,F) = (Modify SW INPUT->downward OR NOT)     
27     &                  ,FlagSV_SBC                                    &   ! FLAG  for SISVAT: (T,F) = (INPUT of Soil & Vege DATA OR NOT in SISVAT)
28     &                  ,FlagSV_UBC                                    &   ! FLAG  for SISVAT: (T,F) = (pkt UpperBC is Von Neuman OR NOT in SISVAT)
29     &                  ,FlagAT                                        &   ! FLAG  for Atm_AT: (T,F) = (Turbulent Transfer active OR NOT)
30     &                  ,TypeAT                                        &   ! TYPE  of  Atm_AT: (e= Ee Duynkerke, K= Ee Kitada, L= EL, H= Ee Huan-R)
31     &                  ,FlagAT_TKE                                    &   ! FLAG  for genTKE: (T,F) = (TKE-e     Model    active OR NOT)
32     &                  ,FlagCM                                        &   ! FLAG  for CMiPhy: (T,F) = (Cloud Microphysics active OR NOT)
33     &                  ,FlagCM_UpD                                    &   ! FLAG  for CMiPhy: (T,F) = (qv & hydrometeors updated OR NOT IN CMiPhy)
34     &                  ,FlagCP                                        &   ! FLAG  for Convection Paramet.
35     &                  ,FlagRT                                        &   ! FLAG  for Radiative Transfer
36     &                  ,FlagS0_SLO                                    &   ! FLAG  for Insolation, Surfac.Slope                 Impact  included NEW
37     &                  ,FlagS0_MtM                                    &   ! FLAG  for Insolation, Surfac.Slope & Mountain Mask Impacts included NEW
38     &                  ,Flag_O                                        &   ! FLAG  for OUTPUT
39     &                  ,FlagVR                                        &   ! FLAG  for OUTPUT for VERIFICATION
40     &                  ,dt0DYn                                        &   ! Time STEP between 2 CALLs of PHY_MAR                     [s] I, fix
41     &                  ,dt0_SV                                        &   ! Time STEP between 2 CALLs of SISVAT                      [s] I, fix
42     &                  ,dt0_AT                                        &   ! Time STEP between 2 CALLs of Atm_AT                      [s] I, fix
43     &                  ,dt0_CM                                        &   ! Time STEP between 2 CALLs of CMiPhy                      [s] I, fix
44     &                  ,dt0_CP                                        &   ! Time STEP between 2 CALLs of CVamnh                      [s] I, fix
45     &                  ,dt0_RT                                        &   ! Time STEP between 2 CALLs of radCEP                      [s] I, fix
46     &                  ,dx                                            &   ! Grid  Mesh size (Horizontal)                             [m] I, fix
47     &                  ,DD_AxX                                        &   ! Grid  x-Axis Direction                              [degree] I, fix
48     &                  ,s_HOST                                        &   ! Grid (Vertical)   of HOST (NORMALIZED PRESSURE assumed)  [-] I, fix
49     &                  ,sh___HOST                                     &   ! Topography                                               [m] I, fix
50     &                  ,sh_a_HOST                                     &   ! Topography Anomaly                                       [m] I, fix NEW
51     &                  ,slopxHOST                                     &   ! Slope, x-direction                                       [-] I, fix NEW
52     &                  ,slopyHOST                                     &   ! Slope, y-direction                                       [-] I, fix NEW
53     &                  ,slopeHOST                                     &   ! Slope                                                    [-] I, fix NEW
54     &                  ,MMaskHOST                                     &   ! Mountain Mask                                            [-] I, fix NEW
55     &                  ,lonh_HOST                                     &   ! Longitude                                             [hour] I, fix
56     &                  ,latr_HOST                                     &   ! Latitude                                            [radian] I, fix
57     &                  ,pkta_HOST                                     &   ! Reduced Potential Temperature                           [XK] I, O
58     &                  ,ptop_HOST                                     &   ! Pressure, Model Top                                    [kPa] I, fix
59     &                  ,psa__HOST                                     &   ! Pressure  Thickness                                    [kPa] I
60     &                  ,gZa__HOST                                     &   ! Geopotential Height                                  [m2/s2] I
61     &                  ,gZam_HOST                                     &   ! Geopotential Height, mid-level                       [m2/s2] I
62     &                  ,Ua___HOST                                     &   ! Wind ,  x-Direction                                    [m/s] I
63     &                  ,Va___HOST                                     &   ! Wind ,  y-Direction                                    [m/s] I
64     &                  ,Wa___HOST                                     &   ! Wind ,  z-Direction                                    [m/s] I
65     &                  ,qv___HOST                                     &   ! Specific  Humidity                                   [kg/kg] I, O
66     &                  ,qw___HOST                                     &   ! Cloud Droplets Concentration                         [kg/kg] I, O
67! #cw&                  ,CCN__HOST                                     &   ! CCN            Concentration                          [-/kg]
68     &                  ,qi___HOST                                     &   ! Cloud Crystals Concentration                         [kg/kg] I, O
69     &                  ,CIN__HOST                                     &   ! CIN            Concentration                          [-/kg] I, O
70     &                  ,CF___HOST                                     &   ! Cloud Fraction                                           [-] I, O
71     &                  ,qs___HOST                                     &   ! Snow Particles Concentration                         [kg/kg] I, O
72     &                  ,qr___HOST                                     &   ! Rain Drops     Concentration                         [kg/kg] I, O
73     &                  ,TKE__HOST                                     &   ! Turbulent Kinetic Energy                             [m2/s2] I, O
74     &                  ,eps__HOST                                     &   ! Turbulent Kinetic Energy Dissipation                 [m2/s3] I, O
75     &                  ,dpkt___dt                                     &   ! Reduced Potential Temperature TENDENCY, ALL Contribut.[KX/s]    O
76     &                  ,dua____dt                                     &   ! Wind Speed       (x-direc.)   TENDENCY, ALL Contribut.[m/s2]    O
77     &                  ,dva____dt                                     &   ! Wind Speed       (y-direc.)   TENDENCY, ALL Contribut.[m/s2]    O
78     &                  ,dqv____dt                                     &   ! Specific          Humidity    TENDENCY, ALL Contr. [kg/kg/s]    O
79     &                  ,dqw____dt                                     &   ! Cloud Droplets Concentration  TENDENCY, ALL Contr. [kg/kg/s]    O
80! #cw&                  ,dCw____dt                                     &   ! CCN            Concentration  TENDENCY, ALL Contr.     [1/s]
81     &                  ,dqi____dt                                     &   ! Cloud Crystals Concentration  TENDENCY, ALL Contr. [kg/kg/s]    O
82     &                  ,dCi____dt                                     &   ! CIN            Concentration  TENDENCY, ALL Contr.     [1/s]    O
83     &                  ,dCF____dt                                     &   ! Cloud Fraction                TENDENCY, ALL Contr.     [1/s]    O
84     &                  ,dqs____dt                                     &   ! Snow Particles Concentration  TENDENCY, ALL Contr. [kg/kg/s]    O
85     &                  ,dqr____dt                                     &   ! Rain Drops     Concentration  TENDENCY, ALL Contr. [kg/kg/s]   O
86! #dT&                  ,dpktSV_dt                                     &   ! Reduced Potential Temperature TENDENCY, SISVAT        [KX/s]  (O)
87! #dT&                  ,dpktAT_dt                                     &   ! Reduced Potential Temperature TENDENCY, Atm_AT        [KX/s]  (O)
88! #dT&                  ,dqv_AT_dt                                     &   ! Specific          Humidity    TENDENCY, Atm_AT     [kg/kg/s]  (O)
89! #dT&                  ,dqw_AT_dt                                     &   ! Cloud Droplets Concentration  TENDENCY, Atm_AT     [kg/kg/s]  (O)
90! #dT&                  ,dqi_AT_dt                                     &   ! Cloud Crystals Concentration  TENDENCY, Atm_AT     [kg/kg/s]  (O)
91! #dT&                  ,dqs_AT_dt                                     &   ! Snow Particles Concentration  TENDENCY, Atm_AT     [kg/kg/s]  (O)
92! #dT&                  ,dqr_AT_dt                                     &   ! Rain Drops     Concentration  TENDENCY, Atm_AT     [kg/kg/s]  (O)
93! #cw&                  ,dCw_AT_dt                                     &   ! CCN            Concentration  TENDENCY, Atm_AT         [1/s]  (O)
94! #dT&                  ,dCi_AT_dt                                     &   ! CIN            Concentration  TENDENCY, Atm_AT         [1/s]  (O)
95! #dT&                  ,dpktCM_dt                                     &   ! Reduced Potential Temperature TENDENCY, CMiPhy        [KX/s]  (O)
96! #dT&                  ,dqv_CM_dt                                     &   ! Specific          Humidity    TENDENCY, CMiPhy     [kg/kg/s]  (O)
97! #dT&                  ,dqw_CM_dt                                     &   ! Cloud Droplets Concentration  TENDENCY, CMiPhy     [kg/kg/s]  (O)
98! #dT&                  ,dCF_CM_dt                                     &   ! Cloud Fraction                TENDENCY, CMiPhy         [1/s]  (O)
99! #dT&                  ,dqi_CM_dt                                     &   ! Cloud Crystals Concentration  TENDENCY, CMiPhy     [kg/kg/s]  (O)
100! #dT&                  ,dqs_CM_dt                                     &   ! Snow Particles Concentration  TENDENCY, CMiPhy     [kg/kg/s]  (O)
101! #dT&                  ,dqr_CM_dt                                     &   ! Rain Drops     Concentration  TENDENCY, CMiPhy     [kg/kg/s]  (O)
102! #cw&                  ,dCw_CM_dt                                     &   ! CCN            Concentration  TENDENCY, CMiPhy         [1/s]  (O)
103! #dT&                  ,dCi_CM_dt                                     &   ! CIN            Concentration  TENDENCY, CMiPhy         [1/s]  (O)
104! #dT&                  ,dpktCP_dt                                     &   ! Reduced Potential Temperature TENDENCY, CVAmnh        [KX/s]  (O)
105! #dT&                  ,dqv_CP_dt                                     &   ! Specific          Humidity    TENDENCY, CVAmnh     [kg/kg/s]  (O)
106! #dT&                  ,dqw_CP_dt                                     &   ! Cloud Droplets Concentration  TENDENCY, CVAmnh     [kg/kg/s]  (O)
107! #dT&                  ,dqi_CP_dt                                     &   ! Cloud Crystals Concentration  TENDENCY, CVAmnh     [kg/kg/s]  (O)
108! #dT&                  ,dpktRT_dt                                     &   ! Reduced Potential Temperature TENDENCY, radCEP        [KX/s]  (O)
109     &                  ,sst__HOST                                     &   ! Ocean     FORCING (SST)                                  [K] I
110! #IP&                  ,sif__HOST                                     &   ! Ocean     FORCING (Sea-Ice Fraction )                    [-] I
111! #AO&                  ,s_T__HOST                                     &   ! Ocean    COUPLING (Surface Temperat.)  n=1: Open Ocean   [-] I,NEMO
112! #AO&                  ,Alb__HOST                                     &   ! Ocean    COUPLING (Surface Albedo   )  n=2: Sea  Ice     [-] I,NEMO
113! #AO&                  ,dSdT2HOST                                     &   ! Ocean    COUPLING ( d(SH Flux) / dT )               [W/m2/K]   O
114! #AO&                  ,dLdT2HOST                                     &   ! Ocean    COUPLING ( d(LH Flux) / dT )               [W/m2/K]   O
115!dead&                  ,it0EXP,it0RUN                                 &   ! Iteration
116     &                  ,Year_H,Mon__H,Day__H,Hour_H,minu_H,sec__H     &   ! Time
117     &                  ,ixq1  ,i0x0  ,mxqq                            &   ! Domain  Dimension: x
118     &                  ,jyq1  ,j0y0  ,myqq                            &   ! Domain  Dimension: y
119     &                  ,mzq   ,mzqq                                   &   ! Domain  Dimension: z
120     &                  ,mwq                                           &   ! Domain  Dimension: mosaic
121     &                  ,kcolq                                         &   ! Domain  Dimension: x * y
122     &                  ,kcolw                                         &   ! Domain  Dimension: x * y * mosaic
123     &                  ,m_azim                                        &   ! Mountain Mask, nb of directions taken into account       [-]
124     &                  ,IOi0SV,IOj0SV,n0pt)                               ! Indices of OUTPUT Grid Point
125
126!------------------------------------------------------------------------------+
127!                                                         Sat 29-Jun-2013  MAR |
128!     subroutine PHY_MAR is the MAR PHYsics Driver                             |
129!     interfaces HOST   variables                                              |
130!            and SISVAT variables                                              |
131!                                                                              |
132!     Applied to: MARthusalem               (variables in MAR***.inc files)    |
133!                                                                              |
134!                                                                              |
135! # OPTIONS: #dT  Distinction among  Tendencies of MAR  Physical Parametr.     |
136! # ^^^^^^^^ #cw  Cloud Condensation Nuclei (CCNw) Microphysics Activation     |
137!                                                                              |
138!                                                                              |
139!     version 3.p.4.1 created by H. Gallee,               Tue 12-Mar-2013      |
140!           Last Modification by H. Gallee,               Sat 29-Jun-2013      |
141!                                                                              |
142!------------------------------------------------------------------------------+
143
144      use Mod_Real
145      use Mod_PHY____dat
146      use Mod_PHY____grd
147      use Mod_PHY____kkl
148      use Mod_PHY_CM_ctr
149      use Mod_PHY_S0_ctr
150      use Mod_SISVAT_ctr
151      use Mod_PHY_CM_dat
152      use Mod_PHY_AT_grd
153      use Mod_PHY_CM_grd
154      use Mod_PHY_CP_grd
155      use Mod_PHY_RT_grd
156      use Mod_PHY_S0_grd
157      use Mod_SISVAT_grd
158      use Mod_PHY_DY_kkl
159      use Mod_PHY_AT_kkl
160      use Mod_PHY_CM_kkl
161      use Mod_PHY_CP_kkl
162      use Mod_PHY_RT_kkl
163      use Mod_PHY_S0_kkl
164      use Mod_SISVAT_kkl
165      use Mod_SISVAT_gpt
166
167      IMPLICIT NONE
168
169      logical                                               ::  FlagSV             !  Flag         (SISVAT)
170      logical                                               ::  FlagSV_Veg         !  Flag         (SISVAT: Vegetation)
171      logical                                               ::  FlagSV_SNo         !  Flag         (SISVAT: Surface * )
172      logical                                               ::  FlagSV_BSn         !  Flag         (SISVAT: Blowing * )
173      logical                                               ::  FlagSV_KzT         !  Flag         (SISVAT: d(KdT/dz)/dz)
174      logical                                               ::  FlagSV_SWD         !  Flag: T/F :  (SISVAT: SW=Down/Abs.)
175      logical                                               ::  FlagSV_SBC         !  Flag: T/F :  (SISVAT: SBC=INP/FIX.)
176      logical                                               ::  FlagSV_UBC         !  Flag: T/F :  (SISVAT: UBC=VonN/Dr.)
177      logical                                               ::  FlagAT             !  Flag         (Turbulent Transfer)
178      character(len=1)                                      ::  TypeAT             !  Type         (Turbulent Transfer)
179      logical                                               ::  FlagAT_TKE         !  Flag         (Turbulent Transfer, TKE-e Model: ON / OFF)
180      logical                                               ::  FlagCM             !  Flag         (Cloud Microphysics)
181      logical                                               ::  FlagCM_UpD         !  Flag: T/F :  (Cloud Microphysics: Update in/out PHY_MAR)
182      logical                                               ::  FlagCP             !  Flag         (Convection Param. )
183      logical                                               ::  FlagRT             !  Flag         (Radiative Transfer)
184      logical                                               ::  FlagS0_SLO         !  FLAG         (Insolation, Surfac.Slope                )
185      logical                                               ::  FlagS0_MtM         !  FLAG         (Insolation, Surfac.Slope & Mountain Mask)
186
187      logical                                               ::  Flag_O             !  Flag         (OUTPUT)
188      logical                                               ::  FlagVR             !  Flag         (OUTPUT for VERIFICATION)
189
190      real                                                  ::  dt0DYn             !  Time   Step  (DYnamics, the shortest)                     [s]
191      real                                                  ::  dt0_SV             !  Time   Step  (SISVAT)                                     [s]
192      real                                                  ::  dt0_AT             !  Time   Step  (Atmo Turb.)                                 [s]
193      real                                                  ::  dt0_CM             !  Time   Step  (Cloud Mic.)                                 [s]
194      real                                                  ::  dt0_CP             !  Time   Step  (Convection)                                 [s]
195      real                                                  ::  dt0_RT             !  Time   Step  (Radiat.Tr.)                                 [s]
196      real                                                  ::  dx                 !  Grid   Size                                               [m]
197      real                                                  ::  DD_AxX             !  x-Axis Direction                                     [degree]
198      real, dimension(mzqq)                                 ::  s_HOST             !  Vertical Coordinate                                       [-]
199      real            , dimension(kcolq)                    ::  sh___HOST          !  Topography                                                [m]
200      real(kind=real8), dimension(kcolq)                    ::  sh_a_HOST          !  Topography Anomaly                                        [m]
201      real(kind=real8), dimension(kcolq)                    ::  slopxHOST          !  Slope, x-direction                                        [-]
202      real(kind=real8), dimension(kcolq)                    ::  slopyHOST          !  Slope, y-direction                                        [-]
203      real(kind=real8), dimension(kcolq)                    ::  slopeHOST          !  Slope                                                     [-]
204      real(kind=real8), dimension(kcolq,m_azim)             ::  MMaskHOST          !  Mountain Mask                                             [-]
205      real, dimension(kcolq)                                ::  lonh_HOST          !  Longitude                                              [hour]
206      real, dimension(kcolq)                                ::  latr_HOST          !   Latitude                                            [radian]
207      real                                                  ::  ptop_HOST          !  Pressure Model Top                                      [kPa]
208
209      real, dimension(kcolq,mzqq)                           ::  pkta_HOST          !  Reduced  Potential Temperature                         [KX/s]
210      real, dimension(kcolq)                                ::  psa__HOST          !  Pressure Thickness                                      [kPa]
211      real, dimension(kcolq,mzqq)                           ::  gZa__HOST          !  Geopotential Height                                   [m2/s2]
212      real, dimension(kcolq,mzqq)                           ::  gZam_HOST          !  Geopotential Height, mid-level                        [m2/s2]
213      real, dimension(kcolq,mzq)                            ::  Ua___HOST          !  Wind Speed, x-direction                                 [m/s]
214      real, dimension(kcolq,mzq)                            ::  Va___HOST          !  Wind Speed, y-direction                                 [m/s]
215      real, dimension(kcolq,mzq) ,INTENT(IN)                ::  Wa___HOST          !  Wind Speed, z-direction                                 [m/s]
216      real, dimension(kcolq,mzqq)                           ::  qv___HOST          !  Specific Humidity                                     [kg/kg]
217      real, dimension(kcolq,mzq)                            ::  qw___HOST          !  Cloud Droplets Concentration                          [kg/kg]
218      real, dimension(kcolq,mzq)                            ::  CCN__HOST          !  CCN            Concentration                           [-/kg]
219      real, dimension(kcolq,mzq)                            ::  qi___HOST          !  Cloud Crystals Concentration                          [kg/kg]
220      real, dimension(kcolq,mzq)                            ::  CIN__HOST          !  CIN            Concentration                           [-/kg]
221      real, dimension(kcolq,mzq)                            ::  CF___HOST          !  Cloud Fraction                                         [-/kg]
222      real, dimension(kcolq,mzq)                            ::  qs___HOST          !  Snow Particles Concentration                          [kg/kg]
223      real, dimension(kcolq,mzq)                            ::  qr___HOST          !  Rain Drops     Concentration                          [kg/kg]
224      real, dimension(kcolq,mzq)                            ::  TKE__HOST          !  Turbulent Kinetic Energy                              [m2/s2]
225      real, dimension(kcolq,mzq)                            ::  eps__HOST          !  Turbulent Kinetic Energy Dissipation                  [m2/s3]
226
227      real, dimension(kcolq,mzq)                            ::  dpkt___dt          !  Reduced  Potential Temperature TENDENCY, ALL Contribut.[KX/s]
228      real, dimension(kcolq,mzq)                            ::  dua____dt          !  Wind Speed       (x-direc.)    TENDENCY, ALL Contribut.[m/s2]
229      real, dimension(kcolq,mzq)                            ::  dva____dt          !  Wind Speed       (y-direc.)    TENDENCY, ALL Contribut.[m/s2]
230      real, dimension(kcolq,mzq)                            ::  dqv____dt          !  Specific           Humidity    TENDENCY, ALL Contr. [kg/kg/s]
231      real, dimension(kcolq,mzq)                            ::  dqw____dt          !  Cloud Droplets Concentration   TENDENCY, ALL Contr. [kg/kg/s]
232! #cw real, dimension(kcolq,mzq)                            ::  dCw____dt          !  CCN            Concentration   TENDENCY, ALL Contr.  [1/kg/s]
233      real, dimension(kcolq,mzq)                            ::  dqi____dt          !  Cloud Crystals Concentration   TENDENCY, ALL Contr. [kg/kg/s]
234      real, dimension(kcolq,mzq)                            ::  dCi____dt          !  CIN            Concentration   TENDENCY, ALL Contr.  [1/kg/s]
235      real, dimension(kcolq,mzq)                            ::  dCF____dt          !  Cloud Fraction                 TENDENCY, ALL Contr. [kg/kg/s]
236      real, dimension(kcolq,mzq)                            ::  dqs____dt          !  Snow Particles Concentration   TENDENCY, ALL Contr. [kg/kg/s]
237      real, dimension(kcolq,mzq)                            ::  dqr____dt          !  Rain Drops     Concentration   TENDENCY, ALL Contr. [kg/kg/s]
238
239      real, dimension(kcolq,mzq)                            ::  dpktSV_dt          !  Reduced  Potential Temperature Tendency, SISVAT        [KX/s]
240      real, dimension(kcolq,mzq)                            ::  dpktAT_dt          !  Reduced  Potential Temperature Tendency, Atm_AT        [KX/s]
241      real, dimension(kcolq,mzq)                            ::  dqv_AT_dt          !  Specific           Humidity    TENDENCY, Atm_AT     [kg/kg/s]
242      real, dimension(kcolq,mzq)                            ::  dqw_AT_dt          !  Cloud Droplets Concentration   TENDENCY, Atm_AT     [kg/kg/s]
243! #cw real, dimension(kcolq,mzq)                            ::  dCw_AT_dt          !  CCN            Concentration   TENDENCY, Atm_AT         [1/s]
244      real, dimension(kcolq,mzq)                            ::  dqi_AT_dt          !  Cloud Crystals Concentration   TENDENCY, Atm_AT     [kg/kg/s]
245      real, dimension(kcolq,mzq)                            ::  dCi_AT_dt          !  CIN            Concentration   TENDENCY, Atm_AT         [1/s]
246      real, dimension(kcolq,mzq)                            ::  dqs_AT_dt          !  Snow Particles Concentration   TENDENCY, Atm_AT     [kg/kg/s]
247      real, dimension(kcolq,mzq)                            ::  dqr_AT_dt          !  Rain Drops     Concentration   TENDENCY, Atm_AT     [kg/kg/s]
248      real, dimension(kcolq,mzq)                            ::  dpktCM_dt          !  Reduced  Potential Temperature Tendency, CMiPhy        [KX/s]
249      real, dimension(kcolq,mzq)                            ::  dqv_CM_dt          !  Specific           Humidity    TENDENCY, CMiPhy     [kg/kg/s]
250      real, dimension(kcolq,mzq)                            ::  dqw_CM_dt          !  Cloud Droplets Concentration   TENDENCY, CMiPhy     [kg/kg/s]
251! #cw real, dimension(kcolq,mzq)                            ::  dCw_CM_dt          !  CCN            Concentration   TENDENCY, CMiPhy         [1/s]
252      real, dimension(kcolq,mzq)                            ::  dqi_CM_dt          !  Cloud Crystals Concentration   TENDENCY, CMiPhy     [kg/kg/s]
253      real, dimension(kcolq,mzq)                            ::  dCi_CM_dt          !  CIN            Concentration   TENDENCY, CMiPhy         [1/s]
254      real, dimension(kcolq,mzq)                            ::  dCF_CM_dt          !  Cloud Fraction                 TENDENCY, CMiPhy         [1/s]
255      real, dimension(kcolq,mzq)                            ::  dqs_CM_dt          !  Snow Particles Concentration   TENDENCY, CMiPhy     [kg/kg/s]
256      real, dimension(kcolq,mzq)                            ::  dqr_CM_dt          !  Rain Drops     Concentration   TENDENCY, CMiPhy     [kg/kg/s]
257      real, dimension(kcolq,mzq)                            ::  dpktCP_dt          !  Reduced  Potential Temperature TENDENCY, CVAmnh        [KX/s]
258      real, dimension(kcolq,mzq)                            ::  dqv_CP_dt          !  Specific           Humidity    TENDENCY, CVAmnh     [kg/kg/s]
259      real, dimension(kcolq,mzq)                            ::  dqw_CP_dt          !  Cloud Droplets Concentration   TENDENCY, CVAmnh     [kg/kg/s]
260      real, dimension(kcolq,mzq)                            ::  dqi_CP_dt          !  Cloud Crystals Concentration   TENDENCY, CVAmnh     [kg/kg/s]
261      real, dimension(kcolq,mzq)                            ::  dpktRT_dt          !  Reduced  Potential Temperature TENDENCY, radCEP        [KX/s]
262
263!dead integer                                               ::  it0EXP             !
264!dead integer                                               ::  it0RUN             !
265      integer                                               ::  Year_H             !  Time                                                   [year]
266      integer                                               ::  Mon__H             !  Time                                                  [month]
267      integer                                               ::  Day__H             !  Time                                                    [Day]
268      integer                                               ::  Hour_H             !  Time                                                   [hour]
269      integer                                               ::  minu_H             !  Time                                                 [minute]
270      integer                                               ::  sec__H             !  Time                                                      [s]
271      integer                                               ::  ixq1,i0x0,mxqq     !  Domain  Dimension: x                                      [-]
272      integer                                               ::  jyq1,j0y0,myqq     !  Domain  Dimension: y                                      [-]
273      integer                                               ::  mzq                !  Domain  Dimension: z                                      [-]
274      integer                                               ::  mzqq               !  Domain  Dimension: z                                      [-]
275      integer                                               ::  mwq                !  Domain  Dimension: mosaic                                 [-]
276      integer                                               ::  kcolq              !  Domain  Dimension: x * y                                  [-]
277      integer                                               ::  kcolw              !  Domain  Dimension: x * y * mosaic                         [-]
278      integer                                               ::  m_azim             !  Mountain Mask, nb of directions taken into account        [-]
279      integer, dimension(n0pt)                              ::  IOi0SV             !
280      integer, dimension(n0pt)                              ::  IOj0SV             !
281      integer                                               ::  n0pt               !
282
283      real, dimension(kcolq)                                ::  sst__HOST          !  Ocean     FORCING (SST)                                   [K]
284! #IP real, dimension(kcolq)                                ::  sif__HOST          !  Ocean     FORCING (Sea-Ice Fraction )                     [-]
285! #AO real, dimension(ixq1:mxqq,jyq1:myqq,mwq)              ::  s_T__HOST          !  A - O    COUPLING                      n=1: Open Ocean    [K]
286! #AO real, dimension(ixq1:mxqq,jyq1:myqq,mwq)              ::  Alb__HOST          !  A - O    COUPLING (Surface Albedo   )  n=2: Sea  Ice      [-]
287! #AO real, dimension(ixq1:mxqq,jyq1:myqq,mwq)              ::  dSdT2HOST          !  A - O    COUPLING ( d(SH Flux) / dT )                [W/m2/K]
288! #AO real, dimension(ixq1:mxqq,jyq1:myqq,mwq)              ::  dLdT2HOST          !  A - O    COUPLING ( d(SH Flux) / dT )                [W/m2/K]
289
290! #TC integer, parameter                                    ::   ntrac  =  28      !
291! #TC real, dimension(ixq1:mxqq,jyq1:myqq,mzq,ntrac)        ::    qxTC             !  Aerosols: Atmospheric  Contentration
292! #TC real, dimension(ixq1:mxqq,jyq1:myqq    ,ntrac)        ::    qsTC             !  Aerosols: Near Surface Contentration
293! #TC real, dimension(ixq1:mxqq,jyq1:myqq    ,ntrac)        ::    uqTC             !  Aerosols: Surf.Flux
294! ---------------------------------------------------------------------------------!
295
296
297
298
299! LOCAL VARIABLES
300! ===============
301
302      real(kind=real8), dimension(kcolq,mzq)                 ::  Wind_HOST            !  Wind Speed, Horizontal                  [m/s]
303      real(kind=real8)                                       ::  dTdz   = 0.0065      ! -d(T) / dz   Lapse Rate                  [K/m]
304      real(kind=real8)                                       ::  dTimAT               !  d(Time) between 2 calls of Atmos.Turbul.  [s]
305
306      integer                                                ::  i   ,j   ,ikl   ,ikp !
307      integer                                                ::  k   ,mn              !
308      integer                                                ::  n   ,ipt ,iwr   ,kk  !
309
310
311
312
313
314!dead     it_RUN = it0RUN
315!dead     it_EXP = it0EXP
316
317
318
319
320! INITIALIZATION
321! ==============
322
323! Initialization of local Variables (used each Time Step)
324! -------------------------------------------------------
325
326      DO k = 1,mzq
327      DO ikl=1,kcolq 
328        Wind_HOST(ikl,k) = sqrt(Ua___HOST(ikl,k)*Ua___HOST(ikl,k)+Va___HOST(ikl,k)*Va___HOST(ikl,k))
329      ENDDO
330      ENDDO
331
332
333
334! Initialization of the run
335! -------------------------
336
337! Martin Control
338     PRINT*, 'Dans PHY_MAR:'
339     PRINT*, 'it_RUN=',it_RUN
340! Martin Control
341
342      IF (it_RUN.LE.1)                                              THEN
343
344
345
346! Initialization of Mod_SISVAT_grd
347! --------------------------------
348
349              jt__SV =              max(1,int(dt0_SV /dt0DYn))
350              dt__SV =                   real(jt__SV)*dt0DYn
351          IF (dt__SV .NE. dt0_SV) write(6,61) dt__SV ,dt0_SV
352 61         format    ('dt__SV =',f9.3,' differs from dt0_SV =',f9.3)
353
354
355! Initialization of Mod_PHY_AT_grd
356! --------------------------------
357
358              jt__AT =              max(1,int(dt0_AT /dt0DYn))
359              dt__AT =                   real(jt__AT)*dt0DYn
360          IF (dt__AT .NE. dt0_AT) write(6,62) dt__AT ,dt0_AT
361 62         format    ('dt__AT =',f9.3,' differs from dt0_AT =',f9.3)
362
363
364! Initialization of Mod_PHY_CM_grd
365! --------------------------------
366
367              jt__CM =              max(1,int(dt0_CM /dt0DYn))
368              dt__CM =                   real(jt__CM)*dt0DYn
369          IF (dt__CM .NE. dt0_CM) write(6,63) dt__CM ,dt0_CM
370 63         format    ('dt__CM =',f9.3,' differs from dt0_CM =',f9.3)
371
372
373! Initialization of Mod_PHY_CP_grd
374! --------------------------------
375
376              jt__CP =              max(1,int(dt0_CP /dt0DYn))
377              dt__CP =                   real(jt__CP)*dt0DYn
378          IF (dt__CP .NE. dt0_CP) write(6,64) dt__CP ,dt0_CP
379 64         format    ('dt__CP =',f9.3,' differs from dt0_CP =',f9.3)
380
381
382! Initialization of Mod_PHY_RT_grd
383! --------------------------------
384
385              jt__RT =              max(1,int(dt0_RT /dt0DYn))
386              dt__RT =                   real(jt__RT)*dt0DYn
387          IF (dt__RT .NE. dt0_RT) write(6,65) dt__RT ,dt0_RT
388 65         format    ('dt__RT =',f9.3,' differs from dt0_RT =',f9.3)
389
390
391! Initialization
392! --------------
393
394! Initialization of 1-D Axes Variables: Vertical   Axis (Atmosphere)
395! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
396              pt__DY     = ptop_HOST
397          DO k=1,mzpp
398              sigma(k)   = s_HOST(k)
399          ENDDO
400
401! Martin control
402!PRINT*,'s_HOST=',s_HOST
403! Martin control
404
405          DO k=1,mzp
406                k1m(k)   = max(k-1, 1)
407                k1p(k)   = min(k+1,mzp)
408                k2m(k)   = max(k-2, 1)
409             dsigma(k)   =  sigma(k+1) - sigma(k)
410              sigmi(k+1) = (sigma(k+1) + sigma(k)) * 0.5
411          END DO
412              sigmi(1)   = 0.0
413              sigmi(mzpp)= 1.0
414
415          DO k=1,mzp
416             dsigmi(k)   =  sigmi(k+1) - sigmi(k)
417
418! Guess of sigma-levels Height
419! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~
420             hsigma(k)   =-(288./.0065)*(sigma(k) ** (287.*.0065/9.81)-1.)
421          END DO
422
423             write(6,501) (sigma(k),k=1,mzpp)
424  501          format(/,'  sigma: ',10f8.5                             &
425     &             ,5(/,'         ',10f8.5))
426
427             write(6,502)(hsigma(k),k=1,mzp)
428  502          format(/,' hsigma: ',10f8.1                             &
429     &             ,5(/,'         ',10f8.1))
430
431
432
433! Initialization of Mod_PHY____dat & Mod_PHY____kkl
434! -------------------------------------------------
435
436! Topography
437! ~~~~~~~~~~
438             sh_MAX      = 0.
439        DO ikl = 1,kcolp
440             sh__AP(ikl) =            sh___HOST(ikl)
441             sha_AP(ikl) =            sh_a_HOST(ikl)
442
443! Surface Slope
444! ~~~~~~~~~~~~~
445             sloxAP(ikl) =            slopxHOST(ikl)
446             sloyAP(ikl) =            slopyHOST(ikl)
447             slopAP(ikl) =            slopeHOST(ikl)
448             sh_MAX      = max(sh_MAX,sh__AP   (ikl))
449        END DO
450             dzaMIN      = hsigma(mzp) *(hsigma(1) -sh_MAX) /hsigma(1)
451
452
453        DO ikl = 1,kcolp
454
455! Geographic Coordinates
456! ~~~~~~~~~~~~~~~~~~~~~~
457
458! Martin rearrangement pour que RADACA ne plante plus:
459
460!           lon__r(ikl)   =     lonh_HOST(ikl) * 2.0 * piNmbr / 24.0
461           IF ((lonh_HOST(ikl) ) .LT. 0) THEN
462             lon__r(ikl) = 360 + (lonh_HOST(ikl) * 2.0 * piNmbr / 24.0)
463            ELSE
464             lon__r(ikl) = (lonh_HOST(ikl) * 2.0 * piNmbr / 24.0)
465           ENDIF
466           lon__h(ikl)   =     lonh_HOST(ikl)
467           lat__r(ikl)   =     latr_HOST(ikl)
468           sinLat(ikl)   = sin(lat__r(ikl))
469           cosLat(ikl)   = cos(lat__r(ikl))
470
471        ENDDO
472
473
474
475! -----------------------------------------------------------------------------!
476! Initialization of  Atm_DY (Counterpart of dynamical variables in PHY_MAR)
477! -------------------------
478
479! Martin control
480!PRINT*,'Avant PHY_Atm_DY_INI'
481!PRINT*,'size(psa_DY)=',size(psa_DY)
482! Martin control
483
484!                            **************
485                       CALL  PHY_Atm_DY_INI
486!                            **************
487
488! Initialization of  Atm_DY: needs plausible Atmospheric Conditions
489! ~~~~~~~~~~~~~~~~~~~~~~~~~ (here  idealized for Temperature Vertical Gradient)
490! Martin control
491!PRINT*,'Apres PHY_Atm_DY_INI'
492!PRINT*,'size(psa_DY)=',size(psa_DY)
493!PRINT*,'Dans PHY_MAR, calcul de Ta__DY:'
494!PRINT*,'mzpp=',mzpp
495!PRINT*,'pt__DY=',pt__DY
496!PRINT*,'Dtdz=',dTdz
497!PRINT*,'RCp=',RCp
498!PRINT*,'minval(pkta(:,mzpp))=',minval(pkta_HOST(:,mzpp))
499!PRINT*,'minval(psa__HOST(:))=',minval(psa__HOST(:))
500!PRINT*,'minval(Z___DY(:,:))=',minval(Z___DY(:,:))
501!PRINT*,'minval(Ta__DY(:,mzpp))=',minval(Ta__DY(:,mzpp))
502! Martin control
503
504        DO ikl = 1,kcolp
505           i                    = ii__AP   (ikl)
506           j                    = jj__AP   (ikl)
507           psa_DY    (ikl     ) = psa__HOST(ikl)
508           Ta__DY    (ikl,mzpp) = pkta_HOST(ikl,mzpp)*(psa_DY(ikl)+pt__DY)**RCp
509!gilles : init
510           Z___DY    (ikl,mzpp) = 0.
511           qv__DY    (ikl,mzpp) = 0.001
512         DO k =   mzp,1,-1
513           Z___DY    (ikl,k   ) =                                      &
514     &                 Ta__DY(ikl,mzpp)                                &
515     &      * (1.0 - ((psa_DY(ikl)*sigma (k)+pt__DY)                   &
516     &               /(psa_DY(ikl)          +pt__DY))                  &
517     &              **(R_DAir     * dTdz    /Grav_F))       / dTdz
518           Ta__DY    (ikl,k) =                                         &
519     &                 Ta__DY(ikl,mzpp)     -Z___DY(ikl,k)  * dTdz
520
521           qv__DY    (ikl,k) =   0.001
522
523           WindDY    (ikl,k) =   Wind_HOST(ikl,k)
524           ua__DY    (ikl,k) =   Ua___HOST(ikl,k)
525           va__DY    (ikl,k) =   va___HOST(ikl,k)
526           wa__DY    (ikl,k) =   Wa___HOST(ikl,k)
527         ENDDO
528        ENDDO
529
530! Martin Control
531!PRINT*,'Avant PHY_Atm_S0_INI'
532!PRINT*,'minval(Ta__DY(:,:))=',minval(Ta__DY(:,:))
533
534
535! -----------------------------------------------------------------------------!
536! Initialization of  Atm_S0 (Insolation and cos of Sun Zenithal Distance)
537! -------------------------
538
539!                            **************
540                       CALL  PHY_Atm_S0_INI
541!                            **************
542
543! Initialization of Mod_PHY_S0_ctr and Mod_PHY_S0_kkl
544! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
545              FaceS0 = FlagS0_SLO
546              MMskS0 = FlagS0_MtM
547          IF (FlagS0_SLO .AND. FlagS0_MtM)                          THEN
548              DO k   = 1,m_azim
549              DO ikl = 1,kcolp
550              cszkS0(ikl,k) = MMaskHOST(ikl,k)
551              ENDDO
552              ENDDO
553          ENDIF
554
555
556! -----------------------------------------------------------------------------!
557! Initialization of  SISVAT
558! -------------------------
559
560! Initialization of Mod_SISVAT_ctr  (SISVAT Switches
561! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~        & Time/Space control Variables)
562        VegMod = FlagSV_Veg
563        SnoMod = FlagSV_SNo
564        BloMod = FlagSV_BSn
565        InpSWD = FlagSV_SWD
566        InpSBC = FlagSV_SBC
567        SVaKzT = FlagSV_KzT
568        SVaUBC = FlagSV_UBC
569
570
571! Initialization of Mod_SISVAT_kkl  (from INPUT from HOST Model)
572! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
573        DO ikp = 1,kcolp
574           i   = ii__AP(ikp)
575           j   = jj__AP(ikp)
576        DO mn  = 1,mwp
577        DO k   = 1,mzp
578           kk  =   mzpp - k
579           WindSV(ikp,mn,k) = Wind_HOST(ikp,k)
580           pkt0SV(ikp,mn,kk)= pkta_HOST(ikp,k)
581        END DO
582           Ua__SV(ikp,mn)   = Ua___HOST(ikp,mzp)
583           Va__SV(ikp,mn)   = Va___HOST(ikp,mzp)
584        END DO
585        END DO
586
587
588! Initialization of SISVAT Variables
589! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
590!                            **************
591                       CALL  PHY_SISVAT_INI
592!                            **************
593
594
595! Initialization of Mod_SISVAT_kkl  (SISVAT OUTPUT Grid Points)
596! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
597           iwr = 0
598        DO ipt = 1,NbPts
599           IF (IOi0SV(ipt).EQ.0)                                    THEN
600               IOi_SV(ipt)=i_x0
601           ELSE
602               IOi_SV(ipt)=IOi0SV(ipt)
603           END IF
604           IF (IOj0SV(ipt).EQ.0)                                    THEN
605               IOj_SV(ipt)=j_y0
606           ELSE
607               IOj_SV(ipt)=IOj0SV(ipt)
608           END IF
609           write(6,*) 'ipt   , IOi_SV, IOj_SV              = '               &
610     &                ,ipt   , IOi_SV(ipt), IOj_SV(ipt)
611        DO n   = 1,mwp
612           iwr = 1+iwr
613          IF (iwr.LE.nbwri)                                         THEN
614           no__SV(iwr) = 0
615           i___SV(iwr) = IOi_SV(ipt)
616           j___SV(iwr) = IOj_SV(ipt)
617           n___SV(iwr) = n
618           write(6,*) 'n     , i___SV, j___SV, n___SV, iwr = '               &
619     &                ,n     , i___SV(iwr), j___SV(iwr), n___SV(iwr), iwr
620          END IF
621        END DO
622        END DO
623
624
625
626! -----------------------------------------------------------------------------!
627! Initialization of  Atm_RT (Radiative Transfert through the Atmosphere)
628! -------------------------
629
630!gilles: PHY_Atm_RT_INI requires correct dates
631          YearTU = Year_H
632          Mon_TU = Mon__H
633          Day_TU = Day__H
634          HourTU = Hour_H
635          minuTU = minu_H
636          sec_TU = sec__H
637
638!                            **************
639                       CALL  PHY_Atm_RT_INI
640!                            **************
641
642! -----------------------------------------------------------------------------!
643! Initialization of  Atm_AT (Turbulent Transfert through the Atmosphere)
644! -------------------------
645
646!                            **************
647                       CALL  PHY_Atm_AT_INI(FlagAT_TKE,TypeAT)
648!                            **************
649
650
651
652! -----------------------------------------------------------------------------!
653! Initialization of  Atm_CP (Convectiv Transfert through the Atmosphere)
654! -------------------------
655
656!                            **************
657                       CALL  PHY_Atm_CP_INI(mzp,kcolp)
658!                            **************
659
660
661
662! -----------------------------------------------------------------------------!
663! Initialization of  CMiPhy (Cloud Microphysical Scheme)
664! -------------------------
665
666! Initialization of Mod_PHY_CM_ctr  (CMiPhy Switches
667! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~        & Time/Space control Variables)
668        CM_UpD = FlagCM_UpD
669
670!                            **************
671                       CALL  PHY_Atm_CM_INI
672!                            **************
673
674      END IF
675
676
677! Martin CONTROL
678!PRINT*, 'Impressions control'
679!call iophys_ecrit('TA__DY_surf',1,'surface temperature','K',Ta__DY(:,61))
680!call iophys_ecrit('TA__DY_air',60,'air temperature','K',Ta__DY(:,1:60))
681
682!PRINT*,'TA__DY(:,61)=',TA__DY(kcolq/2,61)
683
684! Interface: From HOST Model Variables to MARp Physics Variables
685! ==============================================================
686
687
688! Time
689! ----
690
691          YearTU = Year_H
692          Mon_TU = Mon__H
693          Day_TU = Day__H
694          HourTU = Hour_H
695          minuTU = minu_H
696          sec_TU = sec__H
697
698          TimeTU = (float(351)+(float(YearTU) -float(1902)) *float(365)&! Nb Days before YearTU
699     &                        +(float(YearTU) -float(1901)) /float(  4)&! Nb Leap Years
700     &                 + float(njYear(Mon_TU))                         &! Nb Days before Mon_TU
701     &                 + float(njLeap(Mon_TU))                         &! (including Leap Day)
702     &       *max(zer0,un_1-mod(float(YearTU),float(4)))               &!
703     &                        + float(Day_TU)-float(1)  )   *float( 24)&!
704     &             +float(HourTU)                                      &!
705     &           + (float(minuTU) *float(60) +float(sec_TU))/3600.      !
706
707
708
709! -----------------------------------------------------------------------------!
710! Assignation    of  Mod_PHY_AT_grd
711! ---------------------------------
712          IF (it_RUN .EQ. 1) THEN
713            TimeAT = TimeTU-dt0_AT/3600. ! Initialisation à la première itération
714          END IF
715          IF (FlagAT  .AND.      mod(it_RUN-1,jt__AT).EQ.0)         THEN
716!gilles: leapfrog scheme stops time integration regularly
717!              dTimAT =              (TimeTU-TimeAT) * 3600.
718              dTimAT =    dt0_AT
719              TimeAT =               TimeTU
720          END IF
721
722! Martin CONTROL
723PRINT*,'CONTROL PHY_MAR temps'
724PRINT*,'  Year_H =', Year_H
725PRINT*,'  Mon__H =', Mon__H
726PRINT*,'  Day__H =', Day__H
727PRINT*,'  Hour_H =', Hour_H
728PRINT*,'  minu_H =', minu_H
729PRINT*,'  sec__H =', sec__H
730PRINT*,'jt__AT=',jt__AT
731PRINT*,'TimeTU=',TimeTU
732PRINT*,'TimeAT=',TimeAT
733PRINT*,'dTimAT=',dTimAT
734! Martin CONTROL
735
736
737! -----------------------------------------------------------------------------!
738! Assignation    of  Mod_PHY_DY_kkl
739! ---------------------------------
740
741          DO ikl=1,kcolp
742             i  =  ii__AP(ikl)
743             j  =  jj__AP(ikl)
744             k  =  mzpp
745            ExnrDY    (ikl,k) = exp(RCp *log(psa__HOST(ikl)*sigma(k)+pt__DY))
746            pkt_DY    (ikl,k) =              pkta_HOST(ikl,k)
747            Ta__DY    (ikl,k) =              pkta_HOST(ikl,k) * ExnrDY(ikl,k)
748!  CAUTION: Tas_SV_xy is not allowed to be changed by data coming from outside this routine
749            Z___DY    (ikl,k) =              gZa__HOST(ikl,k)         *Grav_I
750            ZmidDY    (ikl,k) =              gZam_HOST(ikl,k)         *Grav_I
751            qv__DY    (ikl,k) =              qv___HOST(ikl,k)
752
753          DO  k               =   1,mzp
754            ExnrDY    (ikl,k) = exp(RCp *log(psa__HOST(ikl)*sigma(k  )+pt__DY))
755            pkt_DY    (ikl,k) =              pkta_HOST(ikl,k)
756            Ta__DY    (ikl,k) =              pkta_HOST(ikl,k) * ExnrDY(ikl,k)
757            roa_DY    (ikl,k) =             (psa__HOST(ikl)*sigma(k  )+pt__DY)     &
758     &                                     /(Ta__DY   (ikl,k)         *R_DAir)
759            roamDY    (ikl,k) =             (psa__HOST(ikl)*sigmi(k+1)+pt__DY)     &
760     &                                     /(Ta__DY   (ikl,k)         *R_DAir)
761            Z___DY    (ikl,k) =              gZa__HOST(ikl,k)         *Grav_I
762            ZmidDY    (ikl,k) =              gZam_HOST(ikl,k)         *Grav_I
763            qv__DY    (ikl,k) =              qv___HOST(ikl,k)
764
765
766
767! -----------------------------------------------------------------------------!
768! Assignation    of  Mod_PHY_AT_kkl
769! ---------------------------------
770
771
772          IF (FlagAT.AND.it_EXP.GT.1.AND.mod(it_RUN-1,jt__AT).EQ.0) THEN
773            TKE_AT    (ikl,k) =              TKE__HOST(ikl,k)
774            eps_AT    (ikl,k) =              eps__HOST(ikl,k)
775            TrT_AT    (ikl,k) =            (TKE_AT    (ikl,k) - TrT_AT    (ikl,k)) &
776     &                                    / dTimAT
777          END IF
778
779
780
781! -----------------------------------------------------------------------------!
782! Assignation    of  Mod_PHY_CM_kkl
783! ---------------------------------
784
785          IF (FlagCM.AND.it_EXP.GT.1.AND.mod(it_RUN-1,jt__CM).EQ.0) THEN
786
787!           IF (qw___HOST(ikl,k).LT.qh_MIN)                         THEN
788!               qv__DY(ikl,k) =             qv__DY(ikl,k) + qw__CM(ikl,k)
789!               qw__CM(ikl,k) =             0.
790! #cw           CCNwCM(ikl,k) =             0.
791!           ELSE
792                qw__CM(ikl,k) =             qw___HOST(ikl,k)
793! #cw           CCNwCM(ikl,k) =             CCN__HOST(ikl,k)
794!           END IF
795
796!           IF (qi___HOST(ikl,k).LT.qh_MIN)                         THEN
797!               qv__DY(ikl,k) =             qv__DY(ikl,k) + qi__CM(ikl,k)
798!               qi__CM(ikl,k) =             0.
799!               CCNiCM(ikl,k) =             0.
800!           ELSE
801                qi__CM(ikl,k) =             qi___HOST(ikl,k)
802                CCNiCM(ikl,k) =             CIN__HOST(ikl,k)
803!           END IF
804
805! Gilles: CF___HOST non sauve & CFraCM reinitialise ici
806!               CFraCM(ikl,k) =             CF___HOST(ikl,k)
807
808!           IF (qw__CM(ikl,k).LT.qh_MIN  .AND.                         &
809!    &          qi__CM(ikl,k).LT.qh_MIN)                            THEN
810!               CFraCM(ikl,k) =             0.
811!           ELSE
812!               CFraCM(ikl,k) =  max(CFrMIN,CF___HOST(ikl,k))
813!           END IF
814
815                qs__CM(ikl,k) =             qs___HOST(ikl,k)
816                qr__CM(ikl,k) =             qr___HOST(ikl,k)
817          END IF
818
819          ENDDO
820
821
822
823! -----------------------------------------------------------------------------!
824! Assignation    of  Mod_SISVAT_gpt                  (A-O FORCING OR COUPLING)
825! ---------------------------------                 
826
827            sst_SB    (ikl)   =              sst__HOST(ikl) 
828! #IP       sif_SB    (ikl)   =              sif__HOST(ikl) 
829! #AO     DO k = 1,mwp
830! #AO       s_T_AO_xyn(i,j,k) =              s_T__HOST(ikl,k)
831! #AO       Alb_AO_xyn(i,j,k) =              Alb__HOST(ikl,k)
832! #AO     ENDDO
833
834          ENDDO
835
836! -----------------------------------------------------------------------------!
837! Assignation    of  Mod_PHY_DY_kkl
838! ---------------------------------
839          DO ikl = 1,kcolp
840             i   = ii__AP(ikl)
841             j   = jj__AP(ikl)
842             psa_DY(ikl)   = psa__HOST(ikl)
843          DO k   = 1,mzp
844             WindDY(ikl,k) = Wind_HOST(ikl,k)
845             ua__DY(ikl,k) = Ua___HOST(ikl,k)
846             va__DY(ikl,k) = va___HOST(ikl,k)
847             wa__DY(ikl,k) = Wa___HOST(ikl,k)  +  sqrt(2.*max(eps6,TKE_AT(ikl,k))/3.)
848          END DO
849          END DO
850
851!                            **************
852                       CALL  PHY_Atm_DY_RUN                                    ! Assignation of MAR Dyn. Variables
853!                            **************
854
855           IF(FlagVR)                                               THEN
856!                            **************
857                       CALL  PHY________OUT('After  PHY_Atm_DY_RUN                             ')
858!                            **************  12345678901234567890123456789012345678901234567890
859!                                                     1         2         3         4         5
860           END IF
861
862
863
864! -----------------------------------------------------------------------------!
865! Saturation Specific Humidity
866! ----------------------------
867
868          IF (FlagCM .OR. &! ***************
869     &        FlagSV)  CALL  PHY_Atm_CM_QSat
870!                            ***************
871
872
873
874! -----------------------------------------------------------------------------!
875! Execution      of  Atm_S0 (Insolation and cos of Sun Zenithal Distance)
876! -------------------------
877
878!                            **************
879                       CALL  PHY_Atm_S0_RUN
880!                            **************
881
882           IF(FlagVR)                                               THEN
883!                            **************
884                       CALL  PHY________OUT('After  PHY_Atm_S0_RUN                             ')
885!                            **************  12345678901234567890123456789012345678901234567890
886!                                                     1         2         3         4         5
887           END IF
888
889
890
891! -----------------------------------------------------------------------------!
892! Execution      of  Atm_RT (Radiative Transfer through the Atmosphere)
893! -------------------------
894
895          IF (FlagRT  .AND.      mod(it_RUN-1,jt__RT).EQ.0)         THEN
896
897!                            **************
898                       CALL  PHY_Atm_RT_RUN(kcolp,ikl0)
899!                            **************
900
901           IF(FlagVR)                                               THEN
902!                            **************
903                       CALL  PHY________OUT('After  PHY_Atm_RT_RUN                             ')
904!                            **************  12345678901234567890123456789012345678901234567890
905!                                                     1         2         3         4         5
906           END IF
907          END IF
908
909
910
911! -----------------------------------------------------------------------------!
912! Execution      of  SISVAT (Soil-Ice-Snow-Vegetation-Atmosphere-Transfer Scheme)
913! -------------------------
914
915          IF (FlagSV)                                               THEN
916
917! Assignation    of Mod_SISVAT_kkl  (from INPUT from HOST Model)
918! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
919            DO ikp = 1,kcolp
920            DO mn  = 1,mwp
921               pkPaSV(ikp,mn)   =              psa__HOST(ikp)           +pt__DY
922               ExnrSV(ikp,mn)   = exp(RCp *log(psa__HOST(ikp)*sigma(mzp)+pt__DY))
923            DO k   = 1,mzp
924               kk  =   mzpp - k
925               WindSV(ikp,mn,k) =              Wind_HOST(ikp,k)
926               pkt0SV(ikp,mn,kk)=              pkta_HOST(ikp,k)
927               qv__SV(ikp,mn,kk)=              qv___HOST(ikp,k)
928               zza_SV(ikp,mn,kk)=             (gZa__HOST(ikp,k)                &  !
929      &                                       -gZa__HOST(ikp,mzpp))     *Grav_I   !
930               roa_SV(ikp,mn,kk)=              roa_DY   (ikp,k)         * 1.e3    ! [kg/m3]
931            END DO
932               Ua__SV(ikp,mn)   =              Ua___HOST(ikp,mzp)
933               Va__SV(ikp,mn)   =              Va___HOST(ikp,mzp)
934               TaT_SV(ikp,mn)   =              pkta_HOST(ikp,mzp) * ExnrSV(ikp,mn)
935            END DO
936            END DO
937
938
939!                            **************
940                       CALL  PHY_SISVAT_RUN(                           &
941!                            **************
942! -----------------------------------------------------------------------------!
943! #TC&                                        qxTC,uqTC  ,qsTC  ,      &       ! Aerosols: Atm.Conc., Surf.Flux
944! -----------------------------------------------------------------------------!
945     &                                     )
946
947           IF(FlagVR)                                               THEN
948!                            **************
949                       CALL  PHY________OUT('After  PHY_SISVAT_RUN                             ')
950!                            **************  12345678901234567890123456789012345678901234567890
951!                                                     1         2         3         4         5
952           END IF
953          END IF
954
955
956
957! -----------------------------------------------------------------------------!
958! Execution      of  Atm_AT (Atmospheric Turbulence Contribution)
959! -------------------------
960
961          IF (FlagAT  .AND.      mod(it_RUN-1,jt__AT).EQ.0)         THEN
962
963!                            **************
964                       CALL  PHY_Atm_AT_RUN(FlagSV_KzT,FlagCM)
965!                            **************
966
967           IF(FlagVR)                                               THEN
968!                            **************
969                       CALL  PHY________OUT('After  PHY_Atm_AT_RUN                             ')
970!                            **************  12345678901234567890123456789012345678901234567890
971!                                                     1         2         3         4         5
972           END IF
973          END IF
974
975
976
977! -----------------------------------------------------------------------------!
978! Execution      of  Atm_CP (Convection Parameterisation - Mass Flux Scheme)
979! -------------------------
980
981          IF (FlagCP  .AND.      mod(it_RUN-1,jt__CP).EQ.0)         THEN
982
983!                            **************
984                       CALL  PHY_Atm_CP_RUN(mzp,kcolp)
985                       ! CALL  PHY_Atm_CP_RUN(mxp,kcolp) ! ancien bug
986!                            **************
987
988           IF(FlagVR)                                               THEN
989!                            **************
990                       CALL  PHY________OUT('After  PHY_Atm_CP_RUN                             ')
991!                            **************  12345678901234567890123456789012345678901234567890
992!                                                     1         2         3         4         5
993           END IF
994          END IF
995
996
997
998! -----------------------------------------------------------------------------!
999! Execution      of  Atm_CM (Cloud Microphysics)
1000! -------------------------
1001
1002          IF (FlagCM  .AND.      mod(it_RUN-1,jt__CM).EQ.0)         THEN
1003
1004!                            **************
1005                       CALL  PHY_Atm_CM_RUN
1006!                            **************
1007
1008
1009           IF(FlagVR)                                               THEN
1010!                            **************
1011                       CALL  PHY________OUT('After  PHY_Atm_CM_RUN                             ')
1012!                            **************  12345678901234567890123456789012345678901234567890
1013!                                                     1         2         3         4         5
1014           END IF
1015          END IF
1016
1017
1018
1019
1020! Assignation of the Tendencies to transfer outside the physical Parameterizations Package
1021! ========================================================================================
1022
1023          DO ikl=1,kcolp
1024             i   = ii__AP(ikl)
1025             j   = jj__AP(ikl)
1026
1027
1028
1029! -----------------------------------------------------------------------------!
1030! Reinitialization of the Tendencies
1031! ----------------------------------
1032
1033          DO k=   1,mzp
1034
1035              dpkt___dt(ikl,k) = 0.0                                    !
1036              dua____dt(ikl,k) = 0.0                                    !
1037              dva____dt(ikl,k) = 0.0                                    !
1038              dqv____dt(ikl,k) = 0.0                                    !
1039              dqw____dt(ikl,k) = 0.0                                    !
1040              dCF____dt(ikl,k) = 0.0                                    !
1041              dqi____dt(ikl,k) = 0.0                                    !
1042              dCi____dt(ikl,k) = 0.0                                    !
1043              dqs____dt(ikl,k) = 0.0                                    !
1044              dqr____dt(ikl,k) = 0.0                                    !
1045
1046
1047
1048! -----------------------------------------------------------------------------!
1049! Tendencies from SISVAT
1050! ----------------------
1051
1052          IF (FlagSV  .AND.      mod(it_RUN-1,jt__SV).EQ.0)         THEN!
1053              dpktSV_dt(ikl,k) =                    dpktSV_gpt(ikl,k)   ! Reduced Potential Temperature TENDENCY, SISVAT        [KX/s]
1054          END IF
1055
1056! Update           of the tendencies
1057! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1058          IF (FlagSV)                                               THEN!
1059              dpkt___dt(ikl,k) = dpkt___dt(ikl,k) + dpktSV_gpt(ikl,k)   !
1060          END IF
1061
1062
1063
1064! -----------------------------------------------------------------------------!
1065! Tendencies from Atm_AT
1066! ----------------------
1067
1068          IF (FlagAT  .AND.      mod(it_RUN-1,jt__AT).EQ.0)         THEN!
1069              dua____dt(ikl,k) =                    dua_AT(ikl,k)       ! Wind Speed       (x-direc.)   TENDENCY, ALL Contribut.[m/s2]
1070              dva____dt(ikl,k) =                    dva_AT(ikl,k)       ! Wind Speed       (y-direc.)   TENDENCY, ALL Contribut.[m/s2]
1071              dpktAT_dt(ikl,k) =                    dpktAT(ikl,k)       ! Reduced Potential Temperature TENDENCY, Atm_AT        [KX/s]
1072              dqv_AT_dt(ikl,k) =                    dqv_AT(ikl,k)       ! Specific          Humidity    TENDENCY, Atm_AT     [kg/kg/s]
1073              dqw_AT_dt(ikl,k) =                    dqw_AT(ikl,k)       ! Cloud Droplets Concentration  TENDENCY, Atm_AT     [kg/kg/s]
1074              dqi_AT_dt(ikl,k) =                    dqi_AT(ikl,k)       ! Cloud Crystals Concentration  TENDENCY, Atm_AT     [kg/kg/s]
1075              dqs_AT_dt(ikl,k) =                    dqs_AT(ikl,k)       ! Snow Particles Concentration  TENDENCY, Atm_AT     [kg/kg/s]
1076              dqr_AT_dt(ikl,k) =                    dqr_AT(ikl,k)       ! Rain Drops     Concentration  TENDENCY, Atm_AT     [kg/kg/s]
1077! #cw         dCw_AT_dt(ikl,k) =                    dCW_AT(ikl,k)       ! CCN            Concentration  TENDENCY, Atm_AT         [1/s]
1078              dCi_AT_dt(ikl,k) =                    dCi_AT(ikl,k)       ! CIN            Concentration  TENDENCY, Atm_AT         [1/s]
1079 
1080! NO Tendencies on TKE and its dissipation ==> TKE__HOST and eps__HOST are updated
1081! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1082              TKE__HOST(ikl,k) =                    TKE_AT(ikl,k)       !
1083              eps__HOST(ikl,k) =                    eps_AT(ikl,k)       !
1084             
1085! Re-Initialization of the TKE Transport Rate
1086! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1087              TrT_AT   (ikl,k) =                    TKE_AT(ikl,k)       !
1088          END IF
1089
1090! Update           of the tendencies
1091! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1092          IF (FlagAT)                                               THEN!
1093              dpkt___dt(ikl,k) = dpkt___dt(ikl,k) + dpktAT(ikl,k)       !
1094              dqv____dt(ikl,k) = dqv____dt(ikl,k) + dqv_AT(ikl,k)       !
1095              dqw____dt(ikl,k) = dqw____dt(ikl,k) + dqw_AT(ikl,k)       !
1096              dqi____dt(ikl,k) = dqi____dt(ikl,k) + dqi_AT(ikl,k)       !
1097              dqs____dt(ikl,k) = dqs____dt(ikl,k) + dqs_AT(ikl,k)       !
1098              dqr____dt(ikl,k) = dqr____dt(ikl,k) + dqr_AT(ikl,k)       !
1099! #cw         dCw____dt(ikl,k) = dCw____dt(ikl,k) + dCw_AT(ikl,k)       !
1100              dCi____dt(ikl,k) = dCi____dt(ikl,k) + dCi_AT(ikl,k)       !
1101          END IF
1102
1103
1104
1105! -----------------------------------------------------------------------------!
1106! Tendencies from CMiPhy (dqw_CM, dCw_CM, dqi_CM, dCi_CM
1107! ----------------------  dqs_CM, dqr_CM  are consumed in CMiPhy then reset to 0
1108!                                                      if FlagCM_UpD  = .TRUE. )
1109
1110          IF (FlagCM  .AND.      mod(it_RUN-1,jt__CM).EQ.0)         THEN!
1111              dpktCM_dt(ikl,k) =                    dpktCM(ikl,k)       ! Reduced Potential Temperature TENDENCY, CMiPhy        [KX/s]
1112              dqv_CM_dt(ikl,k) =                    dqv_CM(ikl,k)       ! Specific          Humidity    TENDENCY, CMiPhy     [kg/kg/s]
1113              dqw_CM_dt(ikl,k) =                    dqw_CM(ikl,k)       ! Cloud Droplets Concentration  TENDENCY, CMiPhy     [kg/kg/s]
1114              dCF_CM_dt(ikl,k) =                    dCF_CM(ikl,k)       ! Cloud Fraction                TENDENCY, CMiPhy     [kg/kg/s]
1115              dqi_CM_dt(ikl,k) =                    dqi_CM(ikl,k)       ! Cloud Crystals Concentration  TENDENCY, CMiPhy     [kg/kg/s]
1116              dqs_CM_dt(ikl,k) =                    dqs_CM(ikl,k)       ! Snow Particles Concentration  TENDENCY, CMiPhy     [kg/kg/s]
1117              dqr_CM_dt(ikl,k) =                    dqr_CM(ikl,k)       ! Rain Drops     Concentration  TENDENCY, CMiPhy     [kg/kg/s]
1118! #cw         dCw_CM_dt(ikl,k) =                    dCW_CM(ikl,k)       ! CCN            Concentration  TENDENCY, CMiPhy         [1/s]
1119              dCi_CM_dt(ikl,k) =                    dCi_CM(ikl,k)       ! CIN            Concentration  TENDENCY, CMiPhy         [1/s]
1120          END IF
1121
1122! Update           of the tendencies
1123! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1124          IF (FlagCM)                                               THEN!
1125              dpkt___dt(ikl,k) = dpkt___dt(ikl,k) + dpktCM(ikl,k)       !
1126              dqv____dt(ikl,k) = dqv____dt(ikl,k) + dqv_CM(ikl,k)       !
1127              dqw____dt(ikl,k) = dqw____dt(ikl,k) + dqw_CM(ikl,k)       !
1128              dCF____dt(ikl,k) = dCF____dt(ikl,k) + dCF_CM(ikl,k)       !
1129              dqi____dt(ikl,k) = dqi____dt(ikl,k) + dqi_CM(ikl,k)       !
1130              dqs____dt(ikl,k) = dqs____dt(ikl,k) + dqs_CM(ikl,k)       !
1131              dqr____dt(ikl,k) = dqr____dt(ikl,k) + dqr_CM(ikl,k)       !
1132! #cw         dCw____dt(ikl,k) = dCw____dt(ikl,k) + dCw_CM(ikl,k)       !
1133              dCi____dt(ikl,k) = dCi____dt(ikl,k) + dCi_CM(ikl,k)       !
1134          END IF
1135
1136
1137
1138! -----------------------------------------------------------------------------!
1139! Tendencies from CVAmnh
1140! ----------------------
1141
1142          IF (FlagCP  .AND.      mod(it_RUN-1,jt__CP).EQ.0)         THEN!
1143              dpktCP_dt(ikl,k) =                    dpktCP(ikl,k)       ! Reduced Potential Temperature TENDENCY, CVAmnh        [KX/s]
1144              dqv_CP_dt(ikl,k) =                    dqv_CP(ikl,k)       ! Specific          Humidity    TENDENCY, CVAmnh     [kg/kg/s]
1145              dqw_CP_dt(ikl,k) =                    dqw_CP(ikl,k)       ! Cloud Droplets Concentration  TENDENCY, CVAmnh     [kg/kg/s]
1146              dqi_CP_dt(ikl,k) =                    dqi_CP(ikl,k)       ! Cloud Crystals Concentration  TENDENCY, CVAmnh     [kg/kg/s]
1147          END IF
1148
1149! Update           of the tendencies
1150! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1151          IF (FlagCP)                                               THEN!
1152              dpkt___dt(ikl,k) = dpkt___dt(ikl,k) + dpktCP(ikl,k)       !
1153              dqv____dt(ikl,k) = dqv____dt(ikl,k) + dqv_CP(ikl,k)       !
1154              dqw____dt(ikl,k) = dqw____dt(ikl,k) + dqw_CP(ikl,k)       !
1155              dqi____dt(ikl,k) = dqi____dt(ikl,k) + dqi_CP(ikl,k)       !
1156          END IF
1157
1158
1159
1160! -----------------------------------------------------------------------------!
1161! Tendencies from radCEP
1162! ----------------------
1163
1164          IF (FlagRT  .AND.      mod(it_RUN-1,jt__RT).EQ.0)         THEN!
1165              dpktRT_dt(ikl,k) =                    dpktRT(ikl,k)       ! Reduced Potential Temperature TENDENCY, radCEP        [KX/s]
1166          END IF
1167
1168! Update           of the tendencies
1169! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1170          IF (FlagRT)                                               THEN!
1171              dpkt___dt(ikl,k) = dpkt___dt(ikl,k) + dpktRT(ikl,k)       !
1172          END IF
1173
1174          ENDDO
1175
1176
1177
1178
1179
1180! Assignation of the Variables  to keep      inside the Physical Parameterizations Package
1181! ========================================================================================
1182
1183!         ...  NOTHING up to now
1184
1185
1186
1187! Assignation of the Variables  to transfer outside the Physical Parameterizations Package
1188! ========================================================================================
1189
1190! -----------------------------------------------------------------------------!
1191! Update of  pkta_HOST, qv___HOST
1192! -------------------------------
1193
1194          DO k=1,mzpp
1195            pkta_HOST(ikl,k) = pkt_DY(ikl,k)                                   !  Always on k=mzpp, possible from 1 to mzp (with dpkt = 0.)
1196            qv___HOST(ikl,k) = qv__DY(ikl,k)                                   !
1197          ENDDO
1198
1199
1200
1201! -----------------------------------------------------------------------------!
1202! Update of  qw___HOST, CF___HOST, qi___HOST, CIN__HOST, qs___HOST, qr___HOST
1203! ---------------------------------------------------------------------------
1204
1205      IF (FlagCM)                                                   THEN
1206          DO k=1,mzp
1207            qw___HOST(ikl,k) = qw__CM(ikl,k)
1208! #cw       CCN__HOST(ikl,k) = CCNwCM(ikl,k)
1209            CF___HOST(ikl,k) = CFraCM(ikl,k)
1210            qi___HOST(ikl,k) = qi__CM(ikl,k)
1211            CIN__HOST(ikl,k) = CCNiCM(ikl,k)
1212            qs___HOST(ikl,k) = qs__CM(ikl,k)
1213            qr___HOST(ikl,k) = qr__CM(ikl,k)
1214          ENDDO
1215      END IF
1216
1217
1218
1219! -----------------------------------------------------------------------------!
1220! Update of  d(S,LH)/dT (needed in NEMO)
1221! ---------------------
1222
1223! #AO     DO k=1,mwpp
1224! #AO       dSdT2HOST(ikl,k) = dSdTAO_xyn(i,j,k)
1225! #AO       dLdT2HOST(ikl,k) = dLdTAO_xyn(i,j,k)
1226! #AO     ENDDO
1227
1228
1229
1230          ENDDO
1231
1232
1233
1234
1235! OUTPUT
1236! ======
1237
1238! OUTPUT of Tendencies
1239! --------------------
1240
1241          IF (FlagVR .OR.                                              &
1242     &       (FLAG_O .AND. ((minuTU.EQ.0 .AND. sec_TU.EQ.0) .OR.       &
1243     &                       it_RUN.EQ.1                        ))) THEN
1244             ikl = ikl0
1245             i   = ii__AP(ikl)
1246             j   = jj__AP(ikl)
1247
1248! pkt TENDENCIES
1249! ~~~~~~~~~~~~~~
1250                write(4,400)
1251 400            format(//,'   pkt TENDENCIES',/,'   **************'/,1x)
1252                write(4,403) Day_TU,Mon_TU,YearTU,HourTU,MinuTU,Sec_TU,it_EXP
1253 403            format(3x,2(i2,'-'),i4,4x,3(i2,'-'),'  Simulation Iteration No ',i6,/,1x)
1254                write(4,404)
1255                write(4,401)
1256 401            format('    |   SISVAT    |   Atm_AT    |   CMiPhy    |   CVAmnh    |   radCEP    |')
1257                write(4,402)
1258 402            format('    |      [K/d]  |      [K/d]  |      [K/d]  |      [K/d]  |      [K/d]  |')
1259                write(4,404)
1260 404            format(4('-'),'+',5(13('-'),'+'))
1261             DO k=1,mzp
1262                write(4,405)        k                                  &
1263     &                     , ExnrDY(ikl,k)*dpktSV_gpt(ikl,k)*86400.    &
1264     &                     , ExnrDY(ikl,k)*dpktAT    (ikl,k)*86400.    &
1265     &                     , ExnrDY(ikl,k)*dpktCM    (ikl,k)*86400.    &
1266     &                     , ExnrDY(ikl,k)*dpktCP    (ikl,k)*86400.    &
1267     &                     , ExnrDY(ikl,k)*dpktRT    (ikl,k)*86400.
1268 405            format(i3,' |',5(f12.6,' |'))
1269
1270                IF (mod(k,20).EQ.0)                                 THEN
1271                write(4,404)
1272                write(4,401)
1273                write(4,402)
1274                write(4,404)
1275                END IF
1276             ENDDO
1277                write(4,404)
1278
1279! pkt TENDENCIES
1280! ~~~~~~~~~~~~~~
1281                write(4,410)
1282 410            format(//,'   Qv  TENDENCIES',/,'   **************'/,1x)
1283                write(4,403) Day_TU,Mon_TU,YearTU,HourTU,MinuTU,Sec_TU,it_EXP
1284!403            format(3x,2(i2,'-'),i4,4x,3(i2,'-'),'  Simulation Iteration No ',i6,/,1x)
1285                write(4,404)
1286                write(4,411)
1287 411            format('    |   SISVAT    |   Atm_AT    |   CMiPhy    |   CVAmnh    |    TOTAL    |')
1288                write(4,412)
1289 412            format('    | [g/kg/min]  | [g/kg/min]  | [g/kg/min]  | [g/kg/min]  | [g/kg/min]  |')
1290                write(4,404)
1291!404            format(4('-'),'+',5(13('-'),'+'))
1292             DO k=1,mzp
1293                write(4,415)        k                                  &
1294     &                     ,               dqv_AT   (ikl,k)*60000.     &
1295     &                     ,               dqv_CM   (ikl,k)*60000.     &
1296     &                     ,               dqv_CP   (ikl,k)*60000.     &
1297     &                     ,               dqv____dt(ikl,k)*60000.
1298 415            format(i3,' |',12x,' |',4(f12.6,' |'))
1299
1300                IF (mod(k,20).EQ.0)                                 THEN
1301                write(4,404)
1302                write(4,411)
1303                write(4,412)
1304                write(4,404)
1305                END IF
1306             ENDDO
1307                write(4,404)
1308          END IF
1309
1310          IF (FLAG_O .AND. ((minuTU.EQ.0 .AND. sec_TU.EQ.0) .OR.       &
1311     &                       it_RUN.EQ.1                        ))  THEN
1312             
1313!                            **************
1314                       CALL  PHY________OUT('After  PHY_MAR                                    ')
1315!                            **************  12345678901234567890123456789012345678901234567890
1316!                                                     1         2         3         4         5
1317          END IF 
1318
1319
1320
1321      return
1322      end
Note: See TracBrowser for help on using the repository browser.