Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

VARxSV.F90 @ 5209

Last change on this file since 5209 was 3331, checked in by acozic, 6 years ago
Add modification for isotopes
Property svn:executable set to ``*
File size: 21.2 KB

Rev	Line
[3331]	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

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

Download in other formats:

Original Format