source: LMDZ4/branches/LMDZ4V5.0-dev/libf/phylmd/write_histrac.h @ 1373

Last change on this file since 1373 was 1279, checked in by Laurent Fairhead, 15 years ago

Merged LMDZ4-dev branch changes r1241:1278 into the trunk
Running trunk and LMDZ4-dev in LMDZOR configuration on local
machine (sequential) and SX8 (4-proc) yields identical results
(restart and restartphy are identical binarily)
Log history from r1241 to r1278 is available by switching to
source:LMDZ4/branches/LMDZ4-dev-20091210

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 2.6 KB
Line 
1!$Id $
2!***************************************
3!  ECRITURE DU FICHIER :  histrac.nc
4!***************************************
5  IF (ecrit_tra > 0. .AND. config_inca == 'none') THEN
6     
7     itau_w = itau_phy + nstep
8     
9     CALL histwrite_phy(nid_tra,"phis",itau_w,pphis)
10     CALL histwrite_phy(nid_tra,"aire",itau_w,airephy)
11
12!TRACEURS
13!----------------
14     DO it=1,nbtr
15        iiq=niadv(it+2)
16
17! CONCENTRATIONS
18        CALL histwrite_phy(nid_tra,tname(iiq),itau_w,tr_seri(:,:,it))
19
20! TD LESSIVAGE       
21        IF (lessivage .AND. aerosol(it)) THEN
22           CALL histwrite_phy(nid_tra,"fl"//tname(iiq),itau_w,flestottr(:,:,it))
23        ENDIF
24
25! TD THERMIQUES
26        IF (iflag_thermals.gt.0) THEN
27           CALL histwrite_phy(nid_tra,"d_tr_th_"//tname(iiq),itau_w,d_tr_th(:,:,it))
28        ENDIF
29
30! TD CONVECTION
31        IF (iflag_con.GE.2) THEN
32           CALL histwrite_phy(nid_tra,"d_tr_cv_"//tname(iiq),itau_w,d_tr_cv(:,:,it))
33        ENDIF
34
35! TD COUCHE-LIMITE
36        CALL histwrite_phy(nid_tra,"d_tr_cl_"//tname(iiq),itau_w,d_tr_cl(:,:,it))
37     ENDDO
38!---------------
39!
40!
41! VENT (niveau 1)   
42     CALL histwrite_phy(nid_tra,"pyu1",itau_w,yu1)
43     CALL histwrite_phy(nid_tra,"pyv1",itau_w,yv1)
44!
45! TEMPERATURE DU SOL
46     zx_tmp_fi2d(:)=ftsol(:,1)         
47     CALL histwrite_phy(nid_tra,"ftsol1",itau_w,zx_tmp_fi2d)
48     zx_tmp_fi2d(:)=ftsol(:,2)
49     CALL histwrite_phy(nid_tra,"ftsol2",itau_w,zx_tmp_fi2d)
50     zx_tmp_fi2d(:)=ftsol(:,3)
51     CALL histwrite_phy(nid_tra,"ftsol3",itau_w,zx_tmp_fi2d)
52     zx_tmp_fi2d(:)=ftsol(:,4)
53     CALL histwrite_phy(nid_tra,"ftsol4",itau_w,zx_tmp_fi2d)
54!     
55! NATURE DU SOL
56     zx_tmp_fi2d(:)=pctsrf(:,1)
57     CALL histwrite_phy(nid_tra,"psrf1",itau_w,zx_tmp_fi2d)
58     zx_tmp_fi2d(:)=pctsrf(:,2)
59     CALL histwrite_phy(nid_tra,"psrf2",itau_w,zx_tmp_fi2d)
60     zx_tmp_fi2d(:)=pctsrf(:,3)
61     CALL histwrite_phy(nid_tra,"psrf3",itau_w,zx_tmp_fi2d)
62     zx_tmp_fi2d(:)=pctsrf(:,4)
63     CALL histwrite_phy(nid_tra,"psrf4",itau_w,zx_tmp_fi2d)
64 
65! DIVERS   
66     CALL histwrite_phy(nid_tra,"pplay",itau_w,pplay)     
67     CALL histwrite_phy(nid_tra,"t",itau_w,t_seri)     
68     CALL histwrite_phy(nid_tra,"mfu",itau_w,pmfu)
69     CALL histwrite_phy(nid_tra,"mfd",itau_w,pmfd)
70     CALL histwrite_phy(nid_tra,"en_u",itau_w,pen_u)
71     CALL histwrite_phy(nid_tra,"en_d",itau_w,pen_d)
72     CALL histwrite_phy(nid_tra,"de_d",itau_w,pde_d)
73     CALL histwrite_phy(nid_tra,"de_u",itau_w,pde_u)
74     CALL histwrite_phy(nid_tra,"coefh",itau_w,coefh)
75
76     IF (ok_sync) THEN
77!$OMP MASTER
78        CALL histsync(nid_tra)
79!$OMP END MASTER
80     ENDIF
81
82  ENDIF !ecrit_tra>0. .AND. config_inca == 'none'
83
Note: See TracBrowser for help on using the repository browser.