source: lmdz_wrf/branches/LMDZ_WRFmeas/WRFV3/lmdz/write_histhfNMC.h @ 25

Last change on this file since 25 was 1, checked in by lfita, 10 years ago
  • -- --- Opening of the WRF+LMDZ coupling repository --- -- -

WRF: version v3.3
LMDZ: version v1818

More details in:

  • Property svn:executable set to *
File size: 2.8 KB
Line 
1!
2! $Header$
3!
4      IF (ok_histNMC(3)) THEN
5!c
6       ndex3d = 0
7       itau_w = itau_phy + itap + start_time * day_step / iphysiq
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!$OMP MASTER
86        call histsync(nid_hfnmc)
87!$OMP END MASTER
88      endif
89!c
90      ENDIF !      (ok_histNMC(3)) THEN
Note: See TracBrowser for help on using the repository browser.