libf/phylmd/coefkz2.f90

SUBROUTINE coefkz2(nsrf, knon, paprs, pplay, t, pcfm, pcfh)

  ! J'introduit un peu de diffusion sauf dans les endroits
  ! ou une forte inversion est presente
  ! On peut dire qu'il represente la convection peu profonde

  use dimens_m
  use indicesol
  use dimphy
  use iniprint
  use SUPHEC_M

  IMPLICIT none

  ! Arguments:
  ! nsrf-----input-I- indicateur de la nature du sol
  ! knon-----input-I- nombre de points a traiter
  ! paprs----input-R- pression a chaque intercouche (en Pa)
  ! pplay----input-R- pression au milieu de chaque couche (en Pa)
  ! t--------input-R- temperature (K)

  ! pcfm-----output-R- coefficients a calculer (vitesse)
  ! pcfh-----output-R- coefficients a calculer (chaleur et humidite)

  ! Arguments:

  INTEGER knon, nsrf
  REAL paprs(klon, klev+1), pplay(klon, klev)
  REAL t(klon, klev)

  REAL pcfm(klon, klev), pcfh(klon, klev)

  ! Quelques constantes et options:

  REAL prandtl
  PARAMETER (prandtl=0.4)
  REAL kstable
  PARAMETER (kstable=0.002)
  REAL mixlen ! constante controlant longueur de melange
  PARAMETER (mixlen=35.0)
  REAL seuil ! au-dela l'inversion est consideree trop faible
  PARAMETER (seuil=-0.02)

  ! Variables locales:

  INTEGER i, k, invb(knon)
  REAL zl2(knon)
  REAL zdthmin(knon), zdthdp

  !----------------------------------------------------------

  ! Initialiser les sorties
  DO k = 1, klev
     DO i = 1, knon
        pcfm(i, k) = 0.0
        pcfh(i, k) = 0.0
     ENDDO
  ENDDO

  ! Chercher la zone d'inversion forte

  DO i = 1, knon
     invb(i) = klev
     zdthmin(i)=0.0
  ENDDO
  DO k = 2, klev/ 2 - 1
     DO i = 1, knon
        zdthdp = (t(i, k) - t(i, k + 1)) / (pplay(i, k) - pplay(i, k + 1)) &
             - RD * 0.5 * (t(i, k) + t(i, k + 1)) / RCPD / paprs(i, k + 1)
        zdthdp = zdthdp * 100.
        IF (pplay(i, k) > 0.8 * paprs(i, 1) .AND. zdthdp < zdthmin(i)) THEN
           zdthmin(i) = zdthdp
           invb(i) = k
        ENDIF
     ENDDO
  ENDDO

  ! Introduire une diffusion:
  DO k = 2, klev
     DO i = 1, knon
        ! si on est sur ocean et s'il n'y a pas d'inversion ou si
        ! l'inversion est trop faible:
        IF ((nsrf.EQ.is_oce) .AND.  &  
             ((invb(i).EQ.klev) .OR. (zdthmin(i) > seuil))) THEN
           zl2(i)=(mixlen*MAX(0.0, (paprs(i, k)-paprs(i, klev+1)) &
                /(paprs(i, 2)-paprs(i, klev+1))))**2
           pcfm(i, k)= zl2(i)* kstable
           pcfh(i, k) = pcfm(i, k) /prandtl ! h et m different
        ENDIF
     ENDDO
  ENDDO

END SUBROUTINE coefkz2
1	SUBROUTINE coefkz2(nsrf, knon, paprs, pplay, t, pcfm, pcfh)
2
3	! J'introduit un peu de diffusion sauf dans les endroits
4	! ou une forte inversion est presente
5	! On peut dire qu'il represente la convection peu profonde
6
7	use dimens_m
8	use indicesol
9	use dimphy
10	use iniprint
11	use SUPHEC_M
12
13	IMPLICIT none
14
15	! Arguments:
16	! nsrf-----input-I- indicateur de la nature du sol
17	! knon-----input-I- nombre de points a traiter
18	! paprs----input-R- pression a chaque intercouche (en Pa)
19	! pplay----input-R- pression au milieu de chaque couche (en Pa)
20	! t--------input-R- temperature (K)
21
22	! pcfm-----output-R- coefficients a calculer (vitesse)
23	! pcfh-----output-R- coefficients a calculer (chaleur et humidite)
24
25	! Arguments:
26
27	INTEGER knon, nsrf
28	REAL paprs(klon, klev+1), pplay(klon, klev)
29	REAL t(klon, klev)
30
31	REAL pcfm(klon, klev), pcfh(klon, klev)
32
33	! Quelques constantes et options:
34
35	REAL prandtl
36	PARAMETER (prandtl=0.4)
37	REAL kstable
38	PARAMETER (kstable=0.002)
39	REAL mixlen ! constante controlant longueur de melange
40	PARAMETER (mixlen=35.0)
41	REAL seuil ! au-dela l'inversion est consideree trop faible
42	PARAMETER (seuil=-0.02)
43
44	! Variables locales:
45
46	INTEGER i, k, invb(knon)
47	REAL zl2(knon)
48	REAL zdthmin(knon), zdthdp
49
50	!----------------------------------------------------------
51
52	! Initialiser les sorties
53	DO k = 1, klev
54	DO i = 1, knon
55	pcfm(i, k) = 0.0
56	pcfh(i, k) = 0.0
57	ENDDO
58	ENDDO
59
60	! Chercher la zone d'inversion forte
61
62	DO i = 1, knon
63	invb(i) = klev
64	zdthmin(i)=0.0
65	ENDDO
66	DO k = 2, klev/ 2 - 1
67	DO i = 1, knon
68	zdthdp = (t(i, k) - t(i, k + 1)) / (pplay(i, k) - pplay(i, k + 1)) &
69	- RD * 0.5 * (t(i, k) + t(i, k + 1)) / RCPD / paprs(i, k + 1)
70	zdthdp = zdthdp * 100.
71	IF (pplay(i, k) > 0.8 * paprs(i, 1) .AND. zdthdp < zdthmin(i)) THEN
72	zdthmin(i) = zdthdp
73	invb(i) = k
74	ENDIF
75	ENDDO
76	ENDDO
77
78	! Introduire une diffusion:
79	DO k = 2, klev
80	DO i = 1, knon
81	! si on est sur ocean et s'il n'y a pas d'inversion ou si
82	! l'inversion est trop faible:
83	IF ((nsrf.EQ.is_oce) .AND. &
84	((invb(i).EQ.klev) .OR. (zdthmin(i) > seuil))) THEN
85	zl2(i)=(mixlen*MAX(0.0, (paprs(i, k)-paprs(i, klev+1)) &
86	/(paprs(i, 2)-paprs(i, klev+1))))**2
87	pcfm(i, k)= zl2(i)* kstable
88	pcfh(i, k) = pcfm(i, k) /prandtl ! h et m different
89	ENDIF
90	ENDDO
91	ENDDO
92
93	END SUBROUTINE coefkz2