/[lmdze]/trunk/phylmd/alboc.f
ViewVC logotype

Diff of /trunk/phylmd/alboc.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/phylmd/albedo.f revision 82 by guez, Wed Mar 5 14:57:53 2014 UTC trunk/phylmd/alboc.f revision 117 by guez, Thu Dec 4 16:48:53 2014 UTC
# Line 1  Line 1 
   
 ! $Header: /home/cvsroot/LMDZ4/libf/phylmd/albedo.F,v 1.2 2005/02/07 15:00:52  
 ! fairhead Exp $  
   
   
   
1  SUBROUTINE alboc(rjour, rlat, albedo)  SUBROUTINE alboc(rjour, rlat, albedo)
2      ! From LMDZ4/libf/phylmd/albedo.F,v 1.2 2005/02/07 15:00:52
3    USE dimens_m    USE dimens_m
4    USE dimphy    USE dimphy
5    USE yomcst    USE yomcst
# Line 126  SUBROUTINE alboc(rjour, rlat, albedo) Line 121  SUBROUTINE alboc(rjour, rlat, albedo)
121    END IF    END IF
122    RETURN    RETURN
123  END SUBROUTINE alboc  END SUBROUTINE alboc
 ! =====================================================================  
 SUBROUTINE alboc_cd(rmu0, albedo)  
   USE dimens_m  
   USE dimphy  
   IMPLICIT NONE  
   ! ======================================================================  
   ! Auteur(s): Z.X. Li (LMD/CNRS)  
   ! date: 19940624  
   ! Calculer l'albedo sur l'ocean en fonction de l'angle zenithal moyen  
   ! Formule due a Larson and Barkstrom (1977) Proc. of the symposium  
   ! on radiation in the atmosphere, 19-28 August 1976, science Press,  
   ! 1977 pp 451-453, ou These de 3eme cycle de Sylvie Joussaume.  
   
   ! Arguments  
   ! rmu0    (in): cosinus de l'angle solaire zenithal  
   ! albedo (out): albedo de surface de l'ocean  
   ! ======================================================================  
   REAL rmu0(klon), albedo(klon)  
   
   REAL fmagic ! un facteur magique pour regler l'albedo  
   ! cc      PARAMETER (fmagic=0.7)  
   ! ccIM => a remplacer  
   ! PARAMETER (fmagic=1.32)  
   PARAMETER (fmagic=1.0)  
   ! PARAMETER (fmagic=0.7)  
   
   REAL fauxo  
   INTEGER i  
   ! ccIM  
   LOGICAL ancien_albedo  
   PARAMETER (ancien_albedo=.FALSE.)  
   ! SAVE albedo  
   
   IF (ancien_albedo) THEN  
   
     DO i = 1, klon  
   
       rmu0(i) = max(rmu0(i), 0.0)  
   
       fauxo = (1.47-acos(rmu0(i)))/0.15  
       albedo(i) = fmagic*(.03+.630/(1.+fauxo*fauxo))  
       albedo(i) = max(min(albedo(i),0.60), 0.04)  
     END DO  
   
     ! nouvel albedo  
   
   ELSE  
   
     DO i = 1, klon  
       rmu0(i) = max(rmu0(i), 0.0)  
       ! IM:orig albedo(i) = 0.058/(rmu0(i) + 0.30)  
       albedo(i) = fmagic*0.058/(rmu0(i)+0.30)  
       albedo(i) = max(min(albedo(i),0.60), 0.04)  
     END DO  
   
   END IF  
   
   RETURN  
 END SUBROUTINE alboc_cd  
 ! ========================================================================  
Legend:
Removed from v.82  
changed lines
  Added in v.117
  ViewVC Help
Powered by ViewVC 1.1.21