Changeset 353 for LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F
- Timestamp:
- Mar 21, 2002, 4:09:32 PM (22 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F
r340 r353 360 360 c 361 361 INTEGER julien 362 INTEGER idayvrai363 SAVE idayvrai364 362 c 365 363 INTEGER lmt_pas … … 555 553 SAVE ecrit_reg ! frequence d'ecriture 556 554 c 555 integer itau_w ! pas de temps ecriture = itap + itau_phy 557 556 c 558 557 c … … 589 588 logical ok_sync 590 589 real date0 590 integer idayref 591 591 592 592 C essai writephys … … 761 761 c Gestion calendrier 762 762 763 CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)764 zjulian = zjulian + day_ini765 766 763 c 767 764 IF (ok_journe) THEN 768 765 c 769 CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)770 zjulian = zjulian + day_ini766 idayref = day_ref 767 CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) 771 768 c 772 769 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) … … 779 776 ENDDO 780 777 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 778 write(*,*)'zx_lon = ',zx_lon(:,1) 779 write(*,*)'zx_lat = ',zx_lat(1,:) 781 780 CALL histbeg("histday", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 782 . 1,iim,1,jjmp1, 0, zjulian, dtime,781 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 783 782 . nhori, nid_day) 783 write(*,*)'Journee ', itau_phy, zjulian 784 784 CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", 785 785 . klev, presnivs, nvert) … … 1002 1002 IF (ok_mensuel) THEN 1003 1003 c 1004 CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)1005 zjulian = zjulian + day_ini1004 idayref = day_ref 1005 CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) 1006 1006 c 1007 1007 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) … … 1015 1015 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 1016 1016 CALL histbeg("histmth.nc", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 1017 . 1,iim,1,jjmp1, 0, zjulian, dtime,1017 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 1018 1018 . nhori, nid_mth) 1019 write(*,*)'Mensuel ', itau_phy, zjulian 1019 1020 CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb", 1020 1021 . klev, presnivs, nvert) … … 1446 1447 IF (ok_instan) THEN 1447 1448 c 1448 CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)1449 zjulian = zjulian + day_ini1449 idayref = day_ref 1450 CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) 1450 1451 c 1451 1452 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) … … 1458 1459 ENDDO 1459 1460 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 1460 CALL histbeg("histins", iim,zx_lon , jjmp1,zx_lat,1461 . 1,iim,1,jjmp1, 0, zjulian, dtime,1461 CALL histbeg("histins", iim,zx_lon(:,1), jjmp1,zx_lat(1,:), 1462 . 1,iim,1,jjmp1, itau_phy, zjulian, dtime, 1462 1463 . nhori, nid_ins) 1464 write(*,*)'Inst ', itau_phy, zjulian 1463 1465 CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", 1464 1466 . klev, presnivs, nvert) … … 1763 1765 c 1764 1766 IF (MOD(itap-1,lmt_pas) .EQ. 0) THEN 1765 idayvrai = NINT(xjour) 1766 PRINT *,' PHYS cond julien ',julien,idayvrai 1767 PRINT *,' PHYS cond julien ',julien 1767 1768 CALL ozonecm( FLOAT(julien), rlat, paprs, wo) 1768 1769 ENDIF … … 2398 2399 zsto = dtime 2399 2400 zout = dtime * FLOAT(ecrit_day) 2401 itau_w = itau_phy + itap 2400 2402 2401 2403 i = NINT(zout/zsto) … … 2418 2420 C 2419 2421 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 2420 CALL histwrite(nid_day,"tsol",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2422 CALL histwrite(nid_day,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2421 2423 c call writephy(fid_day,prof2d_av,'tsol',zxtsol, 2422 2424 c . 'Surface Temperature','K') … … 2425 2427 zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_ter) 2426 2428 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d) 2427 CALL histwrite(nid_day,"tter",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2429 CALL histwrite(nid_day,"tter",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2428 2430 c call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter), 2429 2431 c . 'Surface Temperature','K') … … 2431 2433 zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic) 2432 2434 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 2433 CALL histwrite(nid_day,"tlic",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2435 CALL histwrite(nid_day,"tlic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2434 2436 c call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic), 2435 2437 c . 'Surface Temperature','K') … … 2437 2439 zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce) 2438 2440 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 2439 CALL histwrite(nid_day,"toce",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2441 CALL histwrite(nid_day,"toce",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2440 2442 c call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce), 2441 2443 c . 'Surface Temperature','K') … … 2443 2445 zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic) 2444 2446 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 2445 CALL histwrite(nid_day,"tsic",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2447 CALL histwrite(nid_day,"tsic",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2446 2448 c call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic), 2447 2449 c . 'Surface Temperature','K') … … 2451 2453 ENDDO 2452 2454 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2453 CALL histwrite(nid_day,"psol",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2455 CALL histwrite(nid_day,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2454 2456 c Essai writephys 2455 2457 varname = 'psol' … … 2463 2465 ENDDO 2464 2466 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2465 CALL histwrite(nid_day,"rain",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2467 CALL histwrite(nid_day,"rain",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2466 2468 c call writephy(fid_day,prof2d_av,'rain',zx_tmp_fi2d, 2467 2469 c . 'Precipitation','mm/day') … … 2470 2472 c 2471 2473 CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) 2472 CALL histwrite(nid_day,"snow",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2474 CALL histwrite(nid_day,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2473 2475 c call writephy(fid_day,prof2d_av,'snow',snow_fall, 2474 2476 c . 'Snow','mm/day') 2475 2477 c 2476 2478 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d) 2477 CALL histwrite(nid_day,"snow_cov",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2479 CALL histwrite(nid_day,"snow_cov",itau_w,zx_tmp_2d,iim*jjmp1, 2480 . ndex2d) 2478 2481 c call writephy(fid_day,prof2d_av,'snow_cov',zxsnow, 2479 2482 c . 'Snow cover','mm') 2480 2483 c 2481 2484 CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) 2482 CALL histwrite(nid_day,"evap",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2485 CALL histwrite(nid_day,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2483 2486 c call writephy(fid_day,prof2d_av,'evap',evap, 2484 2487 c . 'Evaporation','mm/day') 2485 2488 c 2486 2489 CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) 2487 CALL histwrite(nid_day,"tops",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2490 CALL histwrite(nid_day,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2488 2491 c call writephy(fid_day,prof2d_av,'tops',topsw, 2489 2492 c . 'Solar rad. at TOA','W/m2') 2490 2493 c 2491 2494 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) 2492 CALL histwrite(nid_day,"topl",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2495 CALL histwrite(nid_day,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2493 2496 c call writephy(fid_day,prof2d_av,'topl',toplw, 2494 2497 c . 'IR rad. at TOA','W/m2') 2495 2498 c 2496 2499 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) 2497 CALL histwrite(nid_day,"sols",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2500 CALL histwrite(nid_day,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2498 2501 c call writephy(fid_day,prof2d_av,'sols',solsw, 2499 2502 c . 'Solar rad. at surf.','W/m2') 2500 2503 c 2501 2504 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) 2502 CALL histwrite(nid_day,"soll",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2505 CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2503 2506 c call writephy(fid_day,prof2d_av,'soll',sollw, 2504 2507 c . 'IR rad. at surface','W/m2') 2505 2508 c 2506 2509 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d) 2507 CALL histwrite(nid_day,"solldown",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2510 CALL histwrite(nid_day,"solldown",itau_w,zx_tmp_2d,iim*jjmp1, 2511 . ndex2d) 2508 2512 c call writephy(fid_day,prof2d_av,'solldown',sollwdown, 2509 2513 c . 'Down. IR rad. at surface','W/m2') 2510 2514 c 2511 2515 CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) 2512 CALL histwrite(nid_day,"bils",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2516 CALL histwrite(nid_day,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2513 2517 c call writephy(fid_day,prof2d_av,'bils',bils, 2514 2518 c . 'Surf. total heat flux','W/m2') 2515 2519 c 2516 2520 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) 2517 CALL histwrite(nid_day,"sens",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2521 CALL histwrite(nid_day,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2518 2522 c call writephy(fid_day,prof2d_av,'sens',sens, 2519 2523 c . 'Sensible heat flux','W/m2') 2520 2524 c 2521 2525 CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) 2522 CALL histwrite(nid_day,"fder",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2526 CALL histwrite(nid_day,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2523 2527 c call writephy(fid_day,prof2d_av,'fder',fder, 2524 2528 c . 'Heat flux derivation','W/m2') … … 2529 2533 zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf) 2530 2534 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2531 CALL histwrite(nid_day,"pourc_"//clnsurf(nsrf),ita p,2535 CALL histwrite(nid_day,"pourc_"//clnsurf(nsrf),itau_w, 2532 2536 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2533 2537 c call writephy(fid_day,prof2d_av,'pourc_'//clnsurf(nsrf), … … 2537 2541 zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf) 2538 2542 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2539 CALL histwrite(nid_day,"tsol_"//clnsurf(nsrf),ita p,2543 CALL histwrite(nid_day,"tsol_"//clnsurf(nsrf),itau_w, 2540 2544 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2541 2545 c call writephy(fid_day,prof2d_av,'tsol_'//clnsurf(nsrf), … … 2545 2549 zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf) 2546 2550 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2547 CALL histwrite(nid_day,"sens_"//clnsurf(nsrf),ita p,2551 CALL histwrite(nid_day,"sens_"//clnsurf(nsrf),itau_w, 2548 2552 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2549 2553 c call writephy(fid_day,prof2d_av,'sens_'//clnsurf(nsrf), … … 2553 2557 zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf) 2554 2558 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2555 CALL histwrite(nid_day,"lat_"//clnsurf(nsrf),ita p,2559 CALL histwrite(nid_day,"lat_"//clnsurf(nsrf),itau_w, 2556 2560 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2557 2561 c call writephy(fid_day,prof2d_av,'lat_'//clnsurf(nsrf), … … 2561 2565 zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf) 2562 2566 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2563 CALL histwrite(nid_day,"taux_"//clnsurf(nsrf),ita p,2567 CALL histwrite(nid_day,"taux_"//clnsurf(nsrf),itau_w, 2564 2568 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2565 2569 c call writephy(fid_day,prof2d_av,'taux_'//clnsurf(nsrf), … … 2569 2573 zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf) 2570 2574 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2571 CALL histwrite(nid_day,"tauy_"//clnsurf(nsrf),ita p,2575 CALL histwrite(nid_day,"tauy_"//clnsurf(nsrf),itau_w, 2572 2576 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2573 2577 c call writephy(fid_day,prof2d_av,'tauy_'//clnsurf(nsrf), … … 2577 2581 zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf) 2578 2582 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2579 CALL histwrite(nid_day,"albe_"//clnsurf(nsrf),ita p,2583 CALL histwrite(nid_day,"albe_"//clnsurf(nsrf),itau_w, 2580 2584 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2581 2585 c call writephy(fid_day,prof2d_av,'albe_'//clnsurf(nsrf), … … 2588 2592 zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf) 2589 2593 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2590 CALL histwrite(nid_day,"rugs_"//clnsurf(nsrf),ita p,2594 CALL histwrite(nid_day,"rugs_"//clnsurf(nsrf),itau_w, 2591 2595 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2592 2596 c call writephy(fid_day,prof2d_av,'rugs_'//clnsurf(nsrf), … … 2600 2604 c$$$ ENDDO 2601 2605 c$$$ CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2602 c$$$ CALL histwrite(nid_day,"sicf",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2606 c$$$ CALL histwrite(nid_day,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2603 2607 c 2604 2608 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) 2605 CALL histwrite(nid_day,"cldl",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2609 CALL histwrite(nid_day,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2606 2610 c call writephy(fid_day,prof2d_av,'cldl',cldl, 2607 2611 c . 'Low-level cloudiness','-') 2608 2612 c 2609 2613 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) 2610 CALL histwrite(nid_day,"cldm",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2614 CALL histwrite(nid_day,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2611 2615 c call writephy(fid_day,prof2d_av,'cldm',cldm, 2612 2616 c . 'Mid-level cloudiness','-') 2613 2617 c 2614 2618 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) 2615 CALL histwrite(nid_day,"cldh",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2619 CALL histwrite(nid_day,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2616 2620 c call writephy(fid_day,prof2d_av,'cldh',cldh, 2617 2621 c . 'High-level cloudiness','-') 2618 2622 c 2619 2623 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d) 2620 CALL histwrite(nid_day,"cldt",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2624 CALL histwrite(nid_day,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2621 2625 c call writephy(fid_day,prof2d_av,'cldt',cldt, 2622 2626 c . 'Total cloudiness','-') 2623 2627 c 2624 2628 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) 2625 CALL histwrite(nid_day,"cldq",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2629 CALL histwrite(nid_day,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2626 2630 c call writephy(fid_day,prof2d_av,'cldq',cldq, 2627 2631 c . 'Cloud liquid water path','-') … … 2630 2634 c 2631 2635 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) 2632 CALL histwrite(nid_day,"temp",ita p,zx_tmp_3d,2636 CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d, 2633 2637 . iim*jjmp1*klev,ndex3d) 2634 2638 c Essai writephys … … 2639 2643 c 2640 2644 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) 2641 CALL histwrite(nid_day,"ovap",ita p,zx_tmp_3d,2645 CALL histwrite(nid_day,"ovap",itau_w,zx_tmp_3d, 2642 2646 . iim*jjmp1*klev,ndex3d) 2643 2647 c call writephy(fid_day,prof3d_av,'ovap',qx(1,1,ivap), … … 2645 2649 c 2646 2650 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) 2647 CALL histwrite(nid_day,"geop",ita p,zx_tmp_3d,2651 CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d, 2648 2652 . iim*jjmp1*klev,ndex3d) 2649 2653 c call writephy(fid_day,prof3d_av,'geop',zphi, … … 2651 2655 c 2652 2656 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 2653 CALL histwrite(nid_day,"vitu",ita p,zx_tmp_3d,2657 CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d, 2654 2658 . iim*jjmp1*klev,ndex3d) 2655 2659 c call writephy(fid_day,prof3d_av,'vitu',u_seri, … … 2657 2661 c 2658 2662 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 2659 CALL histwrite(nid_day,"vitv",ita p,zx_tmp_3d,2663 CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d, 2660 2664 . iim*jjmp1*klev,ndex3d) 2661 2665 c call writephy(fid_day,prof3d_av,'vitv',v_seri, … … 2663 2667 c 2664 2668 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) 2665 CALL histwrite(nid_day,"vitw",ita p,zx_tmp_3d,2669 CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d, 2666 2670 . iim*jjmp1*klev,ndex3d) 2667 2671 c call writephy(fid_day,prof3d_av,'vitw',omega, … … 2669 2673 c 2670 2674 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) 2671 CALL histwrite(nid_day,"pres",ita p,zx_tmp_3d,2675 CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d, 2672 2676 . iim*jjmp1*klev,ndex3d) 2673 2677 c call writephy(fid_day,prof3d_av,'pres',pplay, … … 2690 2694 zsto = dtime 2691 2695 zout = dtime * ecrit_mth 2696 itau_w = itau_phy + itap 2692 2697 2693 2698 i = NINT(zout/zsto) … … 2700 2705 2701 2706 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 2702 CALL histwrite(nid_mth,"tsol",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2707 CALL histwrite(nid_mth,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2703 2708 c 2704 2709 DO i = 1, klon … … 2706 2711 ENDDO 2707 2712 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2708 CALL histwrite(nid_mth,"psol",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2713 CALL histwrite(nid_mth,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2709 2714 c 2710 2715 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsol,zx_tmp_2d) 2711 CALL histwrite(nid_mth,"qsol",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2716 CALL histwrite(nid_mth,"qsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2712 2717 c 2713 2718 DO i = 1, klon … … 2715 2720 ENDDO 2716 2721 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2717 CALL histwrite(nid_mth,"rain",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2722 CALL histwrite(nid_mth,"rain",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2718 2723 c 2719 2724 DO i = 1, klon … … 2721 2726 ENDDO 2722 2727 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2723 CALL histwrite(nid_mth,"plul",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2728 CALL histwrite(nid_mth,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2724 2729 c 2725 2730 DO i = 1, klon … … 2727 2732 ENDDO 2728 2733 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2729 CALL histwrite(nid_mth,"pluc",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2734 CALL histwrite(nid_mth,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2730 2735 c 2731 2736 CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) 2732 CALL histwrite(nid_mth,"snow",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2737 CALL histwrite(nid_mth,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2733 2738 c 2734 2739 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d) 2735 CALL histwrite(nid_mth,"snow_cov",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2740 CALL histwrite(nid_mth,"snow_cov",itau_w,zx_tmp_2d,iim*jjmp1, 2741 . ndex2d) 2736 2742 c 2737 2743 CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) 2738 CALL histwrite(nid_mth,"evap",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2744 CALL histwrite(nid_mth,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2739 2745 c 2740 2746 CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) 2741 CALL histwrite(nid_mth,"tops",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2747 CALL histwrite(nid_mth,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2742 2748 c 2743 2749 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) 2744 CALL histwrite(nid_mth,"topl",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2750 CALL histwrite(nid_mth,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2745 2751 c 2746 2752 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) 2747 CALL histwrite(nid_mth,"sols",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2753 CALL histwrite(nid_mth,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2748 2754 c 2749 2755 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) 2750 CALL histwrite(nid_mth,"soll",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2756 CALL histwrite(nid_mth,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2751 2757 c 2752 2758 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d) 2753 CALL histwrite(nid_mth,"solldown",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2759 CALL histwrite(nid_mth,"solldown",itau_w,zx_tmp_2d,iim*jjmp1, 2760 . ndex2d) 2754 2761 c 2755 2762 CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d) 2756 CALL histwrite(nid_mth,"tops0",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2763 CALL histwrite(nid_mth,"tops0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2757 2764 c 2758 2765 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d) 2759 CALL histwrite(nid_mth,"topl0",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2766 CALL histwrite(nid_mth,"topl0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2760 2767 c 2761 2768 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d) 2762 CALL histwrite(nid_mth,"sols0",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2769 CALL histwrite(nid_mth,"sols0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2763 2770 c 2764 2771 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d) 2765 CALL histwrite(nid_mth,"soll0",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2772 CALL histwrite(nid_mth,"soll0",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2766 2773 c 2767 2774 CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) 2768 CALL histwrite(nid_mth,"bils",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2775 CALL histwrite(nid_mth,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2769 2776 c 2770 2777 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) 2771 CALL histwrite(nid_mth,"sens",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2778 CALL histwrite(nid_mth,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2772 2779 c 2773 2780 CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) 2774 CALL histwrite(nid_mth,"fder",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2781 CALL histwrite(nid_mth,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2775 2782 c 2776 2783 c … … 2779 2786 c ENDDO 2780 2787 c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2781 c CALL histwrite(nid_mth,"frtu",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2788 c CALL histwrite(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2782 2789 c 2783 2790 c DO i = 1, klon … … 2785 2792 c ENDDO 2786 2793 c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2787 c CALL histwrite(nid_mth,"frtv",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2794 c CALL histwrite(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2788 2795 c 2789 2796 DO nsrf = 1, nbsrf … … 2791 2798 zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf) 2792 2799 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2793 CALL histwrite(nid_mth,"pourc_"//clnsurf(nsrf),ita p,2800 CALL histwrite(nid_mth,"pourc_"//clnsurf(nsrf),itau_w, 2794 2801 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2795 2802 C 2796 2803 zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf) 2797 2804 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2798 CALL histwrite(nid_mth,"tsol_"//clnsurf(nsrf),ita p,2805 CALL histwrite(nid_mth,"tsol_"//clnsurf(nsrf),itau_w, 2799 2806 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2800 2807 C 2801 2808 zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf) 2802 2809 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2803 CALL histwrite(nid_mth,"sens_"//clnsurf(nsrf),ita p,2810 CALL histwrite(nid_mth,"sens_"//clnsurf(nsrf),itau_w, 2804 2811 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2805 2812 C 2806 2813 zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf) 2807 2814 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2808 CALL histwrite(nid_mth,"lat_"//clnsurf(nsrf),ita p,2815 CALL histwrite(nid_mth,"lat_"//clnsurf(nsrf),itau_w, 2809 2816 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2810 2817 C 2811 2818 zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf) 2812 2819 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2813 CALL histwrite(nid_mth,"taux_"//clnsurf(nsrf),ita p,2820 CALL histwrite(nid_mth,"taux_"//clnsurf(nsrf),itau_w, 2814 2821 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2815 2822 C 2816 2823 zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf) 2817 2824 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2818 CALL histwrite(nid_mth,"tauy_"//clnsurf(nsrf),ita p,2825 CALL histwrite(nid_mth,"tauy_"//clnsurf(nsrf),itau_w, 2819 2826 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2820 2827 C 2821 2828 zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf) 2822 2829 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2823 CALL histwrite(nid_mth,"albe_"//clnsurf(nsrf),ita p,2830 CALL histwrite(nid_mth,"albe_"//clnsurf(nsrf),itau_w, 2824 2831 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2825 2832 C 2826 2833 zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf) 2827 2834 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 2828 CALL histwrite(nid_mth,"rugs_"//clnsurf(nsrf),ita p,2835 CALL histwrite(nid_mth,"rugs_"//clnsurf(nsrf),itau_w, 2829 2836 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2830 2837 c 2831 2838 zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf) 2832 2839 CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d) 2833 CALL histwrite(nid_mth,"ages_"//clnsurf(nsrf),ita p2840 CALL histwrite(nid_mth,"ages_"//clnsurf(nsrf),itau_w 2834 2841 $ ,zx_tmp_2d,iim*jjmp1,ndex2d) 2835 2842 … … 2839 2846 c$$$ ENDDO 2840 2847 c$$$ CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2841 c$$$ CALL histwrite(nid_mth,"sicf",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2848 c$$$ CALL histwrite(nid_mth,"sicf",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2842 2849 c 2843 2850 CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d) 2844 CALL histwrite(nid_mth,"albs",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2851 CALL histwrite(nid_mth,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2845 2852 CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d) 2846 CALL histwrite(nid_mth,"albslw",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2853 CALL histwrite(nid_mth,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2847 2854 c 2848 2855 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d) 2849 CALL histwrite(nid_mth,"cdrm",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2856 CALL histwrite(nid_mth,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2850 2857 c 2851 2858 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d) 2852 CALL histwrite(nid_mth,"cdrh",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2859 CALL histwrite(nid_mth,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2853 2860 c 2854 2861 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) 2855 CALL histwrite(nid_mth,"cldl",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2862 CALL histwrite(nid_mth,"cldl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2856 2863 c 2857 2864 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) 2858 CALL histwrite(nid_mth,"cldm",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2865 CALL histwrite(nid_mth,"cldm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2859 2866 c 2860 2867 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) 2861 CALL histwrite(nid_mth,"cldh",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2868 CALL histwrite(nid_mth,"cldh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2862 2869 c 2863 2870 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d) 2864 CALL histwrite(nid_mth,"cldt",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2871 CALL histwrite(nid_mth,"cldt",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2865 2872 c 2866 2873 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) 2867 CALL histwrite(nid_mth,"cldq",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2874 CALL histwrite(nid_mth,"cldq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2868 2875 c 2869 2876 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d) 2870 CALL histwrite(nid_mth,"ue",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2877 CALL histwrite(nid_mth,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2871 2878 c 2872 2879 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d) 2873 CALL histwrite(nid_mth,"ve",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2880 CALL histwrite(nid_mth,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2874 2881 c 2875 2882 CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d) 2876 CALL histwrite(nid_mth,"uq",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2883 CALL histwrite(nid_mth,"uq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2877 2884 c 2878 2885 CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d) 2879 CALL histwrite(nid_mth,"vq",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2886 CALL histwrite(nid_mth,"vq",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2880 2887 cKE43 2881 2888 IF (iflag_con .GE. 3) THEN ! sb 2882 2889 c 2883 2890 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d) 2884 CALL histwrite(nid_mth,"cape",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2891 CALL histwrite(nid_mth,"cape",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2885 2892 c 2886 2893 CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d) 2887 CALL histwrite(nid_mth,"pbase",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2894 CALL histwrite(nid_mth,"pbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2888 2895 c 2889 2896 CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d) 2890 CALL histwrite(nid_mth,"ptop",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2897 CALL histwrite(nid_mth,"ptop",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2891 2898 c 2892 2899 CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d) 2893 CALL histwrite(nid_mth,"fbase",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)2900 CALL histwrite(nid_mth,"fbase",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 2894 2901 c 2895 2902 c … … 2900 2907 C 2901 2908 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) 2902 CALL histwrite(nid_mth,"temp",ita p,zx_tmp_3d,2909 CALL histwrite(nid_mth,"temp",itau_w,zx_tmp_3d, 2903 2910 . iim*jjmp1*klev,ndex3d) 2904 2911 c 2905 2912 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) 2906 CALL histwrite(nid_mth,"ovap",ita p,zx_tmp_3d,2913 CALL histwrite(nid_mth,"ovap",itau_w,zx_tmp_3d, 2907 2914 . iim*jjmp1*klev,ndex3d) 2908 2915 c 2909 2916 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) 2910 CALL histwrite(nid_mth,"geop",ita p,zx_tmp_3d,2917 CALL histwrite(nid_mth,"geop",itau_w,zx_tmp_3d, 2911 2918 . iim*jjmp1*klev,ndex3d) 2912 2919 c 2913 2920 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 2914 CALL histwrite(nid_mth,"vitu",ita p,zx_tmp_3d,2921 CALL histwrite(nid_mth,"vitu",itau_w,zx_tmp_3d, 2915 2922 . iim*jjmp1*klev,ndex3d) 2916 2923 c 2917 2924 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 2918 CALL histwrite(nid_mth,"vitv",ita p,zx_tmp_3d,2925 CALL histwrite(nid_mth,"vitv",itau_w,zx_tmp_3d, 2919 2926 . iim*jjmp1*klev,ndex3d) 2920 2927 c 2921 2928 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) 2922 CALL histwrite(nid_mth,"vitw",ita p,zx_tmp_3d,2929 CALL histwrite(nid_mth,"vitw",itau_w,zx_tmp_3d, 2923 2930 . iim*jjmp1*klev,ndex3d) 2924 2931 c 2925 2932 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) 2926 CALL histwrite(nid_mth,"pres",ita p,zx_tmp_3d,2933 CALL histwrite(nid_mth,"pres",itau_w,zx_tmp_3d, 2927 2934 . iim*jjmp1*klev,ndex3d) 2928 2935 c 2929 2936 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d) 2930 CALL histwrite(nid_mth,"rneb",ita p,zx_tmp_3d,2937 CALL histwrite(nid_mth,"rneb",itau_w,zx_tmp_3d, 2931 2938 . iim*jjmp1*klev,ndex3d) 2932 2939 c 2933 2940 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d) 2934 CALL histwrite(nid_mth,"rhum",ita p,zx_tmp_3d,2941 CALL histwrite(nid_mth,"rhum",itau_w,zx_tmp_3d, 2935 2942 . iim*jjmp1*klev,ndex3d) 2936 2943 c 2937 2944 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d) 2938 CALL histwrite(nid_mth,"oliq",ita p,zx_tmp_3d,2945 CALL histwrite(nid_mth,"oliq",itau_w,zx_tmp_3d, 2939 2946 . iim*jjmp1*klev,ndex3d) 2940 2947 c 2941 2948 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d) 2942 CALL histwrite(nid_mth,"dtdyn",ita p,zx_tmp_3d,2949 CALL histwrite(nid_mth,"dtdyn",itau_w,zx_tmp_3d, 2943 2950 . iim*jjmp1*klev,ndex3d) 2944 2951 c 2945 2952 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d) 2946 CALL histwrite(nid_mth,"dqdyn",ita p,zx_tmp_3d,2953 CALL histwrite(nid_mth,"dqdyn",itau_w,zx_tmp_3d, 2947 2954 . iim*jjmp1*klev,ndex3d) 2948 2955 c 2949 2956 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_con, zx_tmp_3d) 2950 CALL histwrite(nid_mth,"dtcon",ita p,zx_tmp_3d,2957 CALL histwrite(nid_mth,"dtcon",itau_w,zx_tmp_3d, 2951 2958 . iim*jjmp1*klev,ndex3d) 2952 2959 c 2953 2960 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d) 2954 CALL histwrite(nid_mth,"dqcon",ita p,zx_tmp_3d,2961 CALL histwrite(nid_mth,"dqcon",itau_w,zx_tmp_3d, 2955 2962 . iim*jjmp1*klev,ndex3d) 2956 2963 c 2957 2964 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_lsc, zx_tmp_3d) 2958 CALL histwrite(nid_mth,"dtlsc",ita p,zx_tmp_3d,2965 CALL histwrite(nid_mth,"dtlsc",itau_w,zx_tmp_3d, 2959 2966 . iim*jjmp1*klev,ndex3d) 2960 2967 c 2961 2968 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d) 2962 CALL histwrite(nid_mth,"dqlsc",ita p,zx_tmp_3d,2969 CALL histwrite(nid_mth,"dqlsc",itau_w,zx_tmp_3d, 2963 2970 . iim*jjmp1*klev,ndex3d) 2964 2971 c 2965 2972 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d) 2966 CALL histwrite(nid_mth,"dtvdf",ita p,zx_tmp_3d,2973 CALL histwrite(nid_mth,"dtvdf",itau_w,zx_tmp_3d, 2967 2974 . iim*jjmp1*klev,ndex3d) 2968 2975 c 2969 2976 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d) 2970 CALL histwrite(nid_mth,"dqvdf",ita p,zx_tmp_3d,2977 CALL histwrite(nid_mth,"dqvdf",itau_w,zx_tmp_3d, 2971 2978 . iim*jjmp1*klev,ndex3d) 2972 2979 c 2973 2980 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_eva, zx_tmp_3d) 2974 CALL histwrite(nid_mth,"dteva",ita p,zx_tmp_3d,2981 CALL histwrite(nid_mth,"dteva",itau_w,zx_tmp_3d, 2975 2982 . iim*jjmp1*klev,ndex3d) 2976 2983 c 2977 2984 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d) 2978 CALL histwrite(nid_mth,"dqeva",ita p,zx_tmp_3d,2985 CALL histwrite(nid_mth,"dqeva",itau_w,zx_tmp_3d, 2979 2986 . iim*jjmp1*klev,ndex3d) 2980 2987 c 2981 2988 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d) 2982 CALL histwrite(nid_mth,"ptconv",ita p,zx_tmp_3d,2989 CALL histwrite(nid_mth,"ptconv",itau_w,zx_tmp_3d, 2983 2990 . iim*(jjmp1)*klev,ndex3d) 2984 2991 c 2985 2992 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d) 2986 CALL histwrite(nid_mth,"ratqs",ita p,zx_tmp_3d,2993 CALL histwrite(nid_mth,"ratqs",itau_w,zx_tmp_3d, 2987 2994 . iim*(jjmp1)*klev,ndex3d) 2988 2995 c 2989 2996 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ajs, zx_tmp_3d) 2990 CALL histwrite(nid_mth,"dtajs",ita p,zx_tmp_3d,2997 CALL histwrite(nid_mth,"dtajs",itau_w,zx_tmp_3d, 2991 2998 . iim*jjmp1*klev,ndex3d) 2992 2999 c 2993 3000 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d) 2994 CALL histwrite(nid_mth,"dqajs",ita p,zx_tmp_3d,3001 CALL histwrite(nid_mth,"dqajs",itau_w,zx_tmp_3d, 2995 3002 . iim*jjmp1*klev,ndex3d) 2996 3003 c 2997 3004 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat, zx_tmp_3d) 2998 CALL histwrite(nid_mth,"dtswr",ita p,zx_tmp_3d,3005 CALL histwrite(nid_mth,"dtswr",itau_w,zx_tmp_3d, 2999 3006 . iim*jjmp1*klev,ndex3d) 3000 3007 c 3001 3008 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat0, zx_tmp_3d) 3002 CALL histwrite(nid_mth,"dtsw0",ita p,zx_tmp_3d,3009 CALL histwrite(nid_mth,"dtsw0",itau_w,zx_tmp_3d, 3003 3010 . iim*jjmp1*klev,ndex3d) 3004 3011 c 3005 3012 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool, zx_tmp_3d) 3006 CALL histwrite(nid_mth,"dtlwr",ita p,zx_tmp_3d,3013 CALL histwrite(nid_mth,"dtlwr",itau_w,zx_tmp_3d, 3007 3014 . iim*jjmp1*klev,ndex3d) 3008 3015 c 3009 3016 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool0, zx_tmp_3d) 3010 CALL histwrite(nid_mth,"dtlw0",ita p,zx_tmp_3d,3017 CALL histwrite(nid_mth,"dtlw0",itau_w,zx_tmp_3d, 3011 3018 . iim*jjmp1*klev,ndex3d) 3012 3019 c 3013 3020 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d) 3014 CALL histwrite(nid_mth,"duvdf",ita p,zx_tmp_3d,3021 CALL histwrite(nid_mth,"duvdf",itau_w,zx_tmp_3d, 3015 3022 . iim*jjmp1*klev,ndex3d) 3016 3023 c 3017 3024 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d) 3018 CALL histwrite(nid_mth,"dvvdf",ita p,zx_tmp_3d,3025 CALL histwrite(nid_mth,"dvvdf",itau_w,zx_tmp_3d, 3019 3026 . iim*jjmp1*klev,ndex3d) 3020 3027 c 3021 3028 IF (ok_orodr) THEN 3022 3029 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d) 3023 CALL histwrite(nid_mth,"duoro",ita p,zx_tmp_3d,3030 CALL histwrite(nid_mth,"duoro",itau_w,zx_tmp_3d, 3024 3031 . iim*jjmp1*klev,ndex3d) 3025 3032 c 3026 3033 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oro, zx_tmp_3d) 3027 CALL histwrite(nid_mth,"dvoro",ita p,zx_tmp_3d,3034 CALL histwrite(nid_mth,"dvoro",itau_w,zx_tmp_3d, 3028 3035 . iim*jjmp1*klev,ndex3d) 3029 3036 c … … 3032 3039 IF (ok_orolf) THEN 3033 3040 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_lif, zx_tmp_3d) 3034 CALL histwrite(nid_mth,"dulif",ita p,zx_tmp_3d,3041 CALL histwrite(nid_mth,"dulif",itau_w,zx_tmp_3d, 3035 3042 . iim*jjmp1*klev,ndex3d) 3036 3043 c 3037 3044 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_lif, zx_tmp_3d) 3038 CALL histwrite(nid_mth,"dvlif",ita p,zx_tmp_3d,3045 CALL histwrite(nid_mth,"dvlif",itau_w,zx_tmp_3d, 3039 3046 . iim*jjmp1*klev,ndex3d) 3040 3047 ENDIF 3041 3048 C 3042 3049 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d) 3043 CALL histwrite(nid_mth,"ozone",ita p,zx_tmp_3d,3050 CALL histwrite(nid_mth,"ozone",itau_w,zx_tmp_3d, 3044 3051 . iim*jjmp1*klev,ndex3d) 3045 3052 c … … 3049 3056 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d) 3050 3057 WRITE(str2,'(i2.2)') iq 3051 CALL histwrite(nid_mth,"trac"//str2,ita p,zx_tmp_3d,3058 CALL histwrite(nid_mth,"trac"//str2,itau_w,zx_tmp_3d, 3052 3059 . iim*jjmp1*klev,ndex3d) 3053 3060 ELSE … … 3061 3068 c 3062 3069 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d) 3063 CALL histwrite(nid_mth,"upwd",ita p,zx_tmp_3d,3070 CALL histwrite(nid_mth,"upwd",itau_w,zx_tmp_3d, 3064 3071 . iim*jjmp1*klev,ndex3d) 3065 3072 c 3066 3073 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d) 3067 CALL histwrite(nid_mth,"dnwd",ita p,zx_tmp_3d,3074 CALL histwrite(nid_mth,"dnwd",itau_w,zx_tmp_3d, 3068 3075 . iim*jjmp1*klev,ndex3d) 3069 3076 c 3070 3077 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d) 3071 CALL histwrite(nid_mth,"dnwd0",ita p,zx_tmp_3d,3078 CALL histwrite(nid_mth,"dnwd0",itau_w,zx_tmp_3d, 3072 3079 . iim*jjmp1*klev,ndex3d) 3073 3080 c 3074 3081 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d) 3075 CALL histwrite(nid_mth,"Ma",ita p,zx_tmp_3d,3082 CALL histwrite(nid_mth,"Ma",itau_w,zx_tmp_3d, 3076 3083 . iim*jjmp1*klev,ndex3d) 3077 3084 c … … 3094 3101 zsto = dtime * ecrit_ins 3095 3102 zout = dtime * ecrit_ins 3103 itau_w = itau_phy + itap 3096 3104 3097 3105 i = NINT(zout/zsto) … … 3107 3115 ENDDO 3108 3116 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 3109 CALL histwrite(nid_ins,"psol",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3117 CALL histwrite(nid_ins,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3110 3118 c 3111 3119 DO i = 1, klon … … 3113 3121 ENDDO 3114 3122 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 3115 CALL histwrite(nid_ins,"rain",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3123 CALL histwrite(nid_ins,"rain",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3116 3124 c 3117 3125 DO i = 1, klon … … 3119 3127 ENDDO 3120 3128 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 3121 CALL histwrite(nid_ins,"plul",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3129 CALL histwrite(nid_ins,"plul",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3122 3130 c 3123 3131 DO i = 1, klon … … 3125 3133 ENDDO 3126 3134 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 3127 CALL histwrite(nid_ins,"pluc",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3135 CALL histwrite(nid_ins,"pluc",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3128 3136 3129 3137 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 3130 CALL histwrite(nid_ins,"tsol",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3138 CALL histwrite(nid_ins,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3131 3139 c 3132 3140 CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) 3133 CALL histwrite(nid_ins,"snow",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3141 CALL histwrite(nid_ins,"snow",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3134 3142 3135 3143 c 3136 3144 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) 3137 CALL histwrite(nid_ins,"topl",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3145 CALL histwrite(nid_ins,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3138 3146 c 3139 3147 CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) 3140 CALL histwrite(nid_ins,"evap",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3148 CALL histwrite(nid_ins,"evap",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3141 3149 c 3142 3150 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) 3143 CALL histwrite(nid_ins,"sols",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3151 CALL histwrite(nid_ins,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3144 3152 c 3145 3153 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) 3146 CALL histwrite(nid_ins,"soll",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3154 CALL histwrite(nid_ins,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3147 3155 c 3148 3156 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d) 3149 CALL histwrite(nid_ins,"solldown",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 3157 CALL histwrite(nid_ins,"solldown",itau_w,zx_tmp_2d,iim*jjmp1, 3158 . ndex2d) 3150 3159 c 3151 3160 CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) 3152 CALL histwrite(nid_ins,"bils",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3161 CALL histwrite(nid_ins,"bils",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3153 3162 c 3154 3163 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) 3155 CALL histwrite(nid_ins,"sens",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3164 CALL histwrite(nid_ins,"sens",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3156 3165 c 3157 3166 CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) 3158 CALL histwrite(nid_ins,"fder",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3167 CALL histwrite(nid_ins,"fder",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3159 3168 c 3160 3169 CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_oce),zx_tmp_2d) 3161 CALL histwrite(nid_ins,"dtsvdfo",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 3170 CALL histwrite(nid_ins,"dtsvdfo",itau_w,zx_tmp_2d,iim*jjmp1, 3171 . ndex2d) 3162 3172 c 3163 3173 CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_ter),zx_tmp_2d) 3164 CALL histwrite(nid_ins,"dtsvdft",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 3174 CALL histwrite(nid_ins,"dtsvdft",itau_w,zx_tmp_2d,iim*jjmp1, 3175 . ndex2d) 3165 3176 c 3166 3177 CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_lic),zx_tmp_2d) 3167 CALL histwrite(nid_ins,"dtsvdfg",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 3178 CALL histwrite(nid_ins,"dtsvdfg",itau_w,zx_tmp_2d,iim*jjmp1, 3179 . ndex2d) 3168 3180 c 3169 3181 CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_sic),zx_tmp_2d) 3170 CALL histwrite(nid_ins,"dtsvdfi",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 3182 CALL histwrite(nid_ins,"dtsvdfi",itau_w,zx_tmp_2d,iim*jjmp1, 3183 . ndex2d) 3171 3184 3172 3185 DO nsrf = 1, nbsrf … … 3174 3187 zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf) 3175 3188 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 3176 CALL histwrite(nid_ins,"pourc_"//clnsurf(nsrf),ita p,3189 CALL histwrite(nid_ins,"pourc_"//clnsurf(nsrf),itau_w, 3177 3190 $ zx_tmp_2d,iim*jjmp1,ndex2d) 3178 3191 C 3179 3192 zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf) 3180 3193 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 3181 CALL histwrite(nid_ins,"sens_"//clnsurf(nsrf),ita p,3194 CALL histwrite(nid_ins,"sens_"//clnsurf(nsrf),itau_w, 3182 3195 $ zx_tmp_2d,iim*jjmp1,ndex2d) 3183 3196 C 3184 3197 zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf) 3185 3198 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 3186 CALL histwrite(nid_ins,"lat_"//clnsurf(nsrf),ita p,3199 CALL histwrite(nid_ins,"lat_"//clnsurf(nsrf),itau_w, 3187 3200 $ zx_tmp_2d,iim*jjmp1,ndex2d) 3188 3201 C 3189 3202 zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf) 3190 3203 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 3191 CALL histwrite(nid_ins,"tsol_"//clnsurf(nsrf),ita p,3204 CALL histwrite(nid_ins,"tsol_"//clnsurf(nsrf),itau_w, 3192 3205 $ zx_tmp_2d,iim*jjmp1,ndex2d) 3193 3206 C 3194 3207 zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf) 3195 3208 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 3196 CALL histwrite(nid_ins,"taux_"//clnsurf(nsrf),ita p,3209 CALL histwrite(nid_ins,"taux_"//clnsurf(nsrf),itau_w, 3197 3210 $ zx_tmp_2d,iim*jjmp1,ndex2d) 3198 3211 C 3199 3212 zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf) 3200 3213 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 3201 CALL histwrite(nid_ins,"tauy_"//clnsurf(nsrf),ita p,3214 CALL histwrite(nid_ins,"tauy_"//clnsurf(nsrf),itau_w, 3202 3215 $ zx_tmp_2d,iim*jjmp1,ndex2d) 3203 3216 C 3204 3217 zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf) 3205 3218 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 3206 CALL histwrite(nid_ins,"rugs_"//clnsurf(nsrf),ita p,3219 CALL histwrite(nid_ins,"rugs_"//clnsurf(nsrf),itau_w, 3207 3220 $ zx_tmp_2d,iim*jjmp1,ndex2d) 3208 3221 C 3209 3222 zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf) 3210 3223 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) 3211 CALL histwrite(nid_ins,"albe_"//clnsurf(nsrf),ita p,3224 CALL histwrite(nid_ins,"albe_"//clnsurf(nsrf),itau_w, 3212 3225 $ zx_tmp_2d,iim*jjmp1,ndex2d) 3213 3226 C 3214 3227 END DO 3215 3228 CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d) 3216 CALL histwrite(nid_ins,"albs",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3229 CALL histwrite(nid_ins,"albs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3217 3230 CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d) 3218 CALL histwrite(nid_ins,"albslw",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3231 CALL histwrite(nid_ins,"albslw",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3219 3232 c 3220 3233 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d) 3221 CALL histwrite(nid_ins,"snow_cov",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 3234 CALL histwrite(nid_ins,"snow_cov",itau_w,zx_tmp_2d,iim*jjmp1, 3235 . ndex2d) 3222 3236 c 3223 3237 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxrugs,zx_tmp_2d) 3224 CALL histwrite(nid_ins,"rugs",ita p,zx_tmp_2d,iim*jjmp1,ndex2d)3238 CALL histwrite(nid_ins,"rugs",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 3225 3239 c 3226 3240 c Champs 3D: 3227 3241 c 3228 3242 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) 3229 CALL histwrite(nid_ins,"temp",ita p,zx_tmp_3d,3243 CALL histwrite(nid_ins,"temp",itau_w,zx_tmp_3d, 3230 3244 . iim*jjmp1*klev,ndex3d) 3231 3245 c 3232 3246 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 3233 CALL histwrite(nid_ins,"vitu",ita p,zx_tmp_3d,3247 CALL histwrite(nid_ins,"vitu",itau_w,zx_tmp_3d, 3234 3248 . iim*jjmp1*klev,ndex3d) 3235 3249 c 3236 3250 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 3237 CALL histwrite(nid_ins,"vitv",ita p,zx_tmp_3d,3251 CALL histwrite(nid_ins,"vitv",itau_w,zx_tmp_3d, 3238 3252 . iim*jjmp1*klev,ndex3d) 3239 3253 c 3240 3254 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) 3241 CALL histwrite(nid_ins,"geop",ita p,zx_tmp_3d,3255 CALL histwrite(nid_ins,"geop",itau_w,zx_tmp_3d, 3242 3256 . iim*jjmp1*klev,ndex3d) 3243 3257 c 3244 3258 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) 3245 CALL histwrite(nid_ins,"pres",ita p,zx_tmp_3d,3259 CALL histwrite(nid_ins,"pres",itau_w,zx_tmp_3d, 3246 3260 . iim*jjmp1*klev,ndex3d) 3247 3261 c 3248 3262 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d) 3249 CALL histwrite(nid_ins,"dtvdf",ita p,zx_tmp_3d,3263 CALL histwrite(nid_ins,"dtvdf",itau_w,zx_tmp_3d, 3250 3264 . iim*jjmp1*klev,ndex3d) 3251 3265 c 3252 3266 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d) 3253 CALL histwrite(nid_ins,"dqvdf",ita p,zx_tmp_3d,3267 CALL histwrite(nid_ins,"dqvdf",itau_w,zx_tmp_3d, 3254 3268 . iim*jjmp1*klev,ndex3d) 3255 3269 … … 3339 3353 c 3340 3354 IF (lafin) THEN 3355 itau_phy = itau_phy + itap 3341 3356 ccc IF (ok_oasis) CALL quitcpl 3342 3357 CALL phyredem ("restartphy.nc",dtime,radpas,co2_ppm,solaire,
Note: See TracChangeset
for help on using the changeset viewer.