| 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 |
|---|
| 723 | PRINT*,'CONTROL PHY_MAR temps' |
|---|
| 724 | PRINT*,' Year_H =', Year_H |
|---|
| 725 | PRINT*,' Mon__H =', Mon__H |
|---|
| 726 | PRINT*,' Day__H =', Day__H |
|---|
| 727 | PRINT*,' Hour_H =', Hour_H |
|---|
| 728 | PRINT*,' minu_H =', minu_H |
|---|
| 729 | PRINT*,' sec__H =', sec__H |
|---|
| 730 | PRINT*,'jt__AT=',jt__AT |
|---|
| 731 | PRINT*,'TimeTU=',TimeTU |
|---|
| 732 | PRINT*,'TimeAT=',TimeAT |
|---|
| 733 | PRINT*,'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 |
|---|
