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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 12 - (hide 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 guez 3 !
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 guez 12 CALL advect( ang, vcov, teta, w, massebx, masseby, du, dv,dteta,
99     $ conser )
100 guez 3
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