source: LMDZ5/branches/LMDZ5-DOFOCO/libf/phylmd/calbeta_clim.F90 @ 5005

Last change on this file since 5005 was 1724, checked in by musat, 12 years ago

Ajout coefficient beta (=Evap/Evap? potentielle) fonction de la latitude
et de la saison (voir these AbdulKhadre? Traore).

Frederic Hourdin
IM

File size: 2.0 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!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.