source: LMDZ4/trunk/libf/phylmd/surf_landice_mod.F90 @ 894

Last change on this file since 894 was 888, checked in by Laurent Fairhead, 17 years ago

Modifications sur l'albedo JG
LF
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 7.0 KB
Line 
1!
2! $Header$
3!
4MODULE surf_landice_mod
5 
6  USE dimphy
7  USE surface_data,     ONLY : ocean, calice, calsno
8  USE fonte_neige_mod,  ONLY : fonte_neige, run_off_lic
9  USE cpl_mod,          ONLY : cpl_send_landice_fields
10  USE calcul_fluxs_mod, ONLY : calcul_fluxs
11
12  IMPLICIT NONE
13
14CONTAINS
15!
16!****************************************************************************************
17!
18  SUBROUTINE surf_landice(itime, dtime, knon, knindex, &
19       swnet, lwnet, tsurf, p1lay, &
20       tq_cdrag, precip_rain, precip_snow, temp_air, spechum, &
21       petAcoef, peqAcoef, petBcoef, peqBcoef, &
22       ps, u1_lay, v1_lay, rugoro, pctsrf, &
23       snow, qsurf, qsol, agesno, &
24       tsoil, z0_new, alb1, alb2, evap, fluxsens, fluxlat, &
25       tsurf_new, dflux_s, dflux_l, pctsrf_lic)
26
27    INCLUDE "indicesol.h"
28    INCLUDE "dimsoil.h"
29    INCLUDE "YOMCST.h"
30    INCLUDE "clesphys.h"
31
32! Input variables
33!****************************************************************************************
34    INTEGER, INTENT(IN)                           :: itime, knon
35    INTEGER, DIMENSION(klon), INTENT(in)          :: knindex
36    REAL, INTENT(in)                              :: dtime
37    REAL, DIMENSION(klon), INTENT(IN)             :: swnet ! net shortwave radiance
38    REAL, DIMENSION(klon), INTENT(IN)             :: lwnet ! net longwave radiance
39    REAL, DIMENSION(klon), INTENT(IN)             :: tsurf
40    REAL, DIMENSION(klon), INTENT(IN)             :: p1lay
41    REAL, DIMENSION(klon), INTENT(IN)             :: tq_cdrag
42    REAL, DIMENSION(klon), INTENT(IN)             :: precip_rain, precip_snow
43    REAL, DIMENSION(klon), INTENT(IN)             :: temp_air, spechum
44    REAL, DIMENSION(klon), INTENT(IN)             :: petAcoef, peqAcoef
45    REAL, DIMENSION(klon), INTENT(IN)             :: petBcoef, peqBcoef
46    REAL, DIMENSION(klon), INTENT(IN)             :: ps
47    REAL, DIMENSION(klon), INTENT(IN)             :: u1_lay, v1_lay
48    REAL, DIMENSION(klon), INTENT(IN)             :: rugoro
49    REAL, DIMENSION(klon,nbsrf), INTENT(IN)       :: pctsrf
50
51! In/Output variables
52!****************************************************************************************
53    REAL, DIMENSION(klon), INTENT(INOUT)          :: snow, qsol
54    REAL, DIMENSION(klon), INTENT(INOUT)          :: agesno
55    REAL, DIMENSION(klon, nsoilmx), INTENT(INOUT) :: tsoil
56
57! Output variables
58!****************************************************************************************
59    REAL, DIMENSION(klon), INTENT(OUT)            :: qsurf
60    REAL, DIMENSION(klon), INTENT(OUT)            :: z0_new
61    REAL, DIMENSION(klon), INTENT(OUT)            :: alb1  ! new albedo in visible SW interval
62    REAL, DIMENSION(klon), INTENT(OUT)            :: alb2  ! new albedo in near IR interval
63    REAL, DIMENSION(klon), INTENT(OUT)            :: evap, fluxsens, fluxlat
64    REAL, DIMENSION(klon), INTENT(OUT)            :: tsurf_new
65    REAL, DIMENSION(klon), INTENT(OUT)            :: dflux_s, dflux_l     
66    REAL, DIMENSION(klon), INTENT(OUT)            :: pctsrf_lic
67
68! Local variables
69!****************************************************************************************
70    REAL, DIMENSION(klon)    :: soilcap, soilflux
71    REAL, DIMENSION(klon)    :: cal, beta, dif_grnd
72    REAL, DIMENSION(klon)    :: zfra, alb_neig
73    REAL, DIMENSION(klon)    :: radsol
74
75! End definition
76!****************************************************************************************
77!
78! Initialize output variables
79    alb2(:) = 999999.
80    alb1(:) = 999999.
81
82!****************************************************************************************
83! Calculate total absorbed radiance at surface
84!
85!****************************************************************************************
86    radsol(:) = 0.0
87    radsol(1:knon) = swnet(1:knon) + lwnet(1:knon)
88
89!****************************************************************************************
90! Soil calculations
91!
92!****************************************************************************************
93    IF (soil_model) THEN
94       CALL soil(dtime, is_lic, knon, snow, tsurf, tsoil, soilcap, soilflux)
95       cal(1:knon) = RCPD / soilcap(1:knon)
96       radsol(1:knon)  = radsol(1:knon) + soilflux(1:knon)
97    ELSE
98       cal = RCPD * calice
99       WHERE (snow > 0.0) cal = RCPD * calsno
100    ENDIF
101
102
103!****************************************************************************************
104! Calulate fluxes
105!
106!****************************************************************************************
107    beta(:) = 1.0
108    dif_grnd(:) = 0.0
109
110    CALL calcul_fluxs(knon, is_lic, dtime, &
111         tsurf, p1lay, cal, beta, tq_cdrag, ps, &
112         precip_rain, precip_snow, snow, qsurf,  &
113         radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, &
114         petAcoef, peqAcoef, petBcoef, peqBcoef, &
115         tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
116
117
118!****************************************************************************************
119! Calculate snow height, age, run-off,..
120!   
121!****************************************************************************************
122    CALL fonte_neige( knon, is_lic, knindex, dtime, &
123         tsurf, precip_rain, precip_snow, &
124         snow, qsol, tsurf_new, evap)
125
126
127!****************************************************************************************
128! Calculate albedo
129!
130!****************************************************************************************
131    CALL albsno(klon,knon,dtime,agesno(:),alb_neig(:), precip_snow(:)) 
132    WHERE (snow(1 : knon) .LT. 0.0001) agesno(1 : knon) = 0.
133    zfra(1:knon) = MAX(0.0,MIN(1.0,snow(1:knon)/(snow(1:knon)+10.0)))
134    alb1(1:knon) = alb_neig(1:knon)*zfra(1:knon) + &
135         0.6 * (1.0-zfra(1:knon))
136!
137!IM: plusieurs choix/tests sur l'albedo des "glaciers continentaux"
138!       alb1(1 : knon)  = 0.6 !IM cf FH/GK
139!       alb1(1 : knon)  = 0.82
140!       alb1(1 : knon)  = 0.77 !211003 Ksta0.77
141!       alb1(1 : knon)  = 0.8 !KstaTER0.8 & LMD_ARMIP5
142!IM: KstaTER0.77 & LMD_ARMIP6   
143
144! Attantion: alb1 and alb2 are the same!
145    alb1(1:knon)  = 0.77
146    alb2(1:knon)  = alb1(1:knon)
147
148
149!****************************************************************************************
150! Rugosity
151!
152!****************************************************************************************
153    z0_new(:) = rugoro(:)
154
155
156!****************************************************************************************
157! Return the pourcentage for this sub-surface
158!
159!****************************************************************************************
160    pctsrf_lic(:) = pctsrf(:,is_lic)
161   
162
163!****************************************************************************************
164! Send run-off on land-ice to coupler if coupled ocean.
165! run_off_lic has been calculated in fonte_neige
166!
167!****************************************************************************************
168    IF (ocean=='couple') THEN
169       CALL cpl_send_landice_fields(itime, knon, knindex, run_off_lic)
170    ENDIF
171
172
173  END SUBROUTINE surf_landice
174!
175!****************************************************************************************
176!
177END MODULE surf_landice_mod
178
179
180
Note: See TracBrowser for help on using the repository browser.