source: LMDZ6/branches/Amaury_dev/libf/phylmd/write_histday_seri.h @ 5151

Last change on this file since 5151 was 5117, checked in by abarral, 4 months ago

rename modules properly lmdz_*
move some unused files to obsolete/
(lint) uppercase fortran keywords
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.9 KB
Line 
1
2! $Header$
3
4      IF (is_sequential) THEN
5     
6      IF (type_run.EQ."AMIP") THEN
7
8      ndex2d = 0
9      itau_w = itau_phy + itap + start_time * day_step_phy
10
11! Champs 2D:
12
13      pi = ACOS(-1.)
14      pir = 4.0*ATAN(1.0) / 180.0
15
16      DO i=1, klon
17       zx_tmp_fi2d(i)=(topsw(i)-toplw(i))
18      ENDDO
19
20      ok_msk=.FALSE.
21      msk(1:klon)=pctsrf(1:klon,is_ter)
22      CALL moyglo_pondaire(klon, zx_tmp_fi2d, cell_area,  &
23           ok_msk, msk, moyglo)
24      zx_tmp_fi2d(1:klon)=moyglo
25
26      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
27      CALL histwrite(nid_day_seri,"bilTOA",itau_w, &
28                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
29
30      ok_msk=.FALSE.
31      CALL moyglo_pondaire(klon, bils, cell_area,  &
32           ok_msk, msk, moyglo)
33      zx_tmp_fi2d(1:klon)=moyglo
34
35      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
36      CALL histwrite(nid_day_seri,"bils",itau_w, &
37                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
38
39      DO k=1, klev
40      DO i=1, klon
41!IM 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
45
46      CALL moyglo_pondaima(klon, klev, zx_tmp_fi3d,  &
47           cell_area, paprs, moyglo)
48      zx_tmp_fi2d(1:klon)=moyglo
49
50      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
51      CALL histwrite(nid_day_seri,"ecin",itau_w, &
52                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d) 
53
54
55
56      CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,cell_area,zx_tmp_2d)
57      airetot=0.
58      DO i=1, klon
59       airetot=airetot+cell_area(i)
60      ENDDO
61!     IF(itap.EQ.1) PRINT*,'airetotphy=',airetot
62
63      airetot=0.
64      DO j=1, nbp_lat
65       DO i=1, nbp_lon
66        airetot=airetot+zx_tmp_2d(i,j)
67       ENDDO
68      ENDDO
69
70!     IF(itap.EQ.1) PRINT*,'airetotij=',airetot,
71!    $ '4piR2',4.*pi*RA*RA
72
73      zx_tmp_fi2d(1:klon)=aam/airetot
74      CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
75      CALL histwrite(nid_day_seri,"momang",itau_w,zx_tmp_2d, &
76                     nbp_lon*nbp_lat,ndex2d)
77
78      zx_tmp_fi2d(1:klon)=torsfc/airetot
79      CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
80      CALL histwrite(nid_day_seri,"torsfc",itau_w,zx_tmp_2d, &
81                     nbp_lon*nbp_lat,ndex2d)
82
83!IM 151004 END
84
85      CALL moyglo_pondmass(klon, klev, t_seri, &
86           cell_area, paprs, moyglo)
87      zx_tmp_fi2d(1:klon)=moyglo
88
89      CALL gr_fi_ecrit(1,klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
90      CALL histwrite(nid_day_seri,"tamv",itau_w, &
91                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
92
93      ok_msk=.FALSE.
94      CALL moyglo_pondaire(klon, paprs(:,1), cell_area,  &
95           ok_msk, msk, moyglo)
96      zx_tmp_fi2d(1:klon)=moyglo
97
98      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
99      CALL histwrite(nid_day_seri,"psol",itau_w, &
100                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
101
102      ok_msk=.FALSE.
103      CALL moyglo_pondaire(klon, evap, cell_area,  &
104           ok_msk, msk, moyglo)
105      zx_tmp_fi2d(1:klon)=moyglo
106
107      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d,zx_tmp_2d)
108      CALL histwrite(nid_day_seri,"evap",itau_w, &
109                     zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
110
111!     DO i=1, klon
112!      zx_tmp_fi2d(i)=SnowFrac(i,is_ter)
113!     ENDDO
114
115!     ok_msk=.TRUE.
116!     msk(1:klon)=pctsrf(1:klon,is_ter)
117!     CALL moyglo_pondaire(klon, zx_tmp_fi2d, cell_area, 
118!    .                     ok_msk, msk, moyglo)
119!     zx_tmp_fi2d(1:klon)=moyglo
120
121!     CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
122!     CALL histwrite(nid_day_seri,"SnowFrac",
123!    .               itau_w,zx_tmp_2d,nbp_lon*nbp_lat,ndex2d) 
124
125!     DO i=1, klon
126!IM 080904    zx_tmp_fi2d(i)=zsnow_mass(i)/330.*rowl
127!      zx_tmp_fi2d(i)=zsnow_mass(i)
128!     ENDDO
129
130!IM 140904   ok_msk=.FALSE.
131!     ok_msk=.TRUE.
132!     msk(1:klon)=pctsrf(1:klon,is_ter)
133!     CALL moyglo_pondaire(klon, zx_tmp_fi2d, cell_area, 
134!    .     ok_msk, msk, moyglo)
135!     zx_tmp_fi2d(1:klon)=moyglo
136
137!     CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat,zx_tmp_fi2d,zx_tmp_2d)
138!     CALL histwrite(nid_day_seri,"snow_depth",itau_w,
139!    .               zx_tmp_2d,nbp_lon*nbp_lat,ndex2d)
140
141      DO i=1, klon
142       zx_tmp_fi2d(i)=ftsol(i,is_oce)
143      ENDDO
144
145      ok_msk=.TRUE.
146      msk(1:klon)=pctsrf(1:klon,is_oce)
147      CALL moyglo_pondaire(klon, zx_tmp_fi2d, cell_area,  &
148           ok_msk, msk, moyglo)
149      zx_tmp_fi2d(1:klon)=moyglo
150
151      CALL gr_fi_ecrit(1, klon,nbp_lon,nbp_lat, zx_tmp_fi2d, zx_tmp_2d)
152      CALL histwrite(nid_day_seri,"tsol_"//clnsurf(is_oce), &
153                     itau_w,zx_tmp_2d,nbp_lon*nbp_lat,ndex2d) 
154
155!=================================================================
156!=================================================================
157!=================================================================
158
159      IF (ok_sync) THEN
160        CALL histsync(nid_day_seri)
161      endif
162
163      ENDIF !fin test sur type_run.EQ."AMIP"
164     
165      ENDIF  ! mono_cpu
Note: See TracBrowser for help on using the repository browser.