/[lmdze]/trunk/Sources/phylmd/calltherm.f
ViewVC logotype

Annotation of /trunk/Sources/phylmd/calltherm.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
Original Path: trunk/libf/phylmd/calltherm.f
File size: 3431 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/phylmd/calltherm.F,v 1.2 2004/12/10 11:27:46 lmdzadmin Exp $
3     !
4     subroutine calltherm(dtime
5     s ,pplay,paprs,pphi
6     s ,u_seri,v_seri,t_seri,q_seri
7     s ,d_u_ajs,d_v_ajs,d_t_ajs,d_q_ajs
8     s ,fm_therm,entr_therm)
9    
10     use dimens_m
11     use dimphy
12     use ctherm
13     implicit none
14    
15 guez 12 REAL, intent(in):: dtime
16 guez 3
17     REAL u_seri(klon,klev),v_seri(klon,klev)
18     REAL t_seri(klon,klev),q_seri(klon,klev)
19     REAL, intent(in):: paprs(klon,klev+1)
20 guez 10 REAL, intent(in):: pplay(klon,klev)
21 guez 3 REAL pphi(klon,klev)
22    
23     CFH Update Thermiques
24     REAL d_t_ajs(klon,klev), d_q_ajs(klon,klev)
25     REAL d_u_ajs(klon,klev),d_v_ajs(klon,klev)
26     real fm_therm(klon,klev+1),entr_therm(klon,klev)
27    
28    
29     c variables locales
30     REAL d_t_the(klon,klev), d_q_the(klon,klev)
31     REAL d_u_the(klon,klev),d_v_the(klon,klev)
32     c
33     real zfm_therm(klon,klev+1),zentr_therm(klon,klev),zdt
34     save zentr_therm,zfm_therm
35    
36     integer i,k, isplit
37    
38     *********************************************************
39    
40     c Modele du thermique
41     c ===================
42     c print*,'thermiques: WARNING on passe t au lieu de t_seri'
43     print*,'avant isplit ',nsplit_thermals
44    
45    
46     fm_therm(:,:)=0.
47     entr_therm(:,:)=0.
48    
49     c tests sur les valeurs negatives de l'eau
50     do k=1,klev
51     do i=1,klon
52     if (.not.q_seri(i,k).ge.0.) then
53     print*,'WARN eau<0 avant therm i=',i,' k=',k
54     s ,' dq,q',d_q_the(i,k),q_seri(i,k)
55     q_seri(i,k)=1.e-15
56     endif
57     enddo
58     enddo
59    
60    
61     zdt=dtime/float(nsplit_thermals)
62     do isplit=1,nsplit_thermals
63    
64     CALL thermcell(klon,klev,zdt
65     s ,pplay,paprs,pphi
66     s ,u_seri,v_seri,t_seri,q_seri
67     s ,d_u_the,d_v_the,d_t_the,d_q_the
68     s ,zfm_therm,zentr_therm
69     s ,r_aspect_thermals,l_mix_thermals,w2di_thermals
70     s ,tho_thermals)
71    
72     c transformation de la derivee en tendance
73     d_t_the(:,:)=d_t_the(:,:)*dtime/float(nsplit_thermals)
74     d_u_the(:,:)=d_u_the(:,:)*dtime/float(nsplit_thermals)
75     d_v_the(:,:)=d_v_the(:,:)*dtime/float(nsplit_thermals)
76     d_q_the(:,:)=d_q_the(:,:)*dtime/float(nsplit_thermals)
77     fm_therm(:,:)=fm_therm(:,:)
78     s +zfm_therm(:,:)/float(nsplit_thermals)
79     entr_therm(:,:)=entr_therm(:,:)
80     s +zentr_therm(:,:)/float(nsplit_thermals)
81     fm_therm(:,klev+1)=0.
82    
83    
84    
85     c accumulation de la tendance
86     d_t_ajs(:,:)=d_t_ajs(:,:)+d_t_the(:,:)
87     d_u_ajs(:,:)=d_u_ajs(:,:)+d_u_the(:,:)
88     d_v_ajs(:,:)=d_v_ajs(:,:)+d_v_the(:,:)
89     d_q_ajs(:,:)=d_q_ajs(:,:)+d_q_the(:,:)
90    
91     c incrementation des variables meteo
92     t_seri(:,:) = t_seri(:,:) + d_t_the(:,:)
93     u_seri(:,:) = u_seri(:,:) + d_u_the(:,:)
94     v_seri(:,:) = v_seri(:,:) + d_v_the(:,:)
95     q_seri(:,:) = q_seri(:,:) + d_q_the(:,:)
96    
97     c tests sur les valeurs negatives de l'eau
98     DO k = 1, klev
99     DO i = 1, klon
100     if (.not.q_seri(i,k).ge.0.) then
101     print*,'WARN eau<0 apres therm i=',i,' k=',k
102     s ,' dq,q',d_q_the(i,k),q_seri(i,k)
103     q_seri(i,k)=1.e-15
104     endif
105     ENDDO
106     ENDDO
107    
108     enddo ! isplit
109    
110     return
111    
112     end
  ViewVC Help
Powered by ViewVC 1.1.21