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

Annotation of /trunk/dyn3d/bernoui.f

Parent Directory Parent Directory | Revision Log Revision Log


Revision 68 - (hide 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 guez 3 !
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 guez 57 use conf_gcm_m
8 guez 27 use filtreg_m, only: filtreg
9 guez 3 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 guez 68 REAL, intent(in):: pphi(ngrid*nlay),pecin(ngrid*nlay)
37     real pbern(ngrid*nlay)
38 guez 3 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 guez 64 CALL filtreg( pbern, jjp1, llm, 2,1, .true.)
57 guez 3 c
58     c-----------------------------------------------------------------------
59     RETURN
60     END

  ViewVC Help
Powered by ViewVC 1.1.21