source: LMDZ4/branches/LMDZ4_V2_patch/libf/phylmd/write_histmthNMC.h @ 5023

Last change on this file since 5023 was 723, checked in by lmdzadmin, 18 years ago

On passe a des ecrit_ins, ecrit_day, etc en nombre de jours (REAL)
On lit frequence ecriture traceurs ecrit_trac dans physiq.def
Correction petits pbs ini_histrac.h
IM
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.8 KB
Line 
1!
2! $Header$
3!
4c     IF (ok_mensuelNMC) THEN
5      IF (ok_mensuel) THEN
6c
7       ndex3d = 0
8       itau_w = itau_phy + itap
9ccc
10c  Champs interpolles sur des niveaux de pression du NMC
11c
12c     PARAMETER(nout=3) !nout=1 : day; =2 : mth; =3 : NMC
13ccc
14      IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN
15ccc
16       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,tsumSTD(:,:,2),
17     $      zx_tmp_3d)
18       CALL histwrite(nid_nmc,"temp",itau_w,zx_tmp_3d,
19     $       iim*jjmp1*nlevSTD,ndex3d)
20c
21       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,phisumSTD(:,:,2),
22     $     zx_tmp_3d)
23       CALL histwrite(nid_nmc,"phi",itau_w,zx_tmp_3d,
24     $       iim*jjmp1*nlevSTD,ndex3d)
25c
26       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,qsumSTD(:,:,2),
27     $     zx_tmp_3d)
28       CALL histwrite(nid_nmc,"q",itau_w,zx_tmp_3d,
29     $       iim*jjmp1*nlevSTD,ndex3d)
30c
31       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,rhsumSTD(:,:,2),
32     $     zx_tmp_3d)
33       CALL histwrite(nid_nmc,"rh",itau_w,zx_tmp_3d,
34     $       iim*jjmp1*nlevSTD,ndex3d)
35c
36       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,usumSTD(:,:,2),
37     $     zx_tmp_3d)
38       CALL histwrite(nid_nmc,"u",itau_w,zx_tmp_3d,
39     $       iim*jjmp1*nlevSTD,ndex3d)
40c
41       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vsumSTD(:,:,2),
42     $     zx_tmp_3d)
43       CALL histwrite(nid_nmc,"v",itau_w,zx_tmp_3d,
44     $       iim*jjmp1*nlevSTD,ndex3d)
45ccc
46      ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN
47ccc
48       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,tsumSTD(:,:,3),
49     $     zx_tmp_3d)
50       CALL histwrite(nid_nmc,"temp",itau_w,zx_tmp_3d,
51     $       iim*jjmp1*nlevSTD,ndex3d)
52c
53       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,phisumSTD(:,:,3),
54     $     zx_tmp_3d)
55       CALL histwrite(nid_nmc,"phi",itau_w,zx_tmp_3d,
56     $       iim*jjmp1*nlevSTD,ndex3d)
57c
58       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,qsumSTD(:,:,3),
59     $     zx_tmp_3d)
60       CALL histwrite(nid_nmc,"q",itau_w,zx_tmp_3d,
61     $       iim*jjmp1*nlevSTD,ndex3d)
62c
63       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,rhsumSTD(:,:,3),
64     $     zx_tmp_3d)
65       CALL histwrite(nid_nmc,"rh",itau_w,zx_tmp_3d,
66     $       iim*jjmp1*nlevSTD,ndex3d)
67c
68       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,usumSTD(:,:,3),
69     $     zx_tmp_3d)
70       CALL histwrite(nid_nmc,"u",itau_w,zx_tmp_3d,
71     $       iim*jjmp1*nlevSTD,ndex3d)
72c
73       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vsumSTD(:,:,3),
74     $     zx_tmp_3d)
75       CALL histwrite(nid_nmc,"v",itau_w,zx_tmp_3d,
76     $       iim*jjmp1*nlevSTD,ndex3d)
77c
78       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,wsumSTD(:,:,3),
79     $     zx_tmp_3d)
80       CALL histwrite(nid_nmc,"w",itau_w,zx_tmp_3d,
81     $       iim*jjmp1*nlevSTD,ndex3d)
82c
83       DO k=1, nlevSTD
84        DO i=1, klon
85         IF(tnondef(i,k,3).NE.1.E+20) THEN
86          zx_tmp_fi3d(i,k) = (100.*tnondef(i,k,3))/ecrit_hf2mth
87         ELSE
88          zx_tmp_fi3d(i,k) = 1.E+20
89         ENDIF
90        ENDDO
91       ENDDO !k=1, nlevSTD
92c
93       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d)
94       CALL histwrite(nid_nmc,"psbg",itau_w,zx_tmp_3d,
95     $       iim*jjmp1*nlevSTD,ndex3d)
96c
97       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,uvsumSTD(:,:,3),
98     $     zx_tmp_3d)
99       CALL histwrite(nid_nmc,"uv",itau_w,zx_tmp_3d,
100     $       iim*jjmp1*nlevSTD,ndex3d)
101c
102       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vqsumSTD(:,:,3),
103     $     zx_tmp_3d)
104       CALL histwrite(nid_nmc,"vq",itau_w,zx_tmp_3d,
105     $       iim*jjmp1*nlevSTD,ndex3d)
106c
107       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vTsumSTD(:,:,3),
108     $     zx_tmp_3d)
109       CALL histwrite(nid_nmc,"vT",itau_w,zx_tmp_3d,
110     $       iim*jjmp1*nlevSTD,ndex3d)
111c
112       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1, wqsumSTD(:,:,3),
113     $     zx_tmp_3d)
114       CALL histwrite(nid_nmc,"wq",itau_w,zx_tmp_3d,
115     $       iim*jjmp1*nlevSTD,ndex3d)
116c
117       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,vphisumSTD(:,:,3),
118     $     zx_tmp_3d)
119       CALL histwrite(nid_nmc,"vphi",itau_w,zx_tmp_3d,
120     $       iim*jjmp1*nlevSTD,ndex3d)
121c
122       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,wTsumSTD(:,:,3),
123     $     zx_tmp_3d)
124       CALL histwrite(nid_nmc,"wT",itau_w,zx_tmp_3d,
125     $       iim*jjmp1*nlevSTD,ndex3d)
126c
127       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,u2sumSTD(:,:,3),
128     $     zx_tmp_3d)
129       CALL histwrite(nid_nmc,"uxu",itau_w,zx_tmp_3d,
130     $       iim*jjmp1*nlevSTD,ndex3d)
131c
132       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,v2sumSTD(:,:,3),
133     $     zx_tmp_3d)
134       CALL histwrite(nid_nmc,"vxv",itau_w,zx_tmp_3d,
135     $       iim*jjmp1*nlevSTD,ndex3d)
136c
137       CALL gr_fi_ecrit(nlevSTD, klon,iim,jjmp1,T2sumSTD(:,:,3),
138     $     zx_tmp_3d)
139       CALL histwrite(nid_nmc,"TxT",itau_w,zx_tmp_3d,
140     $       iim*jjmp1*nlevSTD,ndex3d)
141c
142      ENDIF !type_run
143c
144      if (ok_sync) then
145        call histsync(nid_nmc)
146      endif
147
148      ENDIF
Note: See TracBrowser for help on using the repository browser.