source: LMDZ4/branches/LMDZ4-dev/libf/phylmd/ini_histrac.h @ 1191

Last change on this file since 1191 was 1191, checked in by jghattas, 15 years ago

Reecriture de phytrac et les routines concernes (Anthony Jamelot)

  • les suffix change de F -> F90 (nflxtr.F90,cltracrn.F90,initrrnpb.F90,cvltr.F90,minmaxqfi.F90,cltrac.F90,phytrac.F90)

Traitement d'un nouveau traceur berelium (optionel, toujours pour des
tests)(Anthony Jamelot)

  • radiornpb.F change du nom pour radio_decay.F90 car il traite maintenant tout les traceurs radioactives
  • ajoute init_be.F90

Nouveau interface dans phytrac pour serparer les calculs et appels
specifique a INCA avec les traitements des traceurs specifiques au LMDZ
(JG)

  • ajoute tracinca_mod.F90 pour les appeles a INCA
  • ajoute traclmdz_mod.F90 pour les calculs des traceurs specifiques a LMDZ
  • enleve fichier restartrac et ajoute la variable trs dans restartphy.nc

La convergence numerique a etait rompue uniquement pour les traceurs
LMDZ RN et PB.
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 5.0 KB
Line 
1!
2! $Id $
3!
4#ifdef CPP_IOIPSL
5
6  IF (ecrit_tra>0. .AND. config_inca == 'none') THEN
7!$OMP MASTER
8     CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)
9     CALL histbeg_phy("histrac", itau_phy, zjulian, pdtphys,nhori, nid_tra)
10     CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb",klev, presnivs, nvert)
11
12     zsto = pdtphys
13     zout = ecrit_tra
14     CALL histdef(nid_tra, "phis", "Surface geop. height", "-",   &
15          iim,jj_nb,nhori, 1,1,1, -99, 32,"once",  zsto,zout)
16     CALL histdef(nid_tra, "aire", "Grid area", "-",              &
17          iim,jj_nb,nhori, 1,1,1, -99, 32,"once",  zsto,zout)
18
19!TRACEURS
20!----------------
21     DO it = 1,nbtr
22        iiq = niadv(it+2)
23
24! CONCENTRATIONS
25        CALL histdef(nid_tra, tname(iiq), ttext(iiq), "U/kga",    &
26             iim,jj_nb,nhori, klev,1,klev,nvert, 32,"ave(X)", zsto,zout)
27
28! TD LESSIVAGE
29        IF (lessivage .AND. aerosol(it)) THEN
30           CALL histdef(nid_tra, "fl"//tname(iiq),"Flux "//ttext(iiq), &
31                "U/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
32                "ave(X)", zsto,zout)
33        END IF
34
35! TD THERMIQUES
36        IF (iflag_thermals.gt.0) THEN
37           CALL histdef(nid_tra, "d_tr_th_"//tname(iiq),      &
38                "tendance thermique"// ttext(iiq), "?",       &
39                iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
40                "ave(X)", zsto,zout)
41        ENDIF
42
43! TD CONVECTION
44        IF (iflag_con.GE.2) THEN
45           CALL histdef(nid_tra, "d_tr_cv_"//tname(iiq),   &
46                "tendance convection"// ttext(iiq), "?",   &
47                iim,jj_nb,nhori, klev,1,klev,nvert, 32,    &
48                "ave(X)", zsto,zout)
49        ENDIF
50
51! TD COUCHE-LIMITE
52        CALL histdef(nid_tra, "d_tr_cl_"//tname(iiq),      &
53             "tendance couche limite"// ttext(iiq), "?",   &
54             iim,jj_nb,nhori, klev,1,klev,nvert, 32,       &
55             "ave(X)", zsto,zout)
56     ENDDO
57!---------------   
58!
59! VENT (niveau 1)
60     CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-",      &
61          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
62          "inst(X)",  zout,zout)     
63     CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-",      &
64          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
65          "inst(X)",  zout,zout)
66
67! TEMPERATURE DU SOL
68     CALL histdef(nid_tra, "ftsol1", "temper sol", "-",    &
69          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
70          "inst(X)",  zout,zout)
71     CALL histdef(nid_tra, "ftsol2", "temper sol", "-",    &
72          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
73          "inst(X)",  zout,zout)
74     CALL histdef(nid_tra, "ftsol3", "temper sol", "-",    &
75          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
76          "inst",  zout,zout)
77     CALL histdef(nid_tra, "ftsol4", "temper sol", "-",    &
78          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
79          "inst(X)",  zout,zout)
80
81! NATURE DU SOL
82     CALL histdef(nid_tra, "psrf1", "nature sol", "-",     &
83          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
84          "inst(X)",  zout,zout)
85     CALL histdef(nid_tra, "psrf2", "nature sol", "-",     &
86          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
87          "inst(X)",  zout,zout)
88     CALL histdef(nid_tra, "psrf3", "nature sol", "-",     &
89          iim,jj_nb,nhori, 1,1,1, -99, 32,                 &
90          "inst(X)",  zout,zout)
91     CALL histdef(nid_tra, "psrf4", "nature sol", "-",     &
92          iim,jj_nb,nhori, 1,1,1, -99, 32,                 & 
93          "inst(X)",  zout,zout)
94! DIVERS
95     CALL histdef(nid_tra, "pplay", "flux u mont","-",     &
96          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
97          "inst(X)", zout,zout)
98     CALL histdef(nid_tra, "t", "flux u mont","-",         &
99          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
100          "inst(X)", zout,zout)
101     CALL histdef(nid_tra, "mfu", "flux u mont","-",       &
102          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
103          "ave(X)", zsto,zout)
104     CALL histdef(nid_tra, "mfd", "flux u decen","-",      &
105          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
106          "ave(X)", zsto,zout)
107     CALL histdef(nid_tra, "en_u", "flux u mont","-",      &
108          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
109          "ave(X)", zsto,zout)
110     CALL histdef(nid_tra, "en_d", "flux u mont","-",      &
111          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
112          "ave(X)", zsto,zout)
113     CALL histdef(nid_tra, "de_d", "flux u mont","-",      &
114          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
115          "ave(X)", zsto,zout)
116     CALL histdef(nid_tra, "de_u", "flux u decen","-",     &
117          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
118          "ave(X)", zsto,zout)
119     CALL histdef(nid_tra, "coefh", "turbulent coef","-",  &
120          iim,jj_nb,nhori, klev,1,klev,nvert, 32,          &
121          "ave(X)", zsto,zout)   
122     
123     CALL histend(nid_tra)
124!$OMP END MASTER
125  END IF ! ecrit_tra>0. .AND. config_inca == 'none'
126
127#endif
128 
Note: See TracBrowser for help on using the repository browser.