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

Last change on this file since 2031 was 1907, checked in by lguez, 11 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
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

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.
  • 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: 2.0 KB
RevLine 
[1724]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!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
15implicit none
16integer klon,nt,j,it
17real logbeta(klon),pi
18real lat(klon),lat_radian(klon)
19integer time
20real time_radian
21real lat_sahel,beta(klon)
22real lat_nord,lat_sud
23!==============================================
24
25pi=2.*asin(1.)
26beta=0.
27   
28            !calcul des cordonnees
29
30!  print*,'LATITUDES BETA ',lat_radian
31   time_radian=(time+15.)*2.*pi/360.
32
33   print*,'BETA time_radian',time_radian,time
34
35   lat(:)=180.*lat_radian(:)/pi !lat(:)=lat_radian(:)
36       
37   lat_sahel=-5*sin(time_radian)+13
38   lat_nord=lat_sahel+25.
39   lat_sud=lat_sahel-25.
40do j=1,klon
41    !===========
42        if (lat(j) < 5. ) then
43               
44                logbeta(j)=0.2*(lat(j)-lat_sud)-1.6
45                beta(j)=10**(logbeta(j))
46                beta(j)=max(beta(j),0.03)
47                beta(j)=min(beta(j),0.22)
48!              print*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
49    !===========
50        elseif (lat(j) < 22.) then !lat(j)<22.
51               
52                logbeta(j)=-0.25*(lat(j)-lat_sahel)-1.6
53                beta(j)=10**(logbeta(j))
54                beta(j)=max(beta(j),1.e-2)
55                beta(j)=min(beta(j),0.22)
56!              print*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
57    !===========
58        else
59              logbeta(j)=0.25*(lat(j)-lat_nord)-1.
60              beta(j)=10**(logbeta(j))
61              beta(j)=max(beta(j),1.e-2)
62              beta(j)=min(beta(j),0.25)
63!              print*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
64        endif
65    !===========
66enddo
67end
Note: See TracBrowser for help on using the repository browser.