Changeset 2160 for LMDZ5/branches/testing/libf/phylmd/calbeta_clim.F90

Timestamp:

Nov 28, 2014, 4:36:29 PM (10 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes -r2070:2158 into testing branch. Compilation problems introduced by revision r2155 have been corrected by hand

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/calbeta_clim.F90 (modified) (1 diff)

LMDZ5/branches/testing/libf/phylmd/calbeta_clim.F90

@@ -5 +5 @@
 
 SUBROUTINE calbeta_clim(klon,time,lat_radian,beta)
- 
-!======================================================================
-! Auteur(s): A.K. TRAORE
-!======================================================================
-!USE phys_local_var_mod, ONLY : ideal_beta !pour faire la variable dans le
-! physiq.f pour des sorties directes de beta
-!==============================================
 
-implicit none
-integer klon,nt,j,it
-real logbeta(klon),pi
-real lat(klon),lat_radian(klon)
-integer time
-real time_radian
-real lat_sahel,beta(klon)
-real lat_nord,lat_sud
-!==============================================
+  !======================================================================
+  ! Auteur(s): A.K. TRAORE
+  !======================================================================
 
-pi=2.*asin(1.)
-beta=0.
-    
-            !calcul des cordonnees
+  !USE phys_local_var_mod, ONLY : ideal_beta !pour faire la variable dans le
+  ! physiq.f pour des sorties directes de beta
 
-!  print*,'LATITUDES BETA ',lat_radian
-   time_radian=(time+15.)*2.*pi/360.
+  USE phys_cal_mod, only: year_len
 
-   print*,'BETA time_radian',time_radian,time
+  implicit none
+  integer klon,nt,j,it
+  real logbeta(klon),pi
+  real lat(klon),lat_radian(klon)
+  integer time
+  real time_radian
+  real lat_sahel,beta(klon)
+  real lat_nord,lat_sud
 
-   lat(:)=180.*lat_radian(:)/pi !lat(:)=lat_radian(:)
-        
-   lat_sahel=-5*sin(time_radian)+13
-   lat_nord=lat_sahel+25.
-   lat_sud=lat_sahel-25.
-do j=1,klon
-    !===========
-        if (lat(j) < 5. ) then
-                
-                logbeta(j)=0.2*(lat(j)-lat_sud)-1.6
-                beta(j)=10**(logbeta(j))
-                beta(j)=max(beta(j),0.03)
-                beta(j)=min(beta(j),0.22)
-!              print*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
-    !===========
-        elseif (lat(j) < 22.) then !lat(j)<22.
-                
-                logbeta(j)=-0.25*(lat(j)-lat_sahel)-1.6
-                beta(j)=10**(logbeta(j))
-                beta(j)=max(beta(j),1.e-2)
-                beta(j)=min(beta(j),0.22)
-!              print*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
-    !===========
-        else
-              logbeta(j)=0.25*(lat(j)-lat_nord)-1.
-              beta(j)=10**(logbeta(j))
-              beta(j)=max(beta(j),1.e-2)
-              beta(j)=min(beta(j),0.25)
-!              print*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
-        endif
-    !===========
-enddo
-end
+  include "iniprint.h"
+  !==============================================
+
+  pi=2.*asin(1.)
+  beta=0.
+
+  !calcul des cordonnees
+
+  ! print*,'LATITUDES BETA ',lat_radian
+  time_radian=(time+15.)*2.*pi / year_len
+
+  if (prt_level >= 1) print *, 'time_radian time', time_radian, time
+
+  lat(:)=180.*lat_radian(:)/pi !lat(:)=lat_radian(:)
+
+  lat_sahel=-5*sin(time_radian)+13
+  lat_nord=lat_sahel+25.
+  lat_sud=lat_sahel-25.
+  do j=1,klon
+     !===========
+     if (lat(j) < 5. ) then
+
+        logbeta(j)=0.2*(lat(j)-lat_sud)-1.6
+        beta(j)=10**(logbeta(j))
+        beta(j)=max(beta(j),0.03)
+        beta(j)=min(beta(j),0.22)
+        ! print*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
+        !===========
+     elseif (lat(j) < 22.) then !lat(j)<22.
+
+        logbeta(j)=-0.25*(lat(j)-lat_sahel)-1.6
+        beta(j)=10**(logbeta(j))
+        beta(j)=max(beta(j),1.e-2)
+        beta(j)=min(beta(j),0.22)
+        ! print*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
+        !===========
+     else
+        logbeta(j)=0.25*(lat(j)-lat_nord)-1.
+        beta(j)=10**(logbeta(j))
+        beta(j)=max(beta(j),1.e-2)
+        beta(j)=min(beta(j),0.25)
+        ! print*,'j,lat,lat_radian,beta',j,lat(j),lat_radian(j),beta(j)
+     endif
+     !===========
+  enddo
+
+end SUBROUTINE calbeta_clim