source: LMDZ6/trunk/libf/phylmd/surf_land_bucket_mod.F90 @ 3981

Last change on this file since 3981 was 3974, checked in by asima, 3 years ago

In the "bucket" soil model, soil thermal inertia (formerly a constant value) can also be now a function of soil moisture "qsol"
(F Cheruy's idea; results from "Stage L3" Eve Rebouillat, with E Vignon, A Sima, F Cheruy).

The function depends on iflag_inertie, read from physiq.def (via conf_phys_m.F90), as follows :
(only in physiq.def_NPv6.1 for the time being)


# Controle de la definition de l'inertie thermique du sol dans le modele "bucket":
# 0 (default) : constant
# 1 : fonction lineaire de qsol
# 2 : 2 fonctions lineaires de qsol :
# une pour "sable" (boite lat/lon pour Sahara), une pour le reste (limons et argile)
# 3 : fonction lineaire de qsol, valeurs entre 900 pour qsol=0 et 2000 pour qsol=150

iflag_inertie=0


Relevant changes are in soil.F90 and physiq.def_NPv6.1 ; the other files are adjusted accordingly.

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 6.2 KB
Line 
1!
2MODULE surf_land_bucket_mod
3!
4! Surface land bucket module
5!
6! This module is used when no external land model is choosen.
7!
8  IMPLICIT NONE
9
10CONTAINS
11
12  SUBROUTINE surf_land_bucket(itime, jour, knon, knindex, debut, dtime,&
13       tsurf, p1lay, tq_cdrag, precip_rain, precip_snow, temp_air, &
14       spechum, petAcoef, peqAcoef, petBcoef, peqBcoef, pref, &
15       u1, v1, gustiness, rugoro, swnet, lwnet, &
16       snow, qsol, agesno, tsoil, &
17       qsurf, z0_new, alb1_new, alb2_new, evap, &
18       fluxsens, fluxlat, tsurf_new, dflux_s, dflux_l)
19
20    USE limit_read_mod
21    USE surface_data
22    USE fonte_neige_mod
23    USE calcul_fluxs_mod
24    USE cpl_mod
25    USE dimphy
26    USE geometry_mod, ONLY: longitude,latitude
27    USE mod_grid_phy_lmdz
28    USE mod_phys_lmdz_para
29    USE indice_sol_mod
30!****************************************************************************************
31! Bucket calculations for surface.
32!
33    INCLUDE "clesphys.h"
34    INCLUDE "dimsoil.h"
35    INCLUDE "YOMCST.h"
36
37! Input variables 
38!****************************************************************************************
39    INTEGER, INTENT(IN)                     :: itime, jour, knon
40    INTEGER, DIMENSION(klon), INTENT(IN)    :: knindex
41    LOGICAL, INTENT(IN)                     :: debut
42    REAL, INTENT(IN)                        :: dtime
43    REAL, DIMENSION(klon), INTENT(IN)       :: tsurf
44    REAL, DIMENSION(klon), INTENT(IN)       :: p1lay
45    REAL, DIMENSION(klon), INTENT(IN)       :: tq_cdrag
46    REAL, DIMENSION(klon), INTENT(IN)       :: precip_rain, precip_snow
47    REAL, DIMENSION(klon), INTENT(IN)       :: temp_air, spechum
48    REAL, DIMENSION(klon), INTENT(IN)       :: petAcoef, peqAcoef
49    REAL, DIMENSION(klon), INTENT(IN)       :: petBcoef, peqBcoef
50    REAL, DIMENSION(klon), INTENT(IN)       :: pref
51    REAL, DIMENSION(klon), INTENT(IN)       :: u1, v1, gustiness
52    REAL, DIMENSION(klon), INTENT(IN)       :: rugoro
53    REAL, DIMENSION(klon), INTENT(IN)       :: swnet, lwnet
54
55! In/Output variables
56!****************************************************************************************
57    REAL, DIMENSION(klon), INTENT(INOUT)          :: snow, qsol
58    REAL, DIMENSION(klon), INTENT(INOUT)          :: agesno
59    REAL, DIMENSION(klon, nsoilmx), INTENT(INOUT) :: tsoil
60
61! Output variables
62!****************************************************************************************
63    REAL, DIMENSION(klon), INTENT(OUT)       :: qsurf
64    REAL, DIMENSION(klon), INTENT(OUT)       :: z0_new
65    REAL, DIMENSION(klon), INTENT(OUT)       :: alb1_new, alb2_new
66    REAL, DIMENSION(klon), INTENT(OUT)       :: evap, fluxsens, fluxlat
67    REAL, DIMENSION(klon), INTENT(OUT)       :: tsurf_new
68    REAL, DIMENSION(klon), INTENT(OUT)       :: dflux_s, dflux_l     
69
70! Local variables
71!****************************************************************************************
72    REAL, DIMENSION(klon) :: soilcap, soilflux
73    REAL, DIMENSION(klon) :: cal, beta, dif_grnd, capsol
74    REAL, DIMENSION(klon) :: alb_neig, alb_lim
75    REAL, DIMENSION(klon) :: zfra
76    REAL, DIMENSION(klon) :: radsol       ! total net radiance at surface
77    REAL, DIMENSION(klon) :: u0, v0, u1_lay, v1_lay
78    REAL, DIMENSION(klon) :: dummy_riverflow,dummy_coastalflow
79    INTEGER               :: i
80!
81!****************************************************************************************
82
83
84!
85!* Read from limit.nc : albedo(alb_lim), length of rugosity(z0_new)
86!
87    CALL limit_read_rug_alb(itime, dtime, jour,&
88         knon, knindex, &
89         z0_new, alb_lim)
90!
91!* Calcultaion of fluxes
92!
93
94! calculate total absorbed radiance at surface
95       radsol(:) = 0.0
96       radsol(1:knon) = swnet(1:knon) + lwnet(1:knon)
97
98! calculate constants
99    CALL calbeta(dtime, is_ter, knon, snow, qsol, beta, capsol, dif_grnd)
100    if (type_veget=='betaclim') then
101       CALL calbeta_clim(knon,jour,latitude(knindex(1:knon)),beta)
102    endif
103       
104! calculate temperature, heat capacity and conduction flux in soil
105    IF (soil_model) THEN
106       CALL soil(dtime, is_ter, knon, snow, tsurf, qsol,  &
107        & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil, soilcap, soilflux)
108
109       DO i=1, knon
110          cal(i) = RCPD / soilcap(i)
111          radsol(i) = radsol(i)  + soilflux(i)
112       END DO
113    ELSE
114       cal(:) = RCPD * capsol(:)
115       IF (klon_glo .EQ. 1) THEN
116         cal(:) = 0.
117       ENDIF
118    ENDIF
119   
120! Suppose zero surface speed
121    u0(:)=0.0
122    v0(:)=0.0
123    u1_lay(:) = u1(:) - u0(:)
124    v1_lay(:) = v1(:) - v0(:)
125
126    CALL calcul_fluxs(knon, is_ter, dtime, &
127         tsurf, p1lay, cal, beta, tq_cdrag, tq_cdrag, pref, &
128         precip_rain, precip_snow, snow, qsurf,  &
129         radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, gustiness, &
130         1.,petAcoef, peqAcoef, petBcoef, peqBcoef, &
131         tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
132   
133!
134!* Calculate snow height, run_off, age of snow
135!     
136    CALL fonte_neige( knon, is_ter, knindex, dtime, &
137         tsurf, precip_rain, precip_snow, &
138         snow, qsol, tsurf_new, evap)
139!
140!* Calculate the age of snow
141!
142    CALL albsno(klon,knon,dtime,agesno(:),alb_neig(:), precip_snow(:)) 
143   
144    WHERE (snow(1 : knon) .LT. 0.0001) agesno(1 : knon) = 0.
145   
146    DO i=1, knon
147       zfra(i) = MAX(0.0,MIN(1.0, snow(i)/(snow(i)+10.0)))
148       alb_lim(i)  = alb_neig(i) *zfra(i) + alb_lim(i)*(1.0-zfra(i))
149    END DO
150
151!
152!* Return albedo :
153!    alb1_new and alb2_new are here given the same values
154!
155    alb1_new(:) = 0.0
156    alb2_new(:) = 0.0
157    alb1_new(1:knon) = alb_lim(1:knon)
158    alb2_new(1:knon) = alb_lim(1:knon)
159       
160!
161!* Calculate the rugosity
162!
163    DO i = 1, knon
164       z0_new(i) = MAX(1.5e-05,SQRT(z0_new(i)**2 + rugoro(i)**2))
165    END DO
166
167!* Send to coupler
168!  The run-off from river and coast are not calculated in the bucket modele.
169!  For testing purpose of the coupled modele we put the run-off to zero.
170    IF (type_ocean=='couple') THEN
171       dummy_riverflow(:)   = 0.0
172       dummy_coastalflow(:) = 0.0
173       CALL cpl_send_land_fields(itime, knon, knindex, &
174            dummy_riverflow, dummy_coastalflow)
175    ENDIF
176
177!
178!* End
179!
180  END SUBROUTINE surf_land_bucket
181!
182!****************************************************************************************
183!
184END MODULE surf_land_bucket_mod
Note: See TracBrowser for help on using the repository browser.