/[lmdze]/trunk/Sources/dyn3d/Guide/Read_reanalyse/reanalyse2nat.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/Guide/Read_reanalyse/reanalyse2nat.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC revision 171 by guez, Tue Sep 29 19:48:59 2015 UTC
# Line 18  contains Line 18  contains
18      USE exner_hyb_m, ONLY: exner_hyb      USE exner_hyb_m, ONLY: exner_hyb
19      use massbar_m, only: massbar      use massbar_m, only: massbar
20      use massdair_m, only: massdair      use massdair_m, only: massdair
21        use pres2lev_m, only: pres2lev
22    
23      integer nlevnc      integer nlevnc
24      real, intent(in):: psi(iip1, jjp1)      real, intent(in):: psi(iip1, jjp1)
# Line 42  contains Line 43  contains
43    
44      real p(iip1, jjp1, llmp1)      real p(iip1, jjp1, llmp1)
45      real pks(iip1, jjp1)      real pks(iip1, jjp1)
     real pkf(iip1, jjp1, llm)  
46      real pls(iip1, jjp1, llm)      real pls(iip1, jjp1, llm)
47      real prefkap, unskap      real prefkap, unskap
48    
# Line 53  contains Line 53  contains
53      ! calcul de la pression au milieu des couches      ! calcul de la pression au milieu des couches
54      forall (l = 1: llm + 1) p(:, :, l) = ap(l) + bp(l) * psi      forall (l = 1: llm + 1) p(:, :, l) = ap(l) + bp(l) * psi
55      call massdair(p, masse)      call massdair(p, masse)
56      CALL exner_hyb(psi, p, pks, pk, pkf)      CALL exner_hyb(psi, p, pks, pk)
57    
58      ! Calcul de pls, pression au milieu des couches, en Pascals      ! Calcul de pls, pression au milieu des couches, en Pascals
59      unskap=1./kappa      unskap=1./kappa
Legend:
Removed from v.134  
changed lines
  Added in v.171
  ViewVC Help
Powered by ViewVC 1.1.21