/[lmdze]/trunk/libf/dyn3d/caldyn.f
ViewVC logotype

Contents of /trunk/libf/dyn3d/caldyn.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 12 - (show annotations)
Mon Jul 21 16:05:07 2008 UTC (15 years, 10 months ago) by guez
File size: 3878 byte(s) 
-- Minor modification of input/output:

Created procedure "read_logic". Variables of module "logic" are read
by "read_logic" instead of "conf_gcm". Variable "offline" of module
"conf_gcm" is read from namelist instead of "*.def".

Deleted arguments "dtime", "co2_ppm_etat0", "solaire_etat0",
"tabcntr0" and local variables "radpas", "tab_cntrl" of
"phyetat0". "phyetat0" does not read "controle" in "startphy.nc" any
longer. "phyetat0" now reads global attribute "itau_phy" from
"startphy.nc". "phyredem" does not create variable "controle" in
"startphy.nc" any longer. "phyredem" now writes global attribute
"itau_phy" of "startphy.nc". Deleted argument "tabcntr0" of
"printflag". Removed diagnostic messages written by "printflag" for
comparison of the variable "controle" of "startphy.nc" and the
variables read from "*.def" or namelist input.

-- Removing unwanted functionality:

Removed variable "lunout" from module "iniprint", replaced everywhere
by standard output.

Removed case "ocean == 'couple'" in "clmain", "interfsurf_hq" and
"physiq". Removed procedure "interfoce_cpl".

-- Should not change anything at run time:

Automated creation of graphs in documentation. More documentation on
input files.

Converted Fortran files to free format: "phyredem.f90", "printflag.f90".

Split module "clesphy" into "clesphys" and "clesphys2".

Removed variables "conser", "leapf", "forward", "apphys", "apdiss" and
"statcl" from module "logic". Added arguments "conser" to "advect",
"leapf" to "integrd". Added local variables "forward", "leapf",
"apphys", "conser", "apdiss" in "leapfrog".

Added intent attributes.

Deleted arguments "dtime" of "phyredem", "pdtime" of "flxdtdq", "sh"
of "phytrac", "dt" of "yamada".

Deleted local variables "dtime", "co2_ppm_etat0", "solaire_etat0",
"length", "tabcntr0" in "physiq". Replaced all references to "dtime"
by references to "pdtphys".
1 !
2 ! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/caldyn.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
3 !
4 c
5 c
6 SUBROUTINE caldyn
7 $ (itau,ucov,vcov,teta,ps,masse,pk,pkf,phis ,
8 $ phi,conser,du,dv,dteta,dp,w,pbaru,pbarv,time )
9
10 use dimens_m
11 use paramet_m
12 use comconst
13 use comvert
14 use comgeom
15 use pression_m, only: pression
16
17 IMPLICIT NONE
18
19 c=======================================================================
20 c
21 c Auteur : P. Le Van
22 c
23 c Objet:
24 c ------
25 c
26 c Calcul des tendances dynamiques.
27 c
28 c Modif 04/93 F.Forget
29 c=======================================================================
30
31 c-----------------------------------------------------------------------
32 c 0. Declarations:
33 c ----------------
34
35
36 c Arguments:
37 c ----------
38
39 LOGICAL conser
40
41 INTEGER, intent(in):: itau
42 REAL vcov(ip1jm,llm),ucov(ip1jmp1,llm),teta(ip1jmp1,llm)
43 REAL ps(ip1jmp1),phis(ip1jmp1)
44 REAL, intent(in):: pk(iip1,jjp1,llm)
45 real pkf(ip1jmp1,llm)
46 REAL vcont(ip1jm,llm),ucont(ip1jmp1,llm)
47 REAL phi(ip1jmp1,llm),masse(ip1jmp1,llm)
48 REAL dv(ip1jm,llm),du(ip1jmp1,llm)
49 REAL dteta(ip1jmp1,llm),dp(ip1jmp1)
50 REAL pbaru(ip1jmp1,llm),pbarv(ip1jm,llm)
51 REAL time
52
53 c Local:
54 c ------
55
56 REAL ang(ip1jmp1,llm),p(ip1jmp1,llmp1)
57 REAL massebx(ip1jmp1,llm),masseby(ip1jm,llm),psexbarxy(ip1jm)
58 REAL vorpot(ip1jm,llm)
59 REAL w(ip1jmp1,llm),ecin(ip1jmp1,llm),convm(ip1jmp1,llm)
60 REAL bern(ip1jmp1,llm)
61 REAL massebxy(ip1jm,llm)
62
63 INTEGER ij,l
64
65 c-----------------------------------------------------------------------
66 c Calcul des tendances dynamiques:
67 c --------------------------------
68
69 CALL covcont ( llm , ucov , vcov , ucont, vcont )
70 CALL pression ( ip1jmp1, ap , bp , ps , p )
71 CALL psextbar ( ps , psexbarxy )
72 CALL massdair ( p , masse )
73 CALL massbar ( masse, massebx , masseby )
74 call massbarxy( masse, massebxy )
75 CALL flumass ( massebx, masseby , vcont, ucont ,pbaru, pbarv )
76 CALL dteta1 ( teta , pbaru , pbarv, dteta )
77 CALL convmas ( pbaru, pbarv , convm )
78
79 DO ij =1, ip1jmp1
80 dp( ij ) = convm( ij,1 ) / airesurg( ij )
81 ENDDO
82
83 CALL vitvert ( convm , w )
84 CALL tourpot ( vcov , ucov , massebxy , vorpot )
85 CALL dudv1 ( vorpot , pbaru , pbarv , du , dv )
86 CALL enercin ( vcov , ucov , vcont , ucont , ecin )
87 CALL bernoui ( ip1jmp1, llm , phi , ecin , bern )
88 CALL dudv2 ( teta , pkf , bern , du , dv )
89
90
91 DO l=1,llm
92 DO ij=1,ip1jmp1
93 ang(ij,l) = ucov(ij,l) + constang(ij)
94 ENDDO
95 ENDDO
96
97
98 CALL advect( ang, vcov, teta, w, massebx, masseby, du, dv,dteta,
99 $ conser )
100
101 C WARNING probleme de peridocite de dv sur les PC/linux. Pb d'arrondi
102 C probablement. Observe sur le code compile avec pgf90 3.0-1
103
104 DO l = 1, llm
105 DO ij = 1, ip1jm, iip1
106 IF( dv(ij,l).NE.dv(ij+iim,l) ) THEN
107 c PRINT *,'!!!ATTENTION!!! probleme de periodicite sur vcov',
108 c , ' dans caldyn'
109 c PRINT *,' l, ij = ', l, ij, ij+iim,dv(ij+iim,l),dv(ij,l)
110 dv(ij+iim,l) = dv(ij,l)
111 endif
112 enddo
113 enddo
114 c-----------------------------------------------------------------------
115 c Sorties eventuelles des variables de controle:
116 c ----------------------------------------------
117
118 IF( conser ) THEN
119 CALL sortvarc
120 $ ( itau,ucov,teta,ps,masse,pk,phis,vorpot,phi,bern,dp,time,vcov )
121
122 ENDIF
123
124 RETURN
125 END
  ViewVC Help
Powered by ViewVC 1.1.21