| 1 | ! |
|---|
| 2 | ! $Header: /home/cvsroot/LMDZ4/libf/phylmd/plevel.F,v 1.1.1.1.10.1 2006/08/17 15:41:51 fairhead Exp $ |
|---|
| 3 | ! |
|---|
| 4 | c================================================================ |
|---|
| 5 | c================================================================ |
|---|
| 6 | SUBROUTINE plevel_new(ilon,ilev,klevSTD,lnew,pgcm,pres,Qgcm,Qpres) |
|---|
| 7 | c================================================================ |
|---|
| 8 | c================================================================ |
|---|
| 9 | USE dimphy |
|---|
| 10 | IMPLICIT none |
|---|
| 11 | |
|---|
| 12 | cym#include "dimensions.h" |
|---|
| 13 | cy#include "dimphy.h" |
|---|
| 14 | |
|---|
| 15 | c================================================================ |
|---|
| 16 | c |
|---|
| 17 | c Interpoler des champs 3-D u, v et g du modele a un niveau de |
|---|
| 18 | c pression donnee (pres) |
|---|
| 19 | c |
|---|
| 20 | c INPUT: ilon ----- nombre de points |
|---|
| 21 | c ilev ----- nombre de couches |
|---|
| 22 | c lnew ----- true si on doit reinitialiser les poids |
|---|
| 23 | c pgcm ----- pressions modeles |
|---|
| 24 | c pres ----- pression vers laquelle on interpolle |
|---|
| 25 | c Qgcm ----- champ GCM |
|---|
| 26 | c Qpres ---- champ interpolle au niveau pres |
|---|
| 27 | c |
|---|
| 28 | c================================================================ |
|---|
| 29 | c |
|---|
| 30 | c arguments : |
|---|
| 31 | c ----------- |
|---|
| 32 | |
|---|
| 33 | INTEGER ilon, ilev, klevSTD |
|---|
| 34 | logical lnew |
|---|
| 35 | |
|---|
| 36 | REAL pgcm(ilon,ilev) |
|---|
| 37 | REAL Qgcm(ilon,ilev) |
|---|
| 38 | real pres(klevSTD) |
|---|
| 39 | REAL Qpres(ilon, klevSTD) |
|---|
| 40 | |
|---|
| 41 | c local : |
|---|
| 42 | c ------- |
|---|
| 43 | |
|---|
| 44 | cym INTEGER lt(klon), lb(klon) |
|---|
| 45 | cym REAL ptop, pbot, aist(klon), aisb(klon) |
|---|
| 46 | |
|---|
| 47 | cym save lt,lb,ptop,pbot,aist,aisb |
|---|
| 48 | INTEGER,ALLOCATABLE,SAVE,DIMENSION(:) :: lt,lb |
|---|
| 49 | REAL,ALLOCATABLE,SAVE,DIMENSION(:,:) :: aist,aisb |
|---|
| 50 | c$OMP THREADPRIVATE(lt,lb,aist,aisb) |
|---|
| 51 | REAL,SAVE :: ptop, pbot |
|---|
| 52 | c$OMP THREADPRIVATE(ptop, pbot) |
|---|
| 53 | LOGICAL,SAVE :: first = .true. |
|---|
| 54 | INTEGER :: nlev |
|---|
| 55 | c$OMP THREADPRIVATE(first) |
|---|
| 56 | INTEGER i, k |
|---|
| 57 | c |
|---|
| 58 | if (first) then |
|---|
| 59 | allocate(lt(klon),lb(klon)) |
|---|
| 60 | allocate(aist(klon,klevSTD),aisb(klon, klevSTD)) |
|---|
| 61 | first=.false. |
|---|
| 62 | endif |
|---|
| 63 | |
|---|
| 64 | c===================================================================== |
|---|
| 65 | if (lnew) then |
|---|
| 66 | c on reinitialise les reindicages et les poids |
|---|
| 67 | c===================================================================== |
|---|
| 68 | |
|---|
| 69 | |
|---|
| 70 | c Chercher les 2 couches les plus proches du niveau a obtenir |
|---|
| 71 | c |
|---|
| 72 | c Eventuellement, faire l'extrapolation a partir des deux couches |
|---|
| 73 | c les plus basses ou les deux couches les plus hautes: |
|---|
| 74 | c |
|---|
| 75 | c |
|---|
| 76 | DO nlev = 1, klevSTD |
|---|
| 77 | DO i = 1, klon |
|---|
| 78 | IF ( ABS(pres(nlev)-pgcm(i,ilev) ) .LT. |
|---|
| 79 | & ABS(pres(nlev)-pgcm(i,1)) ) THEN |
|---|
| 80 | lt(i) = ilev ! 2 |
|---|
| 81 | lb(i) = ilev-1 ! 1 |
|---|
| 82 | ELSE |
|---|
| 83 | lt(i) = 2 |
|---|
| 84 | lb(i) = 1 |
|---|
| 85 | ENDIF |
|---|
| 86 | ENDDO |
|---|
| 87 | DO k = 1, ilev-1 |
|---|
| 88 | DO i = 1, klon |
|---|
| 89 | pbot = pgcm(i,k) |
|---|
| 90 | ptop = pgcm(i,k+1) |
|---|
| 91 | IF (ptop.LE.pres(nlev) .AND. pbot.GE.pres(nlev)) THEN |
|---|
| 92 | lt(i) = k+1 |
|---|
| 93 | lb(i) = k |
|---|
| 94 | ENDIF |
|---|
| 95 | ENDDO |
|---|
| 96 | ENDDO |
|---|
| 97 | |
|---|
| 98 | c Interpolation lineaire: |
|---|
| 99 | DO i = 1, klon |
|---|
| 100 | c interpolation en logarithme de pression: |
|---|
| 101 | c |
|---|
| 102 | c ... Modif . P. Le Van ( 20/01/98) .... |
|---|
| 103 | c Modif Frederic Hourdin (3/01/02) |
|---|
| 104 | |
|---|
| 105 | aist(i,nlev) = LOG( pgcm(i,lb(i))/ pres(nlev) ) |
|---|
| 106 | & / LOG( pgcm(i,lb(i))/ pgcm(i,lt(i)) ) |
|---|
| 107 | aisb(i,nlev) = LOG( pres(nlev) / pgcm(i,lt(i)) ) |
|---|
| 108 | & / LOG( pgcm(i,lb(i))/ pgcm(i,lt(i))) |
|---|
| 109 | ENDDO |
|---|
| 110 | ENDDO |
|---|
| 111 | |
|---|
| 112 | ENDIF ! lnew |
|---|
| 113 | |
|---|
| 114 | c====================================================================== |
|---|
| 115 | c inteprollation |
|---|
| 116 | c ET je mets les vents a zero quand je rencontre une montagne |
|---|
| 117 | c====================================================================== |
|---|
| 118 | |
|---|
| 119 | DO nlev = 1, klevSTD |
|---|
| 120 | DO i=1,klon |
|---|
| 121 | IF (pgcm(i,1).LT.pres(nlev)) THEN |
|---|
| 122 | c Qpres(i)=1e33 |
|---|
| 123 | Qpres(i,nlev) = 1e+20 |
|---|
| 124 | ELSE |
|---|
| 125 | Qpres(i,nlev) = |
|---|
| 126 | & Qgcm(i,lb(i))*aisb(i,nlev) + |
|---|
| 127 | & Qgcm(i,lt(i))*aist(i,nlev) |
|---|
| 128 | ENDIF |
|---|
| 129 | ENDDO |
|---|
| 130 | ENDDO |
|---|
| 131 | |
|---|
| 132 | c |
|---|
| 133 | RETURN |
|---|
| 134 | END |
|---|
