source: LMDZ5/trunk/libf/phylmd/calbeta_clim.F90 @ 2098

Last change on this file since 2098 was 2098, checked in by lguez, 10 years ago

Replaced 360 in calbeta_clim by length of current year according to
chosen calendar. Length of current year is given by
ioget_year_len. Since we already need this for ozone, moved the call
to ioget_year_len from physiq to phys_cal_mod and created variable
year_len of module phys_cal_mod.

Control the output from minmaxqfi.

Non-ASCII characters in comments are not always rendered properly and
they risk being lost. See revision 1740.
  • 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
File size: 1.9 KB
Line 
1!
2! $Header: /home/cvsroot/LMDZ4/libf/phylmd/calbeta.F90,v 1.2 2007/06/22 12:49:51
3! fairhead Exp $
4!
5
6SUBROUTINE calbeta_clim(klon,time,lat_radian,beta)
7
8  !======================================================================
9  ! Auteur(s): A.K. TRAORE
10  !======================================================================
11
12  !USE phys_local_var_mod, ONLY : ideal_beta !pour faire la variable dans le
13  ! physiq.f pour des sorties directes de beta
14
15  USE phys_cal_mod, only: year_len
16
17  implicit none
18  integer klon,nt,j,it
19  real logbeta(klon),pi
20  real lat(klon),lat_radian(klon)
21  integer time
22  real time_radian
23  real lat_sahel,beta(klon)
24  real lat_nord,lat_sud
25  !==============================================
26
27  pi=2.*asin(1.)
28  beta=0.
29
30  !calcul des cordonnees
31
32  ! print*,'LATITUDES BETA ',lat_radian
33  time_radian=(time+15.)*2.*pi / year_len
34
35  print*,'BETA time_radian',time_radian,time
36
37  lat(:)=180.*lat_radian(:)/pi !lat(:)=lat_radian(:)
38
39  lat_sahel=-5*sin(time_radian)+13
40  lat_nord=lat_sahel+25.
41  lat_sud=lat_sahel-25.
42  do j=1,klon
43     !===========
44     if (lat(j) < 5. ) then
45
46        logbeta(j)=0.2*(lat(j)-lat_sud)-1.6
47        beta(j)=10**(logbeta(j))
48        beta(j)=max(beta(j),0.03)
49        beta(j)=min(beta(j),0.22)
50        ! print*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
51        !===========
52     elseif (lat(j) < 22.) then !lat(j)<22.
53
54        logbeta(j)=-0.25*(lat(j)-lat_sahel)-1.6
55        beta(j)=10**(logbeta(j))
56        beta(j)=max(beta(j),1.e-2)
57        beta(j)=min(beta(j),0.22)
58        ! print*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
59        !===========
60     else
61        logbeta(j)=0.25*(lat(j)-lat_nord)-1.
62        beta(j)=10**(logbeta(j))
63        beta(j)=max(beta(j),1.e-2)
64        beta(j)=min(beta(j),0.25)
65        ! print*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
66     endif
67     !===========
68  enddo
69
70end SUBROUTINE calbeta_clim
Note: See TracBrowser for help on using the repository browser.