/[lmdze]/trunk/libf/phylmd/orografi.f
ViewVC logotype

Diff of /trunk/libf/phylmd/orografi.f

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

revision 12 by guez, Mon Jul 21 16:05:07 2008 UTC revision 13 by guez, Fri Jul 25 19:59:34 2008 UTC
# Line 35  c Line 35  c
35        REAL, intent(in):: dtime        REAL, intent(in):: dtime
36        REAL, intent(in):: paprs(klon,klev+1)        REAL, intent(in):: paprs(klon,klev+1)
37        REAL, intent(in):: pplay(klon,klev)        REAL, intent(in):: pplay(klon,klev)
38        REAL pmea(nlon),pstd(nlon),psig(nlon),pgam(nlon),pthe(nlon)        REAL pmea(nlon)
39          real, intent(in):: pstd(nlon),psig(nlon)
40          real pgam(nlon),pthe(nlon)
41        REAL ppic(nlon),pval(nlon)        REAL ppic(nlon),pval(nlon)
42        REAL pulow(nlon),pvlow(nlon),pustr(nlon),pvstr(nlon)        REAL pulow(nlon),pvlow(nlon),pustr(nlon),pvstr(nlon)
43        REAL t(nlon,nlev), u(nlon,nlev), v(nlon,nlev)        REAL t(nlon,nlev), u(nlon,nlev), v(nlon,nlev)
# Line 203  c Line 205  c
205        real  pum1(nlon,nlev),        real  pum1(nlon,nlev),
206       *      pvm1(nlon,nlev),       *      pvm1(nlon,nlev),
207       *      ptm1(nlon,nlev),       *      ptm1(nlon,nlev),
208       *      pmea(nlon),pstd(nlon),psig(nlon),       *      pmea(nlon)
209       *      pgamma(nlon),ptheta(nlon),ppic(nlon),pval(nlon),        real, intent(in):: pstd(nlon)
210          real, intent(in):: psig(nlon)
211          real pgamma(nlon),ptheta(nlon),ppic(nlon),pval(nlon),
212       *      pgeom1(nlon,nlev),       *      pgeom1(nlon,nlev),
213       *      papm1(nlon,nlev),       *      papm1(nlon,nlev),
214       *      paphm1(nlon,nlev+1)       *      paphm1(nlon,nlev+1)
# Line 494  c Line 498  c
498       *     pzdep(nlon,klev)       *     pzdep(nlon,klev)
499         real pulow(nlon),pvlow(nlon),ptheta(nlon),pgamma(nlon),pnu(nlon),         real pulow(nlon),pvlow(nlon),ptheta(nlon),pgamma(nlon),pnu(nlon),
500       *     pd1(nlon),pd2(nlon),pdmod(nlon)       *     pd1(nlon),pd2(nlon),pdmod(nlon)
501        real pstd(nlon),pmea(nlon),ppic(nlon),pval(nlon)        real, intent(in):: pstd(nlon)
502          real pmea(nlon),ppic(nlon),pval(nlon)
503  c  c
504  c-----------------------------------------------------------------------  c-----------------------------------------------------------------------
505  c  c
# Line 951  c Line 956  c
956  c  c
957        real prho(nlon,nlev+1),pstab(nlon,nlev+1),ptau(nlon,nlev+1),        real prho(nlon,nlev+1),pstab(nlon,nlev+1),ptau(nlon,nlev+1),
958       *     pvph(nlon,nlev+1),       *     pvph(nlon,nlev+1),
959       *     pgeom1(nlon,nlev),pstd(nlon)       *     pgeom1(nlon,nlev)
960          real, intent(in):: pstd(nlon)
961  c  c
962        real psig(nlon)        real, intent(in):: psig(nlon)
963        real pmea(nlon),ppic(nlon)        real pmea(nlon),ppic(nlon)
964        real pdmod(nlon)        real pdmod(nlon)
965  c  c
# Line 1092  C Line 1098  C
1098       *     PRHO  (NLON,NLEV+1), PVPH (NLON,NLEV+1),       *     PRHO  (NLON,NLEV+1), PVPH (NLON,NLEV+1),
1099       *     PRI   (NLON,NLEV+1), PTAU(NLON,NLEV+1)       *     PRI   (NLON,NLEV+1), PTAU(NLON,NLEV+1)
1100            
1101        REAL pdmod (NLON) , psig(NLON),        REAL pdmod (NLON)
1102       *     pvar(NLON)        real, intent(in):: psig(NLON)
1103          real, intent(in):: pvar(NLON)
1104            
1105  C-----------------------------------------------------------------------  C-----------------------------------------------------------------------
1106  C  C
# Line 1328  c Line 1335  c
1335        REAL, intent(in):: pplay(klon,klev)        REAL, intent(in):: pplay(klon,klev)
1336        REAL, intent(in):: plat(nlon)        REAL, intent(in):: plat(nlon)
1337        real pmea(nlon)        real pmea(nlon)
1338        REAL pstd(nlon)        REAL, intent(in):: pstd(nlon)
1339        REAL ppic(nlon)        REAL ppic(nlon)
1340        REAL pulow(nlon),pvlow(nlon),pustr(nlon),pvstr(nlon)        REAL pulow(nlon),pvlow(nlon),pustr(nlon),pvstr(nlon)
1341        REAL t(nlon,nlev), u(nlon,nlev), v(nlon,nlev)        REAL t(nlon,nlev), u(nlon,nlev), v(nlon,nlev)
# Line 1459  C Line 1466  C
1466       *      PVOM(NLON,NLEV),       *      PVOM(NLON,NLEV),
1467       *      PULOW(NLON),       *      PULOW(NLON),
1468       *      PVLOW(NLON)       *      PVLOW(NLON)
1469        REAL  PUM1(NLON,NLEV),        REAL  PUM1(NLON,NLEV), PVM1(NLON,NLEV), PTM1(NLON,NLEV)
      *      PVM1(NLON,NLEV),  
      *      PTM1(NLON,NLEV)  
