LMDZ
surf_sisvat_mod.F90
Go to the documentation of this file.
1 MODULE surf_sisvat_mod
2  USE dimphy
3  IMPLICIT NONE
4 
5  INTEGER, PARAMETER :: nsnowmx=35
6  INTEGER, PARAMETER :: nsismx=46 ! = nsnowmx + nsoilmx
7 
8 CONTAINS
9 
10 
11 
12  SUBROUTINE surf_sisvat(knon,rlon,rlat, ikl2i, itime, dtime, debut, lafin, &
13  rmu0, swdown, lwdown, pexner, ps, p1lay, &
14  precip_rain, precip_snow, precip_snow_adv, snow_adv, &
15  bl_height, wind_velo, temp_air, dens_air, spechum, tsurf, &
16  rugos, snow_cont_air, alb_soil, slope, cloudf, &
17  radsol, qsol, tsoil, snow, snowhgt, qsnow, to_ice, sissnow, agesno, &
18  acoefh, acoefq, bcoefh, bcoefq, cdragh, &
19  runoff_lic, evap, fluxsens, fluxlat, dflux_s, dflux_l, &
20  tsurf_new, alb1, alb2, alb3, &
21  emis_new, z0_new, qsurf)
22 
23 ! +------------------------------------------------------------------------+
24 ! | |
25 ! | SubRoutine surf_sisvat: Interface between LMDZ and landice scheme |
26 ! | of the SISVAT (Soil/Ice Snow Vegetation Atmosphere Transfer Scheme)|
27 ! | |
28 ! | Author: Heinz Juergen Punge, LSCE June 2009 |
29 ! | based on the MAR-SISVAT interface by Hubert Gallee |
30 ! | |
31 ! +------------------------------------------------------------------------+
32 ! |
33 ! | In the current setup, SISVAT is used only to model the land ice |
34 ! | part of the surface; hence it is called with the compressed variables|
35 ! | from pbl_surface, and only by the surf_landice routine. |
36 ! | |
37 ! | In this interface it is assumed that the partitioning of the soil, |
38 ! | and hence the number of grid points is constant during a simulation, |
39 ! | hence eg. snow properties remain stored in the global SISVAT |
40 ! | variables between the calls and don't need to be handed over as |
41 ! | arguments. When the partitioning is supposed to change, make sure to |
42 ! | update the variables. |
43 ! | |
44 ! | INPUT |
45 ! | ^^^^^ VegMod: SISVAT is set up when .T. |
46 ! | SnoMod: Snow Pack is set up when .T. |
47 ! | reaLBC: Update Bound.Condit.when .T. |
48 ! | |
49 ! | INPUT (via MODULES VARxSV, VARySV, VARtSV) |
50 ! | ^^^^^ xxxxSV: SISVAT/LMDZ interfacing variables |
51 ! | |
52 ! | Preprocessing Option: SISVAT PHYSICS |
53 ! | ^^^^^^^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^ |
54 ! | # #HY |
55 ! | # #SN: Snow Model |
56 ! | # #BS: Blowing Snow Parameterization |
57 ! | |
58 ! | # #DS: diffuse radiation differing from direct |
59 ! | (variable RADsod must still be included) |
60 ! | # #CP: SBL, Col de Porte |
61 ! | # #cp Solar Radiation, Col de Porte |
62 ! | # #AG: Snow Ageing, Col de Porte |
63 ! | |
64 ! | Preprocessing Option: SISVAT IO |
65 ! | ^^^^^^^^^^^^^^^^^^^^^ ^^^^^^^^^ |
66 ! | FILE | CONTENT |
67 ! | ~~~~~~~~~~~~~~~~~~~~~+~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
68 ! | # ANI.yyyymmdd.LAB.nc | #NC: OUTPUT on NetCDF File (Stand Alone EXP.) |
69 ! | # SISVAT_iii_jjj_n | #WV: OUTPUT on ASCII File (SISVAT Variables) |
70 ! | # SISVATtroubles | #VF: OUTPUT on ASCII File (SISVAT Troubles) |
71 ! | | |
72 ! | # | #ES: OUTPUT/Verification: Energy Conservation |
73 ! | # | #E2: OUTPUT/Verification: Energy Consrv.2e pt.|
74 ! | | (no premature stop) |
75 ! | # | #MW: OUTPUT/Verification: H2O Conservation |
76 ! | # | #MS: OUTPUT/Verification: * Mass Conservation |
77 ! | # | #MI: OUTPUT/Verification: SeaIce Conservation |
78 ! | | |
79 ! | # SISVAT__zSn.OUT | #as: OUTPUT/Verification: Snow Layers Agrega. |
80 ! | # SISVAT__SnO.OUT | #aw: OUTPUT/Verification: Albedo Parameteriz. |
81 ! | # SISVATe_qSn.OUT | #em: OUTPUT/Verification: Energy/Water Budget |
82 ! | # SISVATu_qSn.OUT | #su: OUTPUT/Verification: Slush Parameteriz. |
83 ! | # SISVATw_qSo.OUT | #mw: OUTPUT/Verif+Detail: H2O Conservation |
84 ! | | |
85 ! | # SISVAT__GSn.OUT | #VP: OUTPUT/Verification: Snow Properties |
86 ! | # SISVAT__wEq.OUT | #EQ: OUTPUT/Verification: Snow/Ice Water Eqv. |
87 ! | | |
88 ! | # | #VR: VERIFICATION OUTPUT |
89 ! | # | #WR: Additional OUTPUT |
90 ! | |
91 ! +------------------------------------------------------------------------+
92 
93  USE var_sv
94  USE vardsv
95  USE varxsv
96  USE varysv
97  USE vartsv
98 
99  USE vardcp
100  USE varphy, ra_earth=>ra
101  USE yomcst_sisvat
102 
103  IMPLICIT NONE
104 
105 ! +--INTERFACE Variables
106 ! + ===================
107 
108  include "dimsoil.h"
109 
110 
111 ! +--Global Variables
112 ! + ================
113 ! Input Variables for SISVAT
114  INTEGER, INTENT(IN) :: knon
115  INTEGER, INTENT(IN) :: itime
116  REAL, INTENT(IN) :: dtime
117  LOGICAL, INTENT(IN) :: debut ! true if first step
118  LOGICAL, INTENT(IN) :: lafin ! true if last step
119 
120  INTEGER, DIMENSION(klon), INTENT(IN) :: ikl2i ! Index Decompression
121  REAL, DIMENSION(klon), INTENT(IN) :: rlon, rlat
122  REAL, DIMENSION(klon), INTENT(IN) :: rmu0 ! cos sol. zenith angle
123  REAL, DIMENSION(klon), INTENT(IN) :: swdown !
124  REAL, DIMENSION(klon), INTENT(IN) :: lwdown !
125  REAL, DIMENSION(klon), INTENT(IN) :: pexner ! Exner potential
126  REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow
127  REAL, DIMENSION(klon), INTENT(IN) :: precip_snow_adv, snow_adv
128  !Snow Drift
129  REAL, DIMENSION(klon), INTENT(IN) :: bl_height, wind_velo
130  REAL, DIMENSION(klon), INTENT(IN) :: temp_air, spechum, ps,p1lay
131  REAL, DIMENSION(klon), INTENT(IN) :: dens_air, tsurf
132  REAL, DIMENSION(klon), INTENT(IN) :: rugos,snow_cont_air
133  REAL, DIMENSION(klon), INTENT(IN) :: alb_soil, slope
134  REAL, DIMENSION(klon), INTENT(IN) :: cloudf
135  REAL, DIMENSION(klon), INTENT(IN) :: AcoefH, AcoefQ
136  REAL, DIMENSION(klon), INTENT(IN) :: BcoefH, BcoefQ
137  REAL, DIMENSION(klon), INTENT(IN) :: cdragh
138 
139 ! Variables exchanged between LMDZ and SISVAT
140  REAL, DIMENSION(klon,nsoilmx), INTENT(OUT) :: tsoil ! Soil Temperature
141  REAL, DIMENSION(klon), INTENT(OUT) :: qsol ! Soil Water Content
142  REAL, DIMENSION(klon), INTENT(INOUT) :: snow ! Tot snow mass [kg/m2]
143  REAL, DIMENSION(klon), INTENT(IN) :: radsol ! Surface absorbed rad.
144 
145 ! Output Variables from SISVAT
146  REAL, DIMENSION(klon), INTENT(OUT) :: alb1 ! Albedo SW
147  REAL, DIMENSION(klon), INTENT(OUT) :: alb2,alb3 ! Albedo LW
148  REAL, DIMENSION(klon), INTENT(OUT) :: emis_new ! Surface Emissivity
149  REAL, DIMENSION(klon), INTENT(OUT) :: z0_new ! Momentum Roughn Lgt
150  REAL, DIMENSION(klon), INTENT(OUT) :: runoff_lic ! Runoff
151  REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s ! d/dT sens. ht flux
152  REAL, DIMENSION(klon), INTENT(OUT) :: dflux_l ! d/dT latent ht flux
153  REAL, DIMENSION(klon), INTENT(OUT) :: fluxsens ! Sensible ht flux
154  REAL, DIMENSION(klon), INTENT(OUT) :: fluxlat ! Latent heat flux
155  REAL, DIMENSION(klon), INTENT(OUT) :: evap ! Evaporation
156  REAL, DIMENSION(klon), INTENT(OUT) :: agesno ! Snow age (top layer)
157  REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new ! Surface Temperature
158  REAL, DIMENSION(klon), INTENT(OUT) :: qsurf ! Surface Humidity
159  REAL, DIMENSION(klon), INTENT(OUT) :: qsnow ! Total H2O snow[kg/m2]
160  REAL, DIMENSION(klon), INTENT(OUT) :: snowhgt ! Snow height (m)
161  REAL, DIMENSION(klon), INTENT(OUT) :: to_ice ! Snow passed to ice
162  REAL, DIMENSION(klon), INTENT(OUT) :: sissnow ! Snow in model (kg/m2)
163 
164 ! +--OUTPUT for NetCDF File
165 ! + ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
166 ! REAL, DIMENSION(klonv) :: SOsoKL ! Abs Solar Radiation
167 ! REAL, DIMENSION(klonv) :: IRsoKL ! Abs IR Radiation
168 ! REAL, DIMENSION(klonv) :: HSsoKL ! Abs Sensible Ht Flux
169 ! REAL, DIMENSION(klonv) :: HLsoKL ! Abs Latent Heat Flux
170 ! REAL, DIMENSION(klonv) :: HLs_KL ! Evaporation
171 ! REAL, DIMENSION(klonv) :: HLv_KL ! Transpiration
172 !
173 ! REAL, DIMENSION(klonv) :: SOsoNC ! Abs Solar Radiation
174 ! REAL, DIMENSION(klonv) :: IRsoNC ! Abs IR Radiation
175 ! REAL, DIMENSION(klonv) :: HSsoNC ! Abs Sensible Ht Flux
176 ! REAL, DIMENSION(klonv) :: HLsoNC ! Abs Latent Heat Flux
177 ! REAL, DIMENSION(klonv) :: HLs_NC ! Evaporation
178 ! REAL, DIMENSION(klonv) :: HLv_NC ! Transpiration
179 !
180 ! REAL, DIMENSION(klonv,nsoilmx) :: eta_NC ! nsoilmx=nsol+1
181 ! REAL, DIMENSION(klonv,nsoilmx) :: tsolNC !
182 ! REAL, DIMENSION(klonv) :: snowNC !
183 ! INTEGER, DIMENSION(klonv) :: isnoNC !
184 ! INTEGER, DIMENSION(klonv) :: ispiNC !
185 ! INTEGER, DIMENSION(klonv) :: iiceNC !
186 ! REAL, DIMENSION(klonv) :: swaNC !
187 ! REAL, DIMENSION(klonv) :: swsNC !
188 ! INTEGER, DIMENSION(klonv,nsno) :: istoNC !
189 ! REAL, DIMENSION(klonv,nsno) :: dzsnNC !
190 ! REAL, DIMENSION(klonv,nsno) :: rhosnNC !
191 ! REAL, DIMENSION(klonv,nsno) :: etasnNC !
192 ! REAL, DIMENSION(klonv,nsno) :: tsnNC !
193 ! REAL, DIMENSION(klonv,nsno) :: g1snNC !
194 ! REAL, DIMENSION(klonv,nsno) :: g2snNC !
195 ! REAL, DIMENSION(klonv,nsno) :: agsnNC !
196 ! REAL, DIMENSION(klonv) :: meltNC !
197 ! REAL, DIMENSION(klonv) :: refrNC !
198 ! REAL, DIMENSION(klonv) :: alb1NC !
199 ! REAL, DIMENSION(klonv) :: alb2NC !
200 ! REAL, DIMENSION(klonv) :: rnofNC !
201 
202 ! + Optional Variables:
203 ! +--V, dT(a-s) Time Moving Averages
204 ! + ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
205 
206 ! #AW real V__mem(klon,ntaver) ! ntaver defined in LMDZ_SL.inc
207 ! #AW real VVmmem(klon) !
208 ! #AW common/SVeSBLmem/V__mem,VVmmem !
209 ! #AH real T__mem(klon,ntaver) !
210 ! #AH real dTmmem(klon) !
211 ! #AH common/STeSBLmem/T__mem,dTmmem !
212 ! +--u*, u*T*, u*s* Time Moving Averages
213 ! + ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
214 ! #AM real u__mem(klon,ntaver) ! ntaver defined in LMDZ_SL.inc
215 ! #AT real uT_mem(klon,ntaver) !
216 ! #AS real us_mem(klon,ntaver) !
217 ! #AM common/S_eSBLmem/u__mem !
218 ! #AT. ,uT_mem !
219 ! #AS. ,us_mem !
220 !! #AH REAL, DIMENSION(klonv,ntaver) :: T__mem !
221 !! #AH REAL, DIMENSION(klonv) :: dTmmem !
222 !! #AW REAL, DIMENSION(klonv,ntaver) :: V__mem !
223 !! #AM REAL, DIMENSION(klonv) :: u__mem !
224 !! #AT REAL, DIMENSION(klonv) :: uT_mem !
225 !! #AS REAL, DIMENSION(klonv) :: us_mem !
226 !! #AW REAL, DIMENSION(klonv) :: VVmmem !
227 
228 
229 ! +--Internal Variables
230 ! + ===================
231 
232  CHARACTER(len=20) :: fichnom, fn_outfor ! Name for output file
233  INTEGER :: i, ig, ikl, isl, isn, nt
234  INTEGER :: gp_outfor, un_outfor
235 
236  REAL, PARAMETER :: f1=0.5
237  REAL, PARAMETER :: sn_upp=5000.,sn_low=500.
238  REAL, PARAMETER :: sn_add=400.,sn_div=2.
239  ! snow mass upper,lower limit,
240  ! added mass/division lowest layer
241  REAL, PARAMETER :: c1_zuo=12.960e+4, c2_zuo=2.160e+6
242  REAL, PARAMETER :: c3_zuo=1.400e+2, czemin=1.e-3
243  ! Parameters for drainage
244 ! c1_zuo/ 2.796e+4/,c2_zuo/2.160e+6/,c3_zuo/1.400e+2/ ! Tuning
245 ! +... Run Off Parameters
246 ! + 86400*1.5 day ...*25 days (Modif. ETH Camp: 86400*0.3day)
247 ! + (Zuo and Oerlemans 1996, J.Glacio. 42, 305--317)
248 
249  REAL, DIMENSION(klon) :: eps0SL ! surface Emissivity
250  REAL :: zsigma, Ua_min, Us_min
251  REAL :: lambda ! Par. soil discret.
252  REAL, DIMENSION(nsoilmx), SAVE :: dz1,dz2 ! Soil layer thicknesses
253 !$OMP THREADPRIVATE(dz1)
254  LOGICAL, SAVE :: firstcall=.true.,snomod, ok_outfor=.false.!$OMP THREADPRIVATE(firstcall)
255 
256 
257 
258 ! + Optional:
259 !c #BW INTEGER :: noUNIT
260 !c #BW REAL :: BlowST,SnowSB
261 
262 
263 ! +--Internal Variables
264 ! + ==================
265 
266  INTEGER :: ivg,iso
267 
268 !========================================================================
269 
270  snomod=.true.
271  zsigma=1000.
272  dt__sv=dtime
273  IF (ok_outfor) THEN
274  un_outfor=51 ! unit number for point output
275  gp_outfor=79! 633 !79 ! grid point number for point output
276  fn_outfor='outfor_SV.dat'
277  END IF
278 
279 ! write(*,*)'Start of simulation? ',debut !hj
280  IF (debut) THEN
281  firstcall=.true.
282  ini_sv=.false.
283  ELSE
284  firstcall=.false.
285  ini_sv=.true.
286  END IF
287 
288 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
289 
290 
291 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
292 ! + INITIALISATION: BEGIN +++
293 ! + -------------------------
294 ! +
295 ! + Compute soil discretization (as for LMDZ)
296 ! + -----------------------------------------
297  IF (firstcall) THEN
298 
299 ! +--Array size
300  klonv=klon
301  knonv=knon
302  write(*,*)'klon',klon,'klonv',klonv,'knon',knon
303 
304 
305  CALL init_vartsv
306  CALL init_varxsv
307  CALL init_varysv
308 
309  eps0sl(:)=0.
310 ! +--Soil layer thickness
311 ! + -----------------------
312 ! write(*,'(/a)') 'Start SISVAT init: soil discretization ', nsoilmx
313  CALL get_soil_levels(dz1,dz2,lambda)
314  lambsv=lambda
315  dz1_sv(1:knon,1:) = 0.
316  dz2_sv(1:knon,1:) = 0.
317 
318  DO isl = -nsol,0
319  dz_dsv(isl) = 0.5e-3*dz2(1-isl) ! Soil layer thickness
320  DO ikl=1,knon
321  dz1_sv(ikl,isl) = dz1(1-isl) !1.e-3*
322  dz2_sv(ikl,isl) = dz2(1-isl) !1.e-3*
323  END DO
324 ! IF (knon > 0) THEN
325 ! write(*,*)'level:',dz_dSV(isl),dz1_SV(1,isl),dz2_SV(1,isl)
326 ! END IF
327  END DO
328 
329  DO ikl=1,knon
330  eps0sl(ikl )= 1.
331  alb0sv(ikl) = alb_soil(ikl) ! Soil Albedo
332 
333 ! + Soil Upward IR Flux, Water Fluxes, roughness length
334  irs_sv(ikl) = &
335  -eps0sl(ikl)* rsigma*temp_air(ikl) & ! Upward IR Flux
336  *temp_air(ikl)*temp_air(ikl)*temp_air(ikl)
337  tvegsv(ikl) = temp_air(ikl)
338  ! + Soil
339  DO isl = -nsol,0
340  tsissv(ikl,isl) = temp_air(ikl) !tsoil(ikl,1-isl) Soil Temperature
341  eta_sv(ikl,isl) = 0.0001 !etasoil(ikl,1-isl)Soil Water[m3/m3]
342  ro__sv(ikl,isl) = 50. !rosoil(ikl,1-isl) soil water volumic mass
343  END DO
344  END DO
345 
346 ! +--Surface Fall Line Slope
347 ! + -----------------------
348  IF (snomod) THEN
349  DO ikl=1,knon
350  slopsv(ikl) = slope(ikl)
351  swf_sv(ikl) = & ! Normalized Decay of the
352  exp(-dt__sv & ! Surficial Water Content
353  /(c1_zuo & !(Zuo and Oerlemans 1996,
354  +c2_zuo*exp(-c3_zuo*abs(slopsv(ikl))))) ! J.Glacio. 42, 305--317)
355  END DO
356  END IF
357 
358 ! +--SBL Characteristics history, not stored at restart.
359  DO ikl=1,knon
360  DO nt=1,ntaver ! #AW
361  v__mem(ikl,nt)=wind_velo(ikl) ! #AH
362  t__mem(ikl,nt)=temp_air(ikl)-tsurf(ikl)
363  END DO
364  END DO
365 
366 
367 ! + SISVAT_ini (as for use with MAR, but not computing soil layers)
368 ! + -------------------------------------------------------------
369 ! write(*,'(/a)') 'Start SISVAT initialization: SISVAT_ini'
370  CALL sisvat_ini(knon)
371 
372 ! open output file
373  IF (ok_outfor) THEN
374  open(unit=un_outfor,status='new',file=fn_outfor)
375  rewind un_outfor
376  ikl=gp_outfor ! index sur la grille land ice
377  ig=611 ! index sur la grille globale
378  write(un_outfor,501) fn_outfor, ikl, rlon(ig),rlat(ig)
379 501 format(/,a18,/,'Grid point ',i4,' Long',f9.4,' Lat ',f9.4 &
380  & ,/,'++++++++++++++++++++++++++++++++++++++++++++++', &
381  & '++++++++++++++++++++++++++++++++++++++++++++++', &
382  & '++++++++++++++++++++++++++++++++++++++++++++++', &
383  & /,' SWdown + IRdown + Wind + Temp. + Humid. ', &
384  & '+ Press +Precip_l+Precip_s+ Tsrf + Clouds +' &
385  & '+ Zenith + BLhgt + Densair+ Exner +' &
386  & ,/,' sol_SV + IRd_SV + VV__SV + TaT_SV + QaT_SV ', &
387  & '+ ps__SV + drr_SV + dsn_SV + Tsf_SV + cld_SV +' &
388  & '+ coszSV + za__SV + rhT_SV + ExnrSV+' &
389  & ,/,' W/m2 + W/m2 + m/s + K + kg/kg ', &
390  & '+ Pa + kg/m2/s+ kg/m2/s+ K + /1 +' &
391  & '+ - + m + kg/m3 + +' &
392  & ,/,'++++++++++++++++++++++++++++++++++++++++++++++', &
393  & '++++++++++++++++++++++++++++++++++++++++++++++', &
394  & '++++++++++++++++++++++++++++++++++++++++++++++')
395  END IF
396 
397 ! +--Read restart file
398 ! + =================================================
399 
400 ! Martin
401  print*, 'On debranche sisvatetat0'
402 ! Martin
403 
404  ! CALL sisvatetat0("startsis.nc",ikl2i)
405 
406 
407  END IF ! firstcall
408 ! +
409 ! + +++ INITIALISATION: END +++
410 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
411 
412 
413 
414 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
415 ! + READ FORCINGS
416 ! + ------------------------
417 ! + Update Forcings for SISVAT given by the LMDZ model.
418 ! +
419  DO ikl=1,knon
420 
421 ! +--Atmospheric Forcing (INPUT)
422 ! + ^^^^^^^^^^^^^^^^^^^ ^^^^^
423  zsblsv = 1000. ! [m]
424  za__sv(ikl) = bl_height(ikl) ! Height boundary layer [m]
425  ua_min = epsi !
426 ! #VM Ua_min = 0.2 * sqrt(za__SV(ikl) ) !
427  ua_min = 0.2 * sqrt(za__sv(ikl) ) !
428  vv__sv(ikl) = max(ua_min, wind_velo(ikl)) ! Wind velocity [m/s]
429  us_min = 0.01
430  us__sv(ikl) = max(us_min, us__sv(ikl) )
431  tat_sv(ikl) = temp_air(ikl) ! BL top Temperature [K]
432  exnrsv(ikl) = pexner(ikl) ! Exner potential
433  rht_sv(ikl) = dens_air(ikl) ! Air density
434  qat_sv(ikl) = spechum(ikl) ! Specific humidity
435 ! #VX dQa_SV(ikl) = 0. !hj dtDiff/zsigma(mz) ! Water Vapor Flux
436  ps__sv(ikl) = ps(ikl) ! surface pressure [Pa]
437  p1l_sv(ikl) = p1lay(ikl) ! lowest atm. layer press[Pa]
438 
439 ! +--Energy Fluxes (INPUT)
440 ! + ^^^^^^^^^^^^^ ^^^^^
441  coszsv(ikl) = max(czemin,rmu0(ikl)) ! cos(zenith.Dist.)
442  sol_sv(ikl) = swdown(ikl) ! downward Solar
443  ird_sv(ikl) = lwdown(ikl) ! downward IR
444  rsolsv(ikl) = radsol(ikl) ! surface absorbed rad.
445 !hj 110511
446 !hj 121011 IRs_SV(ikl) = min(IRs_SV(ikl),-epsi) ! check upward IR
447 !hj
448 ! +--Water Fluxes (INPUT)
449 ! + ^^^^^^^^^^^^^ ^^^^^
450  drr_sv(ikl) = precip_rain(ikl) ! Rain fall rate [kg/m2/s]
451  dsn_sv(ikl) = precip_snow(ikl) ! Snow fall rate [kg/m2/s]
452 !c #BS dbsnow = -SLussl(i,j,n) ! Erosion
453 !c #BS. *dtPhys *rhT_SV(ikl) /ro_Wat
454 !c #BS dsnbSV(ikl) = snow_adv(ikl) ! min(max(zero,dbsnow)
455 !c #BS. / max(epsi,d_snow),unun)
456 !c #BS dbs_SV(ikl) = snow_cont_air(ikl)
457 !c #BS blowSN(i,j,n) ! [kg/m2]
458 
459 ! +--Soil/Canopy (INPUT)
460 ! + ^^^^^^^^^^^ ^^^^^
461  alb0sv(ikl) = alb_soil(ikl) ! Soil background Albedo
462  acohsv(ikl) = acoefh(ikl)
463  bcohsv(ikl) = bcoefh(ikl)
464  acoqsv(ikl) = acoefq(ikl)
465  bcoqsv(ikl) = bcoefq(ikl)
466  cdh_sv(ikl) = cdragh(ikl)
467  tsf_sv(ikl) = tsurf(ikl) !hj 12 03 2010
468 
469 ! +--Energy Fluxes (INPUT/OUTPUT)
470 ! + ^^^^^^^^^^^^^ ^^^^^^^^^^^^
471  IF (.not.firstcall) THEN
472 ! active hj 110411
473  cld_sv(ikl) = cloudf(ikl) ! Cloudiness
474  END IF
475 
476 ! +--Time Averages of wind and surface-atmosphere temperature difference
477 ! +--for turbulence calculations ! #AA
478 ! Update stored arrays
479  DO nt=1,ntaver-1 ! #AW
480  v__mem(ikl,nt ) = v__mem(ikl,nt+1) ! #AH
481  t__mem(ikl,nt ) = t__mem(ikl,nt+1) ! #AA
482  ENDDO ! #AW
483  v__mem(ikl,ntaver)=wind_velo(ikl) ! #AH
484  t__mem(ikl,ntaver)=temp_air(ikl)-tsurf(ikl) ! #AW
485 ! Calculate averages
486  vvmmem(ikl) = 0.0 ! #AH
487  dtmmem(ikl) = 0.0 ! #AA
488  DO nt=1,ntaver ! #AW
489  vvmmem(ikl) = vvmmem(ikl)+v__mem(ikl,nt) ! #AH
490  dtmmem(ikl) = dtmmem(ikl)+t__mem(ikl,nt) ! #AA
491  ENDDO ! #AW
492  vvmmem(ikl) = vvmmem(ikl)/ntaver ! #AH
493  dtmmem(ikl) = dtmmem(ikl)/ntaver
494  END DO
495 
496 !
497 ! + +++ READ FORCINGS: END +++
498 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
499 
500 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
501 ! + OUTPUT FORCING
502 !
503  IF (ok_outfor) THEN
504  ikl=gp_outfor
505  write(un_outfor,5000) sol_sv(ikl), ird_sv(ikl), vv__sv(ikl), &
506  & tat_sv(ikl), qat_sv(ikl), ps__sv(ikl), &
507  & drr_sv(ikl), dsn_sv(ikl), tsf_sv(ikl), &
508  & cld_sv(ikl), coszsv(ikl), za__sv(ikl), &
509  & rht_sv(ikl), exnrsv(ikl)
510 
511  5000 format(f8.3,' ',f8.3,' ',f8.4,' ',f8.4,' ',f10.8,' ',f8.1,' ', &
512  & f8.6,' ',f8.6,' ',f8.4,' ',f8.6,' ',f8.6,' ',f8.3,' ', &
513  & f8.4,' ',f8.6)
514 
515  ENDIF
516 !
517 ! + OUTPUT FORCINGS: END +++
518 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
519 
520 
521 
522 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
523 ! +--SISVAT EXECUTION
524 ! + ----------------
525 ! +
526 ! write(*,*) ' Start SISVAT execution!'
527 
528  call sisvat(snomod,.false.,1)
529 ! BloMod,jjtime)
530 
531 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
532 ! + RETURN RESULTS
533 ! + --------------
534 ! + Return (compressed) SISVAT variables to LMDZ
535 ! +
536  DO ikl=1,knon ! use only 1:knon (actual ice sheet..)
537  runoff_lic(ikl) = rnofsv(ikl)*dtime ! RunOFF: intensity* time step
538  dflux_s(ikl) = dsdtsv(ikl) ! Sens.H.Flux T-Der.
539  dflux_l(ikl) = dldtsv(ikl) ! Latn.H.Flux T-Der.
540  fluxsens(ikl) = hss_sv(ikl) ! HS
541  fluxlat(ikl) = hls_sv(ikl) ! HL
542  evap(ikl) = hls_sv(ikl)/rlvtt ! Evaporation RLVTT=LhvH2O
543  z0_new(ikl) = z0h_sv(ikl) ! Moment.Roughn.L.
544 
545 
546  snow(ikl) = 0.
547  snowhgt(ikl) = 0.
548  qsnow(ikl) = 0.
549  qsol(ikl) = 0.
550 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
551 ! + Check snow thickness, add if too thin, substract if too thick
552 
553  sissnow(ikl) = 0. !()
554  DO isn = 1,isnosv(ikl)
555  sissnow(ikl) = sissnow(ikl)+dzsnsv(ikl,isn)* ro__sv(ikl,isn)
556  END DO
557 
558  IF (sissnow(ikl) .LE. sn_low) THEN !add snow
559  IF (isnosv(ikl).GE.1) THEN
560  dzsnsv(ikl,1) = dzsnsv(ikl,1) + sn_add/max(ro__sv(ikl,1),epsi)
561  toicsv(ikl) = toicsv(ikl) - sn_add
562  ELSE
563  write(*,*) 'Attention, bare ice... point ',ikl
564  isnosv(ikl) = 1
565  istosv(ikl,1) = 100
566 !ym istoSV(ikl) = 100
567  ro__sv(ikl,1) = 400.
568  dzsnsv(ikl,1) = sn_add/max(ro__sv(ikl,1),epsi) ! 1.
569  eta_sv(ikl,1) = epsi
570  tsissv(ikl,1) = min(tsissv(ikl,0),tfsnow-0.2)
571  g1snsv(ikl,1) = 99.
572  g2snsv(ikl,1) = 0.3
573  agsnsv(ikl,1) = 10.
574  toicsv(ikl) = toicsv(ikl) - sn_add
575  END IF
576  END IF
577 
578  IF (sissnow(ikl) .ge. sn_upp) THEN !thinnen snow layer below
579  dzsnsv(ikl,1) = dzsnsv(ikl,1)/sn_div
580  toicsv(ikl) = toicsv(ikl)+dzsnsv(ikl,1)*ro__sv(ikl,1)/sn_div
581  END IF
582 
583  sissnow(ikl) = 0. !()
584  DO isn = 1,isnosv(ikl)
585  sissnow(ikl) = sissnow(ikl)+dzsnsv(ikl,isn)* ro__sv(ikl,isn)
586 ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
587  snowhgt(ikl) = snowhgt(ikl)+dzsnsv(ikl,isn)
588  qsnow(ikl) = qsnow(ikl)+1e03*eta_sv(ikl,isn)*dzsnsv(ikl,isn)
589  END DO
590  snow(ikl) = sissnow(ikl)+toicsv(ikl)
591  to_ice(ikl) = toicsv(ikl)
592 ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
593 
594 
595  DO isl = -nsol,0
596  tsoil(ikl,1-isl) = tsissv(ikl,isl) ! Soil Temperature
597  qsol(ikl) = qsol(ikl) &
598  +eta_sv(ikl,isl) * dz_dsv(isl)
599  END DO
600  agesno(ikl) = agsnsv(ikl,isnosv(ikl)) ! [day]
601 
602  alb1(ikl) = alb1sv(ikl) ! Albedo VIS
603  alb2(ikl) = ((so1dsv-f1)*alb1sv(ikl) &
604  & +so2dsv*alb2sv(ikl)+so3dsv*alb3sv(ikl))/f1
605  ! Albedo NIR
606  alb3(ikl) = alb3sv(ikl) ! Albedo FIR
607  tsurf_new(ikl) = tsfnsv(ikl)
608 !hj220711 Tsrfsv(ikl) ! Surf.Temperature
609 ! TsisSV(ikl,0) *(1-min(1,isnoSV(ikl))) &
610 ! +TsisSV(ikl,max(1,isnoSV(ikl))) * min(1,isnoSV(ikl))
611 ! tsurf_new(ikl) = max(Ts_Min,tsurf_new(ikl) )
612 ! tsurf_new(ikl) = min(Ts_Max,tsurf_new(ikl) )
613  qsurf(ikl) = qat_sv(ikl)
614  emis_new(ikl) = eps0sl(ikl)
615 
616 !!!hjp 230611 sorties
617 !! qsnow(ikl)=TsisSV(ikl,isnoSV(ikl))
618 !! sissnow(ikl)=TsisSV(ikl,isnoSV(ikl)-1)
619 !! snowhgt(ikl)=TsisSV(ikl,isnoSV(ikl)-2)
620 !! qsol(ikl)=dzsnSV(ikl,isnoSV(ikl)-1)
621 !! agesno(ikl)=dzsnSV(ikl,isnoSV(ikl))
622 
623 
624  END DO
625 
626 
627 ! + -----------------------------
628 ! + END --- RETURN RESULTS
629 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
630 
631 
632 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
633 ! + Uncompress SISVAT output variables & store
634 ! + -----------------------------
635 ! +
636 ! DO ikl = 1,knon
637 ! i = ikl2i(ikl) ! Compression index
638 ! IF (i.GE.0) THEN
639 ! SOsoNC(i) = SOsoKL(ikl) ! Absorb.Sol.Rad.
640 ! IRsoNC(i) = IRsoKL(ikl) ! Absorb.IR Rad.
641 ! HSsoNC(i) = HSsoKL(ikl) ! HS
642 ! HLsoNC(i) = HLsoKL(ikl) ! HL
643 ! HLs_NC(i) = HLs_KL(ikl) ! Evaporation
644 ! HLv_NC(i) = HLv_KL(ikl) ! Transpiration
645 !
646 ! DO isl = -nsol,0
647 ! eta_NC(i,1-isl) = eta_SV(ikl,isl) ! Soil Humidity
648 ! tsolNC(i,1-isl) = TsisSV(ikl,isl) ! Soil Temperature
649 ! END DO
650 ! snowNC(i) = snow(ikl) ! Snow mass
651 ! isnoNC(i) = isnoSV(ikl) ! Nb Snow/Ice Lay.
652 ! ispiNC(i) = ispiSV(ikl) ! Nb Supr.Ice Lay.
653 ! iiceNC(i) = iiceSV(ikl) ! Nb Ice Lay.
654 ! swaNC(i) = rusnSV(ikl) ! Surficial Water
655 ! swsNC(i) = SWS_SV(ikl) ! Surficial Wat.St.
656 ! DO isn = 1,nsno
657 ! istoNC(i,isn) = istoSV(ikl,isn) ! [-]
658 ! dzsnNC(i,isn) = dzsnSV(ikl,isn) ! [m]
659 ! rhosnNC(i,isn) = ro__SV(ikl,isn) ! [kg/m3]
660 ! etasnNC(i,isn) = eta_SV(ikl,isn) ! [m3/m3]
661 ! tsnNC(i,isn) = TsisSV(ikl,isn) ! [K]
662 ! g1snNC(i,isn) = G1snSV(ikl,isn) ! [-] [-]
663 ! g2snNC(i,isn) = G2snSV(ikl,isn) ! [-] [0.0001 m]
664 ! agsnNC(i,isn) = agsnSV(ikl,isn) ! [day]
665 ! END DO
666 !!?c #IB depsubNC(i) = wes_SV(ikl) ! Depo. / Subli.
667 ! meltNC(i) = wem_SV(ikl) ! Melting
668 ! refrNC(i) = wer_SV(ikl) ! Refreezing
669 ! alb1NC(i) = alb1sv(ikl) ! Albedo SW
670 ! alb2NC(i) = (alb1sv(ikl)+3*alb2sv(ikl)+alb3sv(ikl))/5
671  ! Albedo LW
672 ! rnofNC(i) = RnofSV(ikl) ! Run OFF Intensity
673 !
674 !!? SL_z0(i) = Z0m_SV(ikl) ! Moment.Roughn.L.
675 !!? SL_r0(i,j,n) = Z0h_SV(ikl) ! Heat Roughn.L.
676 !!! zWE_NC(i) = zWE_SV(ikl) ! Current *Thick.
677 !!! zWEcNC(i) = zWEcSV(ikl) ! Non-Erod.*Thick.
678 !!! hSalNC(i) = hSalSV(ikl) ! Salt.Layer Height
679 !!! hsenSL(i) = -SLuts(i,j) * rhAir *cp ! Sensible Heat Flux
680 !!! hlatSL(i) = -SLuqs(i,j) * rhAir *Lv_H2O ! Latent Heat Flux
681 ! END IF
682 ! END DO
683 
684 !hj + Put storage to Output file here!
685  IF (lafin) THEN
686 
687  fichnom = "restartsis.nc"
688  CALL sisvatredem("restartsis.nc",ikl2i,rlon,rlat)
689  IF (ok_outfor) THEN
690  close(unit=un_outfor)
691  END IF
692  END IF
693 ! + -----------------------------
694 ! + END --- RETURN RESULTS
695 ! +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
696 
697  END SUBROUTINE surf_sisvat
698 
699 
700 
701 
702 
703 
704 
705  SUBROUTINE get_soil_levels(dz1, dz2, lambda)
706 ! ======================================================================
707 ! Routine to compute the vertical discretization of the soil in analogy
708 ! to LMDZ. In LMDZ it is done in soil.F, which is not used in the case
709 ! of SISVAT, therefore it's needed here.
710 !
711  USE mod_phys_lmdz_mpi_data, ONLY : is_mpi_root
712  USE mod_phys_lmdz_para
713 
714  include "dimsoil.h"
715 
716  REAL, DIMENSION(nsoilmx), INTENT(OUT) :: dz2, dz1
717  REAL, INTENT(OUT) :: lambda
718 
719 
720 !-----------------------------------------------------------------------
721 ! Depthts:
722 ! --------
723  REAL fz,rk,fz1,rk1,rk2
724  REAL min_period, dalph_soil
725  INTEGER ierr,jk
726 
727  fz(rk)=fz1*(dalph_soil**rk-1.)/(dalph_soil-1.)
728 
729 ! write(*,*)'Start soil level computation'
730 !-----------------------------------------------------------------------
731 ! Calculation of some constants
732 ! NB! These constants do not depend on the sub-surfaces
733 !-----------------------------------------------------------------------
734 !-----------------------------------------------------------------------
735 ! ground levels
736 ! grnd=z/l where l is the skin depth of the diurnal cycle:
737 !-----------------------------------------------------------------------
738 
739  min_period=1800. ! en secondes
740  dalph_soil=2. ! rapport entre les epaisseurs de 2 couches succ.
741 ! !$OMP MASTER
742 ! IF (is_mpi_root) THEN
743 ! OPEN(99,file='soil.def',status='old',form='formatted',iostat=ierr)
744 ! IF (ierr == 0) THEN ! Read file only if it exists
745 ! READ(99,*) min_period
746 ! READ(99,*) dalph_soil
747 ! PRINT*,'Discretization for the soil model'
748 ! PRINT*,'First level e-folding depth',min_period, &
749 ! ' dalph',dalph_soil
750 ! CLOSE(99)
751 ! END IF
752 ! ENDIF
753 ! !$OMP END MASTER
754 ! CALL bcast(min_period)
755 ! CALL bcast(dalph_soil)
756 
757 ! la premiere couche represente un dixieme de cycle diurne
758  fz1=sqrt(min_period/3.14)
759 
760  DO jk=1,nsoilmx
761  rk1=jk
762  rk2=jk-1
763  dz2(jk)=fz(rk1)-fz(rk2)
764  ENDDO
765  DO jk=1,nsoilmx-1
766  rk1=jk+.5
767  rk2=jk-.5
768  dz1(jk)=1./(fz(rk1)-fz(rk2))
769  ENDDO
770  lambda=fz(.5)*dz1(1)
771  print*,'full layers, intermediate layers (seconds)'
772  DO jk=1,nsoilmx
773  rk=jk
774  rk1=jk+.5
775  rk2=jk-.5
776  print *,'fz=', &
777  fz(rk1)*fz(rk2)*3.14,fz(rk)*fz(rk)*3.14
778  ENDDO
779 
780  END SUBROUTINE get_soil_levels
781 
782  SUBROUTINE sisvat_ini(knon)
783 
784 !C +------------------------------------------------------------------------+
785 !C | MAR SISVAT_ini Jd 11-10-2007 MAR |
786 !C | SubRoutine SISVAT_ini generates non time dependant SISVAT parameters |
787 !C +------------------------------------------------------------------------+
788 !C | PARAMETERS: klonv: Total Number of columns = |
789 !C | ^^^^^^^^^^ = Total Number of continental grid boxes |
790 !C | X Number of Mosaic Cell per grid box |
791 !C | |
792 !C | INPUT: dt__SV : Time Step [s] |
793 !C | ^^^^^ dz_dSV : Layer Thickness [m] |
794 !C | |
795 !C | OUTPUT: RF__SV : Root Fraction in Layer isl [-] |
796 !C | ^^^^^^ rocsSV : Soil Contrib. to (ro c)_s exclud.Water [J/kg/K] |
797 !C | etamSV : Soil Minimum Humidity [m3/m3] |
798 !C | (based on a prescribed Soil Relative Humidity) |
799 !C | s1__SV : Factor of eta**( b+2) in Hydraul.Diffusiv. |
800 !C | s2__SV : Factor of eta**( b+2) in Hydraul.Conduct. |
801 !C | aKdtSV : KHyd: Piecewise Linear Profile: a * dt [m] |
802 !C | bKdtSV : KHyd: Piecewise Linear Profile: b * dt [m/s] |
803 !C | dzsnSV(0): Soil first Layer Thickness [m] |
804 !C | dzmiSV : Distance between two contiguous levels [m] |
805 !C | dz78SV : 7/8 (Layer Thickness) [m] |
806 !C | dz34SV : 3/4 (Layer Thickness) [m] |
807 !C | dz_8SV : 1/8 (Layer Thickness) [m] |
808 !C | dzAvSV : 1/8 dz_(i-1) + 3/4 dz_(i) + 1/8 dz_(i+1) [m] |
809 !C | dtz_SV : dt/dz [s/m] |
810 !C | OcndSV : Swab Ocean / Soil Ratio [-] |
811 !C | Implic : Implicit Parameter (0.5: Crank-Nicholson) |
812 !C | Explic : Explicit Parameter = 1.0 - Implic |
813 !C | |
814 !C | # OPTIONS: #ER: Richards Equation is not smoothed |
815 !C | # ^^^^^^^ #kd: De Ridder Discretization |
816 !C | # #SH: Hapex-Sahel Values !
817 !C | |
818 !C +------------------------------------------------------------------------+
819 !
820 !
821 
822 !C +--Global Variables
823 !C + ================
824 
825  USE varphy, ra_earth=>ra
826  USE var_sv
827  USE vardsv
828  USE var0sv
829  USE varxsv
830  USE vartsv
831  USE varxsv
832  USE varysv
833  IMPLICIT NONE
834 
835 
836 
837 !C +--Arguments
838 !C + ==================
839  INTEGER,INTENT(IN) :: knon
840 
841 !C +--Internal Variables
842 !C + ==================
843 
844  INTEGER :: ivt ,ist ,ivg ,ikl ,isl ,isn ,ikh
845  INTEGER :: misl_2,nisl_2
846  REAL :: zDepth
847  REAL :: d__eta,eta__1,eta__2,Khyd_1,Khyd_2
848  REAL,PARAMETER :: RHsMin= 0.001 ! Min.Soil Relative Humidity
849  REAL :: PsiMax ! Max.Soil Water Potential
850  REAL :: a_Khyd,b_Khyd ! Piecewis.Water Conductivity
851 
852 
853 !c #WR REAL :: Khyd_x,Khyd_y
854 
855 
856 
857 !C +--Non Time Dependant SISVAT parameters
858 !C + ====================================
859 
860 !C +--Soil Discretization
861 !C + -------------------
862 
863 !C +--Numerical Scheme Parameters
864 !C + ^^^^^^^^^^^^^^^^^^^^^^^^^^^
865  implic = 0.75 ! 0.5 <==> Crank-Nicholson
866  explic = 1.00 - implic !
867 
868 !C +--Soil/Snow Layers Indices
869 !C + ^^^^^^^^^^^^^^^^^^^^^^^^
870  DO isl=-nsol,0
871  islpsv(isl) = isl+1
872  islpsv(isl) = min( islpsv(isl),0)
873  islmsv(isl) = isl-1
874  islmsv(isl) = max(-nsol,islmsv(isl))
875  END DO
876 
877  DO isn=1,nsno
878  isnpsv(isn) = isn+1
879  isnpsv(isn) = min( isnpsv(isn),nsno)
880  END DO
881 
882 !C +--Soil Layers Thicknesses
883 !C + ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
884 ! Not used here as LMDZ method is applied, see SUBROUTINE get_soil_levels!
885 !c #kd IF (nsol.gt.4) THEN
886 !c #kd DO isl=-5,-nsol,-1
887 !c #kd dz_dSV(isl)= 1.
888 !c #kd END DO
889 !c #kd END IF
890 !
891 ! IF (nsol.ne.4) THEN
892 ! DO isl= 0,-nsol,-1
893 ! misl_2 = -mod(isl,2)
894 ! nisl_2 = -isl/2
895 ! dz_dSV(isl)=(((1-misl_2) * 0.001
896 ! . + misl_2 * 0.003) * 10**(nisl_2)) * 4.
897 !C +... dz_dSV(0) = Hapex-Sahel Calibration: 4 mm
898 !
899 !c +SH dz_dSV(isl)=(((1-misl_2) * 0.001
900 !c +SH. + misl_2 * 0.003) * 10**(nisl_2)) * 1.
901 !
902 !c #05 dz_dSV(isl)=(((1-misl_2) * 0.001
903 !c #05. + misl_2 * 0.008) * 10**(nisl_2)) * 0.5
904 ! END DO
905 ! dz_dSV(0) = 0.001
906 ! dz_dSV(-1) = dz_dSV(-1) - dz_dSV(0) + 0.004
907 ! END IF
908 
909  zz_dsv = 0.
910  DO isl=-nsol,0
911  dzmisv(isl) = 0.500*(dz_dsv(isl) +dz_dsv(islmsv(isl)))
912  dziisv(isl) = 0.500* dz_dsv(isl) /dzmisv(isl)
913  dzi_sv(isl) = 0.500* dz_dsv(islmsv(isl))/dzmisv(isl)
914  dtz_sv(isl) = dt__sv /dz_dsv(isl)
915  dz78sv(isl) = 0.875* dz_dsv(isl)
916  dz34sv(isl) = 0.750* dz_dsv(isl)
917  dz_8sv(isl) = 0.125* dz_dsv(isl)
918  dzavsv(isl) = 0.125* dz_dsv(islmsv(isl)) &
919  & + 0.750* dz_dsv(isl) &
920  & + 0.125* dz_dsv(islpsv(isl))
921 !c #ER dz78SV(isl) = dz_dSV(isl)
922 !c #ER dz34SV(isl) = dz_dSV(isl)
923 !c #ER dz_8SV(isl) = 0.
924 !c #ER dzAvSV(isl) = dz_dSV(isl)
925  zz_dsv = zz_dsv+dz_dsv(isl)
926  END DO
927  DO ikl=1,knon !v
928  dzsnsv(ikl,0) = dz_dsv(0)
929  END DO
930 
931 !C +--Conversion to a 50 m Swab Ocean Discretization
932 !C + ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
933  ocndsv = 0.
934  DO isl=-nsol,0
935  ocndsv = ocndsv +dz_dsv(isl)
936  END DO
937  ocndsv = 50. /ocndsv
938 
939 
940 !C +--Secondary Vegetation Parameters
941 !C + -------------------------------
942 !
943 !C +--Minimum Stomatal Resistance (Hapex Sahel Data)
944 !C + (Taylor et al. 1997, J.Hydrol 188-189, p.1047)
945 !C + ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
946 !c #SH DO ivg=1,3 !
947 !c #SH StodSV(ivg) = 210. ! Millet
948 !c #SH END DO !
949 !c #SH StodSV( 4) = 120. ! Sparse Tiger Bush
950 !c #SH DO ivg=5,6 !
951 !c #SH StodSV(ivg) = 80. ! Dense Tiger Bush
952 !c #SH END DO !
953 !c #SH StodSV( 7) = 80. ! Low Trees (Fallow)
954 !c #SH StodSV( 10) = 80. !
955 !
956 !C +--Minimum Stomatal Resistance (Tropical Forest)
957 !C + ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
958 !c #SH StodSV( 8) = 60. ! Medium Trees
959 !c #SH StodSV( 11) = 60. !
960 !c #SH StodSV( 9) = 40. ! High Trees
961 !c #SH StodSV( 12) = 40. !
962 
963 !C +--Root Fraction
964 !C + ^^^^^^^^^^^^^
965 !C + * GENERAL REFERENCE
966 !C + Jackson et al., 1996: A global analysis of root distributions for
967 !C + terrestrial biomes. In Oecologia, 108, 389-411.
968 
969 !C + * ROOT PROFILE
970 !C + The cumulative root fraction Y is given by
971 !C + Y = 1 - beta**d with d the depth (in cm),
972 !C + beta a coefficient (vegetation dependent).
973 
974 !C + * BETA VALUES (for 11 world biomes)
975 !C + 1 boreal forest 0.943
976 !C + 2 crops 0.961
977 !C + 3 desert 0.975
978 !C + 4 sclerophyllous shrubs 0.964
979 !C + 5 temperate coniferous forest 0.976
980 !C + 6 temperate deciduous forest 0.966
981 !C + 7 temperate grassland 0.943
982 !C + 8 tropical deciduous forest 0.961
983 !C + 9 tropical evergreen forest 0.962
984 !C + 10 tropical grassland savanna 0.972
985 !C + 11 tundra 0.914
986 !
987 !C + * ADVISED BETA VALUES FOR MAR
988 !C + (see 'block data SISVAT_dat', variable rbtdSV)
989 !C +
990 !C + SVAT veg. type default West Africa
991 !C + 0 barren soil 0.000 0.000
992 !C + 1 crops low 0.961 (2) 0.961 (2)
993 !C + 2 crops medium 0.961 (2) 0.961 (2)
994 !C + 3 crops high 0.961 (2) 0.961 (2)
995 !C + 4 grass low 0.943 (7) 0.943 (7)
996 !C + 5 grass medium 0.943 (7) 0.964 (4)
997 !C + 6 grass high 0.943 (7) 0.972 (10)
998 !C + 7 broadleaf low 0.966 (6) 0.968 (4,10)
999 !C + 8 broadleaf medium 0.966 (6) 0.962 (8,9)
1000 !C + 9 broadleaf high 0.966 (6) 0.962 (8,9)
1001 !C + 10 needleleaf low 0.976 (5) 0.971 (5,6)
1002 !C + 11 needleleaf medium 0.976 (5) 0.976 (5)
1003 !C + 12 needleleaf high 0.976 (5) 0.976 (5)
1004 
1005 !C + Numbers between brackets refer to Jackson's biomes. For more details
1006 !C + about some choices, see the correspondance between the IGBP and SVAT
1007 !C + vegetation classes (i.e. in NESTOR).
1008 
1009 !C + * WARNING
1010 !C + Most of the roots are located in the first 2 m of soil. The root
1011 !C + fraction per layer depends on the definition of the soil layer
1012 !C + thickness. It will get wrong if a thick layer is defined around 2 m
1013 !C + deep.
1014 
1015  write(*,'(/a)') 'ROOT PROFILES (Jackson, 1996) :'
1016 
1017  DO ivt = 0, nvgt
1018  zdepth = 0.
1019  DO isl = 0, -nsol, -1
1020  IF (ivt .ne. 0) THEN
1021  rf__sv(ivt,isl) = rbtdsv(ivt)**zdepth * &
1022  & (1. - rbtdsv(ivt)**(dz_dsv(isl)*100) )
1023  zdepth = zdepth + dz_dsv(isl)*100 !in cm
1024  ELSE
1025  rf__sv(ivt,isl) = 0.
1026  END IF
1027  END DO
1028  write(*,'(a,i2,a,i3,a,99f10.5:)') &
1029  & ' RF__SV(', ivt, ',', -nsol, ':0) =', rf__sv(ivt,:)
1030  END DO
1031 ! write(6, format( &
1032 ! & ' NOTE: If root fraction is not close to 0 around 2 m deep,', &
1033 ! &/,' Then you should redefine the soil layer thicknesses.', &
1034 ! &/,' See the code for more details.'))
1035 
1036 
1037 !C +--Secondary Soil Parameters
1038 !C + -------------------------------
1039 
1040  DO ist=0,nsot
1041  rocssv(ist)=(1.0-etadsv(ist))*1.2e+6 ! Soil Contrib. to (ro c)_s
1042  s1__sv(ist)= bchdsv(ist) & ! Factor of (eta)**(b+2)
1043  & *psidsv(ist) *ks_dsv(ist) & ! in DR97, Eqn.(3.36)
1044  & /(etadsv(ist)**( bchdsv(ist)+3.)) !
1045  s2__sv(ist)= ks_dsv(ist) & ! Factor of (eta)**(2b+3)
1046  & /(etadsv(ist)**(2.*bchdsv(ist)+3.)) ! in DR97, Eqn.(3.35)
1047 
1048 !C +--Soil Minimum Humidity (from a prescribed minimum relative Humidity)
1049 !C + ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
1050  psimax = -(log(rhsmin))/7.2e-5 ! DR97, Eqn 3.15 Inversion
1051  etamsv(ist) = etadsv(ist) &
1052  & *(psimax/psidsv(ist))**(-min(10.,1./bchdsv(ist)))
1053  END DO
1054  etamsv(12) = 0.
1055 
1056 !C +--Piecewise Hydraulic Conductivity Profiles
1057 !C + ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
1058  DO ist=0,nsot
1059 !c #WR write(6, format(' Type | etaSat | No | eta__1 | eta__2 |', &
1060 !c #WR& ' Khyd_1 | Khyd_x | Khyd_2 | Khyd_y |' &
1061 !c #WR& /,' -----+-----------+----+-----------+-----------+', &
1062 !c #WR& '-----------+-----------+-----------+-----------+'))
1063 
1064  d__eta = etadsv(ist)/nkhy
1065  eta__1 = 0.
1066  eta__2 = d__eta
1067  DO ikh=0,nkhy
1068  khyd_1 = s2__sv(ist) & ! DR97, Eqn.(3.35)
1069  & *(eta__1 **(2. *bchdsv(ist)+3.)) !
1070  khyd_2 = s2__sv(ist) &!
1071  & *(eta__2 **(2. *bchdsv(ist)+3.)) !
1072 
1073  a_khyd = (khyd_2-khyd_1)/d__eta !
1074  b_khyd = khyd_1-a_khyd *eta__1 !
1075 !c #WR Khyd_x = a_Khyd*eta__1 +b_Khyd !
1076 !c #WR Khyd_y = a_Khyd*eta__2 +b_Khyd !
1077  akdtsv(ist,ikh) = a_khyd * dt__sv !
1078  bkdtsv(ist,ikh) = b_khyd * dt__sv !
1079 !c #WR write(6,format(i5,' |',e10.2,' |',i3,' |', 6(e10.2,' |')) &
1080 !c #WR& ist,etadSV(ist),ikh,eta__1, &
1081 !c #WR& eta__2,Khyd_1,Khyd_x,Khyd_2,Khyd_y
1082 
1083  eta__1 = eta__1 + d__eta
1084  eta__2 = eta__2 + d__eta
1085  END DO
1086  END DO
1087 
1088 
1089 !c
1090 ! ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
1091 ! + INITIALISATION: ARRAYS on grid +++
1092 ! + -----------------------------------
1093 ! +
1094  DO ikl=1,knon
1095 
1096 ! + Water Fluxes, roughness length
1097  us__sv(ikl) = 0.25 ! Frict. Velocity
1098  ! (approx ETH-Camp Lefebre etal CD 2005
1099  uts_sv(ikl) = 1. ! u*T* arbitrary
1100  uqs_sv(ikl) = 0.25 ! turb_vel(ikl,3) ! u*q* "
1101  uss_sv(ikl) = 0.25 ! turb_vel(ikl,4) ! u*s* "
1102 !!c #AE usthSV(ikl) = 0.
1103  z0h_sv(ikl) = 0. ! rlength(ikl,9)
1104  z0m_sv(ikl) = 0. ! Moment.Roughn.L.
1105  z0h_sv(ikl) = 0. ! Heat Roughn.L.
1106 !!c #OR Z0roSV(ikl) = 0. ! Orogr. Roughn.L.
1107  lmo_sv(ikl) = 0.
1108  vvs_sv(ikl) = 0.
1109  rrs_sv(ikl) = 0.
1110  dds_sv(ikl) = 0.
1111 
1112 ! + Snow
1113  isnosv(ikl) = 0 ! # snow layers
1114  bufssv(ikl) = 0.
1115  zzsnsv(ikl,:) = 0. ! Snow pack thickness
1116  qsnosv(ikl) = 0. ! BL snow content
1117  zwecsv(ikl) = 0.
1118  dbs_sv(ikl) = 0.
1119  dsnbsv(ikl) = 0.
1120  esnbsv(ikl) = 0.
1121  brossv(ikl) = 0.
1122  bg1ssv(ikl) = 0.
1123  bg2ssv(ikl) = 0.
1124  sws_sv(ikl) = 0.
1125  rnofsv(ikl) = 0. ! RunOFF Intensity
1126  toicsv(ikl) = 0. !hj preli
1127 ! + Clouds !hj1602
1128  cld_sv(ikl) = 0.
1129 
1130 ! + Vegetation
1131  lsmask(ikl) = 1 ! Land/Sea Mask
1132  isotsv(ikl) = 12 ! Soil Type
1133  iwafsv(ikl) = 1 ! Soil Drainage
1134  rootsv(ikl,:) = 0.
1135  ivgtsv(ikl) = 0 ! Vegetation Type
1136  lai0sv(ikl) = 0 ! LAI
1137  glf0sv(ikl) = 0 ! Green Leaf Frac.
1138 ! TVegSV(ikl) = 0.
1139  sncasv(ikl) = 0.
1140  rrcasv(ikl) = 0.
1141  psivsv(ikl) = 0.
1142 
1143 ! + z0(Sastrugi)
1144 ! #SZ Z0SaSL(ikl) = 0.
1145 ! #ZM DO nt=1,ntavSL
1146 ! #ZM SLn_z0(ikl,nt) = 0.5e-6
1147 ! #ZM SLn_b0(ikl,nt) = 0.5e-6
1148 ! #ZM SLn_r0(ikl,nt) = 0.5e-6
1149 ! #ZM END DO
1150 
1151 
1152 ! +--z0(Orography Roughness)
1153 ! #OR SL_z0 (ikl) = min(SL_z0 (ikl),zsigma/3.)
1154 ! #OR SLzoro(ikl) = min(SLzoro(ikl),zsigma/3.)
1155 
1156 ! +--SBL Characteristics
1157 ! #AA DO nt= 1,ntaver
1158 ! #AW V_0aSL(ikl,nt) = ssvSL(ikl)
1159 ! #AH dT0aSL(ikl,nt) = temp_air(ikl)-tsurf(ikl)
1160 ! #AA END DO
1161 
1162 
1163  END DO
1164 
1165 
1166 
1167  return
1168 
1169  END SUBROUTINE sisvat_ini
1170 
1171 
1172 
1173 
1174 
1175 
1176 
1177 !***************************************************************************
1178 
1179  SUBROUTINE sisvatetat0 (fichnom,ikl2i)
1180 
1181  USE dimphy
1182  USE mod_grid_phy_lmdz
1183  USE mod_phys_lmdz_para
1184 
1185  USE iostart
1186  USE var_sv
1187  USE vardsv
1188  USE varxsv
1189  USE vartsv
1190  USE indice_sol_mod
1191 
1192  IMPLICIT none
1193 !======================================================================
1194 ! Auteur(s) HJ PUNGE (LSCE) date: 07/2009
1195 ! Objet: Ecriture de l'etat de redemarrage pour SISVAT
1196 !======================================================================
1197  include "netcdf.inc"
1198 ! include "indicesol.h"
1199 
1200  include "dimsoil.h"
1201  include "clesphys.h"
1202  include "thermcell.h"
1203  include "compbl.h"
1204 
1205 !======================================================================
1206  CHARACTER(LEN=*) :: fichnom
1207 
1208 
1209  INTEGER, DIMENSION(klon), INTENT(IN) :: ikl2i
1210  REAL, DIMENSION(klon) :: rlon
1211  REAL, DIMENSION(klon) :: rlat
1212 
1213 ! les variables globales ecrites dans le fichier restart
1214  REAL, DIMENSION(klon) :: isno
1215  REAL, DIMENSION(klon) :: ispi
1216  REAL, DIMENSION(klon) :: iice
1217  REAL, DIMENSION(klon) :: rusn
1218  REAL, DIMENSION(klon, nsno) :: isto
1219 
1220  REAL, DIMENSION(klon, nsismx) :: Tsis
1221  REAL, DIMENSION(klon, nsismx) :: eta
1222  REAL, DIMENSION(klon, nsismx) :: ro
1223 
1224  REAL, DIMENSION(klon, nsno) :: dzsn
1225  REAL, DIMENSION(klon, nsno) :: G1sn
1226  REAL, DIMENSION(klon, nsno) :: G2sn
1227  REAL, DIMENSION(klon, nsno) :: agsn
1228 
1229  REAL, DIMENSION(klon) :: toic
1230 
1231 ! REAL, DIMENSION(klon) :: IRs
1232 ! REAL, DIMENSION(klon) :: LMO
1233 ! REAL, DIMENSION(klon) :: Bufs
1234 ! REAL, DIMENSION(klon, 9) :: rlength
1235 ! REAL, DIMENSION(klon, 5) :: turb_vel
1236 
1237  INTEGER :: isl, ikl, i, isn , errT, erreta, errro, errdz, snopts
1238  CHARACTER (len=2) :: str2
1239  LOGICAL :: found
1240 
1241  errt=0
1242  errro=0
1243  erreta=0
1244  errdz=0
1245  snopts=0
1246 ! Ouvrir le fichier contenant l'etat initial:
1247 
1248  CALL open_startphy(fichnom)
1249 
1250 ! Lecture des latitudes, longitudes (coordonnees):
1251 
1252  CALL get_field("latitude",rlat,found)
1253  CALL get_field("longitude",rlon,found)
1254 
1255  CALL get_field("n_snows", isno,found)
1256  IF (.NOT. found) THEN
1257  print*, 'phyetat0: Le champ <n_snows> est absent'
1258  print *, 'fichier startsisvat non compatible avec sisvatetat0'
1259  ENDIF
1260 
1261  CALL get_field("n_ice_top",ispi,found)
1262  CALL get_field("n_ice",iice,found)
1263  CALL get_field("surf_water",rusn,found)
1264 ! IF (.NOT. found) THEN
1265 ! PRINT*, 'phyetat0: Le champ <surf_water> est absent'
1266 ! rusn(:)=0.
1267 ! ENDIF
1268 
1269 
1270  CALL get_field("to_ice",toic,found)
1271  IF (.NOT. found) THEN
1272  print*, 'phyetat0: Le champ <to_ice> est absent'
1273  toic(:)=0.
1274  ENDIF
1275 
1276 
1277 ! CALL get_field("IR_soil",IRs,found)
1278 ! CALL get_field("LMO",LMO,found)
1279 ! CALL get_field("snow_buffer",Bufs,found)
1280 ! DO i = 1, 5
1281 ! WRITE(str2,'(i2.2)') i
1282 ! CALL get_field("turb_veloc"//str2, &
1283 ! turb_vel(:,i),found)
1284 ! ENDDO
1285 ! DO i = 1, 9
1286 ! WRITE(str2,'(i2.2)') i
1287 ! CALL get_field("rough_length"//str2, &
1288 ! rlength(:,i),found)
1289 ! ENDDO
1290  DO isn = 1,nsno
1291  IF (isn.LE.99) THEN
1292  WRITE(str2,'(i2.2)') isn
1293  CALL get_field("AGESNOW"//str2, &
1294  agsn(:,isn),found)
1295  ELSE
1296  print*, "Trop de couches"
1297  CALL abort
1298  ENDIF
1299  ENDDO
1300  DO isn = 1,nsno
1301  IF (isn.LE.99) THEN
1302  WRITE(str2,'(i2.2)') isn
1303  CALL get_field("DZSNOW"//str2, &
1304  dzsn(:,isn),found)
1305  ELSE
1306  print*, "Trop de couches"
1307  CALL abort
1308  ENDIF
1309  ENDDO
1310  DO isn = 1,nsno
1311  IF (isn.LE.99) THEN
1312  WRITE(str2,'(i2.2)') isn
1313  CALL get_field("G2SNOW"//str2, &
1314  g2sn(:,isn),found)
1315  ELSE
1316  print*, "Trop de couches"
1317  CALL abort
1318  ENDIF
1319  ENDDO
1320  DO isn = 1,nsno
1321  IF (isn.LE.99) THEN
1322  WRITE(str2,'(i2.2)') isn
1323  CALL get_field("G1SNOW"//str2, &
1324  g1sn(:,isn),found)
1325  ELSE
1326  print*, "Trop de couches"
1327  CALL abort
1328  ENDIF
1329  ENDDO
1330  DO isn = 1,nsismx
1331  IF (isn.LE.99) THEN
1332  WRITE(str2,'(i2.2)') isn
1333  CALL get_field("ETA"//str2, &
1334  eta(:,isn),found)
1335  ELSE
1336  print*, "Trop de couches"
1337  CALL abort
1338  ENDIF
1339  ENDDO
1340  DO isn = 1,nsismx
1341  IF (isn.LE.99) THEN
1342  WRITE(str2,'(i2.2)') isn
1343  CALL get_field("RO"//str2, &
1344  ro(:,isn),found)
1345  ELSE
1346  print*, "Trop de couches"
1347  CALL abort
1348  ENDIF
1349  ENDDO
1350  DO isn = 1,nsismx
1351  IF (isn.LE.99) THEN
1352  WRITE(str2,'(i2.2)') isn
1353  CALL get_field("TSS"//str2, &
1354  tsis(:,isn),found)
1355  ELSE
1356  print*, "Trop de couches"
1357  CALL abort
1358  ENDIF
1359  ENDDO
1360  DO isn = 1,nsno
1361  IF (isn.LE.99) THEN
1362  WRITE(str2,'(i2.2)') isn
1363  CALL get_field("HISTORY"//str2, &
1364  isto(:,isn),found)
1365  ELSE
1366  print*, "Trop de couches"
1367  CALL abort
1368  ENDIF
1369  ENDDO
1370  write(*,*)'Read ',fichnom,' finished!!'
1371 
1372 !*********************************************************************************
1373 ! Compress restart file variables for SISVAT
1374 
1375 
1376  isnosv(:) = 0
1377  ispisv(:) = 0
1378  iicesv(:) = 0
1379 
1380  eta_sv(:,1:nsno) = 0.
1381  tsissv(:,1:nsno) = 0.
1382  istosv(:,:) = 0
1383  ro__sv(:,1:nsno) = 0.
1384  g1snsv(:,:) = 0.
1385  g2snsv(:,:) = 0.
1386  dzsnsv(:,1:nsno) = 0.
1387  agsnsv(:,:) = 0.
1388  rusnsv(:) = 0.
1389  toicsv(:) = 0.
1390 
1391  DO ikl = 1,klon
1392  i = ikl2i(ikl)
1393  IF (i > 0) THEN
1394  isnosv(ikl) = int(isno(i)) ! Nb Snow/Ice Lay.
1395  ispisv(ikl) = int(ispi(i)) ! Nb Supr.Ice Lay.
1396  iicesv(ikl) = int(iice(i)) ! Nb Ice Lay.
1397 
1398 !! IRs_SV(ikl) = IRs(i)
1399 ! LMO_SV(ikl) = LMO(i) !?
1400 ! us__SV(ikl) = turb_vel(i,1)
1401 ! uts_SV(ikl) = turb_vel(i,2)
1402 ! uqs_SV(ikl) = turb_vel(i,3)
1403 ! uss_SV(ikl) = turb_vel(i,4)
1404 ! usthSV(ikl) = turb_vel(i,5)
1405 ! Z0m_SV(ikl) = rlength(i,1) ! Moment.Roughn.L.
1406 ! Z0mmSV(ikl) = rlength(i,2) ! Moment.Roughn.L.
1407 ! Z0mnSV(ikl) = rlength(i,3) ! Moment.Roughn.L.
1408 ! Z0SaSV(ikl) = rlength(i,4) ! Moment.Roughn.L.
1409 ! Z0e_SV(ikl) = rlength(i,5) ! Moment.Roughn.L.
1410 ! Z0emSV(ikl) = rlength(i,6) ! Moment.Roughn.L.
1411 ! Z0enSV(ikl) = rlength(i,7) ! Moment.Roughn.L.
1412 ! Z0roSV(ikl) = rlength(i,8) ! Moment.Roughn.L.
1413 !! Z0h_SV(ikl) = rlength(i,9) ! Moment.Roughn.L.
1414 
1415  DO isl = -nsol,0
1416  ro__sv(ikl,isl) = ro(i,nsno+1-isl) !
1417  eta_sv(ikl,isl) = eta(i,nsno+1-isl) ! Soil Humidity
1418 !hjp 15/10/2010
1419  IF (eta_sv(ikl,isl) <= 1.e-6) THEN !hj check
1420  eta_sv(ikl,isl) = 1.e-6
1421  ENDIF
1422  tsissv(ikl,isl) = tsis(i,nsno+1-isl) ! Soil Temperature
1423  IF (tsissv(ikl,isl) <= 1.) THEN !hj check
1424 ! errT=errT+1
1425  tsissv(ikl,isl) = 273.15
1426  ENDIF
1427 
1428  END DO
1429  write(*,*)'Copy histo', ikl
1430 
1431 
1432  istosv(ikl,0) = 0 ! Snow History
1433 
1434  g1snsv(ikl,0) = 0.5 ! [-]
1435  g2snsv(ikl,0) = 3. ! [-] [0.0001 m]
1436  dzsnsv(ikl,0) = dz_dsv(0) ! [m]
1437  agsnsv(ikl,0) = 0. ! [day]
1438 
1439  DO isn = 1,isnosv(ikl) !nsno
1440  snopts=snopts+1
1441  IF (isto(i,isn) > 10.) THEN !hj check
1442  write(*,*)'Irregular isto',ikl,i,isn,isto(i,isn)
1443  isto(i,isn) = 1.
1444  ENDIF
1445 
1446  istosv(ikl,isn) = int(isto(i,isn)) ! Snow History
1447  ro__sv(ikl,isn) = ro(i,isn) ! [kg/m3]
1448  eta_sv(ikl,isn) = eta(i,isn) ! [m3/m3]
1449  tsissv(ikl,isn) = tsis(i,isn) ! [K]
1450 
1451  IF (tsissv(ikl,isn) <= 1.) THEN !hj check
1452  errt=errt+1
1453  tsissv(ikl,isn) = tsissv(ikl,0)
1454  ENDIF
1455  IF (tsissv(ikl,isn) <= 1.) THEN !hj check
1456  tsissv(ikl,isn) = 263.15
1457  ENDIF
1458  IF (eta_sv(ikl,isn) < 1.e-9) THEN !hj check
1459  eta_sv(ikl,isn) = 1.e-6
1460  erreta=erreta+1
1461  ENDIF
1462  IF (ro__sv(ikl,isn) <= 10.) THEN !hj check
1463  ro__sv(ikl,isn) = 11.
1464  errro=errro+1
1465  ENDIF
1466  write(*,*)ikl,i,isn,tsis(i,isn),g1sn(i,isn)
1467  g1snsv(ikl,isn) = g1sn(i,isn) ! [-] [-]
1468  g2snsv(ikl,isn) = g2sn(i,isn) ! [-] [0.0001 m]
1469  dzsnsv(ikl,isn) = dzsn(i,isn) ! [m]
1470 ! IF (dzsnSV(ikl,isn) < 5.e-6) THEN !hj check
1471 ! dzsnSV(ikl,isn) = 0.000005
1472 ! errdz=errdz+1
1473 ! ENDIF
1474  agsnsv(ikl,isn) = agsn(i,isn) ! [day]
1475  END DO
1476  rusnsv(ikl) = rusn(i) ! Surficial Water
1477  toicsv(ikl) = toic(i) ! bilan snow to ice
1478 !! BufsSV(ikl) = Bufs(i)
1479  END IF
1480  END DO
1481 ! write(*,*)snopts,'snow pts',errT,' T errors,',errro,'ro,'
1482 ! write(*,*)' ',erreta,'eta,',errdz,'dz'
1483  END SUBROUTINE sisvatetat0
1484 
1485 
1486 
1487 
1488  SUBROUTINE sisvatredem (fichnom,ikl2i,rlon,rlat)
1489 !======================================================================
1490 ! Auteur(s) HJ PUNGE (LSCE) date: 07/2009
1491 ! Objet: Ecriture de l'etat de redemarrage pour SISVAT
1492 !======================================================================
1493  USE mod_grid_phy_lmdz
1494  USE mod_phys_lmdz_para
1495  USE iostart
1496  USE var_sv
1497  USE varxsv
1498  USE varysv !hj tmp 12 03 2010
1499  USE vartsv
1500  USE indice_sol_mod
1501 
1502  IMPLICIT none
1503 
1504  include "netcdf.inc"
1505 ! include "indicesol.h"
1506  include "dimsoil.h"
1507  include "clesphys.h"
1508  include "thermcell.h"
1509  include "compbl.h"
1510 
1511 !======================================================================
1512 
1513  CHARACTER(LEN=*) :: fichnom
1514  INTEGER, DIMENSION(klon), INTENT(IN) :: ikl2i
1515  REAL, DIMENSION(klon), INTENT(IN) :: rlon
1516  REAL, DIMENSION(klon), INTENT(IN) :: rlat
1517 
1518 ! les variables globales ecrites dans le fichier restart
1519  REAL, DIMENSION(klon) :: isno
1520  REAL, DIMENSION(klon) :: ispi
1521  REAL, DIMENSION(klon) :: iice
1522  REAL, DIMENSION(klon, nsnowmx) :: isto
1523 
1524  REAL, DIMENSION(klon, nsismx) :: Tsis
1525  REAL, DIMENSION(klon, nsismx) :: eta
1526  REAL, DIMENSION(klon, nsnowmx) :: dzsn
1527  REAL, DIMENSION(klon, nsismx) :: ro
1528  REAL, DIMENSION(klon, nsnowmx) :: G1sn
1529  REAL, DIMENSION(klon, nsnowmx) :: G2sn
1530  REAL, DIMENSION(klon, nsnowmx) :: agsn
1531  REAL, DIMENSION(klon) :: IRs
1532  REAL, DIMENSION(klon) :: LMO
1533  REAL, DIMENSION(klon) :: rusn
1534  REAL, DIMENSION(klon) :: toic
1535  REAL, DIMENSION(klon) :: Bufs
1536  REAL, DIMENSION(klon) :: alb1,alb2,alb3
1537  REAL, DIMENSION(klon, 9) :: rlength
1538  REAL, DIMENSION(klon, 5) :: turb_vel
1539 
1540  INTEGER isl, ikl, i, isn
1541  CHARACTER (len=2) :: str2
1542 
1543  isno(:) = 0
1544  ispi(:) = 0
1545  iice(:) = 0
1546  irs(:) = 0.
1547  lmo(:) = 0.
1548  turb_vel(:,:) = 0.
1549  rlength(:,:) = 0.
1550  eta(:,:) = 0.
1551  tsis(:,:) = 0.
1552  isto(:,:) = 0
1553  ro(:,:) = 0.
1554  g1sn(:,:) = 0.
1555  g2sn(:,:) = 0.
1556  dzsn(:,:) = 0.
1557  agsn(:,:) = 0.
1558  rusn(:) = 0.
1559  toic(:) = 0.
1560  bufs(:) = 0.
1561  alb1(:) = 0.
1562  alb2(:) = 0.
1563  alb3(:) = 0.
1564 
1565 !***************************************************************************
1566 ! Uncompress SISVAT output variables for storage
1567 
1568  DO ikl = 1,klon
1569  i = ikl2i(ikl)
1570  IF (i > 0) THEN
1571  isno(i) = 1.*isnosv(ikl) ! Nb Snow/Ice Lay.
1572  ispi(i) = 1.*ispisv(ikl) ! Nb Supr.Ice Lay.
1573  iice(i) = 1.*iicesv(ikl) ! Nb Ice Lay.
1574 
1575 ! IRs(i) = IRs_SV(ikl)
1576 ! LMO(i) = LMO_SV(ikl)
1577 ! turb_vel(i,1) = us__SV(ikl)
1578 ! turb_vel(i,2) = uts_SV(ikl)
1579 ! turb_vel(i,3) = uqs_SV(ikl)
1580 ! turb_vel(i,4) = uss_SV(ikl)
1581 ! turb_vel(i,5) = usthSV(ikl)
1582 ! rlength(i,1) = Z0m_SV(ikl) ! Moment.Roughn.L.
1583 ! rlength(i,2) = Z0mmSV(ikl) !
1584 ! rlength(i,3) = Z0mnSV(ikl) !
1585 ! rlength(i,4) = Z0SaSV(ikl) !
1586 ! rlength(i,5) = Z0e_SV(ikl) !
1587 ! rlength(i,6) = Z0emSV(ikl) !
1588 ! rlength(i,7) = Z0enSV(ikl) !
1589 ! rlength(i,8) = Z0roSV(ikl) !
1590 ! rlength(i,9) = Z0h_SV(ikl) !
1591 
1592  DO isl = -nsol,0 !
1593  eta(i,nsno+1-isl) = eta_sv(ikl,isl) ! Soil Humidity
1594  tsis(i,nsno+1-isl) = tsissv(ikl,isl) ! Soil Temperature
1595  ro(i,nsno+1-isl) = ro__sv(ikl,isl) ! [kg/m3]
1596  END DO
1597 
1598 
1599  DO isn = 1,nsno
1600  isto(i,isn) = 1.*istosv(ikl,isn) ! Snow History
1601  ro(i,isn) = ro__sv(ikl,isn) ! [kg/m3]
1602  eta(i,isn) = eta_sv(ikl,isn) ! [m3/m3]
1603  tsis(i,isn) = tsissv(ikl,isn) ! [K]
1604  g1sn(i,isn) = g1snsv(ikl,isn) ! [-] [-]
1605  g2sn(i,isn) = g2snsv(ikl,isn) ! [-] [0.0001 m]
1606  dzsn(i,isn) = dzsnsv(ikl,isn) ! [m]
1607  agsn(i,isn) = agsnsv(ikl,isn) ! [day]
1608  END DO
1609  rusn(i) = rusnsv(ikl) ! Surficial Water
1610  toic(i) = toicsv(ikl) ! Surficial Water
1611  alb1(i) = alb1sv(ikl)
1612  alb2(i) = alb2sv(ikl)
1613  alb3(i) = alb3sv(ikl)
1614 ! Bufs(i) = BufsSV(ikl)
1615  END IF
1616  END DO
1617 
1618  CALL open_restartphy(fichnom)
1619  CALL put_field("longitude", &
1620  "Longitudes de la grille physique",rlon)
1621  CALL put_field("latitude","Latitudes de la grille physique",rlat)
1622 
1623  CALL put_field("n_snows", "number of snow/ice layers",isno)
1624  CALL put_field("n_ice_top", "number of top ice layers",ispi)
1625  CALL put_field("n_ice", "number of ice layers",iice)
1626  CALL put_field("IR_soil", "Soil IR flux",irs)
1627  CALL put_field("LMO", "Monin-Obukhov Scale",lmo)
1628  CALL put_field("surf_water", "Surficial water",rusn)
1629  CALL put_field("snow_buffer", "Snow buffer layer",bufs)
1630  CALL put_field("alb_1", "albedo sw",alb1)
1631  CALL put_field("alb_2", "albedo nIR",alb2)
1632  CALL put_field("alb_3", "albedo fIR",alb3)
1633  CALL put_field("to_ice", "Snow passed to ice",toic)
1634 
1635  ! DO i = 1, 5
1636  ! WRITE(str2,'(i2.2)') i
1637  ! CALL put_field("turb_veloc"//str2, &
1638  ! "various turbulent velocities"//str2, &
1639  ! turb_vel(:,i))
1640  ! ENDDO
1641  ! DO i = 1, 9
1642  ! WRITE(str2,'(i2.2)') i
1643  ! CALL put_field("rough_length"//str2, &
1644  ! "various roughness lengths"//str2, &
1645  ! rlength(:,i))
1646  ! ENDDO
1647  DO isn = 1,nsno
1648  IF (isn.LE.99) THEN
1649  WRITE(str2,'(i2.2)') isn
1650  CALL put_field("AGESNOW"//str2, &
1651  "Age de la neige layer No."//str2, &
1652  agsn(:,isn))
1653  ELSE
1654  print*, "Trop de couches"
1655  CALL abort
1656  ENDIF
1657  ENDDO
1658  DO isn = 1,nsno
1659  IF (isn.LE.99) THEN
1660  WRITE(str2,'(i2.2)') isn
1661  CALL put_field("DZSNOW"//str2, &
1662  "Snow/ice thickness layer No."//str2, &
1663  dzsn(:,isn))
1664  ELSE
1665  print*, "Trop de couches"
1666  CALL abort
1667  ENDIF
1668  ENDDO
1669  DO isn = 1,nsno
1670  IF (isn.LE.99) THEN
1671  WRITE(str2,'(i2.2)') isn
1672  CALL put_field("G2SNOW"//str2, &
1673  "Snow Property 2, layer No."//str2, &
1674  g2sn(:,isn))
1675  ELSE
1676  print*, "Trop de couches"
1677  CALL abort
1678  ENDIF
1679  ENDDO
1680  DO isn = 1,nsno
1681  IF (isn.LE.99) THEN
1682  WRITE(str2,'(i2.2)') isn
1683  CALL put_field("G1SNOW"//str2, &
1684  "Snow Property 1, layer No."//str2, &
1685  g1sn(:,isn))
1686  ELSE
1687  print*, "Trop de couches"
1688  CALL abort
1689  ENDIF
1690  ENDDO
1691  DO isn = 1,nsismx
1692  IF (isn.LE.99) THEN
1693  WRITE(str2,'(i2.2)') isn
1694  CALL put_field("ETA"//str2, &
1695  "Soil/snow water content layer No."//str2, &
1696  eta(:,isn))
1697  ELSE
1698  print*, "Trop de couches"
1699  CALL abort
1700  ENDIF
1701  ENDDO
1702  DO isn = 1,nsismx !nsno
1703  IF (isn.LE.99) THEN
1704  WRITE(str2,'(i2.2)') isn
1705  CALL put_field("RO"//str2, &
1706  "Snow density layer No."//str2, &
1707  ro(:,isn))
1708  ELSE
1709  print*, "Trop de couches"
1710  CALL abort
1711  ENDIF
1712  ENDDO
1713  DO isn = 1,nsismx
1714  IF (isn.LE.99) THEN
1715  WRITE(str2,'(i2.2)') isn
1716  CALL put_field("TSS"//str2, &
1717  "Soil/snow temperature layer No."//str2, &
1718  tsis(:,isn))
1719  ELSE
1720  print*, "Trop de couches"
1721  CALL abort
1722  ENDIF
1723  ENDDO
1724  DO isn = 1,nsno
1725  IF (isn.LE.99) THEN
1726  WRITE(str2,'(i2.2)') isn
1727  CALL put_field("HISTORY"//str2, &
1728  "Snow history layer No."//str2, &
1729  isto(:,isn))
1730  ELSE
1731  print*, "Trop de couches"
1732  CALL abort
1733  ENDIF
1734  ENDDO
1735 
1736  END SUBROUTINE sisvatredem
1737 
1738 END MODULE surf_sisvat_mod
vartsv::t__mem
real, dimension(:,:), allocatable, save t__mem
Definition: VARtSV.F90:40
vartsv::lambsv
real, save lambsv
Definition: VARtSV.F90:32
var0sv::akdtsv
real, dimension(0:nsot, 0:nkhy) akdtsv
Definition: VAR0SV.F90:28
varxsv::brossv
real, dimension(:), allocatable, save brossv
Definition: VARxSV.F90:35
var0sv::explic
real explic
Definition: VAR0SV.F90:13
vartsv::bcohsv
real, dimension(:), allocatable, save bcohsv
Definition: VARtSV.F90:20
vardsv::rbtdsv
real, dimension(0:nvgt), parameter rbtdsv
Definition: VARdSV.F90:65
surf_sisvat_mod::nsnowmx
integer, parameter nsnowmx
Definition: surf_sisvat_mod.F90:5
vartsv::vvmmem
real, dimension(:), allocatable, save vvmmem
Definition: VARtSV.F90:38
varxsv::tvegsv
real, dimension(:), allocatable, save tvegsv
Definition: VARxSV.F90:157
varxsv::dsdtsv
real, dimension(:), allocatable, save dsdtsv
Definition: VARxSV.F90:60
vartsv::dz1_sv
real, dimension(:,:), allocatable, save dz1_sv
Definition: VARtSV.F90:12
vartsv
Definition: VARtSV.F90:1
var0sv::rocssv
real, dimension(0:nsot) rocssv
Definition: VAR0SV.F90:24
varxsv
Definition: VARxSV.F90:1
varysv
Definition: VARySV.F90:1
vartsv::ps__sv
real, dimension(:), allocatable, save ps__sv
Definition: VARtSV.F90:24
vartsv::p1l_sv
real, dimension(:), allocatable, save p1l_sv
Definition: VARtSV.F90:26
vardsv
Definition: VARdSV.F90:1
varphy::ra
real, parameter ra
Definition: VARphy.F90:53
var0sv::dzavsv
real, dimension(-nsol:0) dzavsv
Definition: VAR0SV.F90:21
var0sv::dz_8sv
real, dimension(-nsol:0) dz_8sv
Definition: VAR0SV.F90:20
varxsv::exnrsv
real, dimension(:), allocatable, save exnrsv
Definition: VARxSV.F90:58
var0sv::etamsv
real, dimension(0:nsot) etamsv
Definition: VAR0SV.F90:25
varxsv::bufssv
real, dimension(:), allocatable, save bufssv
Definition: VARxSV.F90:174
varxsv::istosv
integer, dimension(:,:), allocatable, save istosv
Definition: VARxSV.F90:101
varxsv::sol_sv
real, dimension(:), allocatable, save sol_sv
Definition: VARxSV.F90:20
iostart::open_startphy
subroutine, public open_startphy(filename)
Definition: iostart.F90:32
dimphy::klon
integer, save klon
Definition: dimphy.F90:3
varxsv::agsnsv
real, dimension(:,:), allocatable, save agsnsv
Definition: VARxSV.F90:172
varxsv::dzsnsv
real, dimension(:,:), allocatable, save dzsnsv
Definition: VARxSV.F90:170
varysv::rootsv
real, dimension(:,:), allocatable, save rootsv
Definition: VARySV.F90:87
vartsv::rsolsv
real, dimension(:), allocatable, save rsolsv
Definition: VARtSV.F90:30
varxsv::lai0sv
real, dimension(:), allocatable, save lai0sv
Definition: VARxSV.F90:73
varxsv::isnosv
integer, dimension(:), allocatable, save isnosv
Definition: VARxSV.F90:95
varxsv::isotsv
integer, dimension(:), allocatable, save isotsv
Definition: VARxSV.F90:11
varxsv::dsnbsv
real, dimension(:), allocatable, save dsnbsv
Definition: VARxSV.F90:29
vartsv::bcoqsv
real, dimension(:), allocatable, save bcoqsv
Definition: VARtSV.F90:22
surf_sisvat_mod::nsismx
integer, parameter nsismx
Definition: surf_sisvat_mod.F90:6
varxsv::sncasv
real, dimension(:), allocatable, save sncasv
Definition: VARxSV.F90:151
varxsv::sws_sv
real, dimension(:), allocatable, save sws_sv
Definition: VARxSV.F90:180
varxsv::iicesv
integer, dimension(:), allocatable, save iicesv
Definition: VARxSV.F90:99
varxsv::g1snsv
real, dimension(:,:), allocatable, save g1snsv
Definition: VARxSV.F90:166
var0sv::dzmisv
real, dimension(-nsol:0) dzmisv
Definition: VAR0SV.F90:14
varxsv::tsissv
real, dimension(:,:), allocatable, save tsissv
Definition: VARxSV.F90:160
varxsv::uts_sv
real, dimension(:), allocatable, save uts_sv
Definition: VARxSV.F90:114
vartsv::toicsv
real, dimension(:), allocatable, save toicsv
Definition: VARtSV.F90:9
varxsv::init_varxsv
subroutine init_varxsv
Definition: VARxSV.F90:232
varxsv::tat_sv
real, dimension(:), allocatable, save tat_sv
Definition: VARxSV.F90:56
surf_sisvat_mod::get_soil_levels
subroutine get_soil_levels(dz1, dz2, lambda)
Definition: surf_sisvat_mod.F90:706
vartsv::acoqsv
real, dimension(:), allocatable, save acoqsv
Definition: VARtSV.F90:22
surf_sisvat_mod
Definition: surf_sisvat_mod.F90:1
vartsv::dz2_sv
real, dimension(:,:), allocatable, save dz2_sv
Definition: VARtSV.F90:14
var0sv::dzi_sv
real, dimension(-nsol:0) dzi_sv
Definition: VAR0SV.F90:15
varxsv::zwecsv
real, dimension(:), allocatable, save zwecsv
Definition: VARxSV.F90:187
var_sv
Definition: VAR_SV.F90:1
varxsv::psivsv
real, dimension(:), allocatable, save psivsv
Definition: VARxSV.F90:155
varxsv::cld_sv
real, dimension(:), allocatable, save cld_sv
Definition: VARxSV.F90:44
var0sv::bkdtsv
real, dimension(0:nsot, 0:nkhy) bkdtsv
Definition: VAR0SV.F90:29
vardsv::so2dsv
real, parameter so2dsv
Definition: VARdSV.F90:26
varphy::tfsnow
real, parameter tfsnow
Definition: VARphy.F90:93
iostart::get_field
Definition: iostart.F90:10
varphy
Definition: VARphy.F90:1
var0sv
Definition: VAR0SV.F90:1
false
!$Id itapm1 ENDIF!IM on interpole les champs sur les niveaux STD de pression!IM a chaque pas de temps de la physique c!positionnement de l argument logique a false c!pour ne pas recalculer deux fois la meme chose!c!a cet effet un appel a plevel_new a ete deplace c!a la fin de la serie d appels c!la boucle DO nlevSTD a ete internalisee c!dans d ou la creation de cette routine c c!CALL false
Definition: calcul_STDlev.h:26
varxsv::dds_sv
real, dimension(:), allocatable, save dds_sv
Definition: VARxSV.F90:54
varxsv::alb0sv
real, dimension(:), allocatable, save alb0sv
Definition: VARxSV.F90:78
varxsv::rusnsv
real, dimension(:), allocatable, save rusnsv
Definition: VARxSV.F90:176
varxsv::glf0sv
real, dimension(:), allocatable, save glf0sv
Definition: VARxSV.F90:75
vartsv::acohsv
real, dimension(:), allocatable, save acohsv
Definition: VARtSV.F90:20
varxsv::bg2ssv
real, dimension(:), allocatable, save bg2ssv
Definition: VARxSV.F90:39
var0sv::dziisv
real, dimension(-nsol:0) dziisv
Definition: VAR0SV.F90:16
sisvat
subroutine sisvat(SnoMod, BloMod, jjtime)
Definition: sisvat.F:2
surf_sisvat_mod::sisvatredem
subroutine sisvatredem(fichnom, ikl2i, rlon, rlat)
Definition: surf_sisvat_mod.F90:1489
varxsv::z0h_sv
real, dimension(:), allocatable, save z0h_sv
Definition: VARxSV.F90:144
varxsv::iwafsv
integer, dimension(:), allocatable, save iwafsv
Definition: VARxSV.F90:13
varxsv::ivgtsv
integer, dimension(:), allocatable, save ivgtsv
Definition: VARxSV.F90:15
var0sv::rf__sv
real, dimension(0:nvgt,-nsol:0) rf__sv
Definition: VAR0SV.F90:23
var0sv::dz78sv
real, dimension(-nsol:0) dz78sv
Definition: VAR0SV.F90:18
vardsv::ks_dsv
real, dimension(0:nsot), parameter ks_dsv
Definition: VARdSV.F90:106
vardsv::ini_sv
logical ini_sv
Definition: VARdSV.F90:11
var_sv::nsno
integer, parameter nsno
Definition: VAR_SV.F90:11
varysv::hss_sv
real, dimension(:), allocatable, save hss_sv
Definition: VARySV.F90:53
mod_phys_lmdz_para
Definition: mod_phys_lmdz_para.F90:4
varxsv::za__sv
real, dimension(:), allocatable, save za__sv
Definition: VARxSV.F90:46
varxsv::vv__sv
real, dimension(:), allocatable, save vv__sv
Definition: VARxSV.F90:48
vartsv::tsf_sv
real, dimension(:), allocatable, save tsf_sv
Definition: VARtSV.F90:16
varxsv::ispisv
integer, dimension(:), allocatable, save ispisv
Definition: VARxSV.F90:97
varxsv::qsnosv
real, dimension(:), allocatable, save qsnosv
Definition: VARxSV.F90:70
varxsv::rht_sv
real, dimension(:), allocatable, save rht_sv
Definition: VARxSV.F90:64
vartsv::dtmmem
real, dimension(:), allocatable, save dtmmem
Definition: VARtSV.F90:42
var0sv::s2__sv
real, dimension(0:nsot) s2__sv
Definition: VAR0SV.F90:27
var_sv::knonv
integer, save knonv
Definition: VAR_SV.F90:15
var0sv::ocndsv
real ocndsv
Definition: VAR0SV.F90:22
varxsv::rnofsv
real, dimension(:), allocatable, save rnofsv
Definition: VARxSV.F90:224
var_sv::nsol
integer, parameter nsol
Definition: VAR_SV.F90:10
yomcst_sisvat::rlvtt
real, parameter rlvtt
Definition: YOMCST_SISVAT.F90:43
var_sv::klonv
integer, save klonv
Definition: VAR_SV.F90:13
var0sv::islmsv
integer, dimension(-nsol:0) islmsv
Definition: VAR0SV.F90:10
vardsv::nvgt
integer, parameter nvgt
Definition: VARdSV.F90:48
varxsv::ird_sv
real, dimension(:), allocatable, save ird_sv
Definition: VARxSV.F90:22
vardsv::etadsv
real, dimension(0:nsot), parameter etadsv
Definition: VARdSV.F90:98
var0sv::dz34sv
real, dimension(-nsol:0) dz34sv
Definition: VAR0SV.F90:19
varxsv::dbs_sv
real, dimension(:), allocatable, save dbs_sv
Definition: VARxSV.F90:33
true
!$Id itapm1 ENDIF!IM on interpole les champs sur les niveaux STD de pression!IM a chaque pas de temps de la physique c!positionnement de l argument logique a false c!pour ne pas recalculer deux fois la meme chose!c!a cet effet un appel a plevel_new a ete deplace c!a la fin de la serie d appels c!la boucle DO nlevSTD a ete internalisee c!dans d ou la creation de cette routine c c!CALL ulevSTD CALL &zphi philevSTD CALL &zx_rh rhlevSTD!DO klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon du jour ou toutes les read_climoz CALL true
Definition: calcul_STDlev.h:139
vardsv::psidsv
real, dimension(0:nsot), parameter psidsv
Definition: VARdSV.F90:102
varxsv::drr_sv
real, dimension(:), allocatable, save drr_sv
Definition: VARxSV.F90:25
varxsv::ro__sv
real, dimension(:,:), allocatable, save ro__sv
Definition: VARxSV.F90:162
varxsv::irs_sv
real, dimension(:), allocatable, save irs_sv
Definition: VARxSV.F90:108
varxsv::lmo_sv
real, dimension(:), allocatable, save lmo_sv
Definition: VARxSV.F90:110
vartsv::cdh_sv
real, dimension(:), allocatable, save cdh_sv
Definition: VARtSV.F90:28
varxsv::z0m_sv
real, dimension(:), allocatable, save z0m_sv
Definition: VARxSV.F90:128
varysv::alb2sv
real, dimension(:), allocatable, save alb2sv
Definition: VARySV.F90:20
varxsv::rrs_sv
real, dimension(:), allocatable, save rrs_sv
Definition: VARxSV.F90:52
var0sv::islpsv
integer, dimension(-nsol:0) islpsv
Definition: VAR0SV.F90:8
varxsv::dsn_sv
real, dimension(:), allocatable, save dsn_sv
Definition: VARxSV.F90:27
surf_sisvat_mod::sisvat_ini
subroutine sisvat_ini(knon)
Definition: surf_sisvat_mod.F90:783
var0sv::dtz_sv
real, dimension(-nsol:0) dtz_sv
Definition: VAR0SV.F90:17
varxsv::lsmask
integer, dimension(:), allocatable, save lsmask
Definition: VARxSV.F90:9
varxsv::zsblsv
real, save zsblsv
Definition: VARxSV.F90:84
yomcst_sisvat
Definition: YOMCST_SISVAT.F90:1
var0sv::s1__sv
real, dimension(0:nsot) s1__sv
Definition: VAR0SV.F90:26
varysv::init_varysv
subroutine init_varysv
Definition: VARySV.F90:94
varxsv::uqs_sv
real, dimension(:), allocatable, save uqs_sv
Definition: VARxSV.F90:118
varxsv::us__sv
real, dimension(:), allocatable, save us__sv
Definition: VARxSV.F90:112
varxsv::uss_sv
real, dimension(:), allocatable, save uss_sv
Definition: VARxSV.F90:120
varxsv::slopsv
real, dimension(:), allocatable, save slopsv
Definition: VARxSV.F90:80
mod_grid_phy_lmdz
Definition: mod_grid_phy_lmdz.F90:4
vartsv::ntaver
integer, parameter ntaver
Definition: VARtSV.F90:7
iostart::put_field
Definition: iostart.F90:18
var0sv::isnpsv
integer, dimension(nsno) isnpsv
Definition: VAR0SV.F90:9
vartsv::init_vartsv
subroutine init_vartsv
Definition: VARtSV.F90:48
vardsv::bchdsv
real, dimension(0:nsot), parameter bchdsv
Definition: VARdSV.F90:111
mod_phys_lmdz_mpi_data::is_mpi_root
logical, save is_mpi_root
Definition: mod_phys_lmdz_mpi_data.F90:39
vardsv::dz_dsv
real, dimension(-nsol:0) dz_dsv
Definition: VARdSV.F90:92
varxsv::eta_sv
real, dimension(:,:), allocatable, save eta_sv
Definition: VARxSV.F90:164
vardsv::so3dsv
real, parameter so3dsv
Definition: VARdSV.F90:26
yomcst_sisvat::rsigma
real, parameter rsigma
Definition: YOMCST_SISVAT.F90:30
vardsv::nsot
integer, parameter nsot
Definition: VARdSV.F90:97
vardcp
Definition: VARdCP.F90:1
vardsv::zz_dsv
real zz_dsv
Definition: VARdSV.F90:95
varxsv::qat_sv
real, dimension(:), allocatable, save qat_sv
Definition: VARxSV.F90:66
varysv::zzsnsv
real, dimension(:,:), allocatable, save zzsnsv
Definition: VARySV.F90:81
varxsv::dt__sv
real, save dt__sv
Definition: VARxSV.F90:86
vartsv::tsfnsv
real, dimension(:), allocatable, save tsfnsv
Definition: VARtSV.F90:18
varxsv::rrcasv
real, dimension(:), allocatable, save rrcasv
Definition: VARxSV.F90:153
vartsv::v__mem
real, dimension(:,:), allocatable, save v__mem
Definition: VARtSV.F90:36
varysv::alb1sv
real, dimension(:), allocatable, save alb1sv
Definition: VARySV.F90:18
var0sv::nkhy
integer, parameter nkhy
Definition: VAR0SV.F90:12
varysv::alb3sv
real, dimension(:), allocatable, save alb3sv
Definition: VARySV.F90:22
surf_sisvat_mod::surf_sisvat
subroutine surf_sisvat(knon, rlon, rlat, ikl2i, itime, dtime, debut, lafin, rmu0, swdown, lwdown, pexner, ps, p1lay, precip_rain, precip_snow, precip_snow_adv, snow_adv, bl_height, wind_velo, temp_air, dens_air, spechum, tsurf, rugos, snow_cont_air, alb_soil, slope, cloudf, radsol, qsol, tsoil, snow, snowhgt, qsnow, to_ice, sissnow, agesno, AcoefH, AcoefQ, BcoefH, BcoefQ, cdragh, runoff_lic, evap, fluxsens, fluxlat, dflux_s, dflux_l, tsurf_new, alb1, alb2, alb3, emis_new, z0_new, qsurf)
Definition: surf_sisvat_mod.F90:22
var0sv::implic
real implic
Definition: VAR0SV.F90:13
dimphy
Definition: dimphy.F90:1
varxsv::coszsv
real, dimension(:), allocatable, save coszsv
Definition: VARxSV.F90:18
surf_sisvat_mod::sisvatetat0
subroutine sisvatetat0(fichnom, ikl2i)
Definition: surf_sisvat_mod.F90:1180
epsi
!$Id!Thermodynamical constants for t0 real clmci real epsi
Definition: cvthermo.h:6
unit
!$Header!integer nvarmx s s unit
Definition: gradsdef.h:20
vardsv::so1dsv
real, parameter so1dsv
Definition: VARdSV.F90:26
varysv::hls_sv
real, dimension(:), allocatable, save hls_sv
Definition: VARySV.F90:55
varxsv::dldtsv
real, dimension(:), allocatable, save dldtsv
Definition: VARxSV.F90:62
varxsv::swf_sv
real, dimension(:), allocatable, save swf_sv
Definition: VARxSV.F90:178
varxsv::g2snsv
real, dimension(:,:), allocatable, save g2snsv
Definition: VARxSV.F90:168
indice_sol_mod
Definition: indice_sol_mod.F90:1
iostart::open_restartphy
subroutine, public open_restartphy(filename)
Definition: iostart.F90:312
varxsv::esnbsv
real, dimension(:), allocatable, save esnbsv
Definition: VARxSV.F90:31
varxsv::vvs_sv
real, dimension(:), allocatable, save vvs_sv
Definition: VARxSV.F90:50
varxsv::bg1ssv
real, dimension(:), allocatable, save bg1ssv
Definition: VARxSV.F90:37
mod_phys_lmdz_mpi_data
Definition: mod_phys_lmdz_mpi_data.F90:4