Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

VARxSV.F90 @ 3792

Last change on this file since 3792 was 3792, checked in by evignon, 4 years ago
Ajout de INLANDSIS, nouvelle interface entre LMDZ et la neige de SISVAT Etienne, 04/01/2021
File size: 20.7 KB

Rev	Line
[3792]	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 :: dQa_SV ! SBL Flux Limitation of Qa
	70	!$OMP THREADPRIVATE(dQa_SV)
	71	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: qsnoSV ! SBL Mean Snow Content
	72	!$OMP THREADPRIVATE(qsnoSV)
	73
	74	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: alb0SV ! Soil Albedo
	75	!$OMP THREADPRIVATE(alb0SV)
	76	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: slopSV ! Snow/Ice/Soil-Water Surf. Slope
	77	!$OMP THREADPRIVATE(slopSV)
	78
	79
	80	REAL,SAVE :: zSBLSV ! SBL Height (Initial Value)
	81	!$OMP THREADPRIVATE(zSBLSV)
	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, DIMENSION(:),ALLOCATABLE,SAVE :: BufsSV ! Snow Buffer Layer
	163	!$OMP THREADPRIVATE(BufsSV)
	164	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: rusnSV ! Surficial Water
	165	!$OMP THREADPRIVATE(rusnSV)
	166	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: SWf_SV ! Normalized Decay
	167	!$OMP THREADPRIVATE(SWf_SV)
	168	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: SWS_SV ! Surficial Water Status
	169	!$OMP THREADPRIVATE(SWS_SV)
	170	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: HFraSV ! Frazil Thickness
	171	!$OMP THREADPRIVATE(HFraSV)
	172
	173	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: zWE_SV ! Current Snow Thickness [mmWE]
	174	!$OMP THREADPRIVATE(zWE_SV)
	175	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: zWEcSV ! Compacted Snow Thickness [mmWE]
	176	!$OMP THREADPRIVATE(zWEcSV)
	177	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: wem_SV ! Only Melting [mmWE]
	178	!$OMP THREADPRIVATE(wem_SV)
	179	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: wer_SV ! Refreezing [mmWE]
	180	!$OMP THREADPRIVATE(wer_SV)
	181	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: wes_SV ! Sublimation [mmWE]
	182	!$OMP THREADPRIVATE(wes_SV)
	183	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: zn4_SV ! snowheight change [mm]
	184	!$OMP THREADPRIVATE(zn4_SV)
	185	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: zn5_SV ! snowheight change [mm]
	186	!$OMP THREADPRIVATE(zn5_SV)
	187
	188	! +--SISVAT OUTPUT Variables
	189	! + -----------------------------
	190
	191	INTEGER,DIMENSION(nb_wri),SAVE :: no__SV ! OUTPUT file Unit Number
	192	!$OMP THREADPRIVATE(no__SV)
	193	INTEGER,DIMENSION(nb_wri),SAVE :: i___SV ! OUTPUT point i Coordinate
	194	!$OMP THREADPRIVATE(i___SV)
	195	INTEGER,DIMENSION(nb_wri),SAVE :: j___SV ! OUTPUT point j Coordinate
	196	!$OMP THREADPRIVATE(j___SV)
	197	INTEGER,DIMENSION(nb_wri),SAVE :: n___SV ! OUTPUT point n Coordinate
	198	!$OMP THREADPRIVATE(n___SV)
	199	INTEGER,DIMENSION(nb_wri),SAVE :: lwriSV ! OUTPUT point vec Index
	200	!$OMP THREADPRIVATE(lwriSV)
	201	!
	202	INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: ii__SV ! WORK point i Coordinate
	203	!$OMP THREADPRIVATE(ii__SV)
	204	INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: jj__SV ! WORK point j Coordinate
	205	!$OMP THREADPRIVATE(jj__SV)
	206	INTEGER, DIMENSION(:),ALLOCATABLE,SAVE :: nn__SV ! WORK point n Coordinate
	207	!$OMP THREADPRIVATE(nn__SV)
	208
	209	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: IRu_SV ! UPward IR Flux (effective)
	210	!$OMP THREADPRIVATE(IRu_SV)
	211	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: hSalSV ! Saltating Layer Height
	212	!$OMP THREADPRIVATE(hSalSV)
	213	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: qSalSV ! Saltating Snow Concentration
	214	!$OMP THREADPRIVATE(qSalSV)
	215	REAL, DIMENSION(:),ALLOCATABLE,SAVE :: RnofSV ! RunOFF Intensity
	216	!$OMP THREADPRIVATE(RnofSV)
	217	REAL, DIMENSION(:,:),ALLOCATABLE,SAVE :: RuofSV ! RunOFF Intensity
	218	!$OMP THREADPRIVATE(RuofSV)
	219
	220
	221
	222
	223	CONTAINS
	224
	225
	226
	227	SUBROUTINE INIT_VARxSV
	228	IMPLICIT NONE
	229
	230	INTEGER :: ikl
	231
	232	ALLOCATE(LSmask(klonv)) ! Land-Sea Mask
	233	ALLOCATE(isotSV(klonv)) ! Soil Type
	234	ALLOCATE(iWaFSV(klonv)) ! Soil Drainage:(1,0)=(y,n)
	235	ALLOCATE(coszSV(klonv)) ! Cosine of Sun zenithal Angle
	236	ALLOCATE(sol_SV(klonv)) ! Downward Solar Radiation
	237	ALLOCATE(IRd_SV(klonv)) ! Downward Longwave Radiation
	238
	239	ALLOCATE(drr_SV(klonv)) ! Rain Intensity [kg/m2/s]
	240	ALLOCATE(dsn_SV(klonv)) ! Snow Intensity [kg/m2/s]
	241	ALLOCATE(dsnbSV(klonv)) ! Idem, from Drift [kg/m2]
	242	ALLOCATE(esnbSV(klonv)) ! Idem, from Drift [kg/m2]
	243	ALLOCATE(dbs_SV(klonv)) ! Drift Amount [kg/m2]
	244	ALLOCATE(dbs_er(klonv)) ! Drift erosion [kg/m2]
	245	ALLOCATE(dbs_ac(klonv)) ! Drift accumulation [kg/m2]
	246	ALLOCATE(BrosSV(klonv))
	247	ALLOCATE(BG1sSV(klonv))
	248	ALLOCATE(BG2sSV(klonv))
	249	ALLOCATE(dz0_SV(klonv)) ! dz0(Sastrugi dh) [m]
	250
	251	ALLOCATE(cld_SV(klonv)) ! Cloudiness (seen from SBL)
	252	ALLOCATE(za__SV(klonv)) ! SBL Height
	253	ALLOCATE(VV__SV(klonv)) !(SBL Top) Wind Velocity
	254	ALLOCATE(VVs_SV(klonv))
	255	ALLOCATE(RRs_SV(klonv))
	256	ALLOCATE(DDs_SV(klonv))
	257	ALLOCATE(TaT_SV(klonv)) ! SBL Top Temperature
	258	ALLOCATE(ExnrSV(klonv)) ! Exner Potential
	259	ALLOCATE(dSdTSV(klonv)) ! Sensible Heat Flux T Derivat.
	260	ALLOCATE(dLdTSV(klonv)) ! Latent Heat Flux T Derivat.
	261	ALLOCATE(rhT_SV(klonv)) ! SBL Top Air Density
	262	ALLOCATE(QaT_SV(klonv)) ! SBL Top Specific Humidity
	263	ALLOCATE(dQa_SV(klonv)) ! SBL Flux Limitation of Qa
	264	ALLOCATE(qsnoSV(klonv)) ! SBL Mean Snow Content
	265
	266	ALLOCATE(alb0SV(klonv)) ! Soil Albedo
	267	ALLOCATE(slopSV(klonv)) ! Snow/Ice/Soil-Water Surf. Slope
	268
	269
	270
	271	! +--SISVAT INPUT/OUTPUT Variables
	272	! + -----------------------------
	273
	274	ALLOCATE(isnoSV(klonv)) ! Nb of Ice/Snow Layers
	275	ALLOCATE(ispiSV(klonv)) ! Uppermost superimposed ice
	276	ALLOCATE(iiceSV(klonv)) ! Nb of Ice Layers
	277	ALLOCATE(istoSV(klonv,0:nsno)) ! Snow Layer History
	278
	279	ALLOCATE(alb_SV(klonv)) ! Surface-Canopy Albedo
	280	ALLOCATE(emi_SV(klonv)) ! Surface-Canopy Emissivity
	281	ALLOCATE(IRs_SV(klonv)) ! Soil IR Flux
	282	ALLOCATE(LMO_SV(klonv)) ! Monin-Obukhov Scale
	283	ALLOCATE(us__SV(klonv)) ! Friction Velocity
	284	ALLOCATE(uts_SV(klonv)) ! Temperature Turbulent Scale
	285	ALLOCATE(cutsSV(klonv)) ! Temperature Turbulent Scale C.
	286	ALLOCATE(uqs_SV(klonv)) ! Spec.Humid. Turbulent Scale
	287	ALLOCATE(uss_SV(klonv)) ! Blow.Snow Turbulent Scale
	288	ALLOCATE(usthSV(klonv)) ! Blowing Snow Erosion Thresh.
	289	ALLOCATE(rCDmSV(klonv)) ! Square Root Contribut. Drag_m
	290	ALLOCATE(rCDhSV(klonv)) ! Square Root Contribut. Drag_h
	291	ALLOCATE(Z0m_SV(klonv)) ! Momentum Roughness Length
	292	ALLOCATE(Z0mmSV(klonv)) ! z0(Momentum, Time Mean) [m]
	293	ALLOCATE(Z0mnSV(klonv)) ! z0(Momentum, instanta.) [m]
	294	ALLOCATE(Z0roSV(klonv)) ! Subgrid Topo Roughness Length
	295	ALLOCATE(Z0SaSV(klonv)) ! z0(Sastrugi h) [m]
	296	ALLOCATE(Z0e_SV(klonv)) ! z0(Snow eroded) [m]
	297	ALLOCATE(Z0emSV(klonv)) ! z0(Snow eroded, Time Mean) [m]
	298	ALLOCATE(Z0enSV(klonv)) ! z0(Snow eroded, instanta.) [m]
	299	ALLOCATE(Z0h_SV(klonv)) ! Heat Roughness Length
	300	ALLOCATE(Z0hmSV(klonv)) ! z0(Heat, Time Mean) [m]
	301	ALLOCATE(Z0hnSV(klonv)) ! z0(Heat, instanta.) [m]
	302
	303
	304	ALLOCATE(TsisSV(klonv,-nsol:nsno)) ! Snow/Ice/Soil-Water Temperature
	305	ALLOCATE(ro__SV(klonv,-nsol:nsno)) ! Snow/Ice/Soil-Water VolumicMass
	306	ALLOCATE(eta_SV(klonv,-nsol:nsno)) ! Snow/Ice/Soil Water Content
	307	ALLOCATE(G1snSV(klonv, 0:nsno)) ! Snow Dendricity/Sphericity
	308	ALLOCATE(G2snSV(klonv, 0:nsno)) ! Snow Sphericity/Size
	309	ALLOCATE(dzsnSV(klonv, 0:nsno)) ! Snow Layer Thickness
	310	ALLOCATE(agsnSV(klonv, 0:nsno)) ! Snow Age
	311	ALLOCATE(BufsSV(klonv)) ! Snow Buffer Layer
	312	ALLOCATE(rusnSV(klonv)) ! Surficial Water
	313	ALLOCATE(SWf_SV(klonv)) ! Normalized Decay
	314	ALLOCATE(SWS_SV(klonv)) ! Surficial Water Status
	315	ALLOCATE(HFraSV(klonv)) ! Frazil Thickness
	316
	317	ALLOCATE(zWE_SV(klonv)) ! Current Snow Thickness [mmWE]
	318	ALLOCATE(zWEcSV(klonv)) ! Compacted Snow Thickness [mmWE]
	319	ALLOCATE(wem_SV(klonv)) ! Only Melting [mmWE]
	320	ALLOCATE(wer_SV(klonv)) ! Refreezing [mmWE]
	321	ALLOCATE(wes_SV(klonv)) ! Sublimation [mmWE]
	322	ALLOCATE(zn4_SV(klonv)) ! snow height change [mm]
	323	ALLOCATE(zn5_SV(klonv)) ! snow height change [mm]
	324
	325	! +--SISVAT OUTPUT Variables
	326	! + -----------------------------
	327
	328	!
	329	ALLOCATE(ii__SV(klonv)) ! WORK point i Coordinate
	330	ALLOCATE(jj__SV(klonv)) ! WORK point j Coordinate
	331	ALLOCATE(nn__SV(klonv)) ! WORK point n Coordinate
	332
	333	ALLOCATE(IRu_SV(klonv)) ! UPward IR Flux (effective)
	334	ALLOCATE(hSalSV(klonv)) ! Saltating Layer Height
	335	ALLOCATE(qSalSV(klonv)) ! Saltating Snow Concentration
	336	ALLOCATE(RnofSV(klonv)) ! RunOFF Intensity
	337	ALLOCATE(RuofSV(klonv,4)) ! RunOFF Intensity
	338
	339
	340	DO ikl=1,klonv
	341	LSmask(ikl) = 0
	342	isotSV(ikl) = 0
	343	iWaFSV(ikl) = 0
	344	isnoSV(ikl) = 0
	345	ispiSV(ikl) = 0
	346	iiceSV(ikl) = 0
	347	istoSV(ikl,:) = 0
	348	ii__SV(ikl) = 0
	349	jj__SV(ikl) = 0
	350	nn__SV(ikl) = 0
	351	END DO
	352	END SUBROUTINE INIT_VARxSV
	353
	354
	355
	356	END MODULE VARxSV

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

Download in other formats:

Original Format