source: LMDZ4/trunk/libf/phytherm/write_histday_seri.h @ 1080

Last change on this file since 1080 was 814, checked in by Laurent Fairhead, 17 years ago

Rajout de la physique utilisant les thermiques FH
LF
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 6.9 KB
RevLine 
[814]1c
2c $Header$
3c
4      IF (is_sequential) THEN
5     
6      IF (type_run.EQ."AMIP") THEN
7c
8      ndex2d = 0
9      itau_w = itau_phy + itap
10c
11c Champs 2D:
12c
13      pi = ACOS(-1.)
14      pir = 4.0*ATAN(1.0) / 180.0
15c
16      DO i=1, klon
17       zx_tmp_fi2d(i)=(topsw(i)-toplw(i))
18      ENDDO
19c
20      ok_msk=.FALSE.
21      msk(1:klon)=pctsrf(1:klon,is_ter)
22      CALL moyglo_pondaire(klon, zx_tmp_fi2d, airephy, 
23     .     ok_msk, msk, moyglo)
24      zx_tmp_fi2d(1:klon)=moyglo
25c
26      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
27      CALL histwrite(nid_day_seri,"bilTOA",itau_w,
28     .               zx_tmp_2d,iim*jjmp1,ndex2d)
29c
30      ok_msk=.FALSE.
31      CALL moyglo_pondaire(klon, bils, airephy, 
32     .     ok_msk, msk, moyglo)
33      zx_tmp_fi2d(1:klon)=moyglo
34c
35      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
36      CALL histwrite(nid_day_seri,"bils",itau_w,
37     .               zx_tmp_2d,iim*jjmp1,ndex2d)
38c
39      DO k=1, klev
40      DO i=1, klon
41cIM 080904    zx_tmp_fi3d(i,k)=u(i,k)**2+v(i,k)**2
42       zx_tmp_fi3d(i,k)=(u(i,k)**2+v(i,k)**2)/2.
43      ENDDO
44      ENDDO
45c
46      CALL moyglo_pondaima(klon, klev, zx_tmp_fi3d, 
47     .     airephy, paprs, moyglo)
48      zx_tmp_fi2d(1:klon)=moyglo
49c
50      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
51      CALL histwrite(nid_day_seri,"ecin",itau_w,
52     .               zx_tmp_2d,iim*jjmp1,ndex2d) 
53c
54cIM 151004 BEG
55      IF(1.EQ.0) THEN
56c
57      DO k=1, klev
58      DO i=1, klon
59       zx_tmp_fi3d(i,k)=u_seri(i,k)*RA*cos(pir* rlat(i))
60      ENDDO
61      ENDDO
62c
63      CALL moyglo_pondaima(klon, klev, zx_tmp_fi3d, 
64     .     airephy, paprs, moyglo)
65      zx_tmp_fi2d(1:klon)=moyglo
66c
67      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
68      CALL histwrite(nid_day_seri,"momang",itau_w,zx_tmp_2d,
69     .               iim*jjmp1,ndex2d)
70c
71c friction torque
72c
73      DO i=1, klon
74       zx_tmp_fi2d(i)=zxfluxu(i,1)*RA* cos(pir* rlat(i))
75      ENDDO
76c
77      ok_msk=.FALSE.
78      CALL moyglo_pondaire(klon, zx_tmp_fi2d, airephy, 
79     .     ok_msk, msk, moyglo)
80      zx_tmp_fi2d(1:klon)=moyglo
81c
82      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
83      CALL histwrite(nid_day_seri,"frictor",itau_w,zx_tmp_2d,
84     .               iim*jjmp1,ndex2d)
85c
86c mountain torque
87c
88cIM 190504 BEG
89      CALL gr_fi_dyn(1,klon,iim+1,jjm+1,airephy,airedyn)
90      CALL gr_fi_dyn(klev+1,klon,iim+1,jjm+1,paprs,padyn)
91      CALL gr_fi_dyn(1,klon,iim+1,jjm+1,rlat,rlatdyn)
92      mountor=0.
93      airetot=0.
94      DO j = 1, jjmp1
95       DO i = 1, iim+1
96        ij=i+(iim+1)*(j-1)
97        zx_tmp(ij)=0.
98        DO k = 1, klev
99         zx_tmp(ij)=zx_tmp(ij)+dudyn(i,j,k)*airedyn(i,j)*
100     $              (padyn(i,j,k+1)-padyn(i,j,k))/RG
101         airetot=airetot+airedyn(i,j)
102        ENDDO
103cIM 190504 mountor=mountor+zx_tmp(ij)*airedyn(i,j)*RA*
104        mountor=mountor+zx_tmp(ij)*RA*
105     $           cos(pir* rlatdyn(i,j))
106       ENDDO
107      ENDDO
108cIM 151004 BEG
109      IF(itap.EQ.1) PRINT*,'airetot=',airetot,airetot/klev
110cIM 151004 END
111cIM 190504      mountor=mountor/(airetot*airetot)
112      mountor=mountor/airetot
113c
114cIM 190504 END
115      zx_tmp_2d(1:iim,1:jjmp1)=mountor
116      CALL histwrite(nid_day_seri,"mountor",itau_w,zx_tmp_2d,
117     .               iim*jjmp1,ndex2d)
118c
119      ENDIF !(1.EQ.0) THEN
120c
121c
122      CALL gr_fi_dyn(1,klon,iim+1,jjm+1,airephy,airedyn)
123      CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d)
124      airetot=0.
125c     DO j = 1, jjmp1
126c      DO i = 1, iim+1
127c       ij=i+(iim+1)*(j-1)
128c       DO k = 1, klev
129c        airetot=airetot+airedyn(i,j)
130c        airetot=airetot+airedyn(i,j)
131c       ENDDO !k
132c      ENDDO !i
133c     ENDDO !j
134c
135      DO i=1, klon
136       airetot=airetot+airephy(i)
137      ENDDO
138c     IF(itap.EQ.1) PRINT*,'airetotphy=',airetot
139c
140      airetot=0.
141      DO j=1, jjmp1
142       DO i=1, iim
143        airetot=airetot+zx_tmp_2d(i,j)
144       ENDDO
145      ENDDO
146c
147c     IF(itap.EQ.1) PRINT*,'airetotij=',airetot,
148c    $ '4piR2',4.*pi*RA*RA
149c
150      zx_tmp_fi2d(1:klon)=aam/airetot
151      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
152      CALL histwrite(nid_day_seri,"momang",itau_w,zx_tmp_2d,
153     .               iim*jjmp1,ndex2d)
154c
155      zx_tmp_fi2d(1:klon)=torsfc/airetot
156      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
157      CALL histwrite(nid_day_seri,"torsfc",itau_w,zx_tmp_2d,
158     .               iim*jjmp1,ndex2d)
159c
160cIM 151004 END
161c
162      CALL moyglo_pondmass(klon, klev, t_seri,
163     .     airephy, paprs, moyglo)
164      zx_tmp_fi2d(1:klon)=moyglo
165c
166      CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
167      CALL histwrite(nid_day_seri,"tamv",itau_w,
168     .               zx_tmp_2d,iim*jjmp1,ndex2d)
169c
170      ok_msk=.FALSE.
171      CALL moyglo_pondaire(klon, paprs(:,1), airephy, 
172     .     ok_msk, msk, moyglo)
173      zx_tmp_fi2d(1:klon)=moyglo
174c
175      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
176      CALL histwrite(nid_day_seri,"psol",itau_w,
177     .               zx_tmp_2d,iim*jjmp1,ndex2d)
178c
179      ok_msk=.FALSE.
180      CALL moyglo_pondaire(klon, evap, airephy, 
181     .     ok_msk, msk, moyglo)
182      zx_tmp_fi2d(1:klon)=moyglo
183c
184      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d)
185      CALL histwrite(nid_day_seri,"evap",itau_w,
186     .               zx_tmp_2d,iim*jjmp1,ndex2d)
187c
188c     DO i=1, klon
189c      zx_tmp_fi2d(i)=SnowFrac(i,is_ter)
190c     ENDDO
191c
192c     ok_msk=.TRUE.
193c     msk(1:klon)=pctsrf_new(1:klon,is_ter)
194c     CALL moyglo_pondaire(klon, zx_tmp_fi2d, airephy, 
195c    .                     ok_msk, msk, moyglo)
196c     zx_tmp_fi2d(1:klon)=moyglo
197c
198c     CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
199c     CALL histwrite(nid_day_seri,"SnowFrac",
200c    .               itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 
201c
202c     DO i=1, klon
203cIM 080904    zx_tmp_fi2d(i)=zsnow_mass(i)/330.*rowl
204c      zx_tmp_fi2d(i)=zsnow_mass(i)
205c     ENDDO
206c
207cIM 140904   ok_msk=.FALSE.
208c     ok_msk=.TRUE.
209c     msk(1:klon)=pctsrf(1:klon,is_ter)
210c     CALL moyglo_pondaire(klon, zx_tmp_fi2d, airephy, 
211c    .     ok_msk, msk, moyglo)
212c     zx_tmp_fi2d(1:klon)=moyglo
213c
214c     CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d)
215c     CALL histwrite(nid_day_seri,"snow_depth",itau_w,
216c    .               zx_tmp_2d,iim*jjmp1,ndex2d)
217c
218      DO i=1, klon
219       zx_tmp_fi2d(i)=ftsol(i,is_oce)
220      ENDDO
221c
222      ok_msk=.TRUE.
223      msk(1:klon)=pctsrf(1:klon,is_oce)
224      CALL moyglo_pondaire(klon, zx_tmp_fi2d, airephy, 
225     .     ok_msk, msk, moyglo)
226      zx_tmp_fi2d(1:klon)=moyglo
227c
228      CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d)
229      CALL histwrite(nid_day_seri,"tsol_"//clnsurf(is_oce),
230     $               itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 
231c
232c=================================================================
233c=================================================================
234c=================================================================
235c
236      if (ok_sync) then
237        call histsync(nid_day_seri)
238      endif
239c
240      ENDIF !fin test sur type_run.EQ."AMIP"
241     
242      ENDIF  ! mono_cpu
Note: See TracBrowser for help on using the repository browser.