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

Last change on this file since 671 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
RevLine 
[524]1!
2! $Header$
3!
4c     IF (ok_mensuelNMC) THEN
5      IF (ok_mensuel) THEN
6c
[644]7       ndex3d = 0
8       itau_w = itau_phy + itap
9ccc
10c  Champs interpolles sur des niveaux de pression du NMC
[524]11c
[644]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)
[524]20c
[644]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)
[524]25c
[644]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)
[524]30c
[644]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)
[524]35c
[644]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
[524]89         ENDIF
[644]90        ENDDO
91       ENDDO !k=1, nlevSTD
[524]92c
[644]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)
[524]96c
[644]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)
[524]101c
[644]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)
[524]106c
[644]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)
[524]111c
[644]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)
[524]116c
[644]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)
[524]121c
[644]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)
[524]126c
[644]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
[524]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.