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