/[lmdze]/trunk/phylmd/Interface_surf/calbeta.f
ViewVC logotype

Contents of /trunk/phylmd/Interface_surf/calbeta.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 99 - (show annotations)
Wed Jul 2 18:39:15 2014 UTC (9 years, 11 months ago) by guez
File size: 1806 byte(s) 
Created procedure test_disvert (following LMDZ). Added procedures
hybrid and funcd in module disvert_m. Upgraded compute_ab from
internal procedure of disvert to module procedure. Added variables y,
ya in module disvert_m. Upgraded s from local variable of procedure
disvert to module variable.

Renamed allowed value of variable vert_sampling in procedure disvert
from "read" to "read_hybrid". Added possibility to read pressure
values, value "read_pressure". Replaced vertical distribution for
value "param" by the distribution "strato_correct" from LMDZ (but kept
the value "param"). In case "tropo", replaced 1 by dsigmin (following
LMDZ). In case "strato", replaced 0.3 by dsigmin (following LMDZ).

Changed computation of bp in procedure compute_ab.

Removed debugindex case in clmain. Removed useless argument rlon of
procedure clmain. Removed useless variables ytaux, ytauy of procedure
clmain.

Removed intermediary variables tsol, qsol, tsolsrf, tslab in procedure
etat0.

Removed variable ok_veget:. coupling with the model Orchid is not
possible. Removed variable ocean: modeling an ocean slab is not
possible.

Removed useless variables tmp_rriv and tmp_rcoa from module
interface_surf.

Moved initialization of variables da, mp, phi in procedure physiq to
to inside the test iflag_con >= 3.
1 SUBROUTINE calbeta(dtime, indice, knon, snow, qsol, vbeta, vcal, vdif)
2 USE dimens_m
3 USE indicesol
4 USE dimphy
5 USE conf_gcm_m
6 USE suphec_m
7 IMPLICIT NONE
8 ! ======================================================================
9 ! Auteur(s): Z.X. Li (LMD/CNRS) (adaptation du GCM du LMD)
10 ! date: 19940414
11 ! ======================================================================
12
13 ! Calculer quelques parametres pour appliquer la couche limite
14 ! ------------------------------------------------------------
15 REAL tau_gl ! temps de relaxation pour la glace de mer
16 ! cc PARAMETER (tau_gl=86400.0*30.0)
17 PARAMETER (tau_gl=86400.0*5.0)
18 REAL mx_eau_sol
19 PARAMETER (mx_eau_sol=150.0)
20
21 REAL calsol, calsno, calice ! epaisseur du sol: 0.15 m
22 PARAMETER (calsol=1.0/(2.5578E+06*0.15))
23 PARAMETER (calsno=1.0/(2.3867E+06*0.15))
24 PARAMETER (calice=1.0/(5.1444E+06*0.15))
25
26 INTEGER i
27
28 REAL dtime
29 REAL snow(klon), qsol(klon)
30 INTEGER indice, knon
31
32 REAL vbeta(klon)
33 REAL vcal(klon)
34 REAL vdif(klon)
35
36
37 IF (indice==is_oce) THEN
38 DO i = 1, knon
39 vcal(i) = 0.0
40 vbeta(i) = 1.0
41 vdif(i) = 0.0
42 END DO
43 END IF
44
45 IF (indice==is_sic) THEN
46 DO i = 1, knon
47 vcal(i) = calice
48 IF (snow(i)>0.0) vcal(i) = calsno
49 vbeta(i) = 1.0
50 vdif(i) = 1.0/tau_gl
51 ! cc vdif(i) = calice/tau_gl ! c'etait une erreur
52 END DO
53 END IF
54
55 IF (indice==is_ter) THEN
56 DO i = 1, knon
57 vcal(i) = calsol
58 IF (snow(i)>0.0) vcal(i) = calsno
59 vbeta(i) = min(2.0*qsol(i)/mx_eau_sol, 1.0)
60 vdif(i) = 0.0
61 END DO
62 END IF
63
64 IF (indice==is_lic) THEN
65 DO i = 1, knon
66 vcal(i) = calice
67 IF (snow(i)>0.0) vcal(i) = calsno
68 vbeta(i) = 1.0
69 vdif(i) = 0.0
70 END DO
71 END IF
72
73 RETURN
74 END SUBROUTINE calbeta
  ViewVC Help
Powered by ViewVC 1.1.21