Context Navigation

PHY_MAR___.f90 @ 2100

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

Inclusion de la physique de MAR

Integration of MAR physics

File size: 73.2 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	DO ikl = 1,kcolp
504	i = ii__AP (ikl)
505	j = jj__AP (ikl)
506	psa_DY (ikl ) = psa__HOST(ikl)
507	Ta__DY (ikl,mzpp) = pkta_HOST(ikl,mzpp)(psa_DY(ikl)+pt__DY)*RCp
508	DO k = mzp,1,-1
509	Z___DY (ikl,k ) = &
510	& Ta__DY(ikl,mzpp) &
511	& * (1.0 - ((psa_DY(ikl)*sigma (k)+pt__DY) &
512	& /(psa_DY(ikl) +pt__DY)) &
513	& *(R_DAir dTdz /Grav_F)) / dTdz
514	Ta__DY (ikl,k) = &
515	& Ta__DY(ikl,mzpp) -Z___DY(ikl,k) * dTdz
516
517	qv__DY (ikl,k) = 0.001
518
519	WindDY (ikl,k) = Wind_HOST(ikl,k)
520	ua__DY (ikl,k) = Ua___HOST(ikl,k)
521	va__DY (ikl,k) = va___HOST(ikl,k)
522	wa__DY (ikl,k) = Wa___HOST(ikl,k)
523	ENDDO
524	ENDDO
525
526	! Martin Control
527	!PRINT*,'Avant PHY_Atm_S0_INI'
528	!PRINT*,'minval(Ta__DY(:,:))=',minval(Ta__DY(:,:))
529
530
531	! -----------------------------------------------------------------------------!
532	! Initialization of Atm_S0 (Insolation and cos of Sun Zenithal Distance)
533	! -------------------------
534
535	! **************
536	CALL PHY_Atm_S0_INI
537	! **************
538
539	! Initialization of Mod_PHY_S0_ctr and Mod_PHY_S0_kkl
540	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
541	FaceS0 = FlagS0_SLO
542	MMskS0 = FlagS0_MtM
543	IF (FlagS0_SLO .AND. FlagS0_MtM) THEN
544	DO k = 1,m_azim
545	DO ikl = 1,kcolp
546	cszkS0(ikl,k) = MMaskHOST(ikl,k)
547	ENDDO
548	ENDDO
549	ENDIF
550
551
552	! -----------------------------------------------------------------------------!
553	! Initialization of SISVAT
554	! -------------------------
555
556	! Initialization of Mod_SISVAT_ctr (SISVAT Switches
557	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ & Time/Space control Variables)
558	VegMod = FlagSV_Veg
559	SnoMod = FlagSV_SNo
560	BloMod = FlagSV_BSn
561	InpSWD = FlagSV_SWD
562	InpSBC = FlagSV_SBC
563	SVaKzT = FlagSV_KzT
564	SVaUBC = FlagSV_UBC
565
566
567	! Initialization of Mod_SISVAT_kkl (from INPUT from HOST Model)
568	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
569	DO ikp = 1,kcolp
570	i = ii__AP(ikp)
571	j = jj__AP(ikp)
572	DO mn = 1,mwp
573	DO k = 1,mzp
574	kk = mzpp - k
575	WindSV(ikp,mn,k) = Wind_HOST(ikp,k)
576	pkt0SV(ikp,mn,kk)= pkta_HOST(ikp,k)
577	END DO
578	Ua__SV(ikp,mn) = Ua___HOST(ikp,mzp)
579	Va__SV(ikp,mn) = Va___HOST(ikp,mzp)
580	END DO
581	END DO
582
583
584	! Initialization of SISVAT Variables
585	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
586	! **************
587	CALL PHY_SISVAT_INI
588	! **************
589
590
591	! Initialization of Mod_SISVAT_kkl (SISVAT OUTPUT Grid Points)
592	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
593	iwr = 0
594	DO ipt = 1,NbPts
595	IF (IOi0SV(ipt).EQ.0) THEN
596	IOi_SV(ipt)=i_x0
597	ELSE
598	IOi_SV(ipt)=IOi0SV(ipt)
599	END IF
600	IF (IOj0SV(ipt).EQ.0) THEN
601	IOj_SV(ipt)=j_y0
602	ELSE
603	IOj_SV(ipt)=IOj0SV(ipt)
604	END IF
605	write(6,*) 'ipt , IOi_SV, IOj_SV = ' &
606	& ,ipt , IOi_SV(ipt), IOj_SV(ipt)
607	DO n = 1,mwp
608	iwr = 1+iwr
609	IF (iwr.LE.nbwri) THEN
610	no__SV(iwr) = 0
611	i___SV(iwr) = IOi_SV(ipt)
612	j___SV(iwr) = IOj_SV(ipt)
613	n___SV(iwr) = n
614	write(6,*) 'n , i___SV, j___SV, n___SV, iwr = ' &
615	& ,n , i___SV(iwr), j___SV(iwr), n___SV(iwr), iwr
616	END IF
617	END DO
618	END DO
619
620
621
622	! -----------------------------------------------------------------------------!
623	! Initialization of Atm_RT (Radiative Transfert through the Atmosphere)
624	! -------------------------
625
626	! **************
627	CALL PHY_Atm_RT_INI
628	! **************
629
630
631
632	! -----------------------------------------------------------------------------!
633	! Initialization of Atm_AT (Turbulent Transfert through the Atmosphere)
634	! -------------------------
635
636	! **************
637	CALL PHY_Atm_AT_INI(FlagAT_TKE,TypeAT)
638	! **************
639
640
641
642	! -----------------------------------------------------------------------------!
643	! Initialization of Atm_CP (Convectiv Transfert through the Atmosphere)
644	! -------------------------
645
646	! **************
647	CALL PHY_Atm_CP_INI(mzp,kcolp)
648	! **************
649
650
651
652	! -----------------------------------------------------------------------------!
653	! Initialization of CMiPhy (Cloud Microphysical Scheme)
654	! -------------------------
655
656	! Initialization of Mod_PHY_CM_ctr (CMiPhy Switches
657	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ & Time/Space control Variables)
658	CM_UpD = FlagCM_UpD
659
660	! **************
661	CALL PHY_Atm_CM_INI
662	! **************
663
664	END IF
665
666
667	! Martin CONTROL
668	!PRINT*, 'Impressions control'
669	!call iophys_ecrit('TA__DY_surf',1,'surface temperature','K',Ta__DY(:,61))
670	!call iophys_ecrit('TA__DY_air',60,'air temperature','K',Ta__DY(:,1:60))
671
672	!PRINT*,'TA__DY(:,61)=',TA__DY(kcolq/2,61)
673
674	! Interface: From HOST Model Variables to MARp Physics Variables
675	! ==============================================================
676
677
678	! Time
679	! ----
680
681	YearTU = Year_H
682	Mon_TU = Mon__H
683	Day_TU = Day__H
684	HourTU = Hour_H
685	minuTU = minu_H
686	sec_TU = sec__H
687
688	TimeTU = (float(351)+(float(YearTU) -float(1902)) *float(365)&! Nb Days before YearTU
689	& +(float(YearTU) -float(1901)) /float( 4)&! Nb Leap Years
690	& + float(njYear(Mon_TU)) &! Nb Days before Mon_TU
691	& + float(njLeap(Mon_TU)) &! (including Leap Day)
692	& *max(zer0,un_1-mod(float(YearTU),float(4))) &!
693	& + float(Day_TU)-float(1) ) *float( 24)&!
694	& +float(HourTU) &!
695	& + (float(minuTU) *float(60) +float(sec_TU))/3600. !
696
697
698
699	! -----------------------------------------------------------------------------!
700	! Assignation of Mod_PHY_AT_grd
701	! ---------------------------------
702	IF (it_RUN .EQ. 1) THEN
703	TimeAT = TimeTU-dt0_AT/3600. ! Initialisation à la première itération
704	END IF
705	IF (FlagAT .AND. mod(it_RUN-1,jt__AT).EQ.0) THEN
706	dTimAT = (TimeTU-TimeAT) * 3600.
707	TimeAT = TimeTU
708	END IF
709
710	! Martin CONTROL
711	PRINT*,'CONTROL PHY_MAR temps'
712	PRINT*,' Year_H =', Year_H
713	PRINT*,' Mon__H =', Mon__H
714	PRINT*,' Day__H =', Day__H
715	PRINT*,' Hour_H =', Hour_H
716	PRINT*,' minu_H =', minu_H
717	PRINT*,' sec__H =', sec__H
718	PRINT*,'jt__AT=',jt__AT
719	PRINT*,'TimeTU=',TimeTU
720	PRINT*,'TimeAT=',TimeAT
721	PRINT*,'dTimAT=',dTimAT
722	! Martin CONTROL
723
724
725	! -----------------------------------------------------------------------------!
726	! Assignation of Mod_PHY_DY_kkl
727	! ---------------------------------
728
729	DO ikl=1,kcolp
730	i = ii__AP(ikl)
731	j = jj__AP(ikl)
732	k = mzpp
733	ExnrDY (ikl,k) = exp(RCp log(psa__HOST(ikl)sigma(k)+pt__DY))
734	pkt_DY (ikl,k) = pkta_HOST(ikl,k)
735	Ta__DY (ikl,k) = pkta_HOST(ikl,k) * ExnrDY(ikl,k)
736	! CAUTION: Tas_SV_xy is not allowed to be changed by data coming from outside this routine
737	Z___DY (ikl,k) = gZa__HOST(ikl,k) *Grav_I
738	ZmidDY (ikl,k) = gZam_HOST(ikl,k) *Grav_I
739	qv__DY (ikl,k) = qv___HOST(ikl,k)
740
741	DO k = 1,mzp
742	ExnrDY (ikl,k) = exp(RCp log(psa__HOST(ikl)sigma(k )+pt__DY))
743	pkt_DY (ikl,k) = pkta_HOST(ikl,k)
744	Ta__DY (ikl,k) = pkta_HOST(ikl,k) * ExnrDY(ikl,k)
745	roa_DY (ikl,k) = (psa__HOST(ikl)*sigma(k )+pt__DY) &
746	& /(Ta__DY (ikl,k) *R_DAir)
747	roamDY (ikl,k) = (psa__HOST(ikl)*sigmi(k+1)+pt__DY) &
748	& /(Ta__DY (ikl,k) *R_DAir)
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
754
755	! -----------------------------------------------------------------------------!
756	! Assignation of Mod_PHY_AT_kkl
757	! ---------------------------------
758
759
760	IF (FlagAT.AND.it_EXP.GT.1.AND.mod(it_RUN-1,jt__AT).EQ.0) THEN
761	TKE_AT (ikl,k) = TKE__HOST(ikl,k)
762	eps_AT (ikl,k) = eps__HOST(ikl,k)
763	TrT_AT (ikl,k) = (TKE_AT (ikl,k) - TrT_AT (ikl,k)) &
764	& / dTimAT
765	END IF
766
767
768
769	! -----------------------------------------------------------------------------!
770	! Assignation of Mod_PHY_CM_kkl
771	! ---------------------------------
772
773	IF (FlagCM.AND.it_EXP.GT.1.AND.mod(it_RUN-1,jt__CM).EQ.0) THEN
774
775	! IF (qw___HOST(ikl,k).LT.qh_MIN) THEN
776	! qv__DY(ikl,k) = qv__DY(ikl,k) + qw__CM(ikl,k)
777	! qw__CM(ikl,k) = 0.
778	! #cw CCNwCM(ikl,k) = 0.
779	! ELSE
780	qw__CM(ikl,k) = qw___HOST(ikl,k)
781	! #cw CCNwCM(ikl,k) = CCN__HOST(ikl,k)
782	! END IF
783
784	! IF (qi___HOST(ikl,k).LT.qh_MIN) THEN
785	! qv__DY(ikl,k) = qv__DY(ikl,k) + qi__CM(ikl,k)
786	! qi__CM(ikl,k) = 0.
787	! CCNiCM(ikl,k) = 0.
788	! ELSE
789	qi__CM(ikl,k) = qi___HOST(ikl,k)
790	CCNiCM(ikl,k) = CIN__HOST(ikl,k)
791	! END IF
792
793	! Gilles: CF___HOST non sauve & CFraCM reinitialise ici
794	! CFraCM(ikl,k) = CF___HOST(ikl,k)
795
796	! IF (qw__CM(ikl,k).LT.qh_MIN .AND. &
797	! & qi__CM(ikl,k).LT.qh_MIN) THEN
798	! CFraCM(ikl,k) = 0.
799	! ELSE
800	! CFraCM(ikl,k) = max(CFrMIN,CF___HOST(ikl,k))
801	! END IF
802
803	qs__CM(ikl,k) = qs___HOST(ikl,k)
804	qr__CM(ikl,k) = qr___HOST(ikl,k)
805	END IF
806
807	ENDDO
808
809
810
811	! -----------------------------------------------------------------------------!
812	! Assignation of Mod_SISVAT_gpt (A-O FORCING OR COUPLING)
813	! ---------------------------------
814
815	sst_SB (ikl) = sst__HOST(ikl)
816	! #IP sif_SB (ikl) = sif__HOST(ikl)
817	! #AO DO k = 1,mwp
818	! #AO s_T_AO_xyn(i,j,k) = s_T__HOST(ikl,k)
819	! #AO Alb_AO_xyn(i,j,k) = Alb__HOST(ikl,k)
820	! #AO ENDDO
821
822	ENDDO
823
824
825
826	! -----------------------------------------------------------------------------!
827	! Assignation of Mod_PHY_DY_kkl
828	! ---------------------------------
829	DO ikl = 1,kcolp
830	i = ii__AP(ikl)
831	j = jj__AP(ikl)
832	psa_DY(ikl) = psa__HOST(ikl)
833	DO k = 1,mzp
834	WindDY(ikl,k) = Wind_HOST(ikl,k)
835	ua__DY(ikl,k) = Ua___HOST(ikl,k)
836	va__DY(ikl,k) = va___HOST(ikl,k)
837	wa__DY(ikl,k) = Wa___HOST(ikl,k) + sqrt(2.*max(eps6,TKE_AT(ikl,k))/3.)
838	END DO
839	END DO
840
841	! **************
842	CALL PHY_Atm_DY_RUN ! Assignation of MAR Dyn. Variables
843	! **************
844
845	IF(FlagVR) THEN
846	! **************
847	CALL PHY________OUT('After PHY_Atm_DY_RUN ')
848	! ************** 12345678901234567890123456789012345678901234567890
849	! 1 2 3 4 5
850	END IF
851
852
853
854	! -----------------------------------------------------------------------------!
855	! Saturation Specific Humidity
856	! ----------------------------
857	IF (FlagCM .OR. &! ***************
858	& FlagSV) CALL PHY_Atm_CM_QSat
859	! ***************
860
861
862
863	! -----------------------------------------------------------------------------!
864	! Execution of Atm_S0 (Insolation and cos of Sun Zenithal Distance)
865	! -------------------------
866
867	! **************
868	CALL PHY_Atm_S0_RUN
869	! **************
870
871	IF(FlagVR) THEN
872	! **************
873	CALL PHY________OUT('After PHY_Atm_S0_RUN ')
874	! ************** 12345678901234567890123456789012345678901234567890
875	! 1 2 3 4 5
876	END IF
877
878
879
880	! -----------------------------------------------------------------------------!
881	! Execution of Atm_RT (Radiative Transfer through the Atmosphere)
882	! -------------------------
883
884	IF (FlagRT .AND. mod(it_RUN-1,jt__RT).EQ.0) THEN
885
886	! **************
887	CALL PHY_Atm_RT_RUN(kcolp,ikl0)
888	! **************
889
890	IF(FlagVR) THEN
891	! **************
892	CALL PHY________OUT('After PHY_Atm_RT_RUN ')
893	! ************** 12345678901234567890123456789012345678901234567890
894	! 1 2 3 4 5
895	END IF
896	END IF
897
898
899
900	! -----------------------------------------------------------------------------!
901	! Execution of SISVAT (Soil-Ice-Snow-Vegetation-Atmosphere-Transfer Scheme)
902	! -------------------------
903
904	IF (FlagSV) THEN
905
906	! Assignation of Mod_SISVAT_kkl (from INPUT from HOST Model)
907	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
908	DO ikp = 1,kcolp
909	DO mn = 1,mwp
910	pkPaSV(ikp,mn) = psa__HOST(ikp) +pt__DY
911	ExnrSV(ikp,mn) = exp(RCp log(psa__HOST(ikp)sigma(mzp)+pt__DY))
912	DO k = 1,mzp
913	kk = mzpp - k
914	WindSV(ikp,mn,k) = Wind_HOST(ikp,k)
915	pkt0SV(ikp,mn,kk)= pkta_HOST(ikp,k)
916	qv__SV(ikp,mn,kk)= qv___HOST(ikp,k)
917	zza_SV(ikp,mn,kk)= (gZa__HOST(ikp,k) & !
918	& -gZa__HOST(ikp,mzpp)) *Grav_I !
919	roa_SV(ikp,mn,kk)= roa_DY (ikp,k) * 1.e3 ! [kg/m3]
920	END DO
921	Ua__SV(ikp,mn) = Ua___HOST(ikp,mzp)
922	Va__SV(ikp,mn) = Va___HOST(ikp,mzp)
923	TaT_SV(ikp,mn) = pkta_HOST(ikp,mzp) * ExnrSV(ikp,mn)
924	END DO
925	END DO
926
927
928	! **************
929	CALL PHY_SISVAT_RUN( &
930	! **************
931	! -----------------------------------------------------------------------------!
932	! #TC& qxTC,uqTC ,qsTC , & ! Aerosols: Atm.Conc., Surf.Flux
933	! -----------------------------------------------------------------------------!
934	& )
935
936	IF(FlagVR) THEN
937	! **************
938	CALL PHY________OUT('After PHY_SISVAT_RUN ')
939	! ************** 12345678901234567890123456789012345678901234567890
940	! 1 2 3 4 5
941	END IF
942	END IF
943
944
945
946	! -----------------------------------------------------------------------------!
947	! Execution of Atm_AT (Atmospheric Turbulence Contribution)
948	! -------------------------
949
950	IF (FlagAT .AND. mod(it_RUN-1,jt__AT).EQ.0) THEN
951
952	! **************
953	CALL PHY_Atm_AT_RUN(FlagSV_KzT,FlagCM)
954	! **************
955
956	IF(FlagVR) THEN
957	! **************
958	CALL PHY________OUT('After PHY_Atm_AT_RUN ')
959	! ************** 12345678901234567890123456789012345678901234567890
960	! 1 2 3 4 5
961	END IF
962	END IF
963
964
965
966	! -----------------------------------------------------------------------------!
967	! Execution of Atm_CP (Convection Parameterisation - Mass Flux Scheme)
968	! -------------------------
969
970	IF (FlagCP .AND. mod(it_RUN-1,jt__CP).EQ.0) THEN
971
972	! **************
973	CALL PHY_Atm_CP_RUN(mzp,kcolp)
974	! CALL PHY_Atm_CP_RUN(mxp,kcolp) ! ancien bug
975	! **************
976
977	IF(FlagVR) THEN
978	! **************
979	CALL PHY________OUT('After PHY_Atm_CP_RUN ')
980	! ************** 12345678901234567890123456789012345678901234567890
981	! 1 2 3 4 5
982	END IF
983	END IF
984
985
986
987	! -----------------------------------------------------------------------------!
988	! Execution of Atm_CM (Cloud Microphysics)
989	! -------------------------
990
991	IF (FlagCM .AND. mod(it_RUN-1,jt__CM).EQ.0) THEN
992
993	! **************
994	CALL PHY_Atm_CM_RUN
995	! **************
996
997
998	IF(FlagVR) THEN
999	! **************
1000	CALL PHY________OUT('After PHY_Atm_CM_RUN ')
1001	! ************** 12345678901234567890123456789012345678901234567890
1002	! 1 2 3 4 5
1003	END IF
1004	END IF
1005
1006
1007
1008
1009	! Assignation of the Tendencies to transfer outside the physical Parameterizations Package
1010	! ========================================================================================
1011
1012	DO ikl=1,kcolp
1013	i = ii__AP(ikl)
1014	j = jj__AP(ikl)
1015
1016
1017
1018	! -----------------------------------------------------------------------------!
1019	! Reinitialization of the Tendencies
1020	! ----------------------------------
1021
1022	DO k= 1,mzp
1023
1024	dpkt___dt(ikl,k) = 0.0 !
1025	dua____dt(ikl,k) = 0.0 !
1026	dva____dt(ikl,k) = 0.0 !
1027	dqv____dt(ikl,k) = 0.0 !
1028	dqw____dt(ikl,k) = 0.0 !
1029	dCF____dt(ikl,k) = 0.0 !
1030	dqi____dt(ikl,k) = 0.0 !
1031	dCi____dt(ikl,k) = 0.0 !
1032	dqs____dt(ikl,k) = 0.0 !
1033	dqr____dt(ikl,k) = 0.0 !
1034
1035
1036
1037	! -----------------------------------------------------------------------------!
1038	! Tendencies from SISVAT
1039	! ----------------------
1040
1041	IF (FlagSV .AND. mod(it_RUN-1,jt__SV).EQ.0) THEN!
1042	dpktSV_dt(ikl,k) = dpktSV_gpt(ikl,k) ! Reduced Potential Temperature TENDENCY, SISVAT [KX/s]
1043	END IF
1044
1045	! Update of the tendencies
1046	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1047	IF (FlagSV) THEN!
1048	dpkt___dt(ikl,k) = dpkt___dt(ikl,k) + dpktSV_gpt(ikl,k) !
1049	END IF
1050
1051
1052
1053	! -----------------------------------------------------------------------------!
1054	! Tendencies from Atm_AT
1055	! ----------------------
1056
1057	IF (FlagAT .AND. mod(it_RUN-1,jt__AT).EQ.0) THEN!
1058	dua____dt(ikl,k) = dua_AT(ikl,k) ! Wind Speed (x-direc.) TENDENCY, ALL Contribut.[m/s2]
1059	dva____dt(ikl,k) = dva_AT(ikl,k) ! Wind Speed (y-direc.) TENDENCY, ALL Contribut.[m/s2]
1060	dpktAT_dt(ikl,k) = dpktAT(ikl,k) ! Reduced Potential Temperature TENDENCY, Atm_AT [KX/s]
1061	dqv_AT_dt(ikl,k) = dqv_AT(ikl,k) ! Specific Humidity TENDENCY, Atm_AT [kg/kg/s]
1062	dqw_AT_dt(ikl,k) = dqw_AT(ikl,k) ! Cloud Droplets Concentration TENDENCY, Atm_AT [kg/kg/s]
1063	dqi_AT_dt(ikl,k) = dqi_AT(ikl,k) ! Cloud Crystals Concentration TENDENCY, Atm_AT [kg/kg/s]
1064	dqs_AT_dt(ikl,k) = dqs_AT(ikl,k) ! Snow Particles Concentration TENDENCY, Atm_AT [kg/kg/s]
1065	dqr_AT_dt(ikl,k) = dqr_AT(ikl,k) ! Rain Drops Concentration TENDENCY, Atm_AT [kg/kg/s]
1066	! #cw dCw_AT_dt(ikl,k) = dCW_AT(ikl,k) ! CCN Concentration TENDENCY, Atm_AT [1/s]
1067	dCi_AT_dt(ikl,k) = dCi_AT(ikl,k) ! CIN Concentration TENDENCY, Atm_AT [1/s]
1068
1069	! NO Tendencies on TKE and its dissipation ==> TKE__HOST and eps__HOST are updated
1070	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1071	TKE__HOST(ikl,k) = TKE_AT(ikl,k) !
1072	eps__HOST(ikl,k) = eps_AT(ikl,k) !
1073
1074	! Re-Initialization of the TKE Transport Rate
1075	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1076	TrT_AT (ikl,k) = TKE_AT(ikl,k) !
1077	END IF
1078
1079	! Update of the tendencies
1080	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1081	IF (FlagAT) THEN!
1082	dpkt___dt(ikl,k) = dpkt___dt(ikl,k) + dpktAT(ikl,k) !
1083	dqv____dt(ikl,k) = dqv____dt(ikl,k) + dqv_AT(ikl,k) !
1084	dqw____dt(ikl,k) = dqw____dt(ikl,k) + dqw_AT(ikl,k) !
1085	dqi____dt(ikl,k) = dqi____dt(ikl,k) + dqi_AT(ikl,k) !
1086	dqs____dt(ikl,k) = dqs____dt(ikl,k) + dqs_AT(ikl,k) !
1087	dqr____dt(ikl,k) = dqr____dt(ikl,k) + dqr_AT(ikl,k) !
1088	! #cw dCw____dt(ikl,k) = dCw____dt(ikl,k) + dCw_AT(ikl,k) !
1089	dCi____dt(ikl,k) = dCi____dt(ikl,k) + dCi_AT(ikl,k) !
1090	END IF
1091
1092
1093
1094	! -----------------------------------------------------------------------------!
1095	! Tendencies from CMiPhy (dqw_CM, dCw_CM, dqi_CM, dCi_CM
1096	! ---------------------- dqs_CM, dqr_CM are consumed in CMiPhy then reset to 0
1097	! if FlagCM_UpD = .TRUE. )
1098
1099	IF (FlagCM .AND. mod(it_RUN-1,jt__CM).EQ.0) THEN!
1100	dpktCM_dt(ikl,k) = dpktCM(ikl,k) ! Reduced Potential Temperature TENDENCY, CMiPhy [KX/s]
1101	dqv_CM_dt(ikl,k) = dqv_CM(ikl,k) ! Specific Humidity TENDENCY, CMiPhy [kg/kg/s]
1102	dqw_CM_dt(ikl,k) = dqw_CM(ikl,k) ! Cloud Droplets Concentration TENDENCY, CMiPhy [kg/kg/s]
1103	dCF_CM_dt(ikl,k) = dCF_CM(ikl,k) ! Cloud Fraction TENDENCY, CMiPhy [kg/kg/s]
1104	dqi_CM_dt(ikl,k) = dqi_CM(ikl,k) ! Cloud Crystals Concentration TENDENCY, CMiPhy [kg/kg/s]
1105	dqs_CM_dt(ikl,k) = dqs_CM(ikl,k) ! Snow Particles Concentration TENDENCY, CMiPhy [kg/kg/s]
1106	dqr_CM_dt(ikl,k) = dqr_CM(ikl,k) ! Rain Drops Concentration TENDENCY, CMiPhy [kg/kg/s]
1107	! #cw dCw_CM_dt(ikl,k) = dCW_CM(ikl,k) ! CCN Concentration TENDENCY, CMiPhy [1/s]
1108	dCi_CM_dt(ikl,k) = dCi_CM(ikl,k) ! CIN Concentration TENDENCY, CMiPhy [1/s]
1109	END IF
1110
1111	! Update of the tendencies
1112	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1113	IF (FlagCM) THEN!
1114	dpkt___dt(ikl,k) = dpkt___dt(ikl,k) + dpktCM(ikl,k) !
1115	dqv____dt(ikl,k) = dqv____dt(ikl,k) + dqv_CM(ikl,k) !
1116	dqw____dt(ikl,k) = dqw____dt(ikl,k) + dqw_CM(ikl,k) !
1117	dCF____dt(ikl,k) = dCF____dt(ikl,k) + dCF_CM(ikl,k) !
1118	dqi____dt(ikl,k) = dqi____dt(ikl,k) + dqi_CM(ikl,k) !
1119	dqs____dt(ikl,k) = dqs____dt(ikl,k) + dqs_CM(ikl,k) !
1120	dqr____dt(ikl,k) = dqr____dt(ikl,k) + dqr_CM(ikl,k) !
1121	! #cw dCw____dt(ikl,k) = dCw____dt(ikl,k) + dCw_CM(ikl,k) !
1122	dCi____dt(ikl,k) = dCi____dt(ikl,k) + dCi_CM(ikl,k) !
1123	END IF
1124
1125
1126
1127	! -----------------------------------------------------------------------------!
1128	! Tendencies from CVAmnh
1129	! ----------------------
1130
1131	IF (FlagCP .AND. mod(it_RUN-1,jt__CP).EQ.0) THEN!
1132	dpktCP_dt(ikl,k) = dpktCP(ikl,k) ! Reduced Potential Temperature TENDENCY, CVAmnh [KX/s]
1133	dqv_CP_dt(ikl,k) = dqv_CP(ikl,k) ! Specific Humidity TENDENCY, CVAmnh [kg/kg/s]
1134	dqw_CP_dt(ikl,k) = dqw_CP(ikl,k) ! Cloud Droplets Concentration TENDENCY, CVAmnh [kg/kg/s]
1135	dqi_CP_dt(ikl,k) = dqi_CP(ikl,k) ! Cloud Crystals Concentration TENDENCY, CVAmnh [kg/kg/s]
1136	END IF
1137
1138	! Update of the tendencies
1139	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1140	IF (FlagCP) THEN!
1141	dpkt___dt(ikl,k) = dpkt___dt(ikl,k) + dpktCP(ikl,k) !
1142	dqv____dt(ikl,k) = dqv____dt(ikl,k) + dqv_CP(ikl,k) !
1143	dqw____dt(ikl,k) = dqw____dt(ikl,k) + dqw_CP(ikl,k) !
1144	dqi____dt(ikl,k) = dqi____dt(ikl,k) + dqi_CP(ikl,k) !
1145	END IF
1146
1147
1148
1149	! -----------------------------------------------------------------------------!
1150	! Tendencies from radCEP
1151	! ----------------------
1152
1153	IF (FlagRT .AND. mod(it_RUN-1,jt__RT).EQ.0) THEN!
1154	dpktRT_dt(ikl,k) = dpktRT(ikl,k) ! Reduced Potential Temperature TENDENCY, radCEP [KX/s]
1155	END IF
1156
1157	! Update of the tendencies
1158	! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
1159	IF (FlagRT) THEN!
1160	dpkt___dt(ikl,k) = dpkt___dt(ikl,k) + dpktRT(ikl,k) !
1161	END IF
1162
1163	ENDDO
1164
1165
1166
1167
1168
1169	! Assignation of the Variables to keep inside the Physical Parameterizations Package
1170	! ========================================================================================
1171
1172	! ... NOTHING up to now
1173
1174
1175
1176	! Assignation of the Variables to transfer outside the Physical Parameterizations Package
1177	! ========================================================================================
1178
1179	! -----------------------------------------------------------------------------!
1180	! Update of pkta_HOST, qv___HOST
1181	! -------------------------------
1182
1183	DO k=1,mzpp
1184	pkta_HOST(ikl,k) = pkt_DY(ikl,k) ! Always on k=mzpp, possible from 1 to mzp (with dpkt = 0.)
1185	qv___HOST(ikl,k) = qv__DY(ikl,k) !
1186	ENDDO
1187
1188
1189
1190	! -----------------------------------------------------------------------------!
1191	! Update of qw___HOST, CF___HOST, qi___HOST, CIN__HOST, qs___HOST, qr___HOST
1192	! ---------------------------------------------------------------------------
1193
1194	IF (FlagCM) THEN
1195	DO k=1,mzp
1196	qw___HOST(ikl,k) = qw__CM(ikl,k)
1197	! #cw CCN__HOST(ikl,k) = CCNwCM(ikl,k)
1198	CF___HOST(ikl,k) = CFraCM(ikl,k)
1199	qi___HOST(ikl,k) = qi__CM(ikl,k)
1200	CIN__HOST(ikl,k) = CCNiCM(ikl,k)
1201	qs___HOST(ikl,k) = qs__CM(ikl,k)
1202	qr___HOST(ikl,k) = qr__CM(ikl,k)
1203	ENDDO
1204	END IF
1205
1206
1207
1208	! -----------------------------------------------------------------------------!
1209	! Update of d(S,LH)/dT (needed in NEMO)
1210	! ---------------------
1211
1212	! #AO DO k=1,mwpp
1213	! #AO dSdT2HOST(ikl,k) = dSdTAO_xyn(i,j,k)
1214	! #AO dLdT2HOST(ikl,k) = dLdTAO_xyn(i,j,k)
1215	! #AO ENDDO
1216
1217
1218
1219	ENDDO
1220
1221
1222
1223
1224	! OUTPUT
1225	! ======
1226
1227	! OUTPUT of Tendencies
1228	! --------------------
1229
1230	IF (FlagVR .OR. &
1231	& (FLAG_O .AND. ((minuTU.EQ.0 .AND. sec_TU.EQ.0) .OR. &
1232	& it_RUN.EQ.1 ))) THEN
1233	ikl = ikl0
1234	i = ii__AP(ikl)
1235	j = jj__AP(ikl)
1236
1237	! pkt TENDENCIES
1238	! ~~~~~~~~~~~~~~
1239	write(4,400)
1240	400 format(//,' pkt TENDENCIES',/,' **************'/,1x)
1241	write(4,403) Day_TU,Mon_TU,YearTU,HourTU,MinuTU,Sec_TU,it_EXP
1242	403 format(3x,2(i2,'-'),i4,4x,3(i2,'-'),' Simulation Iteration No ',i6,/,1x)
1243	write(4,404)
1244	write(4,401)
1245	401 format(' \| SISVAT \| Atm_AT \| CMiPhy \| CVAmnh \| radCEP \|')
1246	write(4,402)
1247	402 format(' \| [K/d] \| [K/d] \| [K/d] \| [K/d] \| [K/d] \|')
1248	write(4,404)
1249	404 format(4('-'),'+',5(13('-'),'+'))
1250	DO k=1,mzp
1251	write(4,405) k &
1252	& , ExnrDY(ikl,k)dpktSV_gpt(ikl,k)86400. &
1253	& , ExnrDY(ikl,k)dpktAT (ikl,k)86400. &
1254	& , ExnrDY(ikl,k)dpktCM (ikl,k)86400. &
1255	& , ExnrDY(ikl,k)dpktCP (ikl,k)86400. &
1256	& , ExnrDY(ikl,k)dpktRT (ikl,k)86400.
1257	405 format(i3,' \|',5(f12.6,' \|'))
1258
1259	IF (mod(k,20).EQ.0) THEN
1260	write(4,404)
1261	write(4,401)
1262	write(4,402)
1263	write(4,404)
1264	END IF
1265	ENDDO
1266	write(4,404)
1267
1268	! pkt TENDENCIES
1269	! ~~~~~~~~~~~~~~
1270	write(4,410)
1271	410 format(//,' Qv TENDENCIES',/,' **************'/,1x)
1272	write(4,403) Day_TU,Mon_TU,YearTU,HourTU,MinuTU,Sec_TU,it_EXP
1273	!403 format(3x,2(i2,'-'),i4,4x,3(i2,'-'),' Simulation Iteration No ',i6,/,1x)
1274	write(4,404)
1275	write(4,411)
1276	411 format(' \| SISVAT \| Atm_AT \| CMiPhy \| CVAmnh \| TOTAL \|')
1277	write(4,412)
1278	412 format(' \| [g/kg/min] \| [g/kg/min] \| [g/kg/min] \| [g/kg/min] \| [g/kg/min] \|')
1279	write(4,404)
1280	!404 format(4('-'),'+',5(13('-'),'+'))
1281	DO k=1,mzp
1282	write(4,415) k &
1283	& , dqv_AT (ikl,k)*60000. &
1284	& , dqv_CM (ikl,k)*60000. &
1285	& , dqv_CP (ikl,k)*60000. &
1286	& , dqv____dt(ikl,k)*60000.
1287	415 format(i3,' \|',12x,' \|',4(f12.6,' \|'))
1288
1289	IF (mod(k,20).EQ.0) THEN
1290	write(4,404)
1291	write(4,411)
1292	write(4,412)
1293	write(4,404)
1294	END IF
1295	ENDDO
1296	write(4,404)
1297	END IF
1298
1299	IF (FLAG_O .AND. ((minuTU.EQ.0 .AND. sec_TU.EQ.0) .OR. &
1300	& it_RUN.EQ.1 )) THEN
1301
1302	! **************
1303	CALL PHY________OUT('After PHY_MAR ')
1304	! ************** 12345678901234567890123456789012345678901234567890
1305	! 1 2 3 4 5
1306	END IF
1307
1308
1309
1310	return
1311	end

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format