Context Navigation

← Previous Revision
Latest Revision
Next Revision →
Normal
Revision Log

PHY_SISVAT_INP.f90 @ 5456

Last change on this file since 5456 was 2089, checked in by Laurent Fairhead, 10 years ago

Inclusion de la physique de MAR

Integration of MAR physics

File size: 27.4 KB

Rev	Line
[2089]	1	subroutine PHY_SISVAT_INP(FlagSV_SBC,TsHOST,kcolq)
	2
	3	!------------------------------------------------------------------------------+
	4	! Mon 1-Jul-2013 MAR \|
	5	! SubRoutine PHY_SISVAT_INP initializes \|
	6	! SISVAT (Soil \|
	7	! Ice \|
	8	! Snow \|
	9	! Vegetation \|
	10	! Atmosphere \|
	11	! Transfer Scheme) \|
	12	! \|
	13	! USING DATA \|
	14	! \|
	15	! Here DATA are academic \|
	16	! \|
	17	! version 3.p.4.1 created by H. Gallee, Wed 1-May-2013 \|
	18	! Last Modification by H. Gallee, Mon 1-Jul-2013 \|
	19	! \|
	20	!------------------------------------------------------------------------------+
	21
	22
	23	! Global Variables
	24	! =================
	25
	26	use Mod_Real
	27	use Mod_PHY____grd
	28	use Mod_PHY____dat
	29	use Mod_SISVAT_dat
	30	use Mod_SISVAT_grd
	31	use Mod_SISVAT_kkl
	32	use Mod_SISVAT_gpt
	33
	34
	35	IMPLICIT NONE
	36
	37
	38	logical :: FlagSV_SBC ! DATA Surface Flag
	39	integer :: kcolq !
	40	real(kind=real8), dimension(kcolq) :: TsHOST ! (Initial) Surface Temperature [K]
	41
	42
	43
	44
	45	! Local Variables
	46	! =================
	47
	48	! FOR INPUT
	49	! ---------
	50
	51	! integer ,ALLOCATABLE ,dimension(:,:,:) :: Mask___xyn ! Land(1)-Sea(0) Mask (Mosaic) [-]
	52	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: Frac___xyn ! Grid Cell Fraction (Mosaic) [-]
	53	! integer ,ALLOCATABLE ,dimension(:,:,:) :: iVgT___xyn ! Vegetat. Type (Mosaic) [-]
	54	! integer ,ALLOCATABLE ,dimension(:,:,:) :: iSoT___xyn ! Soil Type (Mosaic) [-]
	55	! integer ,ALLOCATABLE ,dimension(:,:,:) :: iDra___xyn ! Soil Drainage Switch (Mosaic) [-]
	56	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: Albs___xyn ! Dry Soil Albedo (Mosaic) [-]
	57	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: LAI____xyn ! Leaf Area Index (Mosaic) [m2/m2]
	58	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: GLF____xyn ! Green Leaf Fraction (Mosaic) [-]
	59
	60
	61
	62	! FOR INPUT/OUTPUT (RESTART)
	63	! ----------------
	64
	65	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: FSbz___xyn ! Fallen Snow Buffer (Mosaic) [mm w.e.]
	66	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: Tveg___xyn ! Vegetation Temperature (Mosaic) [K]
	67	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: Sveg___xyn ! Canopy intercepted Snow (Mosaic) [mm w.e.]
	68	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: Rveg___xyn ! Canopy intercepted Rain (Mosaic) [mm w.e.]
	69	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: pveg___xyn ! Vegetation Water Potential (Mosaic) [m]
	70
	71	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: z0_____xyn ! Roughness Length Momentum (Mosaic) [m]
	72	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: z0ro___xyn ! Orographic Roughness Length (Mosaic) [m]
	73	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: r0_____xyn ! Roughness Length Scalars (Mosaic) [m]
	74	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: uss____xyn ! Blown Snow Flux Turbulent Scale [kg/kg m/s]
	75	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: uds____xyn ! Blown Dust Flux Turbulent Scale [kg/kg m/s]
	76	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: BSbz___xyn ! Eroded Snow Buffer (Mosaic) [mm w.e.]
	77	! #AE real(kind=real8),ALLOCATABLE ,dimension(:,:) :: usth___xyn ! Friction Velocity Saltation Threshold [m/s]
	78
	79	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: Alb____xyn ! Surface Albedo (Mosaic) [-]
	80	! integer ,ALLOCATABLE ,dimension(:,:,:) :: SNnb___xyn ! Nb of Snow Layers (Mosaic) [-]
	81	! integer ,ALLOCATABLE ,dimension(:,:,:) :: SInb___xyn ! Nb of S.Ice Layers (Mosaic) [-]
	82	! integer ,ALLOCATABLE ,dimension(:,:,:) :: ICnb___xyn ! Nb of Ice Layers (Mosaic) [-]
	83
	84	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: SnWE___xyn ! Snow Pack Thickness (Mosaic) [mm w.e.]
	85
	86	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: wes____xyn ! Sublimation Budget (Mosaic) [mm w.e.]
	87	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: wer____xyn ! Melting Budget (Mosaic) [mm w.e.]
	88	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: wem____xyn ! Refreezing Budget (Mosaic) [mm w.e.]
	89	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: wee____xyn ! Evapotranspiration Budget (Mosaic) [mm w.e.]
	90
	91	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:,:):: Tsis__dxyn ! Soil-Ice-Snow Temperature (Mosaic) [K]
	92	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:,:):: Qsis__dxyn ! Soil-Ice-Snow Humidity (Mosaic) [kg/kg]
	93	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:,:):: SNdz__dxyn ! Snow Layer Thickness (Mosaic) [m]
	94	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:,:):: SNro__dxyn ! Snow Layer Density (Mosaic) [kg/m3]
	95	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:,:):: SNag__dxyn ! Snow Layer Age (Mosaic) [day]
	96	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:,:):: SNg1__dxyn ! Snow Layer Property (d./s.)(Mosaic) [-]
	97	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:,:):: SNg2__dxyn ! Snow Layer Property (s./s.)(Mosaic) [-]
	98	! integer ,ALLOCATABLE ,dimension(:,:,:,:):: SNhi__dxyn ! Snow Layer History (Mosaic) [-]
	99
	100	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: SWz____xyn ! Surficial Water Mass (Mosaic) [mm w.e.]
	101	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: SWx____xyn ! Surficial Water Status (Mosaic) [0,1=m,f]
	102	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: SWHy___xyn ! Surficial Water: Normaliz.Decay (Mosaic) [-]
	103	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: ROdt___xyn ! Run-Off Intensity (Mosaic) [mm w.e./s]
	104	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: ROF____xyn ! Cumulative Run-Off (Mosaic) [mm w.e.]
	105
	106
	107	! A-O Coupling
	108	! ------------
	109
	110	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: Alb_AO_xyn ! Ocean COUPLING (Surface Albedo) n=2 sea-ice [-]
	111	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: s_T_AO_xyn ! Ocean COUPLING (SST) n=1 opn-watr [K]
	112	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: dSdTAO_xyn ! Ocean COUPLING (d(SH Flux) / dT) [W/m2/K]
	113	! real(kind=real8),ALLOCATABLE ,dimension(:,:,:) :: dLdTAO_xyn ! Ocean COUPLING (d(LH Flux) / dT) [W/m2/K]
	114
	115
	116	! Dummy
	117	! -----
	118
	119	integer :: i ,j ,n ,l ,nm ,ikp
	120
	121
	122
	123
	124
	125
	126	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	127	! !
	128	! NON ACADEMIC INITIALIZATION !
	129	! !
	130	IF (FlagSV_SBC) THEN
	131	! !
	132	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	133
	134
	135
	136
	137	! +++++++++++++++++++++++++++++++++++++++++++++++++++++
	138	! INITIALIZATION of a non-1D case, by using INPUT Files
	139	! (HERE FROM NESTOR)
	140
	141	! Martin CONTROL
	142	Print*,'mxpp=',mxpp
	143	Print*,'mypp=',mypp
	144	Print*,'ixp1=',ixp1
	145	Print*,'jyp1=',jyp1
	146
	147	IF (mxpp-ixp1.GT.0 .OR. mypp-jyp1.GT.0) THEN
	148	! +++++++++++++++++++++++++++++++++++++++++++++++++++++
	149
	150
	151
	152
	153	! ================================
	154	! FIRST INITIALIZATION (BEGIN)
	155
	156	IF (it_EXP .LE. 1) THEN
	157	! ================================
	158
	159
	160	! ================================
	161	END IF ! (it_EXP .LE. 1)
	162
	163	! FIRST INITIALIZATION (END)
	164	! ================================
	165
	166
	167
	168
	169	! Update Sea-Ice Fraction
	170	! ==========================
	171
	172	! #IP IF ( inpSBC ) THEN
	173
	174	! ******
	175	! #IP call INIsic(newsicSI ) ! ATTENTION 2 premiers Arguments supprimés
	176	! ******
	177
	178	! #IP END IF
	179
	180
	181
	182
	183	! Update LAI and/or Green Leaf Fraction
	184	! =====================================
	185
	186	! #GP glfFIX = .true.
	187	! #CM IF (vegmod .AND. inpSBC .AND. .NOT. glfFIX) THEN
	188
	189	! ******
	190	! #CM call INIglf(newglfSIS) ! ATTENTION 2 premiers Arguments supprimés
	191	! ******
	192
	193	! #CM END IF
	194
	195
	196
	197
	198
	199
	200	! +++++++++++++++++++++++++++++++++++++++++++++++++++++
	201	! INITIALIZATION of a 1D case, by prescribing INPUT
	202
	203	ELSE
	204	! +++++++++++++++++++++++++++++++++++++++++++++++++++++
	205
	206
	207
	208
	209	! ================================
	210	! FIRST INITIALIZATION (BEGIN)
	211
	212	IF (it_EXP .LE. 1) THEN
	213	! ================================
	214
	215
	216
	217	! ALLOCATION
	218	! ==========
	219
	220	! FOR INPUT
	221	! ---------
	222
	223	! allocate ( Mask___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Land(1)-Sea(0) Mask (Mosaic) [-]
	224	! allocate ( Frac___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Grid Cell Fraction (Mosaic) [-]
	225	! allocate ( iVgT___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Vegetat. Type (Mosaic) [-]
	226	! allocate ( iSoT___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Soil Type (Mosaic) [-]
	227	! allocate ( iDra___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Soil Drainage Switch (Mosaic) [-]
	228	! allocate ( Albs___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Dry Soil Albedo (Mosaic) [-]
	229	! allocate ( LAI____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Leaf Area Index (Mosaic) [m2/m2]
	230	! allocate ( GLF____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Green Leaf Fraction (Mosaic) [-]
	231
	232
	233	! FOR INPUT/OUTPUT (RESTART)
	234	! ----------------
	235
	236	! allocate ( FSbz___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Fallen Snow Buffer (Mosaic) [mm w.e.]
	237	! allocate ( Tveg___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Vegetation Temperature (Mosaic) [K]
	238	! allocate ( Sveg___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Canopy intercepted Snow (Mosaic) [mm w.e.]
	239	! allocate ( Rveg___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Canopy intercepted Water (Mosaic) [mm w.e.]
	240	! allocate ( pveg___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Vegetation Water Potential (Mosaic) [m]
	241
	242	! allocate ( z0_____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Roughness Length Momentum (Mosaic) [m]
	243	! allocate ( z0ro___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! z0(Orographic) [m]
	244	! allocate ( r0_____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Roughness Length Scalars (Mosaic) [m]
	245	! allocate ( uss____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Blown Snow Flux Turbulent Scale [kg/kg m/s]
	246	! allocate ( BSbz___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Eroded Snow Buffer (Mosaic) [mm w.e.]
	247	! #AE allocate ( usth___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Friction Velocity Saltation Threshold [m/s]
	248
	249	! allocate ( Alb____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Surface Albedo (Mosaic) [-]
	250	! allocate ( SNnb___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Nb of Snow Layers (Mosaic) [-]
	251	! allocate ( SInb___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Nb of S.Ice Layers (Mosaic) [-]
	252	! allocate ( ICnb___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Nb of Ice Layers (Mosaic) [-]
	253
	254	! allocate ( SnWE___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Snow Pack Thickness (Mosaic) [mm w.e.]
	255
	256	! allocate ( wes____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Sublimation Budget (Mosaic) [mm w.e.]
	257	! allocate ( wer____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Melting Budget (Mosaic) [mm w.e.]
	258	! allocate ( wem____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Refreezing Budget (Mosaic) [mm w.e.]
	259	! allocate ( wee____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Evapotranspiration Budget (Mosaic) [mm w.e.]
	260
	261	! allocate ( Tsis__dxyn(ixp1:mxpp,jyp1:mypp,mwp,-nsoil:nsnow) )
	262	! ! Soil-Ice-Snow Temperature (Mosaic) [K]
	263	! allocate ( Qsis__dxyn(ixp1:mxpp,jyp1:mypp,mwp,-nsoil:nsnow) )
	264	! ! Soil-Ice-Snow Humidity (Mosaic) [kg/kg]
	265	! allocate ( SNdz__dxyn(ixp1:mxpp,jyp1:mypp,mwp,nsnow) ) ! Snow Layer Thickness (Mosaic) [m]
	266	! allocate ( SNro__dxyn(ixp1:mxpp,jyp1:mypp,mwp,nsnow) ) ! Snow Layer Density (Mosaic) [kg/m3]
	267	! allocate ( SNag__dxyn(ixp1:mxpp,jyp1:mypp,mwp,nsnow) ) ! Snow Layer Age (Mosaic) [day]
	268	! allocate ( SNg1__dxyn(ixp1:mxpp,jyp1:mypp,mwp,nsnow) ) ! Snow Layer Property (d./s.)(Mosaic) [-]
	269	! allocate ( SNg2__dxyn(ixp1:mxpp,jyp1:mypp,mwp,nsnow) ) ! Snow Layer Property (s./s.)(Mosaic) [-]
	270	! allocate ( SNhi__dxyn(ixp1:mxpp,jyp1:mypp,mwp,nsnow) ) ! Snow Layer History (Mosaic) [-]
	271
	272	! allocate ( SWz____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Surficial Water Mass (Mosaic) [mm w.e.]
	273	! allocate ( SWx____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Surficial Water Status (Mosaic) [0,1=m,f]
	274	! allocate ( SWHy___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Normalized Decay of Surficial Water Content [-]
	275	! allocate ( ROdt___xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Run-Off Intensity (Mosaic) [mm w.e./s]
	276	! allocate ( ROF____xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Cumulative Run-Off (Mosaic) [mm w.e.]
	277
	278
	279	! A-O Coupling
	280	! ------------
	281
	282	! allocate ( Alb_AO_xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Ocean COUPLING (Surface Albedo) n=2 sea-ice [-]
	283	! allocate ( Alb_AO_xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Ocean COUPLING (Surface Albedo) n=2 sea-ice [-]
	284	! allocate ( dSdTAO_xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Ocean COUPLING (d(SH Flux) / dT) [W/m2/K]
	285	! allocate ( dLdTAO_xyn(ixp1:mxpp,jyp1:mypp,mwp) ) ! Ocean COUPLING (d(LH Flux) / dT) [W/m2/K]
	286
	287
	288	! Martin CONTROL
	289	PRINT*,'On initialise des paramètres en dur dans PHY_SISVAT_INP'
	290
	291	! ------------------------------------------
	292	! From 2-D Horizontal to 1-D Horizontal Grid
	293	! ------------------------------------------
	294
	295	DO ikp=1,kcolp
	296	i = ii__AP(ikp)
	297	j = jj__AP(ikp)
	298
	299
	300	!MaskSV_gpt(ikp) = 1 ! Land(1)-Sea(0) Mask (Grid Cell) [-]
	301	! Martin : aquaplanète
	302	MaskSV_gpt(ikp) = 0 ! Land(1)-Sea(0) Mask (Grid Cell) [-]
	303
	304	DO n=1,mwp
	305	! Mask___xyn(i,j,n) = 1 ! Land(1)-Sea(0) Mask (Mosaic) [-]
	306	!LSmask (ikp,n) = 1 ! Land(1)-Sea(0) Mask (Mosaic) [-]
	307	LSmask (ikp,n) = 0 ! Land(1)-Sea(0) Mask (Mosaic) [-]
	308
	309	SlopSV (ikp,n) = 0.0 ! Assignation of Slope to a Particular Value
	310
	311	! Frac___xyn(i,j,n) = 1.0/float(mwp) ! Grid Cell Fraction (Mosaic) [-]
	312	FracSV (ikp,n) = 1.0/float(mwp) ! Grid Cell Fraction (Mosaic) [-]
	313
	314	! z0_____xyn(i,j,n) = 1.e-6 ! Roughness Length Momentum (Mosaic) [m]
	315	!z0__SV (ikp,n) = 1.e-6 ! Roughness Length Momentum (Mosaic) [m]
	316	z0__SV (ikp,n) = 1.e-3 ! Roughness Length Momentum (Mosaic) [m]
	317	! (Sensitive if non-zero Orographic Roughness)
	318	! r0_____xyn(i,j,n) = 1.e-6 ! Roughness Length Scalars (Mosaic) [m]
	319	!Z0h_SV (ikp,n) = 1.e-6 ! Roughness Length Scalars (Mosaic) [m]
	320	Z0h_SV (ikp,n) = 1.e-3 ! Roughness Length Scalars (Mosaic) [m]
	321
	322
	323	! Soil and Vegetation
	324	! -------------------
	325
	326	! iVgT___xyn(i,j,n) = 5 ! Vegetat. Type (Mosaic) [-] (Sensitive)
	327	iVgTSV (ikp,n) = 0 ! Vegetat. Type (Mosaic) [-] (Sensitive)
	328	! iSoT___xyn(i,j,n) = 5 ! Soil Type (Mosaic) [-]
	329	iSoTSV (ikp,n) = 0 ! Soil Type (Mosaic) [-]
	330	! iDra___xyn(i,j,n) = 1 ! Soil Drainage Switch (Mosaic) [-]
	331	iWaFSV (ikp,n) = 1 ! Soil Drainage Switch (Mosaic) [-]
	332	! Albs___xyn(i,j,n) = 0.15 ! Dry Soil Albedo (Mosaic) [-]
	333	alb0SV (ikp,n) = 0.15 ! Dry Soil Albedo (Mosaic) [-]
	334	! LAI0SV (ikp,n) = LAI____xyn(i,j,n) ! LAI
	335	! glf0SV (ikp,n) = GLF____xyn(i,j,n) ! Green Leaf Frac.
	336
	337
	338	! Snow Pack Initialization: NO Snow Pack
	339	! --------------------------------------
	340
	341	! SNnb___xyn(i,j,n) = 0 ! Nb of Snow Layers (Mosaic) [-] (Sensitive, maybe)
	342	isnoSV (ikp,n) = 0 ! Nb of Snow Layers (Mosaic) [-] (Sensitive, maybe)
	343	! SInb___xyn(i,j,n) = 0 ! Nb of S.Ice Layers (Mosaic) [-]
	344	ispiSV (ikp,n) = 0 ! Nb of S.Ice Layers (Mosaic) [-]
	345	! ICnb___xyn(i,j,n) = 0 ! Nb of Ice Layers (Mosaic) [-] (Sensitive, maybe)
	346	iiceSV (ikp,n) = 0 ! Nb of Ice Layers (Mosaic) [-] (Sensitive, maybe)
	347
	348	DO l=-nsol,0
	349	! Tsis__dxyn(i,j,n,l) = TsHOST(ikp) ! Soil-Ice-Snow Temperature (Mosaic) [K]
	350	! TsisSV (ikp,n,l) = Tsis__dxyn(i,j,n,l) !
	351	TsisSV (ikp,n,l) = TsHOST(ikp) !
	352	! Qsis__dxyn(i,j,n,l) = &! Humidity (Mosaic) [kg/kg]
	353	! & etadSV(iVgTSV(ikp,n)) ! (Saturation is assumed)
	354	! eta_SV (ikp,n,l) = Qsis__dxyn(i,j,n,l) !
	355	eta_SV (ikp,n,l) = &!
	356	& etadSV(iVgTSV(ikp,n)) !
	357	ENDDO
	358
	359	DO l=1,nsno
	360	! Tsis__dxyn(i,j,n,l) = 273.15 ! Soil-Ice-Snow Temperature (Mosaic) [K]
	361	! TsisSV (ikp,n,l) = Tsis__dxyn(i,j,n,l) !
	362	TsisSV (ikp,n,l) = 273.15 !
	363	! Qsis__dxyn(i,j,n,l) = 0.00 ! Soil-Ice-Snow Humidity (Mosaic) [kg/kg]
	364	! eta_SV (ikp,n,l) = Qsis__dxyn(i,j,n,l) !
	365	eta_SV (ikp,n,l) = 0.00 !
	366	! SNdz__dxyn(i,j,n,l) = 0.00 ! Snow Layer Thickness (Mosaic) [m]
	367	! dzsnSV (ikp,n,l) = SNdz__dxyn(i,j,n,l) !
	368	dzsnSV (ikp,n,l) = 0.00 !
	369	! SNro__dxyn(i,j,n,l) = 0.00 ! Snow Layer Density (Mosaic) [kg/m3]
	370	! ro__SV (ikp,n,l) = SNro__dxyn(i,j,n,l) !
	371	ro__SV (ikp,n,l) = 0.00 !
	372	! SNag__dxyn(i,j,n,l) = 0.00 ! Snow Layer Age (Mosaic) [day]
	373	! agsnSV (ikp,n,l) = SNag__dxyn(i,j,n,l) !
	374	agsnSV (ikp,n,l) = 0.00 !
	375	! SNg1__dxyn(i,j,n,l) = 0.00 ! Snow Layer Property (d./s.)(Mosaic) [-]
	376	! G1snSV (ikp,n,l) = SNg1__dxyn(i,j,n,l) !
	377	G1snSV (ikp,n,l) = 0.00 !
	378	! SNg2__dxyn(i,j,n,l) = 0.00 ! Snow Layer Property (s./s.)(Mosaic) [-]
	379	! G2snSV (ikp,n,l) = SNg2__dxyn(i,j,n,l) !
	380	G2snSV (ikp,n,l) = 0.00 !
	381	! SNhi__dxyn(i,j,n,l) = 0 ! Snow Layer History (Mosaic)
	382	! istoSV (ikp,n,l) = SNhi__dxyn(i,j,n,l) !
	383	istoSV (ikp,n,l) = 0 !
	384	ENDDO
	385
	386	ENDDO
	387	ENDDO
	388
	389
	390	! ================================
	391	END IF ! (it_EXP .LE. 1)
	392
	393	! FIRST INITIALIZATION (END)
	394	! ================================
	395
	396
	397
	398
	399	! OUTPUT
	400	! ======
	401
	402	! dSdTAO_xyn(i,j,n) = dSdTSV(ikp,n) ! Sens.H.Flux T-Der. (A-O Coupling) [W/m2/K]
	403	! dLdTAO_xyn(i,j,n) = dLdTSV(ikp,n) ! Latn.H.Flux T-Der. (A-O Coupling) [W/m2/K]
	404	!!#SN SnWE___xyn(i,j,n) = zWE_SV(ikp,n) ! Current *Thick.
	405	! ROdt___xyn(i,j,n) = RnofSV(ikp,n) ! Run OFF Intensity
	406
	407
	408
	409
	410	! RESTART
	411	! =======
	412
	413	!!#SN BufsSV(ikp,n) = FSbz___xyn(i,j,n) ! Snow Buffer Lay.
	414	!!#SN rusnSV(ikp,n) = SWz____xyn(i,j,n) ! Surficial Water
	415	!!#SN SWS_SV(ikp,n) = SWx____xyn(i,j,n) ! Surficial Water Status
	416	!!#SN SWf_SV(ikp,n) = SWHySV_xyn(i,j,n) ! Normalized Decay
	417
	418
	419
	420
	421	! +++++++++++++++++++++++++++++++++++++++++++++++++++++
	422	END IF
	423	! +++++++++++++++++++++++++++++++++++++++++++++++++++++
	424
	425
	426
	427
	428
	429
	430	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	431	! !
	432	! ACADEMIC INITIALIZATION !
	433	! !
	434	ELSE ! IF (.NOT.FlagSV_SBC) !
	435	! !
	436	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	437
	438
	439	! ==============================
	440	! FIRST INITIALIZATION (BEGIN)
	441	PRINT*, 'CONTROL PHY_SISVAT_INP'
	442	PRINT*,'kcolq=',kcolq
	443	PRINT*,'kcolp=',kcolp
	444	PRINT*,'kmwp=',mwp
	445
	446	IF (it_EXP .LE. 1) THEN
	447	! ==============================
	448
	449	DO ikp=1,kcolp
	450
	451	MaskSV_gpt(ikp) = 0 ! Land-Sea Mask
	452
	453	DO nm=1,mwp
	454
	455	!LSMask(ikp,nm) = 1 ! Land-Sea Mask
	456	! Test aquaplanète (Martin):
	457	LSMask(ikp,nm) = 0 ! Land-Sea Mask
	458	SlopSV(ikp,nm) = 0. ! Surface Slope
	459	FracSV(ikp,nm) = 1. / float(mwp) ! Mesh Fraction
	460	z0__SV(ikp,nm) = 1. ! z0 Momentum
	461	z0h_SV(ikp,nm) = 1. ! z0 Scalars
	462	!z0__SV(ikp,nm) = 0.1 ! z0 Momentum
	463	!z0h_SV(ikp,nm) = 0.1 ! z0 Scalars
	464	!z0__SV(ikp,nm) = 0.001 ! z0 Momentum
	465	!z0h_SV(ikp,nm) = 0.00001 ! z0 Scalars
	466	isnoSV(ikp,nm) = 0 ! NO Snow
	467	ispiSV(ikp,nm) = 0 ! NO Ice (Superimposed)
	468	iiceSV(ikp,nm) = 0 ! NO Ice
	469
	470	!iVgTSV(ikp,nm) = 5 ! GRASS MEDIUM Vegetat. Type (Mosaic) [-] (see Mod_SISVAT_dat)
	471	iVgTSV(ikp,nm) = 0 ! GRASS MEDIUM Vegetat. Type (Mosaic) [-] (see Mod_SISVAT_dat)
	472	!iSoTSV(ikp,nm) = 5 ! LOAM Soil Type (Mosaic) [-] (see Mod_SISVAT_dat)
	473	iSoTSV(ikp,nm) = 0 ! LOAM Soil Type (Mosaic) [-] (see Mod_SISVAT_dat)
	474	iWaFSV(ikp,nm) = 1 ! Soil Drainage Switch (Mosaic)(1=free)[-]
	475	alb0SV(ikp,nm) = 0.15 ! Dry Soil Albedo [-]
	476
	477	DO l=-nsol,0
	478	TsisSV(ikp,nm,l) = TsHOST(ikp) ! Soil-Ice-Snow Temperature (Mosaic) [K]
	479	eta_SV(ikp,nm,l) = &! Humidity (Mosaic) [kg/kg]
	480	& etadSV(iVgTSV(ikp,nm)) ! (Saturation is assumed)
	481	ENDDO
	482
	483	DO l=1,nsno
	484	TsisSV(ikp,nm,l) = 273.15 ! Soil-Ice-Snow Temperature (Mosaic) [K]
	485	eta_SV(ikp,nm,l) = 0.00 ! Soil-Ice-Snow Humidity (Mosaic) [kg/kg]
	486	dzsnSV(ikp,nm,l) = 0.00 ! Snow Layer Thickness (Mosaic) [m]
	487	ro__SV(ikp,nm,l) = 0.00 ! Snow Layer Density (Mosaic) [kg/m3]
	488	agsnSV(ikp,nm,l) = 0.00 ! Snow Layer Age (Mosaic) [day]
	489	G1snSV(ikp,nm,l) = 0.00 ! Snow Layer Property (d./s.)(Mosaic) [-]
	490	G2snSV(ikp,nm,l) = 0.00 ! Snow Layer Property (s./s.)(Mosaic) [-]
	491	istoSV(ikp,nm,l) = 0.00 ! Snow Layer History (Mosaic) [-]
	492	ENDDO
	493
	494	ENDDO
	495	ENDDO
	496
	497
	498	! ==============================
	499	END IF ! (it_EXP .LE. 1)
	500
	501	! FIRST INITIALIZATION (END)
	502	! ==============================
	503
	504
	505	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	506	!
	507	END IF ! (.NOT.FlagSV_SBC)
	508	!
	509	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	510
	511
	512	end subroutine PHY_SISVAT_INP

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

Download in other formats:

Original Format