1 | ! |
---|
2 | ! $Header$ |
---|
3 | ! |
---|
4 | IF (ok_histNMC(3)) THEN |
---|
5 | c |
---|
6 | ndex3d = 0 |
---|
7 | itau_w = itau_phy + itap |
---|
8 | ccc |
---|
9 | c Champs interpolles sur des niveaux de pression du NMC |
---|
10 | c |
---|
11 | c PARAMETER(nout=3) |
---|
12 | c nout=1 : in=pdtphys, out=mth |
---|
13 | c nout=2 : in=pdtphys, out=day |
---|
14 | c nout=3 : in=pdtphys, out=hf |
---|
15 | ccc |
---|
16 | CALL histwrite_phy(nid_hfnmc,lNMC,"tnondef",itau_w, |
---|
17 | $tnondef(:,:,3)) |
---|
18 | c |
---|
19 | CALL histwrite_phy(nid_hfnmc,lNMC,"ta",itau_w, |
---|
20 | $twriteSTD3) |
---|
21 | c |
---|
22 | CALL histwrite_phy(nid_hfnmc,lNMC,"zg",itau_w, |
---|
23 | $phiwriteSTD3) |
---|
24 | c |
---|
25 | CALL histwrite_phy(nid_hfnmc,lNMC,"hus",itau_w, |
---|
26 | $qwriteSTD3) |
---|
27 | c |
---|
28 | CALL histwrite_phy(nid_hfnmc,lNMC,"hur",itau_w, |
---|
29 | $rhwriteSTD3) |
---|
30 | c |
---|
31 | CALL histwrite_phy(nid_hfnmc,lNMC,"ua",itau_w, |
---|
32 | $uwriteSTD3) |
---|
33 | c |
---|
34 | CALL histwrite_phy(nid_hfnmc,lNMC,"va",itau_w, |
---|
35 | $vwriteSTD3) |
---|
36 | c |
---|
37 | CALL histwrite_phy(nid_hfnmc,lNMC,"wap",itau_w, |
---|
38 | $wwriteSTD3) |
---|
39 | c |
---|
40 | IF (1.EQ.0) THEN |
---|
41 | c |
---|
42 | DO k=1, nlevSTD |
---|
43 | DO i=1, klon |
---|
44 | IF(tnondef(i,k,3).NE.missing_val) THEN |
---|
45 | zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,3))/freq_moyNMC(3) |
---|
46 | ELSE |
---|
47 | zx_tmp_fiNC(i,k) = missing_val |
---|
48 | ENDIF |
---|
49 | ENDDO |
---|
50 | ENDDO !k=1, nlevSTD |
---|
51 | c |
---|
52 | CALL histwrite_phy(nid_hfnmc,lNMC,"psbg",itau_w, |
---|
53 | $zx_tmp_fiNC) |
---|
54 | c |
---|
55 | CALL histwrite_phy(nid_hfnmc,lNMC,"uv",itau_w, |
---|
56 | $uvsumSTD(:,:,3)) |
---|
57 | c |
---|
58 | CALL histwrite_phy(nid_hfnmc,lNMC,"vq",itau_w, |
---|
59 | $vqsumSTD(:,:,3)) |
---|
60 | c |
---|
61 | CALL histwrite_phy(nid_hfnmc,lNMC,"vT",itau_w, |
---|
62 | $vTsumSTD(:,:,3)) |
---|
63 | c |
---|
64 | CALL histwrite_phy(nid_hfnmc,lNMC,"wq",itau_w, |
---|
65 | $wqsumSTD(:,:,3)) |
---|
66 | c |
---|
67 | CALL histwrite_phy(nid_hfnmc,lNMC,"vphi",itau_w, |
---|
68 | $vphisumSTD(:,:,3)) |
---|
69 | c |
---|
70 | CALL histwrite_phy(nid_hfnmc,lNMC,"wT",itau_w, |
---|
71 | $wTsumSTD(:,:,3)) |
---|
72 | c |
---|
73 | CALL histwrite_phy(nid_hfnmc,lNMC,"uxu",itau_w, |
---|
74 | $u2sumSTD(:,:,3)) |
---|
75 | c |
---|
76 | CALL histwrite_phy(nid_hfnmc,lNMC,"vxv",itau_w, |
---|
77 | $v2sumSTD(:,:,3)) |
---|
78 | c |
---|
79 | CALL histwrite_phy(nid_hfnmc,lNMC,"TxT",itau_w, |
---|
80 | $T2sumSTD(:,:,3)) |
---|
81 | c |
---|
82 | ENDIF !(1.EQ.0) THEN |
---|
83 | c |
---|
84 | if (ok_sync) then |
---|
85 | c$OMP MASTER |
---|
86 | call histsync(nid_hfnmc) |
---|
87 | c$OMP END MASTER |
---|
88 | endif |
---|
89 | c |
---|
90 | ENDIF ! (ok_histNMC(3)) THEN |
---|