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