source: LMDZ4/trunk/libf/phylmd/ocean_forced_mod.F90 @ 1006

Last change on this file since 1006 was 996, checked in by lsce, 16 years ago
  • Modifications liées au calcul des nouveau sous-fractions
  • Nettoyage de ocean slab : il reste uniquement la version avec glace de mer forcé
  • Nouveaux variables pour distiguer la version et type d'ocean : type_ocean=force/slab/couple, version_ocean=opa8/nemo pour couplé ou version_ocean=sicOBS pour slab

JG
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 9.4 KB
RevLine 
[781]1!
2! $Header$
3!
4MODULE ocean_forced_mod
5!
6! This module is used for both the sub-surfaces ocean and sea-ice for the case of a
7! forced ocean,  "ocean=force".
8!
9  USE surface_data,     ONLY : calice, calsno, tau_gl
10  USE fonte_neige_mod,  ONLY : fonte_neige
11  USE calcul_fluxs_mod, ONLY : calcul_fluxs
12  USE dimphy
13
14  IMPLICIT NONE
15
16CONTAINS
17!
18!****************************************************************************************
19!
20  SUBROUTINE ocean_forced_noice(itime, dtime, jour, knon, knindex, &
21       debut, &
22       p1lay, tq_cdrag, precip_rain, precip_snow, &
23       temp_air, spechum, &
24       petAcoef, peqAcoef, petBcoef, peqBcoef, &
25       ps, u1_lay, v1_lay, &
[888]26       radsol, snow, agesno, &
27       qsurf, evap, fluxsens, fluxlat, &
[996]28       tsurf_new, dflux_s, dflux_l)
[781]29!
30! This subroutine treats the "open ocean", all grid points that are not entierly covered
31! by ice.
[996]32! The routine receives data from climatologie file limit.nc and does some calculations at the
[781]33! surface.
34!
[996]35    USE limit_read_mod
[793]36    INCLUDE "indicesol.h"
37    INCLUDE "YOMCST.h"
[781]38
39! Input arguments
40!****************************************************************************************
41    INTEGER, INTENT(IN)                      :: itime, jour, knon
42    INTEGER, DIMENSION(klon), INTENT(IN)     :: knindex
43    LOGICAL, INTENT(IN)                      :: debut
44    REAL, INTENT(IN)                         :: dtime
45    REAL, DIMENSION(klon), INTENT(IN)        :: p1lay
46    REAL, DIMENSION(klon), INTENT(IN)        :: tq_cdrag
47    REAL, DIMENSION(klon), INTENT(IN)        :: precip_rain, precip_snow
48    REAL, DIMENSION(klon), INTENT(IN)        :: temp_air, spechum
49    REAL, DIMENSION(klon), INTENT(IN)        :: petAcoef, peqAcoef
50    REAL, DIMENSION(klon), INTENT(IN)        :: petBcoef, peqBcoef
51    REAL, DIMENSION(klon), INTENT(IN)        :: ps
52    REAL, DIMENSION(klon), INTENT(IN)        :: u1_lay, v1_lay
53
54! In/Output arguments
55!****************************************************************************************
56    REAL, DIMENSION(klon), INTENT(INOUT)     :: radsol
57    REAL, DIMENSION(klon), INTENT(INOUT)     :: snow
58    REAL, DIMENSION(klon), INTENT(INOUT)     :: agesno !? put to 0 in ocean
59 
60! Output arguments
61!****************************************************************************************
62    REAL, DIMENSION(klon), INTENT(OUT)       :: qsurf
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
67! Local variables
68!****************************************************************************************
69    INTEGER                     :: i
70    REAL, DIMENSION(klon)       :: cal, beta, dif_grnd
71    REAL, DIMENSION(klon)       :: alb_neig, tsurf_lim, zx_sl
72    LOGICAL                     :: check=.FALSE.
73
74!****************************************************************************************
75! Start calculation
76!****************************************************************************************
77    IF (check) WRITE(*,*)' Entering ocean_forced_noice'
78
79!****************************************************************************************
80! 1)   
[996]81! Read sea-surface temperature from file limit.nc
[781]82!
83!****************************************************************************************
[996]84    CALL limit_read_sst(knon,knindex,tsurf_lim)
85
[781]86!****************************************************************************************
87! 2)
88! Flux calculation
89!
90!****************************************************************************************
91! Set some variables for calcul_fluxs
92    cal = 0.
93    beta = 1.
94    dif_grnd = 0.
95    alb_neig(:) = 0.
96    agesno(:) = 0.
97 
98! Calcul de tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l and qsurf
99    CALL calcul_fluxs(knon, is_oce, dtime, &
100         tsurf_lim, p1lay, cal, beta, tq_cdrag, ps, &
101         precip_rain, precip_snow, snow, qsurf,  &
102         radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, &
103         petAcoef, peqAcoef, petBcoef, peqBcoef, &
104         tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
105
106
107  END SUBROUTINE ocean_forced_noice
108!
109!****************************************************************************************
110!
111  SUBROUTINE ocean_forced_ice(itime, dtime, jour, knon, knindex, &
112       debut, &
113       tsurf_in, p1lay, tq_cdrag, precip_rain, precip_snow, temp_air, spechum, &
114       petAcoef, peqAcoef, petBcoef, peqBcoef, &
115       ps, u1_lay, v1_lay, &
[888]116       radsol, snow, qsol, agesno, tsoil, &
117       qsurf, alb1_new, alb2_new, evap, fluxsens, fluxlat, &
[996]118       tsurf_new, dflux_s, dflux_l)
[781]119!
120! This subroutine treats the ocean where there is ice.
121! The routine reads data from climatologie file and does flux calculations at the
122! surface.
[996]123!
124    USE limit_read_mod
125
[793]126    INCLUDE "indicesol.h"
[781]127    INCLUDE "dimsoil.h"
[793]128    INCLUDE "YOMCST.h"
129    INCLUDE "clesphys.h"
[781]130
131! Input arguments
132!****************************************************************************************
133    INTEGER, INTENT(IN)                  :: itime, jour, knon
134    INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
135    LOGICAL, INTENT(IN)                  :: debut
136    REAL, INTENT(IN)                     :: dtime
137    REAL, DIMENSION(klon), INTENT(IN)    :: tsurf_in
138    REAL, DIMENSION(klon), INTENT(IN)    :: p1lay
139    REAL, DIMENSION(klon), INTENT(IN)    :: tq_cdrag
140    REAL, DIMENSION(klon), INTENT(IN)    :: precip_rain, precip_snow
141    REAL, DIMENSION(klon), INTENT(IN)    :: temp_air, spechum
142    REAL, DIMENSION(klon), INTENT(IN)    :: petAcoef, peqAcoef
143    REAL, DIMENSION(klon), INTENT(IN)    :: petBcoef, peqBcoef
144    REAL, DIMENSION(klon), INTENT(IN)    :: ps
145    REAL, DIMENSION(klon), INTENT(IN)    :: u1_lay, v1_lay
146
147! In/Output arguments
148!****************************************************************************************
149    REAL, DIMENSION(klon), INTENT(INOUT)          :: radsol
150    REAL, DIMENSION(klon), INTENT(INOUT)          :: snow, qsol
151    REAL, DIMENSION(klon), INTENT(INOUT)          :: agesno
152    REAL, DIMENSION(klon, nsoilmx), INTENT(INOUT) :: tsoil
153
154! Output arguments
155!****************************************************************************************
156    REAL, DIMENSION(klon), INTENT(OUT)            :: qsurf
[888]157    REAL, DIMENSION(klon), INTENT(OUT)            :: alb1_new  ! new albedo in visible SW interval
158    REAL, DIMENSION(klon), INTENT(OUT)            :: alb2_new  ! new albedo in near IR interval
[781]159    REAL, DIMENSION(klon), INTENT(OUT)            :: evap, fluxsens, fluxlat
[888]160    REAL, DIMENSION(klon), INTENT(OUT)            :: tsurf_new
[781]161    REAL, DIMENSION(klon), INTENT(OUT)            :: dflux_s, dflux_l     
162
163! Local variables
164!****************************************************************************************
165    LOGICAL                     :: check=.FALSE.
166    INTEGER                     :: i
167    REAL                        :: zfra
168    REAL, PARAMETER             :: t_grnd=271.35
169    REAL, DIMENSION(klon)       :: cal, beta, dif_grnd, capsol
170    REAL, DIMENSION(klon)       :: alb_neig, tsurf_tmp
171    REAL, DIMENSION(klon)       :: soilcap, soilflux
172
173!****************************************************************************************
174! Start calculation
175!****************************************************************************************
176    IF (check) WRITE(*,*)'Entering surface_seaice, knon=',knon
177
178!****************************************************************************************
[996]179! 1)
[781]180! Flux calculation : tsurf_new, evap, fluxlat, fluxsens,
181!                    dflux_s, dflux_l and qsurf
182!****************************************************************************************
[996]183    tsurf_tmp(:) = tsurf_in(:)
[781]184
185! calculate the parameters cal, beta, capsol and dif_grnd
[996]186    CALL calbeta(dtime, is_sic, knon, snow, qsol, beta, capsol, dif_grnd)
[781]187
188   
189    IF (soil_model) THEN
190! update tsoil and calculate soilcap and soilflux
[996]191       CALL soil(dtime, is_sic, knon, snow, tsurf_tmp, tsoil,soilcap, soilflux)
[781]192       cal(1:knon) = RCPD / soilcap(1:knon)
193       radsol(1:knon) = radsol(1:knon)  + soilflux(1:knon)
194       dif_grnd = 1.0 / tau_gl
195    ELSE
196       dif_grnd = 1.0 / tau_gl
197       cal = RCPD * calice
198       WHERE (snow > 0.0) cal = RCPD * calsno
199    ENDIF
200
201    beta = 1.0
202    CALL calcul_fluxs(knon, is_sic, dtime, &
203         tsurf_tmp, p1lay, cal, beta, tq_cdrag, ps, &
204         precip_rain, precip_snow, snow, qsurf,  &
205         radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, &
206         petAcoef, peqAcoef, petBcoef, peqBcoef, &
207         tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
208
209!****************************************************************************************
[996]210! 2)
[781]211! Calculations due to snow and runoff
212!
213!****************************************************************************************
214    CALL fonte_neige( knon, is_sic, knindex, dtime, &
215         tsurf_tmp, precip_rain, precip_snow, &
216         snow, qsol, tsurf_new, evap)
217   
218! Calculation of albedo at snow (alb_neig) and update the age of snow (agesno)
219!
220    CALL albsno(klon, knon, dtime, agesno(:), alb_neig(:), precip_snow(:)) 
221
222    WHERE (snow(1:knon) .LT. 0.0001) agesno(1:knon) = 0.
223
[888]224    alb1_new(:) = 0.0
[781]225    DO i=1, knon
226       zfra = MAX(0.0,MIN(1.0,snow(i)/(snow(i)+10.0)))
[888]227       alb1_new(i) = alb_neig(i) * zfra +  0.6 * (1.0-zfra)
[781]228    ENDDO
229
[888]230    alb2_new(:) = alb1_new(:)
231
[781]232  END SUBROUTINE ocean_forced_ice
233!
234!****************************************************************************************
235!
236END MODULE ocean_forced_mod
237
238
239
240
241
242
Note: See TracBrowser for help on using the repository browser.