phylmd/Interface_surf/alboc.f

module alboc_m

  IMPLICIT NONE

contains

  SUBROUTINE alboc(rjour, rlat, albedo)

    ! From LMDZ4/libf/phylmd/albedo.F, version 1.2 2005/02/07 15:00:52

    ! Auteur(s): Z.X. Li (LMD/CNRS) (adaptation du GCM du LMD)
    ! Date: le 16 mars 1995
    ! Objet: Calculer l'albedo sur l'ocean
    ! Methode: Integrer numeriquement l'albedo pendant une journee

    USE dimphy, only: klon
    USE yomcst, only: r_incl
    USE orbite_m, ONLY: orbite

    ! Arguments;
    ! rjour (in, R) : jour dans l'annee (a compter du 1 janvier)
    ! rlat (in, R) : latitude en degre
    ! albedo (out, R): albedo obtenu (de 0 a 1)

    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)
    INTEGER npts ! il controle la precision de l'integration
    PARAMETER (npts=120) ! 120 correspond a l'interval 6 minutes

    REAL rlat(klon), rjour, albedo(klon)
    REAL zdist, zlonsun, zpi, zdeclin
    REAL rmu, alb, srmu, salb, fauxo, aa, bb
    INTEGER i, k
    ! ccIM
    LOGICAL ancien_albedo
    PARAMETER (ancien_albedo=.FALSE.)
    ! SAVE albedo

    IF (ancien_albedo) THEN

       zpi = 4.*atan(1.)

       ! Calculer la longitude vraie de l'orbite terrestre:
       CALL orbite(rjour, zlonsun, zdist)

       ! Calculer la declinaison du soleil (qui varie entre + et - R_incl):
       zdeclin = asin(sin(zlonsun*zpi/180.0)*sin(r_incl*zpi/180.0))

       DO i = 1, klon
          aa = sin(rlat(i)*zpi/180.0)*sin(zdeclin)
          bb = cos(rlat(i)*zpi/180.0)*cos(zdeclin)

          ! Midi local (angle du temps = 0.0):
          rmu = aa + bb*cos(0.0)
          rmu = max(0.0, rmu)
          fauxo = (1.47-acos(rmu))/.15
          alb = 0.03 + 0.630/(1.+fauxo*fauxo)
          srmu = rmu
          salb = alb*rmu

          ! Faire l'integration numerique de midi a minuit (le facteur 2
          ! prend en compte l'autre moitie de la journee):
          DO k = 1, npts
             rmu = aa + bb*cos(float(k)/float(npts)*zpi)
             rmu = max(0.0, rmu)
             fauxo = (1.47-acos(rmu))/.15
             alb = 0.03 + 0.630/(1.+fauxo*fauxo)
             srmu = srmu + rmu*2.0
             salb = salb + alb*rmu*2.0
          END DO
          IF (srmu/=0.0) THEN
             albedo(i) = salb/srmu*fmagic
          ELSE ! nuit polaire (on peut prendre une valeur quelconque)
             albedo(i) = fmagic
          END IF
       END DO

       ! nouvel albedo

    ELSE

       zpi = 4.*atan(1.)

       ! Calculer la longitude vraie de l'orbite terrestre:
       CALL orbite(rjour, zlonsun, zdist)

       ! Calculer la declinaison du soleil (qui varie entre + et - R_incl):
       zdeclin = asin(sin(zlonsun*zpi/180.0)*sin(r_incl*zpi/180.0))

       DO i = 1, klon
          aa = sin(rlat(i)*zpi/180.0)*sin(zdeclin)
          bb = cos(rlat(i)*zpi/180.0)*cos(zdeclin)

          ! Midi local (angle du temps = 0.0):
          rmu = aa + bb*cos(0.0)
          rmu = max(0.0, rmu)
          ! IM cf. PB alb = 0.058/(rmu + 0.30)
          ! alb = 0.058/(rmu + 0.30) * 1.5
          alb = 0.058/(rmu+0.30)*1.2
          ! alb = 0.058/(rmu + 0.30) * 1.3
          srmu = rmu
          salb = alb*rmu

          ! Faire l'integration numerique de midi a minuit (le facteur 2
          ! prend en compte l'autre moitie de la journee):
          DO k = 1, npts
             rmu = aa + bb*cos(float(k)/float(npts)*zpi)
             rmu = max(0.0, rmu)
             ! IM cf. PB alb = 0.058/(rmu + 0.30)
             ! alb = 0.058/(rmu + 0.30) * 1.5
             alb = 0.058/(rmu+0.30)*1.2
             ! alb = 0.058/(rmu + 0.30) * 1.3
             srmu = srmu + rmu*2.0
             salb = salb + alb*rmu*2.0
          END DO
          IF (srmu/=0.0) THEN
             albedo(i) = salb/srmu*fmagic
          ELSE ! nuit polaire (on peut prendre une valeur quelconque)
             albedo(i) = fmagic
          END IF
       END DO
    END IF

  END SUBROUTINE alboc

