| 1 | MODULE VARxSV |
|---|
| 2 | |
|---|
| 3 | USE VAR_SV, only : klonv, nsol, nsno, nb_wri |
|---|
| 4 | |
|---|
| 5 | IMPLICIT NONE |
|---|
| 6 | ! +--SISVAT INPUT Variables |
|---|
| 7 | ! + ----------------------------- |
|---|
| 8 | |
|---|
| 9 | INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: LSmask ! Land-Sea Mask |
|---|
| 10 | !$OMP THREADPRIVATE(LSmask) |
|---|
| 11 | INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: isotSV ! Soil Type |
|---|
| 12 | !$OMP THREADPRIVATE(isotSV) |
|---|
| 13 | INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: iWaFSV ! Soil Drainage:(1,0)=(y,n) |
|---|
| 14 | !$OMP THREADPRIVATE(iWaFSV) |
|---|
| 15 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: coszSV ! Cosine of Sun zenithal Angle |
|---|
| 16 | !$OMP THREADPRIVATE(coszSV) |
|---|
| 17 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: sol_SV ! Downward Solar Radiation |
|---|
| 18 | !$OMP THREADPRIVATE(sol_SV) |
|---|
| 19 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: IRd_SV ! Downward Longwave Radiation |
|---|
| 20 | !$OMP THREADPRIVATE(IRd_SV) |
|---|
| 21 | |
|---|
| 22 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: drr_SV ! Rain Intensity [kg/m2/s] |
|---|
| 23 | !$OMP THREADPRIVATE(drr_SV) |
|---|
| 24 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: dsn_SV ! Snow Intensity [kg/m2/s] |
|---|
| 25 | !$OMP THREADPRIVATE(dsn_SV) |
|---|
| 26 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: dsnbSV ! Idem, fraction, from Drift [-] |
|---|
| 27 | !$OMP THREADPRIVATE(dsnbSV) |
|---|
| 28 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: esnbSV ! Idem, fraction, from Drift [-] |
|---|
| 29 | !$OMP THREADPRIVATE(esnbSV) |
|---|
| 30 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: dbs_SV ! Drift Amount [kg/m2] |
|---|
| 31 | !$OMP THREADPRIVATE(dbs_SV) |
|---|
| 32 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: dbs_Er ! BS Erosion [kg/m2] |
|---|
| 33 | !$OMP THREADPRIVATE(dbs_Er) |
|---|
| 34 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: dbs_Ac ! BS Accumulation [kg/m2] |
|---|
| 35 | !$OMP THREADPRIVATE(dbs_Ac) |
|---|
| 36 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: BrosSV ! Buffer Snow Layer Density |
|---|
| 37 | !$OMP THREADPRIVATE(BrosSV) |
|---|
| 38 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: BG1sSV ! Buffer Snow Layer Dendr/Sphe[-] |
|---|
| 39 | !$OMP THREADPRIVATE(BG1sSV) |
|---|
| 40 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: BG2sSV ! Buffer Snow Layer Spher/Size[-][0.0001m] |
|---|
| 41 | !$OMP THREADPRIVATE(BG2sSV) |
|---|
| 42 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: dz0_SV ! dz0(Sastrugi dh) [m] |
|---|
| 43 | !$OMP THREADPRIVATE(dz0_SV) |
|---|
| 44 | |
|---|
| 45 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: cld_SV ! Cloudiness (seen from SBL) |
|---|
| 46 | !$OMP THREADPRIVATE(cld_SV) |
|---|
| 47 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: za__SV ! SBL Height |
|---|
| 48 | !$OMP THREADPRIVATE(za__SV) |
|---|
| 49 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: VV__SV !(SBL Top) Wind Velocity |
|---|
| 50 | !$OMP THREADPRIVATE(VV__SV) |
|---|
| 51 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: VVs_SV !(Sastr,V) Relevance |
|---|
| 52 | !$OMP THREADPRIVATE(VVs_SV) |
|---|
| 53 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: RRs_SV !(Sastr,V) Counter |
|---|
| 54 | !$OMP THREADPRIVATE(RRs_SV) |
|---|
| 55 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: DDs_SV !(Sastr,V) Angle |
|---|
| 56 | !$OMP THREADPRIVATE(DDs_SV) |
|---|
| 57 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: TaT_SV ! SBL Top Temperature |
|---|
| 58 | !$OMP THREADPRIVATE(TaT_SV) |
|---|
| 59 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: ExnrSV ! Exner Potential |
|---|
| 60 | !$OMP THREADPRIVATE(ExnrSV) |
|---|
| 61 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: dSdTSV ! Sensible Heat Flux T Derivat. |
|---|
| 62 | !$OMP THREADPRIVATE(dSdTSV) |
|---|
| 63 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: dLdTSV ! Latent Heat Flux T Derivat. |
|---|
| 64 | !$OMP THREADPRIVATE(dLdTSV) |
|---|
| 65 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: rhT_SV ! SBL Top Air Density |
|---|
| 66 | !$OMP THREADPRIVATE(rhT_SV) |
|---|
| 67 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: QaT_SV ! SBL Top Specific Humidity |
|---|
| 68 | !$OMP THREADPRIVATE(QaT_SV) |
|---|
| 69 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: QsT_SV ! SBL Top Specific Humidity |
|---|
| 70 | !$OMP THREADPRIVATE(QsT_SV) |
|---|
| 71 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: dQa_SV ! SBL Flux Limitation of Qa |
|---|
| 72 | !$OMP THREADPRIVATE(dQa_SV) |
|---|
| 73 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: qsnoSV ! SBL Mean Snow Content |
|---|
| 74 | !$OMP THREADPRIVATE(qsnoSV) |
|---|
| 75 | |
|---|
| 76 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: alb0SV ! Soil Albedo |
|---|
| 77 | !$OMP THREADPRIVATE(alb0SV) |
|---|
| 78 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: slopSV ! Snow/Ice/Soil-Water Surf. Slope |
|---|
| 79 | !$OMP THREADPRIVATE(slopSV) |
|---|
| 80 | |
|---|
| 81 | |
|---|
| 82 | REAL,SAVE :: dt__SV ! Time Step |
|---|
| 83 | !$OMP THREADPRIVATE(dt__SV) |
|---|
| 84 | CHARACTER (len=18),SAVE :: daHost ! Date Host Model |
|---|
| 85 | !$OMP THREADPRIVATE(daHost) |
|---|
| 86 | |
|---|
| 87 | |
|---|
| 88 | ! +--SISVAT INPUT/OUTPUT Variables |
|---|
| 89 | ! + ----------------------------- |
|---|
| 90 | |
|---|
| 91 | INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: isnoSV ! Nb of Ice/Snow Layers |
|---|
| 92 | !$OMP THREADPRIVATE(isnoSV) |
|---|
| 93 | INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: ispiSV ! Uppermost superimposed ice |
|---|
| 94 | !$OMP THREADPRIVATE(ispiSV) |
|---|
| 95 | INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: iiceSV ! Nb of Ice Layers |
|---|
| 96 | !$OMP THREADPRIVATE(iiceSV) |
|---|
| 97 | INTEGER ,ALLOCATABLE,SAVE :: istoSV(:,:) ! Snow Layer History |
|---|
| 98 | !$OMP THREADPRIVATE(istoSV) |
|---|
| 99 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: alb_SV ! Surface Albedo |
|---|
| 100 | !$OMP THREADPRIVATE(alb_SV) |
|---|
| 101 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: emi_SV ! Surface Emissivity |
|---|
| 102 | !$OMP THREADPRIVATE(emi_SV) |
|---|
| 103 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: IRs_SV ! Soil IR Flux |
|---|
| 104 | !$OMP THREADPRIVATE(IRs_SV) |
|---|
| 105 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: LMO_SV ! Monin-Obukhov Scale |
|---|
| 106 | !$OMP THREADPRIVATE(LMO_SV) |
|---|
| 107 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: us__SV ! Friction Velocity |
|---|
| 108 | !$OMP THREADPRIVATE(us__SV) |
|---|
| 109 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: uts_SV ! Temperature Turbulent Scale |
|---|
| 110 | !$OMP THREADPRIVATE(uts_SV) |
|---|
| 111 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: cutsSV ! Temperature Turbulent Scale C. |
|---|
| 112 | !$OMP THREADPRIVATE(cutsSV) |
|---|
| 113 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: uqs_SV ! Spec.Humid. Turbulent Scale |
|---|
| 114 | !$OMP THREADPRIVATE(uqs_SV) |
|---|
| 115 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: uss_SV ! Blow.Snow Turbulent Scale |
|---|
| 116 | !$OMP THREADPRIVATE(uss_SV) |
|---|
| 117 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: usthSV ! Blowing Snow Erosion Thresh. |
|---|
| 118 | !$OMP THREADPRIVATE(usthSV) |
|---|
| 119 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: rCDmSV ! Square Root Contribut. Drag_m |
|---|
| 120 | !$OMP THREADPRIVATE(rCDmSV) |
|---|
| 121 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: rCDhSV ! Square Root Contribut. Drag_h |
|---|
| 122 | !$OMP THREADPRIVATE(rCDhSV) |
|---|
| 123 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: Z0m_SV ! Momentum Roughness Length |
|---|
| 124 | !$OMP THREADPRIVATE(Z0m_SV) |
|---|
| 125 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: Z0mmSV ! z0(Momentum, Time Mean) [m] |
|---|
| 126 | !$OMP THREADPRIVATE(Z0mmSV) |
|---|
| 127 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: Z0mnSV ! z0(Momentum, instanta.) [m] |
|---|
| 128 | !$OMP THREADPRIVATE(Z0mnSV) |
|---|
| 129 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: Z0roSV ! Subgrid Topo Roughness Length |
|---|
| 130 | !$OMP THREADPRIVATE(Z0roSV) |
|---|
| 131 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: Z0SaSV ! z0(Sastrugi h) [m] |
|---|
| 132 | !$OMP THREADPRIVATE(Z0SaSV) |
|---|
| 133 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: Z0e_SV ! z0(Snow eroded) [m] |
|---|
| 134 | !$OMP THREADPRIVATE(Z0e_SV) |
|---|
| 135 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: Z0emSV ! z0(Snow eroded, Time Mean) [m] |
|---|
| 136 | !$OMP THREADPRIVATE(Z0emSV) |
|---|
| 137 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: Z0enSV ! z0(Snow eroded, instanta.) [m] |
|---|
| 138 | !$OMP THREADPRIVATE(Z0enSV) |
|---|
| 139 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: Z0h_SV ! Heat Roughness Length |
|---|
| 140 | !$OMP THREADPRIVATE(Z0h_SV) |
|---|
| 141 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: Z0hmSV ! z0(Heat, Time Mean) [m] |
|---|
| 142 | !$OMP THREADPRIVATE(Z0hmSV) |
|---|
| 143 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: Z0hnSV ! z0(Heat, instanta.) [m] |
|---|
| 144 | !$OMP THREADPRIVATE(Z0hnSV) |
|---|
| 145 | |
|---|
| 146 | |
|---|
| 147 | |
|---|
| 148 | REAL ,ALLOCATABLE,SAVE :: TsisSV(:,:) ! Snow/Ice/Soil-Water Temperature |
|---|
| 149 | !$OMP THREADPRIVATE(TsisSV) |
|---|
| 150 | REAL ,ALLOCATABLE,SAVE :: ro__SV(:,:) ! Snow/Ice/Soil-Water VolumicMass |
|---|
| 151 | !$OMP THREADPRIVATE(ro__SV) |
|---|
| 152 | REAL,ALLOCATABLE,SAVE :: eta_SV(:,:) ! Snow/Ice/Soil Water Content |
|---|
| 153 | !$OMP THREADPRIVATE(eta_SV) |
|---|
| 154 | REAL,ALLOCATABLE,SAVE :: G1snSV(:,:) ! Snow Dendricity/Sphericity |
|---|
| 155 | !$OMP THREADPRIVATE(G1snSV) |
|---|
| 156 | REAL,ALLOCATABLE,SAVE :: G2snSV(:,:) ! Snow Sphericity/Size |
|---|
| 157 | !$OMP THREADPRIVATE(G2snSV) |
|---|
| 158 | REAL,ALLOCATABLE,SAVE :: dzsnSV(:,:) ! Snow Layer Thickness |
|---|
| 159 | !$OMP THREADPRIVATE(dzsnSV) |
|---|
| 160 | REAL,ALLOCATABLE,SAVE :: agsnSV(:,:) ! Snow Age |
|---|
| 161 | !$OMP THREADPRIVATE(agsnSV) |
|---|
| 162 | REAL,ALLOCATABLE,SAVE :: DOPsnSV(:,:) ! Snow optical diameter [m] |
|---|
| 163 | !$OMP THREADPRIVATE(DOPsnSV) |
|---|
| 164 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: BufsSV ! Snow Buffer Layer |
|---|
| 165 | !$OMP THREADPRIVATE(BufsSV) |
|---|
| 166 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: rusnSV ! Surficial Water |
|---|
| 167 | !$OMP THREADPRIVATE(rusnSV) |
|---|
| 168 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: SWf_SV ! Normalized Decay |
|---|
| 169 | !$OMP THREADPRIVATE(SWf_SV) |
|---|
| 170 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: SWS_SV ! Surficial Water Status |
|---|
| 171 | !$OMP THREADPRIVATE(SWS_SV) |
|---|
| 172 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: HFraSV ! Frazil Thickness |
|---|
| 173 | !$OMP THREADPRIVATE(HFraSV) |
|---|
| 174 | |
|---|
| 175 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: zWE_SV ! Current Snow Thickness [mmWE] |
|---|
| 176 | !$OMP THREADPRIVATE(zWE_SV) |
|---|
| 177 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: zWEcSV ! Compacted Snow Thickness [mmWE] |
|---|
| 178 | !$OMP THREADPRIVATE(zWEcSV) |
|---|
| 179 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: wem_SV ! Only Melting [mmWE] |
|---|
| 180 | !$OMP THREADPRIVATE(wem_SV) |
|---|
| 181 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: wer_SV ! Refreezing [mmWE] |
|---|
| 182 | !$OMP THREADPRIVATE(wer_SV) |
|---|
| 183 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: wes_SV ! Sublimation [mmWE] |
|---|
| 184 | !$OMP THREADPRIVATE(wes_SV) |
|---|
| 185 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: zn4_SV ! snowheight change [mm] |
|---|
| 186 | !$OMP THREADPRIVATE(zn4_SV) |
|---|
| 187 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: zn5_SV ! snowheight change [mm] |
|---|
| 188 | !$OMP THREADPRIVATE(zn5_SV) |
|---|
| 189 | |
|---|
| 190 | ! +--SISVAT OUTPUT Variables |
|---|
| 191 | ! + ----------------------------- |
|---|
| 192 | |
|---|
| 193 | INTEGER,DIMENSION(nb_wri),SAVE :: no__SV ! OUTPUT file Unit Number |
|---|
| 194 | !$OMP THREADPRIVATE(no__SV) |
|---|
| 195 | INTEGER,DIMENSION(nb_wri),SAVE :: i___SV ! OUTPUT point i Coordinate |
|---|
| 196 | !$OMP THREADPRIVATE(i___SV) |
|---|
| 197 | INTEGER,DIMENSION(nb_wri),SAVE :: j___SV ! OUTPUT point j Coordinate |
|---|
| 198 | !$OMP THREADPRIVATE(j___SV) |
|---|
| 199 | INTEGER,DIMENSION(nb_wri),SAVE :: n___SV ! OUTPUT point n Coordinate |
|---|
| 200 | !$OMP THREADPRIVATE(n___SV) |
|---|
| 201 | INTEGER,DIMENSION(nb_wri),SAVE :: lwriSV ! OUTPUT point vec Index |
|---|
| 202 | !$OMP THREADPRIVATE(lwriSV) |
|---|
| 203 | ! |
|---|
| 204 | INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: ii__SV ! WORK point i Coordinate |
|---|
| 205 | !$OMP THREADPRIVATE(ii__SV) |
|---|
| 206 | INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: jj__SV ! WORK point j Coordinate |
|---|
| 207 | !$OMP THREADPRIVATE(jj__SV) |
|---|
| 208 | INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: nn__SV ! WORK point n Coordinate |
|---|
| 209 | !$OMP THREADPRIVATE(nn__SV) |
|---|
| 210 | |
|---|
| 211 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: IRu_SV ! UPward IR Flux (effective) |
|---|
| 212 | !$OMP THREADPRIVATE(IRu_SV) |
|---|
| 213 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: hSalSV ! Saltating Layer Height |
|---|
| 214 | !$OMP THREADPRIVATE(hSalSV) |
|---|
| 215 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: qSalSV ! Saltating Snow Concentration |
|---|
| 216 | !$OMP THREADPRIVATE(qSalSV) |
|---|
| 217 | REAL, DIMENSION(:),ALLOCATABLE,SAVE :: RnofSV ! RunOFF Intensity |
|---|
| 218 | !$OMP THREADPRIVATE(RnofSV) |
|---|
| 219 | REAL, DIMENSION(:,:),ALLOCATABLE,SAVE :: RuofSV ! RunOFF Intensity |
|---|
| 220 | !$OMP THREADPRIVATE(RuofSV) |
|---|
| 221 | |
|---|
| 222 | |
|---|
| 223 | |
|---|
| 224 | |
|---|
| 225 | CONTAINS |
|---|
| 226 | |
|---|
| 227 | |
|---|
| 228 | |
|---|
| 229 | SUBROUTINE INIT_VARxSV |
|---|
| 230 | IMPLICIT NONE |
|---|
| 231 | |
|---|
| 232 | INTEGER :: ikl |
|---|
| 233 | |
|---|
| 234 | ALLOCATE(LSmask(klonv)) ! Land-Sea Mask |
|---|
| 235 | ALLOCATE(isotSV(klonv)) ! Soil Type |
|---|
| 236 | ALLOCATE(iWaFSV(klonv)) ! Soil Drainage:(1,0)=(y,n) |
|---|
| 237 | ALLOCATE(coszSV(klonv)) ! Cosine of Sun zenithal Angle |
|---|
| 238 | ALLOCATE(sol_SV(klonv)) ! Downward Solar Radiation |
|---|
| 239 | ALLOCATE(IRd_SV(klonv)) ! Downward Longwave Radiation |
|---|
| 240 | |
|---|
| 241 | ALLOCATE(drr_SV(klonv)) ! Rain Intensity [kg/m2/s] |
|---|
| 242 | ALLOCATE(dsn_SV(klonv)) ! Snow Intensity [kg/m2/s] |
|---|
| 243 | ALLOCATE(dsnbSV(klonv)) ! Idem, from Drift [kg/m2] |
|---|
| 244 | ALLOCATE(esnbSV(klonv)) ! Idem, from Drift [kg/m2] |
|---|
| 245 | ALLOCATE(dbs_SV(klonv)) ! Drift Amount [kg/m2] |
|---|
| 246 | ALLOCATE(dbs_er(klonv)) ! Drift erosion [kg/m2] |
|---|
| 247 | ALLOCATE(dbs_ac(klonv)) ! Drift accumulation [kg/m2] |
|---|
| 248 | ALLOCATE(BrosSV(klonv)) |
|---|
| 249 | ALLOCATE(BG1sSV(klonv)) |
|---|
| 250 | ALLOCATE(BG2sSV(klonv)) |
|---|
| 251 | ALLOCATE(dz0_SV(klonv)) ! dz0(Sastrugi dh) [m] |
|---|
| 252 | |
|---|
| 253 | ALLOCATE(cld_SV(klonv)) ! Cloudiness (seen from SBL) |
|---|
| 254 | ALLOCATE(za__SV(klonv)) ! SBL Height |
|---|
| 255 | ALLOCATE(VV__SV(klonv)) !(SBL Top) Wind Velocity |
|---|
| 256 | ALLOCATE(VVs_SV(klonv)) |
|---|
| 257 | ALLOCATE(RRs_SV(klonv)) |
|---|
| 258 | ALLOCATE(DDs_SV(klonv)) |
|---|
| 259 | ALLOCATE(TaT_SV(klonv)) ! SBL Top Temperature |
|---|
| 260 | ALLOCATE(ExnrSV(klonv)) ! Exner Potential |
|---|
| 261 | ALLOCATE(dSdTSV(klonv)) ! Sensible Heat Flux T Derivat. |
|---|
| 262 | ALLOCATE(dLdTSV(klonv)) ! Latent Heat Flux T Derivat. |
|---|
| 263 | ALLOCATE(rhT_SV(klonv)) ! SBL Top Air Density |
|---|
| 264 | ALLOCATE(QaT_SV(klonv)) ! SBL Top Specific Humidity |
|---|
| 265 | ALLOCATE(QsT_SV(klonv)) ! surface Specific Humidity |
|---|
| 266 | ALLOCATE(dQa_SV(klonv)) ! SBL Flux Limitation of Qa |
|---|
| 267 | ALLOCATE(qsnoSV(klonv)) ! SBL Mean Snow Content |
|---|
| 268 | |
|---|
| 269 | ALLOCATE(alb0SV(klonv)) ! Soil Albedo |
|---|
| 270 | ALLOCATE(slopSV(klonv)) ! Snow/Ice/Soil-Water Surf. Slope |
|---|
| 271 | |
|---|
| 272 | |
|---|
| 273 | |
|---|
| 274 | ! +--SISVAT INPUT/OUTPUT Variables |
|---|
| 275 | ! + ----------------------------- |
|---|
| 276 | |
|---|
| 277 | ALLOCATE(isnoSV(klonv)) ! Nb of Ice/Snow Layers |
|---|
| 278 | ALLOCATE(ispiSV(klonv)) ! Uppermost superimposed ice |
|---|
| 279 | ALLOCATE(iiceSV(klonv)) ! Nb of Ice Layers |
|---|
| 280 | ALLOCATE(istoSV(klonv,0:nsno)) ! Snow Layer History |
|---|
| 281 | |
|---|
| 282 | ALLOCATE(alb_SV(klonv)) ! Surface-Canopy Albedo |
|---|
| 283 | ALLOCATE(emi_SV(klonv)) ! Surface-Canopy Emissivity |
|---|
| 284 | ALLOCATE(IRs_SV(klonv)) ! Soil IR Flux |
|---|
| 285 | ALLOCATE(LMO_SV(klonv)) ! Monin-Obukhov Scale |
|---|
| 286 | ALLOCATE(us__SV(klonv)) ! Friction Velocity |
|---|
| 287 | ALLOCATE(uts_SV(klonv)) ! Temperature Turbulent Scale |
|---|
| 288 | ALLOCATE(cutsSV(klonv)) ! Temperature Turbulent Scale C. |
|---|
| 289 | ALLOCATE(uqs_SV(klonv)) ! Spec.Humid. Turbulent Scale |
|---|
| 290 | ALLOCATE(uss_SV(klonv)) ! Blow.Snow Turbulent Scale |
|---|
| 291 | ALLOCATE(usthSV(klonv)) ! Blowing Snow Erosion Thresh. |
|---|
| 292 | ALLOCATE(rCDmSV(klonv)) ! Square Root Contribut. Drag_m |
|---|
| 293 | ALLOCATE(rCDhSV(klonv)) ! Square Root Contribut. Drag_h |
|---|
| 294 | ALLOCATE(Z0m_SV(klonv)) ! Momentum Roughness Length |
|---|
| 295 | ALLOCATE(Z0mmSV(klonv)) ! z0(Momentum, Time Mean) [m] |
|---|
| 296 | ALLOCATE(Z0mnSV(klonv)) ! z0(Momentum, instanta.) [m] |
|---|
| 297 | ALLOCATE(Z0roSV(klonv)) ! Subgrid Topo Roughness Length |
|---|
| 298 | ALLOCATE(Z0SaSV(klonv)) ! z0(Sastrugi h) [m] |
|---|
| 299 | ALLOCATE(Z0e_SV(klonv)) ! z0(Snow eroded) [m] |
|---|
| 300 | ALLOCATE(Z0emSV(klonv)) ! z0(Snow eroded, Time Mean) [m] |
|---|
| 301 | ALLOCATE(Z0enSV(klonv)) ! z0(Snow eroded, instanta.) [m] |
|---|
| 302 | ALLOCATE(Z0h_SV(klonv)) ! Heat Roughness Length |
|---|
| 303 | ALLOCATE(Z0hmSV(klonv)) ! z0(Heat, Time Mean) [m] |
|---|
| 304 | ALLOCATE(Z0hnSV(klonv)) ! z0(Heat, instanta.) [m] |
|---|
| 305 | |
|---|
| 306 | |
|---|
| 307 | ALLOCATE(TsisSV(klonv,-nsol:nsno)) ! Snow/Ice/Soil-Water Temperature |
|---|
| 308 | ALLOCATE(ro__SV(klonv,-nsol:nsno)) ! Snow/Ice/Soil-Water VolumicMass |
|---|
| 309 | ALLOCATE(eta_SV(klonv,-nsol:nsno)) ! Snow/Ice/Soil Water Content |
|---|
| 310 | ALLOCATE(G1snSV(klonv, 0:nsno)) ! Snow Dendricity/Sphericity |
|---|
| 311 | ALLOCATE(G2snSV(klonv, 0:nsno)) ! Snow Sphericity/Size |
|---|
| 312 | ALLOCATE(dzsnSV(klonv, 0:nsno)) ! Snow Layer Thickness |
|---|
| 313 | ALLOCATE(agsnSV(klonv, 0:nsno)) ! Snow Age |
|---|
| 314 | ALLOCATE(DOPsnSV(klonv, 0:nsno)) ! Snow Optical diameter |
|---|
| 315 | ALLOCATE(BufsSV(klonv)) ! Snow Buffer Layer |
|---|
| 316 | ALLOCATE(rusnSV(klonv)) ! Surficial Water |
|---|
| 317 | ALLOCATE(SWf_SV(klonv)) ! Normalized Decay |
|---|
| 318 | ALLOCATE(SWS_SV(klonv)) ! Surficial Water Status |
|---|
| 319 | ALLOCATE(HFraSV(klonv)) ! Frazil Thickness |
|---|
| 320 | |
|---|
| 321 | ALLOCATE(zWE_SV(klonv)) ! Current Snow Thickness [mmWE] |
|---|
| 322 | ALLOCATE(zWEcSV(klonv)) ! Compacted Snow Thickness [mmWE] |
|---|
| 323 | ALLOCATE(wem_SV(klonv)) ! Only Melting [mmWE] |
|---|
| 324 | ALLOCATE(wer_SV(klonv)) ! Refreezing [mmWE] |
|---|
| 325 | ALLOCATE(wes_SV(klonv)) ! Sublimation [mmWE] |
|---|
| 326 | ALLOCATE(zn4_SV(klonv)) ! snow height change [mm] |
|---|
| 327 | ALLOCATE(zn5_SV(klonv)) ! snow height change [mm] |
|---|
| 328 | |
|---|
| 329 | ! +--SISVAT OUTPUT Variables |
|---|
| 330 | ! + ----------------------------- |
|---|
| 331 | |
|---|
| 332 | ! |
|---|
| 333 | ALLOCATE(ii__SV(klonv)) ! WORK point i Coordinate |
|---|
| 334 | ALLOCATE(jj__SV(klonv)) ! WORK point j Coordinate |
|---|
| 335 | ALLOCATE(nn__SV(klonv)) ! WORK point n Coordinate |
|---|
| 336 | |
|---|
| 337 | ALLOCATE(IRu_SV(klonv)) ! UPward IR Flux (effective) |
|---|
| 338 | ALLOCATE(hSalSV(klonv)) ! Saltating Layer Height |
|---|
| 339 | ALLOCATE(qSalSV(klonv)) ! Saltating Snow Concentration |
|---|
| 340 | ALLOCATE(RnofSV(klonv)) ! RunOFF Intensity |
|---|
| 341 | ALLOCATE(RuofSV(klonv,4)) ! RunOFF Intensity |
|---|
| 342 | |
|---|
| 343 | |
|---|
| 344 | DO ikl=1,klonv |
|---|
| 345 | |
|---|
| 346 | |
|---|
| 347 | isnoSV(ikl) =0. |
|---|
| 348 | ispiSV(ikl) =0. |
|---|
| 349 | iiceSV(ikl) =0. |
|---|
| 350 | istoSV(ikl,:)=0. |
|---|
| 351 | alb_SV(ikl) =0. |
|---|
| 352 | emi_SV(ikl) =0. |
|---|
| 353 | IRs_SV(ikl) =0. |
|---|
| 354 | LMO_SV(ikl) =0. |
|---|
| 355 | us__SV(ikl) =0. |
|---|
| 356 | uts_SV(ikl) =0. |
|---|
| 357 | cutsSV(ikl) =0. |
|---|
| 358 | uqs_SV(ikl) =0. |
|---|
| 359 | uss_SV(ikl) =0. |
|---|
| 360 | usthSV(ikl) =0. |
|---|
| 361 | rCDmSV(ikl) =0. |
|---|
| 362 | rCDhSV(ikl) =0. |
|---|
| 363 | Z0m_SV(ikl) =0. |
|---|
| 364 | Z0mmSV(ikl) =0. |
|---|
| 365 | Z0mnSV(ikl) =0. |
|---|
| 366 | Z0roSV(ikl) =0. |
|---|
| 367 | Z0SaSV(ikl) =0. |
|---|
| 368 | Z0e_SV(ikl) =0. |
|---|
| 369 | Z0emSV(ikl) =0. |
|---|
| 370 | Z0enSV(ikl) =0. |
|---|
| 371 | Z0h_SV(ikl) =0. |
|---|
| 372 | Z0hmSV(ikl) =0. |
|---|
| 373 | Z0hnSV(ikl) =0. |
|---|
| 374 | |
|---|
| 375 | |
|---|
| 376 | TsisSV(ikl,:) =0. |
|---|
| 377 | ro__SV(ikl,:) =0. |
|---|
| 378 | eta_SV(ikl,:) =0. |
|---|
| 379 | G1snSV(ikl,:) =0. |
|---|
| 380 | G2snSV(ikl,:) =0. |
|---|
| 381 | dzsnSV(ikl,:) =0. |
|---|
| 382 | agsnSV(ikl,:) =0. |
|---|
| 383 | DOPsnSV(ikl,:) =0. |
|---|
| 384 | BufsSV(ikl) =0. |
|---|
| 385 | rusnSV(ikl) =0. |
|---|
| 386 | SWf_SV(ikl) =0. |
|---|
| 387 | SWS_SV(ikl) =0. |
|---|
| 388 | HFraSV(ikl) =0. |
|---|
| 389 | |
|---|
| 390 | zWE_SV(ikl) =0. |
|---|
| 391 | zWEcSV(ikl) =0. |
|---|
| 392 | wem_SV(ikl) =0. |
|---|
| 393 | wer_SV(ikl) =0. |
|---|
| 394 | wes_SV(ikl) =0. |
|---|
| 395 | zn4_SV(ikl) =0. |
|---|
| 396 | zn5_SV(ikl) =0. |
|---|
| 397 | |
|---|
| 398 | |
|---|
| 399 | ii__SV(ikl) =0. |
|---|
| 400 | jj__SV(ikl) =0. |
|---|
| 401 | nn__SV(ikl) =0. |
|---|
| 402 | |
|---|
| 403 | IRu_SV(ikl) =0. |
|---|
| 404 | hSalSV(ikl) =0. |
|---|
| 405 | qSalSV(ikl) =0. |
|---|
| 406 | RnofSV(ikl) =0. |
|---|
| 407 | RuofSV(ikl,:) =0. |
|---|
| 408 | |
|---|
| 409 | |
|---|
| 410 | |
|---|
| 411 | |
|---|
| 412 | END DO |
|---|
| 413 | END SUBROUTINE INIT_VARxSV |
|---|
| 414 | |
|---|
| 415 | |
|---|
| 416 | |
|---|
| 417 | END MODULE VARxSV |
|---|
