- Timestamp:
- Feb 7, 2005, 4:47:11 PM (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/IPSL-CM4_IPCC_patches/libf/phylmd/write_histhf.h
r524 r589 98 98 zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter) 99 99 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 100 CALL histwrite(nid_hf,"SWnetOR",itau_w, 101 $ zx_tmp_2d,iim*jjmp1,ndex2d) 102 c 103 zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol(1:klon)) 104 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d) 100 105 CALL histwrite(nid_hf,"SWdownOR",itau_w, 101 106 $ zx_tmp_2d,iim*jjmp1,ndex2d) … … 166 171 167 172 endif 173 if (ok_hf) then 174 175 c Comprendre comment marche el i=nint(zout/zsto) 176 c 177 c print*,'ACRITURE HF !!! ACRITURE HF !!! ACRITURE HF !!! ' 178 ndex2d = 0 179 ndex3d = 0 180 c 181 zsto = dtime 182 zout = dtime * ecrit_hf 183 itau_w = itau_phy + itap 184 c 185 IF(lev_histhf.GE.1) THEN 186 c 187 c i = NINT(zout/zsto) 188 c CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) 189 c CALL histwrite(nid_hf,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d) 190 c 191 c i = NINT(zout/zsto) 192 c CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) 193 c CALL histwrite(nid_hf,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d) 194 C 195 CALL gr_fi_ecrit(1, klon,iim,jjmp1, paire_ter, zx_tmp_2d) 196 CALL histwrite(nid_hf,"aireTER",itau_w, 197 $ zx_tmp_2d,iim*jjmp1,ndex2d) 198 c 199 DO i=1, klon 200 zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic) 201 ENDDO 202 c 203 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d) 204 CALL histwrite(nid_hf,"contfracATM",itau_w, 205 $ zx_tmp_2d,iim*jjmp1,ndex2d) 206 c 207 CALL gr_fi_ecrit(1,klon,iim,jjmp1,pctsrf_new(:,is_ter),zx_tmp_2d) 208 CALL histwrite(nid_hf,"contfracOR",itau_w, 209 $ zx_tmp_2d,iim*jjmp1,ndex2d) 210 c 211 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d) 212 CALL histwrite(nid_hf,"t2m",itau_w,zx_tmp_2d,iim*jjmp1, 213 . ndex2d) 214 c 215 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d) 216 CALL histwrite(nid_hf,"q2m",itau_w,zx_tmp_2d,iim*jjmp1, 217 . ndex2d) 218 c 219 DO i = 1, klon 220 zx_tmp_fi2d(i) = paprs(i,1) 221 ENDDO 222 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 223 CALL histwrite(nid_hf,"psol",itau_w,zx_tmp_2d,iim*jjmp1, 224 . ndex2d) 225 c 226 DO i = 1, klon 227 zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) 228 ENDDO 229 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 230 CALL histwrite(nid_hf,"rain",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 231 c 232 c ENSEMBLES BEG 233 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 234 CALL histwrite(nid_hf,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 235 c 236 CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d) 237 CALL histwrite(nid_hf,"slp",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 238 c 239 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d) 240 CALL histwrite(nid_hf,"u10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 241 c 242 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d) 243 CALL histwrite(nid_hf,"v10m",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 244 c 245 DO i=1, klon 246 zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i)) 247 ENDDO 248 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 249 CALL histwrite(nid_hf,"wind10m",itau_w,zx_tmp_2d, 250 $ iim*jjmp1,ndex2d) 251 c 252 DO k=1, nlevENS 253 IF(clev(k).EQ."500") THEN 254 CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d) 255 CALL histwrite(nid_hf,"phi"//clev(k),itau_w,zx_tmp_2d, 256 $ iim*jjmp1,ndex2d) 257 ENDIF !clev(k).EQ."500" 258 ENDDO 259 c 260 ENDIF !lev_histhf.GE.1 261 c 262 IF(lev_histhf.GE.2) THEN 263 c 264 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d) 265 CALL histwrite(nid_hf,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 266 c 267 c -- LOOP 268 zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter) 269 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 270 CALL histwrite(nid_hf,"SWnetOR",itau_w, 271 $ zx_tmp_2d,iim*jjmp1,ndex2d) 272 c 273 CALL gr_fi_ecrit(1, klon,iim,jjmp1, swdownor , zx_tmp_2d) 274 CALL histwrite(nid_hf,"SWdownOR",itau_w, 275 $ zx_tmp_2d,iim*jjmp1,ndex2d) 276 c -- LOOP 277 c 278 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d) 279 CALL histwrite(nid_hf,"LWdownOR",itau_w,zx_tmp_2d,iim*jjmp1, 280 $ ndex2d) 281 c 282 c 283 ENDIF !lev_histhf.GE.2 284 c 285 IF(lev_histhf.GE.3) THEN 286 c 287 DO k=1, nlevENS 288 c 289 CALL gr_fi_ecrit(1, klon,iim,jjmp1, tlev(:,k),zx_tmp_2d) 290 CALL histwrite(nid_hf,"t"//clev(k),itau_w,zx_tmp_2d, 291 $ iim*jjmp1,ndex2d) 292 c 293 IF(clev(k).NE."500") THEN !clev(k).NE."500" 294 CALL gr_fi_ecrit(1, klon,iim,jjmp1, philev(:,k),zx_tmp_2d) 295 CALL histwrite(nid_hf,"phi"//clev(k),itau_w,zx_tmp_2d, 296 $ iim*jjmp1,ndex2d) 297 ENDIF !clev(k).NE."500" 298 c 299 CALL gr_fi_ecrit(1, klon,iim,jjmp1, qlev(:,k),zx_tmp_2d) 300 CALL histwrite(nid_hf,"q"//clev(k),itau_w,zx_tmp_2d, 301 $ iim*jjmp1,ndex2d) 302 c 303 IF(1.EQ.0) THEN 304 CALL gr_fi_ecrit(1, klon,iim,jjmp1, rhlev(:,k),zx_tmp_2d) 305 CALL histwrite(nid_hf,"rh"//clev(k),itau_w,zx_tmp_2d, 306 $ iim*jjmp1,ndex2d) 307 ENDIF !1.EQ.0 308 c 309 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ulev(:,k),zx_tmp_2d) 310 CALL histwrite(nid_hf,"u"//clev(k),itau_w,zx_tmp_2d, 311 $ iim*jjmp1,ndex2d) 312 c 313 CALL gr_fi_ecrit(1, klon,iim,jjmp1, vlev(:,k),zx_tmp_2d) 314 CALL histwrite(nid_hf,"v"//clev(k),itau_w,zx_tmp_2d, 315 $ iim*jjmp1,ndex2d) 316 c 317 ENDDO !nlevENS 318 c 319 IF(1.EQ.0) THEN 320 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d) 321 CALL histwrite(nid_hf,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 322 c 323 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d) 324 CALL histwrite(nid_hf,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 325 ENDIF !(1.EQ.0) THEN 326 c 327 ENDIF !lev_histhf.GE.3 328 c 329 IF(lev_histhf.GE.4) THEN 330 c 331 #define histhf3d 332 #ifdef histhf3d 333 #include "write_histhf3d.h" 334 #endif 335 c 336 ENDIF !lev_histhf.GE.4 337 c 338 if (ok_sync) then 339 call histsync(nid_hf) 340 endif 341 342 endif
Note: See TracChangeset
for help on using the changeset viewer.