/[lmdze]/trunk/dyn3d/bernoui.f
ViewVC logotype

Contents of /trunk/dyn3d/bernoui.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 68 - (show annotations)
Wed Nov 14 16:59:30 2012 UTC (11 years, 6 months ago) by guez
Original Path: trunk/libf/dyn3d/bernoui.f
File size: 1580 byte(s)
Split "flincom.f90" into "flinclo.f90", "flinfindcood.f90",
"flininfo.f90" and "flinopen_nozoom.f90", in directory
"IOIPSL/Flincom".

Renamed "etat0_lim" to "ce0l", as in LMDZ.

Split "readsulfate.f" into "readsulfate.f90", "readsulfate_preind.f90"
and "getso4fromfile.f90".

In etat0, renamed variable q3d to q, as in "dynredem1". Replaced calls
to Flicom procedures by calls to NetCDF95.

In leapfrog, added call to writehist.

Extracted ASCII art from "grid_noro" into a file
"grid_noro.txt". Transformed explicit-shape local arrays into
automatic arrays, so that test on values of iim and jjm is no longer
needed. Test on weight:
          IF (weight(ii, jj) /= 0.) THEN
is useless. There is already a test before:
    if (any(weight == 0.)) stop "zero weight in grid_noro"

In "aeropt", replaced duplicated lines with different values of inu by
a loop on inu.

Removed arguments of "conf_phys". Corresponding variables are now
defined in "physiq", in a namelist. In "conf_phys", read a namelist
instead of using getin.

1 !
2 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/bernoui.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
3 !
4 SUBROUTINE bernoui (ngrid,nlay,pphi,pecin,pbern)
5 use dimens_m
6 use paramet_m
7 use conf_gcm_m
8 use filtreg_m, only: filtreg
9 IMPLICIT NONE
10
11 c=======================================================================
12 c
13 c Auteur: P. Le Van
14 c -------
15 c
16 c Objet:
17 c ------
18 c calcul de la fonction de Bernouilli aux niveaux s .....
19 c phi et ecin sont des arguments d'entree pour le s-pg .......
20 c bern est un argument de sortie pour le s-pg ......
21 c
22 c fonction de Bernouilli = bern = filtre de( geopotentiel +
23 c energ.cinet.)
24 c
25 c=======================================================================
26 c
27 c-----------------------------------------------------------------------
28 c Decalrations:
29 c -------------
30 c
31 c
32 c Arguments:
33 c ----------
34 c
35 INTEGER nlay,ngrid
36 REAL, intent(in):: pphi(ngrid*nlay),pecin(ngrid*nlay)
37 real pbern(ngrid*nlay)
38 c
39 c Local:
40 c ------
41 c
42 INTEGER ijl
43 c
44 c-----------------------------------------------------------------------
45 c calcul de Bernouilli:
46 c ---------------------
47 c
48 DO 4 ijl = 1,ngrid*nlay
49 pbern( ijl ) = pphi( ijl ) + pecin( ijl )
50 4 CONTINUE
51 c
52 c-----------------------------------------------------------------------
53 c filtre:
54 c -------
55 c
56 CALL filtreg( pbern, jjp1, llm, 2,1, .true.)
57 c
58 c-----------------------------------------------------------------------
59 RETURN
60 END

  ViewVC Help
Powered by ViewVC 1.1.21