| 1 | ! $Id: pres2lev.F 1179 2009-06-11 14:18:47Z jghattas $ |
|---|
| 2 | |
|---|
| 3 | MODULE lmdz_pres2lev |
|---|
| 4 | IMPLICIT NONE; PRIVATE |
|---|
| 5 | PUBLIC pres2lev |
|---|
| 6 | |
|---|
| 7 | CONTAINS |
|---|
| 8 | |
|---|
| 9 | !****************************************************** |
|---|
| 10 | SUBROUTINE pres2lev(varo, varn, lmo, lmn, po, pn, ni, nj, ok_invertp) |
|---|
| 11 | |
|---|
| 12 | ! interpolation lineaire pour passer |
|---|
| 13 | ! a une nouvelle discretisation verticale pour |
|---|
| 14 | ! les variables de GCM |
|---|
| 15 | ! Francois Forget (01/1995) |
|---|
| 16 | ! MOdif remy roca 12/97 pour passer de pres2sig |
|---|
| 17 | ! Modif F.Codron 07/08 po en 3D |
|---|
| 18 | !********************************************************** |
|---|
| 19 | |
|---|
| 20 | IMPLICIT NONE |
|---|
| 21 | |
|---|
| 22 | ! Declarations: |
|---|
| 23 | ! ============== |
|---|
| 24 | |
|---|
| 25 | ! ARGUMENTS |
|---|
| 26 | ! """"""""" |
|---|
| 27 | LOGICAL, INTENT(IN) :: ok_invertp |
|---|
| 28 | INTEGER, INTENT(IN) :: lmo ! dimensions ancienne couches |
|---|
| 29 | INTEGER, INTENT(IN) :: lmn ! dimensions nouvelle couches |
|---|
| 30 | |
|---|
| 31 | INTEGER, INTENT(IN) :: ni, nj ! nombre de point horizontal |
|---|
| 32 | REAL, INTENT(IN) :: po(ni * nj, lmo) ! niveau de pression ancienne grille |
|---|
| 33 | REAL, INTENT(IN) :: pn(ni * nj, lmn) ! niveau de pression nouvelle grille |
|---|
| 34 | |
|---|
| 35 | REAL, INTENT(IN) :: varo(ni * nj, lmo) ! var dans l'ancienne grille |
|---|
| 36 | REAL, INTENT(OUT) :: varn(ni * nj, lmn) ! var dans la nouvelle grille |
|---|
| 37 | |
|---|
| 38 | REAL :: zvaro(ni * nj, lmo), zpo(ni * nj, lmo) |
|---|
| 39 | |
|---|
| 40 | ! Autres variables |
|---|
| 41 | ! """""""""""""""" |
|---|
| 42 | INTEGER :: ln, lo, k |
|---|
| 43 | REAL :: coef |
|---|
| 44 | |
|---|
| 45 | |
|---|
| 46 | ! Inversion de l'ordre des niveaux verticaux |
|---|
| 47 | IF (ok_invertp) THEN |
|---|
| 48 | DO lo = 1, lmo |
|---|
| 49 | DO k = 1, ni * nj |
|---|
| 50 | zpo(k, lo) = po(k, lmo + 1 - lo) |
|---|
| 51 | zvaro(k, lo) = varo(k, lmo + 1 - lo) |
|---|
| 52 | ENDDO |
|---|
| 53 | ENDDO |
|---|
| 54 | ELSE |
|---|
| 55 | DO lo = 1, lmo |
|---|
| 56 | DO k = 1, ni * nj |
|---|
| 57 | zpo(k, lo) = po(k, lo) |
|---|
| 58 | zvaro(k, lo) = varo(k, lo) |
|---|
| 59 | ENDDO |
|---|
| 60 | ENDDO |
|---|
| 61 | ENDIF |
|---|
| 62 | |
|---|
| 63 | DO ln = 1, lmn |
|---|
| 64 | DO lo = 1, lmo - 1 |
|---|
| 65 | DO k = 1, ni * nj |
|---|
| 66 | IF (pn(k, ln) >= zpo(k, 1)) THEN |
|---|
| 67 | varn(k, ln) = zvaro(k, 1) |
|---|
| 68 | ELSE IF (pn(k, ln) <= zpo(k, lmo)) THEN |
|---|
| 69 | varn(k, ln) = zvaro(k, lmo) |
|---|
| 70 | ELSE IF (pn(k, ln) <= zpo(k, lo) .AND. pn(k, ln) > zpo(k, lo + 1)) THEN |
|---|
| 71 | coef = (pn(k, ln) - zpo(k, lo)) / (zpo(k, lo + 1) - zpo(k, lo)) |
|---|
| 72 | varn(k, ln) = zvaro(k, lo) + coef * (zvaro(k, lo + 1) - zvaro(k, lo)) |
|---|
| 73 | ENDIF |
|---|
| 74 | |
|---|
| 75 | ENDDO |
|---|
| 76 | ENDDO |
|---|
| 77 | ENDDO |
|---|
| 78 | |
|---|
| 79 | END SUBROUTINE pres2lev |
|---|
| 80 | |
|---|
| 81 | END MODULE lmdz_pres2lev |
|---|
