Changeset 1279 for LMDZ4/trunk/libf/phylmd/ini_histrac.h
- Timestamp:
- Dec 10, 2009, 10:02:56 AM (15 years ago)
- Location:
- LMDZ4/trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk
- Property svn:mergeinfo changed
/LMDZ4/branches/LMDZ4-dev merged: 1150-1162,1164-1193,1195-1231,1234-1235,1237-1240,1242-1274,1276
- Property svn:mergeinfo changed
-
LMDZ4/trunk/libf/phylmd/ini_histrac.h
r1146 r1279 1 1 ! 2 ! $ Header$2 ! $Id $ 3 3 ! 4 IF (ecrit_tra>0. .AND. config_inca == 'none') THEN 5 c$OMP MASTER 6 CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian) 7 c 8 CALL histbeg_phy("histrac", itau_phy, zjulian, pdtphys, 9 . nhori, nid_tra) 10 CALL histvert(nid_tra, "presnivs", "Vertical levels", "mb", 11 . klev, presnivs, nvert) 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) 12 9 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) 13 16 17 !TRACEURS 18 !---------------- 19 DO it = 1,nbtr 20 iiq = niadv(it+2) 14 21 15 zsto = pdtphys 16 zout = ecrit_tra 17 c 18 CALL histdef(nid_tra, "phis", "Surface geop. height", "-", 19 . iim,jj_nb,nhori, 1,1,1, -99, 32, 20 . "once", zsto,zout) 21 c 22 CALL histdef(nid_tra, "aire", "Grid area", "-", 23 . iim,jj_nb,nhori, 1,1,1, -99, 32, 24 . "once", zsto,zout) 25 DO it=1,nbtr 26 C champ 2D 27 iq=it+2 28 iiq=niadv(iq) 29 CALL histdef(nid_tra, tname(iiq), ttext(iiq), "U/kga", 30 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 31 . "ave(X)", zsto,zout) 32 if (lessivage) THEN 33 CALL histdef(nid_tra, "fl"//tname(iiq),"Flux "//ttext(iiq), 34 . "U/m2/s",iim,jj_nb,nhori, klev,1,klev,nvert, 32, 35 . "ave(X)", zsto,zout) 36 endif 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) 37 25 38 c---Ajout Olivia 39 CALL histdef(nid_tra, "d_tr_th_"//tname(iiq), 40 . "tendance thermique"// ttext(iiq), "?", 41 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 42 . "ave(X)", zsto,zout) 43 c 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 !(iflag_con.GE.2) then 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 c---fin Olivia 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 55 32 56 ENDDO 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 57 40 58 CALL histdef(nid_tra, "pyu1", "Vent niv 1", "-", 59 . iim,jj_nb,nhori, 1,1,1, -99, 32, 60 . "inst(X)", zout,zout) 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 61 48 62 CALL histdef(nid_tra, "pyv1", "Vent niv 1", "-", 63 . iim,jj_nb,nhori, 1,1,1, -99, 32, 64 . "inst(X)", zout,zout) 65 CALL histdef(nid_tra, "psrf1", "nature sol", "-", 66 . iim,jj_nb,nhori, 1,1,1, -99, 32, 67 . "inst(X)", zout,zout) 68 CALL histdef(nid_tra, "psrf2", "nature sol", "-", 69 . iim,jj_nb,nhori, 1,1,1, -99, 32, 70 . "inst(X)", zout,zout) 71 CALL histdef(nid_tra, "psrf3", "nature sol", "-", 72 . iim,jj_nb,nhori, 1,1,1, -99, 32, 73 . "inst(X)", zout,zout) 74 CALL histdef(nid_tra, "psrf4", "nature sol", "-", 75 . iim,jj_nb,nhori, 1,1,1, -99, 32, 76 . "inst(X)", zout,zout) 77 CALL histdef(nid_tra, "ftsol1", "temper sol", "-", 78 . iim,jj_nb,nhori, 1,1,1, -99, 32, 79 . "inst(X)", zout,zout) 80 CALL histdef(nid_tra, "ftsol2", "temper sol", "-", 81 . iim,jj_nb,nhori, 1,1,1, -99, 32, 82 . "inst(X)", zout,zout) 83 CALL histdef(nid_tra, "ftsol3", "temper sol", "-", 84 . iim,jj_nb,nhori, 1,1,1, -99, 32, 85 . "inst", zout,zout) 86 CALL histdef(nid_tra, "ftsol4", "temper sol", "-", 87 . iim,jj_nb,nhori, 1,1,1, -99, 32, 88 . "inst(X)", zout,zout) 89 CALL histdef(nid_tra, "pplay", "flux u mont","-", 90 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 91 . "inst(X)", zout,zout) 92 CALL histdef(nid_tra, "t", "flux u mont","-", 93 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 94 . "inst(X)", zout,zout) 95 CALL histdef(nid_tra, "mfu", "flux u mont","-", 96 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 97 . "ave(X)", zsto,zout) 98 CALL histdef(nid_tra, "mfd", "flux u decen","-", 99 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 100 . "ave(X)", zsto,zout) 101 CALL histdef(nid_tra, "en_u", "flux u mont","-", 102 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 103 . "ave(X)", zsto,zout) 104 CALL histdef(nid_tra, "en_d", "flux u mont","-", 105 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 106 . "ave(X)", zsto,zout) 107 CALL histdef(nid_tra, "de_d", "flux u mont","-", 108 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 109 . "ave(X)", zsto,zout) 110 CALL histdef(nid_tra, "de_u", "flux u decen","-", 111 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 112 . "ave(X)", zsto,zout) 113 CALL histdef(nid_tra, "coefh", "turbulent coef","-", 114 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 115 . "ave(X)", zsto,zout) 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) 116 64 117 c 118 CALL histend(nid_tra) 119 ndex2d = 0 120 ndex3d = 0 121 ndex = 0 122 c$OMP END MASTER 123 END IF ! ecrit_tra>0. .AND. config_inca == 'none' 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
Note: See TracChangeset
for help on using the changeset viewer.