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	module alboc_m
2
3	IMPLICIT NONE
4
5	contains
6
7	SUBROUTINE alboc(rjour, rlat, albedo)
8
9	! From LMDZ4/libf/phylmd/albedo.F, version 1.2 2005/02/07 15:00:52
10
11	! 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
16	USE dimphy, only: klon
17	USE yomcst, only: r_incl
18	USE orbite_m, ONLY: orbite
19
20	! 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
25	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
34	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
43	IF (ancien_albedo) THEN
44
45	zpi = 4.*atan(1.)
46
47	! Calculer la longitude vraie de l'orbite terrestre:
48	CALL orbite(rjour, zlonsun, zdist)
49
50	! 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
53	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
57	! 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
65	! 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
82	! nouvel albedo
83
84	ELSE
85
86	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