end module alboc_m
1	guez	125	module alboc_m
2
3	guez	81	IMPLICIT NONE
4
5	guez	125	contains
6	guez	81
7	guez	125	SUBROUTINE alboc(rjour, rlat, albedo)
8	guez	81
9	guez	125	! From LMDZ4/libf/phylmd/albedo.F, version 1.2 2005/02/07 15:00:52
10	guez	81
11	guez	125	! Auteur(s): Z.X. Li (LMD/CNRS) (adaptation du GCM du LMD)
12			! Date: le 16 mars 1995
13			! Objet: Calculer l'albedo sur l'ocean
14			! Methode: Integrer numeriquement l'albedo pendant une journee
15	guez	81
16	guez	125	USE dimphy, only: klon
17			USE yomcst, only: r_incl
18			USE orbite_m, ONLY: orbite
19	guez	81
20	guez	125	! Arguments;
21			! rjour (in, R) : jour dans l'annee (a compter du 1 janvier)
22			! rlat (in, R) : latitude en degre
23			! albedo (out, R): albedo obtenu (de 0 a 1)
24	guez	81
25	guez	125	REAL fmagic ! un facteur magique pour regler l'albedo
26			! cc PARAMETER (fmagic=0.7)
27			! ccIM => a remplacer
28			! PARAMETER (fmagic=1.32)
29			PARAMETER (fmagic=1.0)
30			! PARAMETER (fmagic=0.7)
31			INTEGER npts ! il controle la precision de l'integration
32			PARAMETER (npts=120) ! 120 correspond a l'interval 6 minutes
33	guez	81
34	guez	125	REAL rlat(klon), rjour, albedo(klon)
35			REAL zdist, zlonsun, zpi, zdeclin
36			REAL rmu, alb, srmu, salb, fauxo, aa, bb
37			INTEGER i, k
38			! ccIM
39			LOGICAL ancien_albedo
40			PARAMETER (ancien_albedo=.FALSE.)
41			! SAVE albedo
42	guez	81
43	guez	125	IF (ancien_albedo) THEN
44	guez	81
45	guez	125	zpi = 4.*atan(1.)
46	guez	81
47	guez	125	! Calculer la longitude vraie de l'orbite terrestre:
48			CALL orbite(rjour, zlonsun, zdist)
49	guez	81
50	guez	125	! Calculer la declinaison du soleil (qui varie entre + et - R_incl):
51			zdeclin = asin(sin(zlonsunzpi/180.0)sin(r_incl*zpi/180.0))
52	guez	81
53	guez	125	DO i = 1, klon
54			aa = sin(rlat(i)zpi/180.0)sin(zdeclin)
55			bb = cos(rlat(i)zpi/180.0)cos(zdeclin)
56	guez	81
57	guez	125	! Midi local (angle du temps = 0.0):
58			rmu = aa + bb*cos(0.0)
59			rmu = max(0.0, rmu)
60			fauxo = (1.47-acos(rmu))/.15
61			alb = 0.03 + 0.630/(1.+fauxo*fauxo)
62			srmu = rmu
63			salb = alb*rmu
64	guez	81
65	guez	125	! Faire l'integration numerique de midi a minuit (le facteur 2
66			! prend en compte l'autre moitie de la journee):
67			DO k = 1, npts
68			rmu = aa + bbcos(float(k)/float(npts)zpi)
69			rmu = max(0.0, rmu)
70			fauxo = (1.47-acos(rmu))/.15
71			alb = 0.03 + 0.630/(1.+fauxo*fauxo)
72			srmu = srmu + rmu*2.0
73			salb = salb + albrmu2.0
74			END DO
75			IF (srmu/=0.0) THEN
76			albedo(i) = salb/srmu*fmagic
77			ELSE ! nuit polaire (on peut prendre une valeur quelconque)
78			albedo(i) = fmagic
79			END IF
80			END DO
81	guez	81
82	guez	125	! nouvel albedo
83	guez	81
84	guez	125	ELSE
85	guez	81
86	guez	125	zpi = 4.*atan(1.)
87
88			! Calculer la longitude vraie de l'orbite terrestre:
89			CALL orbite(rjour, zlonsun, zdist)
90
91			! Calculer la declinaison du soleil (qui varie entre + et - R_incl):
92			zdeclin = asin(sin(zlonsunzpi/180.0)sin(r_incl*zpi/180.0))
93
94			DO i = 1, klon
95			aa = sin(rlat(i)zpi/180.0)sin(zdeclin)
96			bb = cos(rlat(i)zpi/180.0)cos(zdeclin)
97
98			! Midi local (angle du temps = 0.0):
99			rmu = aa + bb*cos(0.0)
100			rmu = max(0.0, rmu)
101			! IM cf. PB alb = 0.058/(rmu + 0.30)
102			! alb = 0.058/(rmu + 0.30) * 1.5
103			alb = 0.058/(rmu+0.30)*1.2
104			! alb = 0.058/(rmu + 0.30) * 1.3
105			srmu = rmu
106			salb = alb*rmu
107
108			! Faire l'integration numerique de midi a minuit (le facteur 2
109			! prend en compte l'autre moitie de la journee):
110			DO k = 1, npts
111			rmu = aa + bbcos(float(k)/float(npts)zpi)
112			rmu = max(0.0, rmu)
113			! IM cf. PB alb = 0.058/(rmu + 0.30)
114			! alb = 0.058/(rmu + 0.30) * 1.5
115			alb = 0.058/(rmu+0.30)*1.2
116			! alb = 0.058/(rmu + 0.30) * 1.3
117			srmu = srmu + rmu*2.0
118			salb = salb + albrmu2.0
119			END DO
120			IF (srmu/=0.0) THEN
121			albedo(i) = salb/srmu*fmagic
122			ELSE ! nuit polaire (on peut prendre une valeur quelconque)
123			albedo(i) = fmagic
124			END IF
125			END DO
126			END IF
127
128			END SUBROUTINE alboc
129
130			end module alboc_m