source: lmdz_wrf/trunk/WRFV3/lmdz/write_histmthNMC.h @ 1465

Last change on this file since 1465 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:

File size: 3.6 KB
Line 
1!
2! $Id: write_histmthNMC.h 1665 2012-10-09 13:35:26Z fairhead $
3!
4      IF (ok_histNMC(1)) 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_mthnmc,lNMC,"tnondef",itau_w,                           &
17       &tnondef(:,:,1))
18!c
19      CALL histwrite_phy(nid_mthnmc,lNMC,"ta",itau_w,                                &
20       &twriteSTD(:,:,1))
21!c
22      CALL histwrite_phy(nid_mthnmc,lNMC,"zg",itau_w,                                &
23       &phiwriteSTD(:,:,1))
24!c
25       CALL histwrite_phy(nid_mthnmc,lNMC,"hus",itau_w,                              &
26       &qwriteSTD(:,:,1))
27!c
28       CALL histwrite_phy(nid_mthnmc,lNMC,"hur",itau_w,                              &
29       &rhwriteSTD(:,:,1))
30!c
31       CALL histwrite_phy(nid_mthnmc,lNMC,"ua",itau_w,                               &
32       &uwriteSTD(:,:,1))
33!c
34       CALL histwrite_phy(nid_mthnmc,lNMC,"va",itau_w,                               &
35       &vwriteSTD(:,:,1))
36!c
37       CALL histwrite_phy(nid_mthnmc,lNMC,"wap",itau_w,                              &
38       &wwriteSTD(:,:,1))
39!c
40       DO k=1, nlevSTD
41        DO i=1, klon
42         IF(tnondef(i,k,1).NE.missing_val) THEN
43          zx_tmp_fiNC(i,k) = (100.*tnondef(i,k,1))/freq_moyNMC(1)
44         ELSE
45          zx_tmp_fiNC(i,k) = missing_val
46         ENDIF
47        ENDDO
48       ENDDO !k=1, nlevSTD
49!c
50      CALL histwrite_phy(nid_mthnmc,lNMC,"psbg",itau_w,                              &
51       &zx_tmp_fiNC)
52!c
53      CALL histwrite_phy(nid_mthnmc,lNMC,"uv",itau_w,                                &
54       &uvsumSTD(:,:,1))
55!c
56      CALL histwrite_phy(nid_mthnmc,lNMC,"vq",itau_w,                                &
57       &vqsumSTD(:,:,1))
58!c
59      CALL histwrite_phy(nid_mthnmc,lNMC,"vT",itau_w,                                &
60       &vTsumSTD(:,:,1))
61!c
62      CALL histwrite_phy(nid_mthnmc,lNMC,"wq",itau_w,                                &
63       &wqsumSTD(:,:,1))
64!c
65      CALL histwrite_phy(nid_mthnmc,lNMC,"vphi",itau_w,                              &
66       &vphisumSTD(:,:,1))
67!c
68       CALL histwrite_phy(nid_mthnmc,lNMC,"wT",itau_w,                               &
69       &wTsumSTD(:,:,1))
70!c
71       CALL histwrite_phy(nid_mthnmc,lNMC,"uxu",itau_w,                              &
72       &u2sumSTD(:,:,1))
73!c
74       CALL histwrite_phy(nid_mthnmc,lNMC,"vxv",itau_w,                              &
75       &v2sumSTD(:,:,1))
76!c
77       CALL histwrite_phy(nid_mthnmc,lNMC,"TxT",itau_w,                              &
78       &T2sumSTD(:,:,1))
79!c
80       DO k=1, nlevSTD
81        DO i=1, klon
82         IF(O3sumSTD(i,k,1).NE.missing_val) THEN
83          zx_tmp_fiNC(i,k) = O3sumSTD(i,k,1) * 1.e+9
84         ELSE
85          zx_tmp_fiNC(i,k) = missing_val
86         ENDIF
87        ENDDO
88       ENDDO !k=1, nlevSTD
89       CALL histwrite_phy(nid_mthnmc,lNMC,"tro3",itau_w,                             &
90       & zx_tmp_fiNC)
91!c
92       if (read_climoz == 2) THEN
93       DO k=1, nlevSTD
94        DO i=1, klon
95         IF(O3daysumSTD(i,k,1).NE.missing_val) THEN
96          zx_tmp_fiNC(i,k) = O3daysumSTD(i,k,1) * 1.e+9
97         ELSE
98          zx_tmp_fiNC(i,k) = missing_val
99         ENDIF
100        ENDDO
101       ENDDO !k=1, nlevSTD
102!c
103        CALL histwrite_phy(nid_mthnmc,lNMC,"tro3_daylight",                          &
104       &itau_w, zx_tmp_fiNC)
105       endif
106!c
107      if (ok_sync) then
108!$OMP MASTER
109        call histsync(nid_mthnmc)
110!$OMP END MASTER
111      endif
112!c
113      ENDIF !(ok_histNMC(1)) THEN
Note: See TracBrowser for help on using the repository browser.