source: LMDZ6/branches/Amaury_dev/libf/phylmd/calbeta_clim.F90 @ 5103

Last change on this file since 5103 was 5103, checked in by abarral, 8 weeks ago

Handle CPP_INLANDSIS in lmdz_cppkeys_wrapper.F90
Remove obsolete key wrgrads_thermcell, _ADV_HALO, _ADV_HALLO, isminmax
Remove redundant uses of CPPKEY_INCA (thanks acozic)
Remove obsolete misc/write_field.F90
Remove unused ioipsl_* wrappers
Remove calls to WriteField_u with wrong signature
Convert .F -> .[fF]90
(lint) uppercase fortran operators
[note: 1d and iso still broken - working on it]
  • 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
5SUBROUTINE calbeta_clim(klon,time,lat_radian,beta)
6
7  !======================================================================
8  ! Auteur(s): A.K. TRAORE
9  !======================================================================
10
11  !USE phys_local_var_mod, ONLY: ideal_beta !pour faire la variable dans le
12  ! physiq.f pour des sorties directes de beta
13
14  USE phys_cal_mod, only: year_len
15  USE print_control_mod, ONLY: prt_level
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
28  pi=2.*asin(1.)
29  beta=0.
30
31  !calcul des cordonnees
32
33  ! PRINT*,'LATITUDES BETA ',lat_radian
34  time_radian=(time+15.)*2.*pi / year_len
35
36  if (prt_level >= 1) print *, 'time_radian time', time_radian, time
37
38  lat(:)=180.*lat_radian(:)/pi !lat(:)=lat_radian(:)
39
40  lat_sahel=-5*sin(time_radian)+13
41  lat_nord=lat_sahel+25.
42  lat_sud=lat_sahel-25.
43  do j=1,klon
44     !===========
45     if (lat(j) < 5. ) then
46
47        logbeta(j)=0.2*(lat(j)-lat_sud)-1.6
48        beta(j)=10**(logbeta(j))
49        beta(j)=max(beta(j),0.03)
50        beta(j)=min(beta(j),0.22)
51        ! PRINT*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
52        !===========
53     elseif (lat(j) < 22.) then !lat(j)<22.
54
55        logbeta(j)=-0.25*(lat(j)-lat_sahel)-1.6
56        beta(j)=10**(logbeta(j))
57        beta(j)=max(beta(j),1.e-2)
58        beta(j)=min(beta(j),0.22)
59        ! PRINT*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
60        !===========
61     else
62        logbeta(j)=0.25*(lat(j)-lat_nord)-1.
63        beta(j)=10**(logbeta(j))
64        beta(j)=max(beta(j),1.e-2)
65        beta(j)=min(beta(j),0.25)
66        ! PRINT*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
67     endif
68     !===========
69  enddo
70
71END SUBROUTINE  calbeta_clim
Note: See TracBrowser for help on using the repository browser.