1470        real, intent(in):: PLAT(NLON)        real, intent(in):: PLAT(NLON)
1471        real PMEA(NLON),        real PMEA(NLON)
1472       *      PVAROR(NLON),        real, intent(in):: PVAROR(NLON)
1473       *      ppic(NLON),        real ppic(NLON), PGEOM1(NLON,NLEV), PAPHM1(NLON,NLEV+1)
      *      PGEOM1(NLON,NLEV),  
      *      PAPHM1(NLON,NLEV+1)  
1474  C  C
1475        INTEGER  KTEST(NLON)        INTEGER  KTEST(NLON)
1476        real, intent(in):: ptsphy        real, intent(in):: ptsphy
# Line 1601  C Line 1604  C
1604        DO 301 JL=1,klon        DO 301 JL=1,klon
1605        IF(KTEST(JL).EQ.1) THEN        IF(KTEST(JL).EQ.1) THEN
1606         ZTAU(JL,KLEV+1)= - GKLIFT*ZRHO(JL,KLEV+1)*2.*ROMEGA*         ZTAU(JL,KLEV+1)= - GKLIFT*ZRHO(JL,KLEV+1)*2.*ROMEGA*
 C    *                 (2*PVAROR(JL)+PMEA(JL))*  
1607       *                 2*PVAROR(JL)*       *                 2*PVAROR(JL)*
1608       *                 SIN(ZPI/180.*PLAT(JL))*PVLOW(JL)       *                 SIN(ZPI/180.*PLAT(JL))*PVLOW(JL)
1609         ZTAV(JL,KLEV+1)=   GKLIFT*ZRHO(JL,KLEV+1)*2.*ROMEGA*         ZTAV(JL,KLEV+1)=   GKLIFT*ZRHO(JL,KLEV+1)*2.*ROMEGA*
 C    *                 (2*PVAROR(JL)+PMEA(JL))*  
1610       *                 2*PVAROR(JL)*       *                 2*PVAROR(JL)*
1611       *                 SIN(ZPI/180.*PLAT(JL))*PULOW(JL)       *                 SIN(ZPI/180.*PLAT(JL))*PULOW(JL)
1612        ELSE        ELSE
Legend:
Removed from v.12  
changed lines
  Added in v.13
  ViewVC Help
Powered by ViewVC 1.1.21