Changeset 902 for trunk/LMDZ.VENUS/libf/phyvenus/ini_histday.h
- Timestamp:
- Mar 13, 2013, 11:06:21 AM (12 years ago)
- File:
-
- 1 edited
-
trunk/LMDZ.VENUS/libf/phyvenus/ini_histday.h (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.VENUS/libf/phyvenus/ini_histday.h
r892 r902 6 6 zsto = dtime 7 7 zout = dtime * REAL(ecrit_day) 8 zsto1= dtime * REAL(ecrit_day)9 8 10 9 idayref = day_ref … … 12 11 13 12 call histbeg_phy("histday.nc",itau_phy, 14 . zjulian, zsto,nhori,nid_day)13 . zjulian,dtime,nhori,nid_day) 15 14 16 15 !$OMP MASTER … … 39 38 . "ave(X)", zsto,zout) 40 39 c 41 CALL histdef(nid_day, "ue", "Zonal energy transport", "-", 42 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 43 . "ave(X)", zsto,zout) 44 c 45 CALL histdef(nid_day, "ve", "Merid energy transport", "-", 46 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 47 . "ave(X)", zsto,zout) 40 c CALL histdef(nid_day, "ue", "Zonal energy transport", "-", 41 c . iim,jj_nb,nhori, 1,1,1, nvert, 32, 42 c . "ave(X)", zsto,zout) 43 c 44 c CALL histdef(nid_day, "ve", "Merid energy transport", "-", 45 c . iim,jj_nb,nhori, 1,1,1, nvert, 32, 46 c . "ave(X)", zsto,zout) 47 c 48 c CALL histdef(nid_day, "cdragh", "Drag coef on T", "-", 49 c . iim,jj_nb,nhori, 1,1,1, nvert, 32, 50 c . "ave(X)", zsto,zout) 51 c 52 c CALL histdef(nid_day, "cdragm", "Drag coef on U", "-", 53 c . iim,jj_nb,nhori, 1,1,1, nvert, 32, 54 c . "ave(X)", zsto,zout) 48 55 c 49 56 ENDIF !lev_histday.GE.1 … … 78 85 . "ave(X)", zsto,zout) 79 86 c 87 CALL histdef(nid_day, "dudyn", "Dynamics dU", "m/s2", 88 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 89 . "ave(X)", zsto,zout) 90 c 91 CALL histdef(nid_day, "duvdf", "Boundary-layer dU", "m/s2", 92 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 93 . "ave(X)", zsto,zout) 94 c 95 c CALL histdef(nid_day, "mang", "Angular momentum", "kg m2/s", 96 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 97 c . "ave(X)", zsto,zout) 98 c 99 c CALL histdef(nid_day, "Kz", "vertical diffusion coef", "m2/s", 100 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 101 c . "ave(X)", zsto,zout) 102 c 103 c plusieurs traceurs 104 if (iflag_trac.eq.1) THEN 105 DO iq=1,nqmax 106 IF (iq.LE.99) THEN 107 WRITE(str2,'(i2.2)') iq 108 CALL histdef(nid_day, tname(iq), ttext(iq), "ppm", 109 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 110 . "ave(X)", zsto,zout) 111 ELSE 112 PRINT*, "Trop de traceurs" 113 CALL abort 114 ENDIF 115 ENDDO 116 endif 117 c 118 CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2", 119 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 120 . "ave(X)", zsto,zout) 121 c 80 122 ENDIF !lev_histday.GE.2 81 123 c … … 87 129 c 2D 88 130 c 89 CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2",90 . iim,jj_nb,nhori, 1,1,1, nvert, 32,91 . "ave(X)", zsto1,zout)92 c93 131 CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2", 94 132 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 95 . "ave(X)", zsto 1,zout)133 . "ave(X)", zsto,zout) 96 134 c 97 135 CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2", 98 136 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 99 . "ave(X)", zsto 1,zout)137 . "ave(X)", zsto,zout) 100 138 c 101 139 CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2", 102 140 . iim,jj_nb,nhori, 1,1,1, nvert, 32, 103 . "ave(X)", zsto 1,zout)141 . "ave(X)", zsto,zout) 104 142 c 105 143 c 3D … … 107 145 CALL histdef(nid_day, "SWnet", "Net SW flux","W/m2", 108 146 . iim,jj_nb,nhori, klev,1,klev,nvert, 109 . 32, "ave(X)", zsto 1,zout)147 . 32, "ave(X)", zsto,zout) 110 148 c 111 149 CALL histdef(nid_day, "LWnet", "Net LW flux","W/m2", 112 150 . iim,jj_nb,nhori, klev,1,klev,nvert, 113 . 32, "ave(X)", zsto1,zout) 151 . 32, "ave(X)", zsto,zout) 152 c 153 CALL histdef(nid_day, "fluxvdf", "PBL net flux","W/m2", 154 . iim,jj_nb,nhori, klev,1,klev,nvert, 155 . 32, "ave(X)", zsto,zout) 156 c 157 CALL histdef(nid_day, "fluxdyn", "Dyn. net flux","W/m2", 158 . iim,jj_nb,nhori, klev,1,klev,nvert, 159 . 32, "ave(X)", zsto,zout) 160 c 161 CALL histdef(nid_day, "fluxajs", "Dry adj. net flux","W/m2", 162 . iim,jj_nb,nhori, klev,1,klev,nvert, 163 . 32, "ave(X)", zsto,zout) 164 c 165 c CALL histdef(nid_day, "fluxec", "Cin. net flux","W/m2", 166 c . iim,jj_nb,nhori, klev,1,klev,nvert, 167 c . 32, "ave(X)", zsto,zout) 114 168 c 115 169 ENDIF !lev_histday.GE.3 … … 122 176 . "ave(X)", zsto,zout) 123 177 c 124 c 125 c 126 c 178 c CALL histdef(nid_day, "dtphy", "Physics dT", "K/s", 179 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 180 c . "ave(X)", zsto,zout) 127 181 c 128 182 CALL histdef(nid_day, "dtvdf", "Boundary-layer dT", "K/s", … … 142 196 . "ave(X)", zsto,zout) 143 197 c 144 c CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s", 145 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 146 c . "ave(X)", zsto,zout) 147 c 148 CALL histdef(nid_mth, "dudyn", "Dynamics dU", "m/s2", 149 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 150 . "ave(X)", zsto,zout) 151 c 152 CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2", 153 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 154 . "ave(X)", zsto,zout) 198 c CALL histdef(nid_day, "dtec", "Cinetic dissip dT", "K/s", 199 c . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 200 c . "ave(X)", zsto,zout) 155 201 c 156 202 CALL histdef(nid_mth, "duajs", "Dry convection dU", "m/s2", … … 174 220 c------------------------------------------------------- 175 221 IF(lev_histday.GE.5) THEN 176 c177 if (nqmax.GE.3) THEN178 DO iq=1,nqmax-2179 IF (iq.LE.99) THEN180 WRITE(str2,'(i2.2)') iq181 CALL histdef(nid_day, "trac"//str2, "Tracer No."//str2, "-",182 . iim,jj_nb,nhori, klev,1,klev,nvert, 32,183 . "ave(X)", zsto,zout)184 ELSE185 PRINT*, "Trop de traceurs"186 CALL abort187 ENDIF188 ENDDO189 endif190 222 c 191 223 c call histdef(nid_day, "taux",
Note: See TracChangeset
for help on using the changeset viewer.
