source: LMDZ4/trunk/libf/phylmd/write_histmthNMC.h @ 644

Last change on this file since 644 was 644, checked in by Laurent Fairhead, 20 years ago

Synchronisation avec tous les diagnostiques de Ionela IM
Inclusion du slab ocean IM
LF
  • 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))/FLOAT(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.