source: LMDZ4/trunk/libf/phylmd/ocean_slab_mod.F90 @ 1645

Last change on this file since 1645 was 1067, checked in by Laurent Fairhead, 16 years ago
  • Modifications lie au premier niveau du modele pour la diffusion turbulent

du vent.

  • Preparation pour un couplage des courrant oceaniques.

JG
  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 6.4 KB
Line 
1!
2MODULE ocean_slab_mod
3!
4! This module is used for both surface ocean and sea-ice when using the slab ocean,
5! "ocean=slab".
6!
7  IMPLICIT NONE
8  PRIVATE
9  PUBLIC :: ocean_slab_frac, ocean_slab_noice
10
11CONTAINS
12!
13!****************************************************************************************
14!
15  SUBROUTINE ocean_slab_frac(itime, dtime, jour, pctsrf, is_modified)
16
17    USE dimphy
18    USE limit_read_mod
19    USE surface_data
20    INCLUDE "indicesol.h"
21!    INCLUDE "clesphys.h"
22
23! Arguments
24!****************************************************************************************
25    INTEGER, INTENT(IN)                        :: itime   ! numero du pas de temps courant
26    INTEGER, INTENT(IN)                        :: jour    ! jour a lire dans l'annee
27    REAL   , INTENT(IN)                        :: dtime   ! pas de temps de la physique (en s)
28    REAL, DIMENSION(klon,nbsrf), INTENT(INOUT) :: pctsrf  ! sub-surface fraction
29    LOGICAL, INTENT(OUT)                       :: is_modified ! true if pctsrf is modified at this time step
30
31! Local variables
32!****************************************************************************************
33    CHARACTER (len = 80)   :: abort_message
34    CHARACTER (len = 20)   :: modname = 'ocean_slab_frac'
35
36
37    IF (version_ocean == 'sicOBS') THEN   
38       CALL limit_read_frac(itime, dtime, jour, pctsrf, is_modified)
39    ELSE
40       abort_message='Ocean slab model without forced sea-ice fractions has to be rewritten!!!'
41       CALL abort_gcm(modname,abort_message,1)
42! Here should sea-ice/ocean fraction either be calculated or returned if saved as a module varaiable
43! (in the case the new fractions are calculated in ocean_slab_ice or ocean_slab_noice subroutines). 
44    END IF
45
46  END SUBROUTINE ocean_slab_frac
47!
48!****************************************************************************************
49!
50  SUBROUTINE ocean_slab_noice( &
51       itime, dtime, jour, knon, knindex, &
52       p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum, &
53       AcoefH, AcoefQ, BcoefH, BcoefQ, &
54       AcoefU, AcoefV, BcoefU, BcoefV, &
55       ps, u1, v1, tsurf_in, &
56       radsol, snow, agesno, &
57       qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
58       tsurf_new, dflux_s, dflux_l, lmt_bils)
59   
60    USE dimphy
61    USE calcul_fluxs_mod
62 
63    INCLUDE "indicesol.h"
64    INCLUDE "iniprint.h"
65
66! Input arguments
67!****************************************************************************************
68    INTEGER, INTENT(IN)                  :: itime
69    INTEGER, INTENT(IN)                  :: jour
70    INTEGER, INTENT(IN)                  :: knon
71    INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
72    REAL, INTENT(IN)                     :: dtime
73    REAL, DIMENSION(klon), INTENT(IN)    :: p1lay
74    REAL, DIMENSION(klon), INTENT(IN)    :: cdragh, cdragm
75    REAL, DIMENSION(klon), INTENT(IN)    :: precip_rain, precip_snow
76    REAL, DIMENSION(klon), INTENT(IN)    :: temp_air, spechum
77    REAL, DIMENSION(klon), INTENT(IN)    :: AcoefH, AcoefQ, BcoefH, BcoefQ
78    REAL, DIMENSION(klon), INTENT(IN)    :: AcoefU, AcoefV, BcoefU, BcoefV
79    REAL, DIMENSION(klon), INTENT(IN)    :: ps
80    REAL, DIMENSION(klon), INTENT(IN)    :: u1, v1
81    REAL, DIMENSION(klon), INTENT(IN)    :: tsurf_in
82
83! In/Output arguments
84!****************************************************************************************
85    REAL, DIMENSION(klon), INTENT(INOUT) :: radsol
86    REAL, DIMENSION(klon), INTENT(INOUT) :: snow
87    REAL, DIMENSION(klon), INTENT(INOUT) :: agesno
88   
89! Output arguments
90!****************************************************************************************
91    REAL, DIMENSION(klon), INTENT(OUT)   :: qsurf
92    REAL, DIMENSION(klon), INTENT(OUT)   :: evap, fluxsens, fluxlat
93    REAL, DIMENSION(klon), INTENT(OUT)   :: flux_u1, flux_v1
94    REAL, DIMENSION(klon), INTENT(OUT)   :: tsurf_new
95    REAL, DIMENSION(klon), INTENT(OUT)   :: dflux_s, dflux_l     
96    REAL, DIMENSION(klon), INTENT(OUT)   :: lmt_bils
97
98! Local variables
99!****************************************************************************************
100    INTEGER               :: i
101    REAL, DIMENSION(klon) :: cal, beta, dif_grnd
102    REAL, DIMENSION(klon) :: lmt_bils_oce, lmt_foce, diff_sst
103    REAL, DIMENSION(klon) :: u0, v0
104    REAL, DIMENSION(klon) :: u1_lay, v1_lay
105    REAL                  :: calc_bils_oce, deltat
106    REAL, PARAMETER       :: cyang=50.0 * 4.228e+06 ! capacite calorifique volumetrique de l'eau J/(m2 K)
107
108!****************************************************************************************
109! 1) Flux calculation
110!
111!****************************************************************************************
112    cal(:)      = 0.
113    beta(:)     = 1.
114    dif_grnd(:) = 0.
115    agesno(:)   = 0.
116   
117! Suppose zero surface speed
118    u0(:)=0.0
119    v0(:)=0.0
120    u1_lay(:) = u1(:) - u0(:)
121    v1_lay(:) = v1(:) - v0(:)
122
123    CALL calcul_fluxs(knon, is_oce, dtime, &
124         tsurf_in, p1lay, cal, beta, cdragh, ps, &
125         precip_rain, precip_snow, snow, qsurf,  &
126         radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, &
127         AcoefH, AcoefQ, BcoefH, BcoefQ, &
128         tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
129
130! - Flux calculation at first modele level for U and V
131    CALL calcul_flux_wind(knon, dtime, &
132         u0, v0, u1, v1, cdragm, &
133         AcoefU, AcoefV, BcoefU, BcoefV, &
134         p1lay, temp_air, &
135         flux_u1, flux_v1) 
136
137!****************************************************************************************
138! 2) Get global variables lmt_bils and lmt_foce from file limit_slab.nc
139!
140!****************************************************************************************
141    CALL limit_slab(itime, dtime, jour, lmt_bils, lmt_foce, diff_sst)  ! global pour un processus
142
143    lmt_bils_oce(:) = 0.
144    WHERE (lmt_foce > 0.)
145       lmt_bils_oce = lmt_bils / lmt_foce ! global
146    END WHERE
147
148!****************************************************************************************
149! 3) Recalculate new temperature
150!
151!****************************************************************************************
152    DO i = 1, knon
153       calc_bils_oce = radsol(i) + fluxsens(i) + fluxlat(i)
154       deltat        = (calc_bils_oce - lmt_bils_oce(knindex(i)))*dtime/cyang +diff_sst(knindex(i))
155       tsurf_new(i)  = tsurf_in(i) + deltat
156    END DO
157
158  END SUBROUTINE ocean_slab_noice
159!
160!****************************************************************************************
161!
162END MODULE ocean_slab_mod
Note: See TracBrowser for help on using the repository browser.