Changeset 80
- Timestamp:
- Apr 18, 2000, 10:37:55 AM (25 years ago)
- File:
-
- 1 edited
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F (modified) (30 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F
r79 r80 50 50 c====================================================================== 51 51 #include "dimensions.h" 52 integer jjmp1 53 parameter (jjmp1=jjm+1-1/jjm) 52 54 #include "dimphy.h" 53 55 #include "regdim.h" … … 72 74 PARAMETER (itimestep=1800) 73 75 EXTERNAL fromcpl, intocpl, inicma 74 REAL cpl_sst(iim,jjm +1), cpl_sic(iim,jjm+1)75 REAL cpl_alb_sst(iim,jjm +1), cpl_alb_sic(iim,jjm+1)76 REAL cpl_sst(iim,jjmp1), cpl_sic(iim,jjmp1) 77 REAL cpl_alb_sst(iim,jjmp1), cpl_alb_sic(iim,jjmp1) 76 78 c====================================================================== 77 79 c ok_ocean indique l'utilisation du modele oceanique "slab ocean", … … 120 122 c 121 123 LOGICAL ok_instan ! sortir le fichier instantane 122 PARAMETER (ok_instan=. FALSE.)124 PARAMETER (ok_instan=.TRUE.) 123 125 c 124 126 LOGICAL ok_region ! sortir le fichier regional … … 152 154 REAL presnivs(klev) 153 155 REAL znivsig(klev) 156 REAL zsurf(nbsrf) 154 157 155 158 REAL u(klon,klev) … … 283 286 REAL yv1(klon) ! vents dans la premiere couche V 284 287 LOGICAL offline ! Controle du stockage ds "physique" 285 PARAMETER (offline=. true.)288 PARAMETER (offline=.false.) 286 289 INTEGER physid 287 290 REAL pfrac_impa(klon,klev)! Produits des coefs lessivage impaction … … 421 424 c 422 425 REAL zphi(klon,klev) 423 REAL zx_tmp_x(iim), zx_tmp_y (jjm+1)424 REAL zx_relief(iim,jjm +1)425 REAL zx_aire(iim,jjm +1)426 REAL zx_tmp_x(iim), zx_tmp_yjjmp1 427 REAL zx_relief(iim,jjmp1) 428 REAL zx_aire(iim,jjmp1) 426 429 c 427 430 c Variables du changement … … 473 476 SAVE ecrit_reg ! frequence d'ecriture 474 477 c 475 REAL oas_sols(klon), z_sols(iim,jjm +1)478 REAL oas_sols(klon), z_sols(iim,jjmp1) 476 479 SAVE oas_sols 477 REAL oas_nsol(klon), z_nsol(iim,jjm +1)480 REAL oas_nsol(klon), z_nsol(iim,jjmp1) 478 481 SAVE oas_nsol 479 REAL oas_rain(klon), z_rain(iim,jjm +1)482 REAL oas_rain(klon), z_rain(iim,jjmp1) 480 483 SAVE oas_rain 481 REAL oas_snow(klon), z_snow(iim,jjm +1)484 REAL oas_snow(klon), z_snow(iim,jjmp1) 482 485 SAVE oas_snow 483 REAL oas_evap(klon), z_evap(iim,jjm +1)486 REAL oas_evap(klon), z_evap(iim,jjmp1) 484 487 SAVE oas_evap 485 REAL oas_ruis(klon), z_ruis(iim,jjm +1)488 REAL oas_ruis(klon), z_ruis(iim,jjmp1) 486 489 SAVE oas_ruis 487 REAL oas_tsol(klon), z_tsol(iim,jjm +1)490 REAL oas_tsol(klon), z_tsol(iim,jjmp1) 488 491 SAVE oas_tsol 489 REAL oas_fder(klon), z_fder(iim,jjm +1)492 REAL oas_fder(klon), z_fder(iim,jjmp1) 490 493 SAVE oas_fder 491 REAL oas_albe(klon), z_albe(iim,jjm +1)494 REAL oas_albe(klon), z_albe(iim,jjmp1) 492 495 SAVE oas_albe 493 REAL oas_taux(klon), z_taux(iim,jjm +1)496 REAL oas_taux(klon), z_taux(iim,jjmp1) 494 497 SAVE oas_taux 495 REAL oas_tauy(klon), z_tauy(iim,jjm +1)498 REAL oas_tauy(klon), z_tauy(iim,jjmp1) 496 499 SAVE oas_tauy 497 REAL oas_ruisoce(klon), z_ruisoce(iim,jjm +1)500 REAL oas_ruisoce(klon), z_ruisoce(iim,jjmp1) 498 501 SAVE oas_ruisoce 499 REAL oas_ruisriv(klon), z_ruisriv(iim,jjm +1)502 REAL oas_ruisriv(klon), z_ruisriv(iim,jjmp1) 500 503 SAVE oas_ruisriv 501 504 c … … 518 521 REAL tabcntr0( length ) 519 522 c 520 INTEGER ndex2d(iim* (jjm+1)),ndex3d(iim*(jjm+1)*klev)523 INTEGER ndex2d(iim*jjmp1),ndex3d(iim*jjmp1*klev) 521 524 REAL zx_tmp_fi2d(klon) 522 REAL zx_tmp_2d(iim,jjm +1), zx_tmp_3d(iim,jjm+1,klev)523 REAL zx_lon(iim,jjm +1), zx_lat(iim,jjm+1)525 REAL zx_tmp_2d(iim,jjmp1), zx_tmp_3d(iim,jjmp1,klev) 526 REAL zx_lon(iim,jjmp1), zx_lat(iim,jjmp1) 524 527 c 525 528 INTEGER nid_day, nid_mth, nid_ins 526 529 SAVE nid_day, nid_mth, nid_ins 527 530 c 528 INTEGER nhori, nvert 531 INTEGER nhori, nvert,nvertsf 529 532 REAL zsto, zout, zjulian 530 533 … … 661 664 ENDIF 662 665 ccc ecrit_ins = NINT(86400./dtime *0.5) ! 2 fois par jour 663 ecrit_ins = NINT(86400./dtime *0.25) ! tous les jours 666 ccc ecrit_ins = NINT(86400./dtime *0.25) ! 4 fois par jour 667 ecrit_ins = NINT(86400./dtime/48.) ! a chaque pas de temps 664 668 IF (ok_instan) THEN 665 669 PRINT*, 'La frequence de sortie instant. est de ', ecrit_ins … … 676 680 zjulian = zjulian + day_ini 677 681 c 678 CALL gr_fi_ecrit(1,klon,iim,jjm +1,rlon,zx_lon)682 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) 679 683 DO i = 1, iim 680 684 zx_lon(i,1) = rlon(i+1) 681 zx_lon(i,jjm +1) = rlon(i+1)685 zx_lon(i,jjmp1) = rlon(i+1) 682 686 ENDDO 683 687 DO ll=1,klev 684 688 znivsig(ll)=float(ll) 685 689 ENDDO 686 CALL gr_fi_ecrit(1,klon,iim,jjm +1,rlat,zx_lat)687 CALL histbeg("histday", iim,zx_lon, jjm +1,zx_lat,688 . 1,iim,1,jjm +1, 0, zjulian, dtime,690 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 691 CALL histbeg("histday", iim,zx_lon, jjmp1,zx_lat, 692 . 1,iim,1,jjmp1, 0, zjulian, dtime, 689 693 . nhori, nid_day) 690 694 c CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", … … 697 701 c 698 702 CALL histdef(nid_day, "phis", "Surface geop. height", "-", 699 . iim,jjm +1,nhori, 1,1,1, -99, 32,703 . iim,jjmp1,nhori, 1,1,1, -99, 32, 700 704 . "once", zsto,zout) 701 705 c 702 706 CALL histdef(nid_day, "aire", "Grid area", "-", 703 . iim,jjm +1,nhori, 1,1,1, -99, 32,707 . iim,jjmp1,nhori, 1,1,1, -99, 32, 704 708 . "once", zsto,zout) 705 709 c … … 707 711 c 708 712 CALL histdef(nid_day, "tsol", "Surface Temperature", "K", 709 . iim,jjm +1,nhori, 1,1,1, -99, 32,713 . iim,jjmp1,nhori, 1,1,1, -99, 32, 710 714 . "ave(X)", zsto,zout) 711 715 c 712 716 CALL histdef(nid_day, "psol", "Surface Pressure", "Pa", 713 . iim,jjm +1,nhori, 1,1,1, -99, 32,717 . iim,jjmp1,nhori, 1,1,1, -99, 32, 714 718 . "ave(X)", zsto,zout) 715 719 c 716 720 CALL histdef(nid_day, "rain", "Precipitation", "mm/day", 717 . iim,jjm +1,nhori, 1,1,1, -99, 32,721 . iim,jjmp1,nhori, 1,1,1, -99, 32, 718 722 . "ave(X)", zsto,zout) 719 723 c 720 724 CALL histdef(nid_day, "snow", "Snow fall", "mm/day", 721 . iim,jjm +1,nhori, 1,1,1, -99, 32,725 . iim,jjmp1,nhori, 1,1,1, -99, 32, 722 726 . "ave(X)", zsto,zout) 723 727 c 724 728 CALL histdef(nid_day, "evap", "Evaporation", "mm/day", 725 . iim,jjm +1,nhori, 1,1,1, -99, 32,729 . iim,jjmp1,nhori, 1,1,1, -99, 32, 726 730 . "ave(X)", zsto,zout) 727 731 c 728 732 CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2", 729 . iim,jjm +1,nhori, 1,1,1, -99, 32,733 . iim,jjmp1,nhori, 1,1,1, -99, 32, 730 734 . "ave(X)", zsto,zout) 731 735 c 732 736 CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2", 733 . iim,jjm +1,nhori, 1,1,1, -99, 32,737 . iim,jjmp1,nhori, 1,1,1, -99, 32, 734 738 . "ave(X)", zsto,zout) 735 739 c 736 740 CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2", 737 . iim,jjm +1,nhori, 1,1,1, -99, 32,741 . iim,jjmp1,nhori, 1,1,1, -99, 32, 738 742 . "ave(X)", zsto,zout) 739 743 c 740 744 CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2", 741 . iim,jjm +1,nhori, 1,1,1, -99, 32,745 . iim,jjmp1,nhori, 1,1,1, -99, 32, 742 746 . "ave(X)", zsto,zout) 743 747 c 744 748 CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2", 745 . iim,jjm +1,nhori, 1,1,1, -99, 32,749 . iim,jjmp1,nhori, 1,1,1, -99, 32, 746 750 . "ave(X)", zsto,zout) 747 751 c 748 752 CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2", 749 . iim,jjm +1,nhori, 1,1,1, -99, 32,753 . iim,jjmp1,nhori, 1,1,1, -99, 32, 750 754 . "ave(X)", zsto,zout) 751 755 c 752 756 CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2", 753 . iim,jjm +1,nhori, 1,1,1, -99, 32,757 . iim,jjmp1,nhori, 1,1,1, -99, 32, 754 758 . "ave(X)", zsto,zout) 755 759 c 756 760 CALL histdef(nid_day, "frtu", "Zonal wind stress", "Pa", 757 . iim,jjm +1,nhori, 1,1,1, -99, 32,761 . iim,jjmp1,nhori, 1,1,1, -99, 32, 758 762 . "ave(X)", zsto,zout) 759 763 c 760 764 CALL histdef(nid_day, "frtv", "Meridional wind stress", "Pa", 761 . iim,jjm +1,nhori, 1,1,1, -99, 32,765 . iim,jjmp1,nhori, 1,1,1, -99, 32, 762 766 . "ave(X)", zsto,zout) 763 767 c 764 768 CALL histdef(nid_day, "ruis", "Runoff", "mm/day", 765 . iim,jjm +1,nhori, 1,1,1, -99, 32,769 . iim,jjmp1,nhori, 1,1,1, -99, 32, 766 770 . "ave(X)", zsto,zout) 767 771 c 768 772 CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-", 769 . iim,jjm +1,nhori, 1,1,1, -99, 32,773 . iim,jjmp1,nhori, 1,1,1, -99, 32, 770 774 . "ave(X)", zsto,zout) 771 775 c 772 776 CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-", 773 . iim,jjm +1,nhori, 1,1,1, -99, 32,777 . iim,jjmp1,nhori, 1,1,1, -99, 32, 774 778 . "ave(X)", zsto,zout) 775 779 c 776 780 CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-", 777 . iim,jjm +1,nhori, 1,1,1, -99, 32,781 . iim,jjmp1,nhori, 1,1,1, -99, 32, 778 782 . "ave(X)", zsto,zout) 779 783 c 780 784 CALL histdef(nid_day, "cldh", "High-level cloudiness", "-", 781 . iim,jjm +1,nhori, 1,1,1, -99, 32,785 . iim,jjmp1,nhori, 1,1,1, -99, 32, 782 786 . "ave(X)", zsto,zout) 783 787 c 784 788 CALL histdef(nid_day, "cldt", "Total cloudiness", "-", 785 . iim,jjm +1,nhori, 1,1,1, -99, 32,789 . iim,jjmp1,nhori, 1,1,1, -99, 32, 786 790 . "ave(X)", zsto,zout) 787 791 c 788 792 CALL histdef(nid_day, "cldq", "Cloud liquid water path", "-", 789 . iim,jjm +1,nhori, 1,1,1, -99, 32,793 . iim,jjmp1,nhori, 1,1,1, -99, 32, 790 794 . "ave(X)", zsto,zout) 791 795 c … … 793 797 c 794 798 CALL histdef(nid_day, "temp", "Air temperature", "K", 795 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,799 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 796 800 . "ave(X)", zsto,zout) 797 801 c 798 802 CALL histdef(nid_day, "ovap", "Specific humidity", "Kg/Kg", 799 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,803 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 800 804 . "ave(X)", zsto,zout) 801 805 c 802 806 CALL histdef(nid_day, "geop", "Geopotential height", "m", 803 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,807 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 804 808 . "ave(X)", zsto,zout) 805 809 c 806 810 CALL histdef(nid_day, "vitu", "Zonal wind", "m/s", 807 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,811 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 808 812 . "ave(X)", zsto,zout) 809 813 c 810 814 CALL histdef(nid_day, "vitv", "Meridional wind", "m/s", 811 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,815 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 812 816 . "ave(X)", zsto,zout) 813 817 c 814 818 CALL histdef(nid_day, "vitw", "Vertical wind", "m/s", 815 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,819 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 816 820 . "ave(X)", zsto,zout) 817 821 c 818 822 CALL histdef(nid_day, "pres", "Air pressure", "Pa", 819 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,823 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 820 824 . "ave(X)", zsto,zout) 821 825 c … … 832 836 zjulian = zjulian + day_ini 833 837 c 834 CALL gr_fi_ecrit(1,klon,iim,jjm +1,rlon,zx_lon)838 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) 835 839 DO i = 1, iim 836 840 zx_lon(i,1) = rlon(i+1) 837 zx_lon(i,jjm +1) = rlon(i+1)841 zx_lon(i,jjmp1) = rlon(i+1) 838 842 ENDDO 839 843 DO ll=1,klev 840 844 znivsig(ll)=float(ll) 841 845 ENDDO 842 CALL gr_fi_ecrit(1,klon,iim,jjm +1,rlat,zx_lat)843 CALL histbeg("histmth", iim,zx_lon, jjm +1,zx_lat,844 . 1,iim,1,jjm +1, 0, zjulian, dtime,846 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 847 CALL histbeg("histmth", iim,zx_lon, jjmp1,zx_lat, 848 . 1,iim,1,jjmp1, 0, zjulian, dtime, 845 849 . nhori, nid_mth) 846 850 c CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb", … … 853 857 c 854 858 CALL histdef(nid_mth, "phis", "Surface geop. height", "-", 855 . iim,jjm +1,nhori, 1,1,1, -99, 32,859 . iim,jjmp1,nhori, 1,1,1, -99, 32, 856 860 . "once", zsto,zout) 857 861 c 858 862 CALL histdef(nid_mth, "aire", "Grid area", "-", 859 . iim,jjm +1,nhori, 1,1,1, -99, 32,863 . iim,jjmp1,nhori, 1,1,1, -99, 32, 860 864 . "once", zsto,zout) 861 865 c … … 863 867 c 864 868 CALL histdef(nid_mth, "tsol", "Surface Temperature", "K", 865 . iim,jjm +1,nhori, 1,1,1, -99, 32,869 . iim,jjmp1,nhori, 1,1,1, -99, 32, 866 870 . "ave(X)", zsto,zout) 867 871 c 868 872 CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa", 869 . iim,jjm +1,nhori, 1,1,1, -99, 32,873 . iim,jjmp1,nhori, 1,1,1, -99, 32, 870 874 . "ave(X)", zsto,zout) 871 875 c 872 876 CALL histdef(nid_mth, "qsol", "Surface humidity", "mm", 873 . iim,jjm +1,nhori, 1,1,1, -99, 32,877 . iim,jjmp1,nhori, 1,1,1, -99, 32, 874 878 . "ave(X)", zsto,zout) 875 879 c 876 880 CALL histdef(nid_mth, "rain", "Precipitation", "mm/day", 877 . iim,jjm +1,nhori, 1,1,1, -99, 32,881 . iim,jjmp1,nhori, 1,1,1, -99, 32, 878 882 . "ave(X)", zsto,zout) 879 883 c 880 884 CALL histdef(nid_mth, "plul", "Large-scale Precip.", "mm/day", 881 . iim,jjm +1,nhori, 1,1,1, -99, 32,885 . iim,jjmp1,nhori, 1,1,1, -99, 32, 882 886 . "ave(X)", zsto,zout) 883 887 c 884 888 CALL histdef(nid_mth, "pluc", "Convective Precip.", "mm/day", 885 . iim,jjm +1,nhori, 1,1,1, -99, 32,889 . iim,jjmp1,nhori, 1,1,1, -99, 32, 886 890 . "ave(X)", zsto,zout) 887 891 c 888 892 CALL histdef(nid_mth, "snow", "Snow fall", "mm/day", 889 . iim,jjm +1,nhori, 1,1,1, -99, 32,893 . iim,jjmp1,nhori, 1,1,1, -99, 32, 890 894 . "ave(X)", zsto,zout) 891 895 c 892 896 CALL histdef(nid_mth, "ages", "Snow age", "day", 893 . iim,jjm +1,nhori, 1,1,1, -99, 32,897 . iim,jjmp1,nhori, 1,1,1, -99, 32, 894 898 . "ave(X)", zsto,zout) 895 899 c 896 900 CALL histdef(nid_mth, "evap", "Evaporation", "mm/day", 897 . iim,jjm +1,nhori, 1,1,1, -99, 32,901 . iim,jjmp1,nhori, 1,1,1, -99, 32, 898 902 . "ave(X)", zsto,zout) 899 903 c 900 904 CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2", 901 . iim,jjm +1,nhori, 1,1,1, -99, 32,905 . iim,jjmp1,nhori, 1,1,1, -99, 32, 902 906 . "ave(X)", zsto,zout) 903 907 c 904 908 CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2", 905 . iim,jjm +1,nhori, 1,1,1, -99, 32,909 . iim,jjmp1,nhori, 1,1,1, -99, 32, 906 910 . "ave(X)", zsto,zout) 907 911 c 908 912 CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2", 909 . iim,jjm +1,nhori, 1,1,1, -99, 32,913 . iim,jjmp1,nhori, 1,1,1, -99, 32, 910 914 . "ave(X)", zsto,zout) 911 915 c 912 916 CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2", 913 . iim,jjm +1,nhori, 1,1,1, -99, 32,917 . iim,jjmp1,nhori, 1,1,1, -99, 32, 914 918 . "ave(X)", zsto,zout) 915 919 c 916 920 CALL histdef(nid_mth, "tops0", "Solar rad. at TOA", "W/m2", 917 . iim,jjm +1,nhori, 1,1,1, -99, 32,921 . iim,jjmp1,nhori, 1,1,1, -99, 32, 918 922 . "ave(X)", zsto,zout) 919 923 c 920 924 CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2", 921 . iim,jjm +1,nhori, 1,1,1, -99, 32,925 . iim,jjmp1,nhori, 1,1,1, -99, 32, 922 926 . "ave(X)", zsto,zout) 923 927 c 924 928 CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2", 925 . iim,jjm +1,nhori, 1,1,1, -99, 32,929 . iim,jjmp1,nhori, 1,1,1, -99, 32, 926 930 . "ave(X)", zsto,zout) 927 931 c 928 932 CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2", 929 . iim,jjm +1,nhori, 1,1,1, -99, 32,933 . iim,jjmp1,nhori, 1,1,1, -99, 32, 930 934 . "ave(X)", zsto,zout) 931 935 c 932 936 CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2", 933 . iim,jjm +1,nhori, 1,1,1, -99, 32,937 . iim,jjmp1,nhori, 1,1,1, -99, 32, 934 938 . "ave(X)", zsto,zout) 935 939 c 936 940 CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2", 937 . iim,jjm +1,nhori, 1,1,1, -99, 32,941 . iim,jjmp1,nhori, 1,1,1, -99, 32, 938 942 . "ave(X)", zsto,zout) 939 943 c 940 944 CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2", 941 . iim,jjm +1,nhori, 1,1,1, -99, 32,945 . iim,jjmp1,nhori, 1,1,1, -99, 32, 942 946 . "ave(X)", zsto,zout) 943 947 c 944 948 CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa", 945 . iim,jjm +1,nhori, 1,1,1, -99, 32,949 . iim,jjmp1,nhori, 1,1,1, -99, 32, 946 950 . "ave(X)", zsto,zout) 947 951 c 948 952 CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa", 949 . iim,jjm +1,nhori, 1,1,1, -99, 32,953 . iim,jjmp1,nhori, 1,1,1, -99, 32, 950 954 . "ave(X)", zsto,zout) 951 955 c 952 956 CALL histdef(nid_mth, "ruis", "Runoff", "mm/day", 953 . iim,jjm +1,nhori, 1,1,1, -99, 32,957 . iim,jjmp1,nhori, 1,1,1, -99, 32, 954 958 . "ave(X)", zsto,zout) 955 959 c 956 960 CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-", 957 . iim,jjm +1,nhori, 1,1,1, -99, 32,961 . iim,jjmp1,nhori, 1,1,1, -99, 32, 958 962 . "ave(X)", zsto,zout) 959 963 c 960 964 CALL histdef(nid_mth, "albs", "Surface albedo", "-", 961 . iim,jjm +1,nhori, 1,1,1, -99, 32,965 . iim,jjmp1,nhori, 1,1,1, -99, 32, 962 966 . "ave(X)", zsto,zout) 963 967 c 964 968 CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-", 965 . iim,jjm +1,nhori, 1,1,1, -99, 32,969 . iim,jjmp1,nhori, 1,1,1, -99, 32, 966 970 . "ave(X)", zsto,zout) 967 971 c 968 972 CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-", 969 . iim,jjm +1,nhori, 1,1,1, -99, 32,973 . iim,jjmp1,nhori, 1,1,1, -99, 32, 970 974 . "ave(X)", zsto,zout) 971 975 c 972 976 CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-", 973 . iim,jjm +1,nhori, 1,1,1, -99, 32,977 . iim,jjmp1,nhori, 1,1,1, -99, 32, 974 978 . "ave(X)", zsto,zout) 975 979 c 976 980 CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-", 977 . iim,jjm +1,nhori, 1,1,1, -99, 32,981 . iim,jjmp1,nhori, 1,1,1, -99, 32, 978 982 . "ave(X)", zsto,zout) 979 983 c 980 984 CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-", 981 . iim,jjm +1,nhori, 1,1,1, -99, 32,985 . iim,jjmp1,nhori, 1,1,1, -99, 32, 982 986 . "ave(X)", zsto,zout) 983 987 c 984 988 CALL histdef(nid_mth, "cldt", "Total cloudiness", "-", 985 . iim,jjm +1,nhori, 1,1,1, -99, 32,989 . iim,jjmp1,nhori, 1,1,1, -99, 32, 986 990 . "ave(X)", zsto,zout) 987 991 c 988 992 CALL histdef(nid_mth, "cldq", "Cloud liquid water path", "-", 989 . iim,jjm +1,nhori, 1,1,1, -99, 32,993 . iim,jjmp1,nhori, 1,1,1, -99, 32, 990 994 . "ave(X)", zsto,zout) 991 995 c 992 996 CALL histdef(nid_mth, "ue", "Zonal energy transport", "-", 993 . iim,jjm +1,nhori, 1,1,1, -99, 32,997 . iim,jjmp1,nhori, 1,1,1, -99, 32, 994 998 . "ave(X)", zsto,zout) 995 999 c 996 1000 CALL histdef(nid_mth, "ve", "Merid energy transport", "-", 997 . iim,jjm +1,nhori, 1,1,1, -99, 32,1001 . iim,jjmp1,nhori, 1,1,1, -99, 32, 998 1002 . "ave(X)", zsto,zout) 999 1003 c 1000 1004 CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-", 1001 . iim,jjm +1,nhori, 1,1,1, -99, 32,1005 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1002 1006 . "ave(X)", zsto,zout) 1003 1007 c 1004 1008 CALL histdef(nid_mth, "vq", "Merid humidity transport", "-", 1005 . iim,jjm +1,nhori, 1,1,1, -99, 32,1009 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1006 1010 . "ave(X)", zsto,zout) 1007 1011 c … … 1009 1013 c 1010 1014 CALL histdef(nid_mth, "temp", "Air temperature", "K", 1011 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1015 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1012 1016 . "ave(X)", zsto,zout) 1013 1017 c 1014 1018 CALL histdef(nid_mth, "ovap", "Specific humidity", "Kg/Kg", 1015 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1019 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1016 1020 . "ave(X)", zsto,zout) 1017 1021 c 1018 1022 CALL histdef(nid_mth, "geop", "Geopotential height", "m", 1019 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1023 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1020 1024 . "ave(X)", zsto,zout) 1021 1025 c 1022 1026 CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s", 1023 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1027 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1024 1028 . "ave(X)", zsto,zout) 1025 1029 c 1026 1030 CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s", 1027 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1031 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1028 1032 . "ave(X)", zsto,zout) 1029 1033 c 1030 1034 CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s", 1031 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1035 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1032 1036 . "ave(X)", zsto,zout) 1033 1037 c 1034 1038 CALL histdef(nid_mth, "pres", "Air pressure", "Pa", 1035 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1039 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1036 1040 . "ave(X)", zsto,zout) 1037 1041 c 1038 1042 CALL histdef(nid_mth, "rneb", "Cloud fraction", "-", 1039 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1043 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1040 1044 . "ave(X)", zsto,zout) 1041 1045 c 1042 1046 CALL histdef(nid_mth, "rhum", "Relative humidity", "-", 1043 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1047 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1044 1048 . "ave(X)", zsto,zout) 1045 1049 c 1046 1050 CALL histdef(nid_mth, "oliq", "Liquid water content", "kg/kg", 1047 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1051 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1048 1052 . "ave(X)", zsto,zout) 1049 1053 c 1050 1054 CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s", 1051 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1055 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1052 1056 . "ave(X)", zsto,zout) 1053 1057 c 1054 1058 CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "Kg/Kg/s", 1055 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1059 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1056 1060 . "ave(X)", zsto,zout) 1057 1061 c 1058 1062 CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s", 1059 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1063 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1060 1064 . "ave(X)", zsto,zout) 1061 1065 c 1062 1066 CALL histdef(nid_mth, "dqcon", "Convection dQ", "Kg/Kg/s", 1063 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1067 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1064 1068 . "ave(X)", zsto,zout) 1065 1069 c 1066 1070 CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s", 1067 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1071 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1068 1072 . "ave(X)", zsto,zout) 1069 1073 c 1070 1074 CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "Kg/Kg/s", 1071 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1075 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1072 1076 . "ave(X)", zsto,zout) 1073 1077 c 1074 1078 CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s", 1075 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1079 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1076 1080 . "ave(X)", zsto,zout) 1077 1081 c 1078 1082 CALL histdef(nid_mth, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s", 1079 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1083 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1080 1084 . "ave(X)", zsto,zout) 1081 1085 c 1082 1086 CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s", 1083 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1087 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1084 1088 . "ave(X)", zsto,zout) 1085 1089 c 1086 1090 CALL histdef(nid_mth, "dqeva", "Reevaporation dQ", "Kg/Kg/s", 1087 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1091 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1088 1092 . "ave(X)", zsto,zout) 1089 1093 c 1090 1094 CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s", 1091 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1095 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1092 1096 . "ave(X)", zsto,zout) 1093 1097 1094 1098 CALL histdef(nid_mth, "dqajs", "Dry adjust. dQ", "Kg/Kg/s", 1095 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1099 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1096 1100 . "ave(X)", zsto,zout) 1097 1101 c 1098 1102 CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s", 1099 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1103 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1100 1104 . "ave(X)", zsto,zout) 1101 1105 c 1102 1106 CALL histdef(nid_mth, "dtsw0", "SW radiation dT", "K/s", 1103 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1107 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1104 1108 . "ave(X)", zsto,zout) 1105 1109 c 1106 1110 CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s", 1107 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1111 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1108 1112 . "ave(X)", zsto,zout) 1109 1113 c 1110 1114 CALL histdef(nid_mth, "dtlw0", "LW radiation dT", "K/s", 1111 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1115 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1112 1116 . "ave(X)", zsto,zout) 1113 1117 c 1114 1118 CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2", 1115 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1119 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1116 1120 . "ave(X)", zsto,zout) 1117 1121 c 1118 1122 CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2", 1119 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1123 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1120 1124 . "ave(X)", zsto,zout) 1121 1125 c 1122 1126 IF (ok_orodr) THEN 1123 1127 CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2", 1124 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1128 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1125 1129 . "ave(X)", zsto,zout) 1126 1130 c 1127 1131 CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2", 1128 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1132 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1129 1133 . "ave(X)", zsto,zout) 1130 1134 c … … 1133 1137 IF (ok_orolf) THEN 1134 1138 CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2", 1135 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1139 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1136 1140 . "ave(X)", zsto,zout) 1137 1141 c 1138 1142 CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2", 1139 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1143 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1140 1144 . "ave(X)", zsto,zout) 1141 1145 ENDIF 1142 1146 C 1143 1147 CALL histdef(nid_mth, "ozone", "Ozone concentration", "-", 1144 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1148 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1145 1149 . "ave(X)", zsto,zout) 1146 1150 c … … 1150 1154 WRITE(str2,'(i2.2)') iq 1151 1155 CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-", 1152 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1156 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1153 1157 . "ave(X)", zsto,zout) 1154 1158 ELSE … … 1172 1176 zjulian = zjulian + day_ini 1173 1177 c 1174 CALL gr_fi_ecrit(1,klon,iim,jjm +1,rlon,zx_lon)1178 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) 1175 1179 DO i = 1, iim 1176 1180 zx_lon(i,1) = rlon(i+1) 1177 zx_lon(i,jjm +1) = rlon(i+1)1181 zx_lon(i,jjmp1) = rlon(i+1) 1178 1182 ENDDO 1179 1183 DO ll=1,klev 1180 1184 znivsig(ll)=float(ll) 1181 1185 ENDDO 1182 CALL gr_fi_ecrit(1,klon,iim,jjm +1,rlat,zx_lat)1183 CALL histbeg("histins", iim,zx_lon, jjm +1,zx_lat,1184 . 1,iim,1,jjm +1, 0, zjulian, dtime,1186 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 1187 CALL histbeg("histins", iim,zx_lon, jjmp1,zx_lat, 1188 . 1,iim,1,jjmp1, 0, zjulian, dtime, 1185 1189 . nhori, nid_ins) 1186 1190 c CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", … … 1189 1193 . klev, znivsig, nvert) 1190 1194 c 1195 c 1191 1196 zsto = dtime * ecrit_ins 1192 1197 zout = dtime * ecrit_ins 1193 1198 C 1194 1199 CALL histdef(nid_ins, "phis", "Surface geop. height", "-", 1195 . iim,jjm +1,nhori, 1,1,1, -99, 32,1200 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1196 1201 . "once", zsto,zout) 1197 1202 c 1198 1203 CALL histdef(nid_ins, "aire", "Grid area", "-", 1199 . iim,jjm +1,nhori, 1,1,1, -99, 32,1204 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1200 1205 . "once", zsto,zout) 1201 1206 c … … 1203 1208 c 1204 1209 CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa", 1205 . iim,jjm +1,nhori, 1,1,1, -99, 32,1210 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1206 1211 . "inst(X)", zsto,zout) 1207 1212 c 1208 1213 CALL histdef(nid_ins, "topl", "OLR", "W/m2", 1209 . iim,jjm +1,nhori, 1,1,1, -99, 32,1214 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1210 1215 . "inst(X)", zsto,zout) 1211 1216 c 1217 CALL histdef(nid_ins, "evap", "Evaporation", "mm/day", 1218 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1219 . "inst(X)", zsto,zout) 1220 c 1221 CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2", 1222 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1223 . "inst(X)", zsto,zout) 1224 c 1225 CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2", 1226 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1227 . "inst(X)", zsto,zout) 1228 c 1229 CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2", 1230 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1231 . "inst(X)", zsto,zout) 1232 c 1233 CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2", 1234 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1235 . "inst(X)", zsto,zout) 1236 c 1237 CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2", 1238 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1239 . "inst(X)", zsto,zout) 1240 c 1241 CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s", 1242 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1243 . "inst(X)", zsto,zout) 1244 c 1245 CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s", 1246 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1247 . "inst(X)", zsto,zout) 1248 c 1249 CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s", 1250 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1251 . "inst(X)", zsto,zout) 1252 c 1253 CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s", 1254 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1255 . "inst(X)", zsto,zout) 1256 c 1212 1257 c Champs 3D: 1213 1258 c 1214 1259 CALL histdef(nid_ins, "temp", "Temperature", "K", 1215 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1260 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1216 1261 . "inst(X)", zsto,zout) 1217 1262 c 1218 1263 CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s", 1219 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1264 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1220 1265 . "inst(X)", zsto,zout) 1221 1266 c 1222 1267 CALL histdef(nid_ins, "vitv", "Merid wind", "m/s", 1223 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1268 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1224 1269 . "inst(X)", zsto,zout) 1225 1270 c 1226 1271 CALL histdef(nid_ins, "geop", "Geopotential height", "m", 1227 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1272 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1228 1273 . "inst(X)", zsto,zout) 1229 1274 c 1230 1275 CALL histdef(nid_ins, "pres", "Air pressure", "Pa", 1231 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1276 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1232 1277 . "inst(X)", zsto,zout) 1233 1278 c 1279 CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s", 1280 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1281 . "inst(X)", zsto,zout) 1282 c 1283 CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s", 1284 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1285 . "inst(X)", zsto,zout) 1286 c 1287 1234 1288 CALL histend(nid_ins) 1235 1289 c … … 1368 1422 IF (ok_oasis .AND. MOD(itap-1,nexca).EQ.0) THEN 1369 1423 C 1370 CALL fromcpl(itap, (jjm+1)*iim,1424 CALL fromcpl(itap,jjmp1*iim, 1371 1425 . cpl_sst,cpl_sic,cpl_alb_sst,cpl_alb_sic) 1372 1426 DO i = 1, iim-1 ! un seul point pour le pole nord … … 1377 1431 ENDDO 1378 1432 DO i = 2, iim ! un seul point pour le pole sud 1379 cpl_sst(i,jjm +1) = cpl_sst(1,jjm+1)1380 cpl_sic(i,jjm +1) = cpl_sic(1,jjm+1)1381 cpl_alb_sst(i,jjm +1) = cpl_alb_sst(1,jjm+1)1382 cpl_alb_sic(i,jjm +1) = cpl_alb_sic(1,jjm+1)1433 cpl_sst(i,jjmp1) = cpl_sst(1,jjmp1) 1434 cpl_sic(i,jjmp1) = cpl_sic(1,jjmp1) 1435 cpl_alb_sst(i,jjmp1) = cpl_alb_sst(1,jjmp1) 1436 cpl_alb_sic(i,jjmp1) = cpl_alb_sic(1,jjmp1) 1383 1437 ENDDO 1384 1438 c … … 1407 1461 . pctsrf(ig,is_sic).GT.epsfra) THEN 1408 1462 pctsrf(ig,is_oce) = pctsrf(ig,is_oce) 1409 . - (cpl_sic(1,jjm +1)-pctsrf(ig,is_sic))1410 pctsrf(ig,is_sic) = cpl_sic(1,jjm +1)1411 lmt_sst(ig) = cpl_sst(1,jjm +1)1463 . - (cpl_sic(1,jjmp1)-pctsrf(ig,is_sic)) 1464 pctsrf(ig,is_sic) = cpl_sic(1,jjmp1) 1465 lmt_sst(ig) = cpl_sst(1,jjmp1) 1412 1466 ENDIF 1413 1467 c … … 1481 1535 fder(i) = dsens(i) + devap(i) 1482 1536 ENDDO 1537 1483 1538 DO k = 1, klev 1484 1539 DO i = 1, klon … … 1975 2030 c 1976 2031 i = NINT(zout/zsto) 1977 CALL gr_fi_ecrit(1,klon,iim,jjm +1,pphis,zx_tmp_2d)1978 CALL histwrite(nid_day,"phis",i,zx_tmp_2d,iim* (jjm+1),ndex2d)2032 CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) 2033 CALL histwrite(nid_day,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d) 1979 2034 c 1980 2035 i = NINT(zout/zsto) 1981 CALL gr_fi_ecrit(1,klon,iim,jjm +1,paire,zx_tmp_2d)1982 CALL histwrite(nid_day,"aire",i,zx_tmp_2d,iim* (jjm+1),ndex2d)2036 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) 2037 CALL histwrite(nid_day,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d) 1983 2038 C 1984 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zxtsol,zx_tmp_2d)1985 CALL histwrite(nid_day,"tsol",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2039 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 2040 CALL histwrite(nid_day,"tsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 1986 2041 c 1987 2042 DO i = 1, klon 1988 2043 zx_tmp_fi2d(i) = paprs(i,1) 1989 2044 ENDDO 1990 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zx_tmp_fi2d,zx_tmp_2d)1991 CALL histwrite(nid_day,"psol",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2045 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2046 CALL histwrite(nid_day,"psol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 1992 2047 c 1993 2048 DO i = 1, klon 1994 2049 zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) 1995 2050 ENDDO 1996 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zx_tmp_fi2d,zx_tmp_2d)1997 CALL histwrite(nid_day,"rain",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)1998 c 1999 CALL gr_fi_ecrit(1, klon,iim,jjm +1, snow_fall,zx_tmp_2d)2000 CALL histwrite(nid_day,"snow",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2001 c 2002 CALL gr_fi_ecrit(1, klon,iim,jjm +1, evap,zx_tmp_2d)2003 CALL histwrite(nid_day,"evap",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2004 c 2005 CALL gr_fi_ecrit(1, klon,iim,jjm +1, topsw,zx_tmp_2d)2006 CALL histwrite(nid_day,"tops",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2007 c 2008 CALL gr_fi_ecrit(1, klon,iim,jjm +1, toplw,zx_tmp_2d)2009 CALL histwrite(nid_day,"topl",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2010 c 2011 CALL gr_fi_ecrit(1, klon,iim,jjm +1, solsw,zx_tmp_2d)2012 CALL histwrite(nid_day,"sols",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2013 c 2014 CALL gr_fi_ecrit(1, klon,iim,jjm +1, sollw,zx_tmp_2d)2015 CALL histwrite(nid_day,"soll",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2016 c 2017 CALL gr_fi_ecrit(1, klon,iim,jjm +1, bils,zx_tmp_2d)2018 CALL histwrite(nid_day,"bils",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2019 c 2020 CALL gr_fi_ecrit(1, klon,iim,jjm +1, sens,zx_tmp_2d)2021 CALL histwrite(nid_day,"sens",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2022 c 2023 CALL gr_fi_ecrit(1, klon,iim,jjm +1, fder,zx_tmp_2d)2024 CALL histwrite(nid_day,"fder",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2025 c 2026 CALL gr_fi_ecrit(1, klon,iim,jjm +1, ruis,zx_tmp_2d)2027 CALL histwrite(nid_day,"ruis",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2051 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2052 CALL histwrite(nid_day,"rain",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2053 c 2054 CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) 2055 CALL histwrite(nid_day,"snow",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2056 c 2057 CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) 2058 CALL histwrite(nid_day,"evap",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2059 c 2060 CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) 2061 CALL histwrite(nid_day,"tops",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2062 c 2063 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) 2064 CALL histwrite(nid_day,"topl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2065 c 2066 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) 2067 CALL histwrite(nid_day,"sols",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2068 c 2069 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) 2070 CALL histwrite(nid_day,"soll",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2071 c 2072 CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) 2073 CALL histwrite(nid_day,"bils",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2074 c 2075 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) 2076 CALL histwrite(nid_day,"sens",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2077 c 2078 CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) 2079 CALL histwrite(nid_day,"fder",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2080 c 2081 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ruis,zx_tmp_2d) 2082 CALL histwrite(nid_day,"ruis",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2028 2083 c 2029 2084 DO i = 1, klon 2030 2085 zx_tmp_fi2d(i) = fluxu(i,1) 2031 2086 ENDDO 2032 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zx_tmp_fi2d,zx_tmp_2d)2033 CALL histwrite(nid_day,"frtu",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2087 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2088 CALL histwrite(nid_day,"frtu",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2034 2089 c 2035 2090 DO i = 1, klon 2036 2091 zx_tmp_fi2d(i) = fluxv(i,1) 2037 2092 ENDDO 2038 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zx_tmp_fi2d,zx_tmp_2d)2039 CALL histwrite(nid_day,"frtv",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2093 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2094 CALL histwrite(nid_day,"frtv",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2040 2095 c 2041 2096 DO i = 1, klon 2042 2097 zx_tmp_fi2d(i) = pctsrf(i,is_sic) 2043 2098 ENDDO 2044 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zx_tmp_fi2d,zx_tmp_2d)2045 CALL histwrite(nid_day,"sicf",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2046 c 2047 CALL gr_fi_ecrit(1, klon,iim,jjm +1, cldl,zx_tmp_2d)2048 CALL histwrite(nid_day,"cldl",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2049 c 2050 CALL gr_fi_ecrit(1, klon,iim,jjm +1, cldm,zx_tmp_2d)2051 CALL histwrite(nid_day,"cldm",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2052 c 2053 CALL gr_fi_ecrit(1, klon,iim,jjm +1, cldh,zx_tmp_2d)2054 CALL histwrite(nid_day,"cldh",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2055 c 2056 CALL gr_fi_ecrit(1, klon,iim,jjm +1, cldt,zx_tmp_2d)2057 CALL histwrite(nid_day,"cldt",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2058 c 2059 CALL gr_fi_ecrit(1, klon,iim,jjm +1, cldq,zx_tmp_2d)2060 CALL histwrite(nid_day,"cldq",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2099 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2100 CALL histwrite(nid_day,"sicf",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2101 c 2102 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) 2103 CALL histwrite(nid_day,"cldl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2104 c 2105 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) 2106 CALL histwrite(nid_day,"cldm",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2107 c 2108 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) 2109 CALL histwrite(nid_day,"cldh",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2110 c 2111 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d) 2112 CALL histwrite(nid_day,"cldt",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2113 c 2114 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) 2115 CALL histwrite(nid_day,"cldq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2061 2116 c 2062 2117 c Champs 3D: 2063 2118 c 2064 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, t_seri, zx_tmp_3d)2119 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) 2065 2120 CALL histwrite(nid_day,"temp",itap,zx_tmp_3d, 2066 . iim* (jjm+1)*klev,ndex3d)2067 c 2068 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, qx(1,1,ivap), zx_tmp_3d)2121 . iim*jjmp1*klev,ndex3d) 2122 c 2123 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) 2069 2124 CALL histwrite(nid_day,"ovap",itap,zx_tmp_3d, 2070 . iim* (jjm+1)*klev,ndex3d)2071 c 2072 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, zphi, zx_tmp_3d)2125 . iim*jjmp1*klev,ndex3d) 2126 c 2127 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) 2073 2128 CALL histwrite(nid_day,"geop",itap,zx_tmp_3d, 2074 . iim* (jjm+1)*klev,ndex3d)2075 c 2076 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, u_seri, zx_tmp_3d)2129 . iim*jjmp1*klev,ndex3d) 2130 c 2131 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 2077 2132 CALL histwrite(nid_day,"vitu",itap,zx_tmp_3d, 2078 . iim* (jjm+1)*klev,ndex3d)2079 c 2080 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, v_seri, zx_tmp_3d)2133 . iim*jjmp1*klev,ndex3d) 2134 c 2135 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 2081 2136 CALL histwrite(nid_day,"vitv",itap,zx_tmp_3d, 2082 . iim* (jjm+1)*klev,ndex3d)2083 c 2084 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, omega, zx_tmp_3d)2137 . iim*jjmp1*klev,ndex3d) 2138 c 2139 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) 2085 2140 CALL histwrite(nid_day,"vitw",itap,zx_tmp_3d, 2086 . iim* (jjm+1)*klev,ndex3d)2087 c 2088 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, pplay, zx_tmp_3d)2141 . iim*jjmp1*klev,ndex3d) 2142 c 2143 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) 2089 2144 CALL histwrite(nid_day,"pres",itap,zx_tmp_3d, 2090 . iim* (jjm+1)*klev,ndex3d)2145 . iim*jjmp1*klev,ndex3d) 2091 2146 c 2092 2147 if (ok_sync) then … … 2103 2158 c 2104 2159 i = NINT(zout/zsto) 2105 CALL gr_fi_ecrit(1,klon,iim,jjm +1,pphis,zx_tmp_2d)2106 CALL histwrite(nid_mth,"phis",i,zx_tmp_2d,iim* (jjm+1),ndex2d)2160 CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) 2161 CALL histwrite(nid_mth,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d) 2107 2162 C 2108 2163 i = NINT(zout/zsto) 2109 CALL gr_fi_ecrit(1,klon,iim,jjm +1,paire,zx_tmp_2d)2110 CALL histwrite(nid_mth,"aire",i,zx_tmp_2d,iim* (jjm+1),ndex2d)2111 2112 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zxtsol,zx_tmp_2d)2113 CALL histwrite(nid_mth,"tsol",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2164 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) 2165 CALL histwrite(nid_mth,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d) 2166 2167 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 2168 CALL histwrite(nid_mth,"tsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2114 2169 c 2115 2170 DO i = 1, klon 2116 2171 zx_tmp_fi2d(i) = paprs(i,1) 2117 2172 ENDDO 2118 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zx_tmp_fi2d,zx_tmp_2d)2119 CALL histwrite(nid_mth,"psol",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2120 c 2121 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zxqsol,zx_tmp_2d)2122 CALL histwrite(nid_mth,"qsol",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2173 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2174 CALL histwrite(nid_mth,"psol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2175 c 2176 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsol,zx_tmp_2d) 2177 CALL histwrite(nid_mth,"qsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2123 2178 c 2124 2179 DO i = 1, klon 2125 2180 zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) 2126 2181 ENDDO 2127 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zx_tmp_fi2d,zx_tmp_2d)2128 CALL histwrite(nid_mth,"rain",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2182 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2183 CALL histwrite(nid_mth,"rain",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2129 2184 c 2130 2185 DO i = 1, klon 2131 2186 zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i) 2132 2187 ENDDO 2133 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zx_tmp_fi2d,zx_tmp_2d)2134 CALL histwrite(nid_mth,"plul",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2188 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2189 CALL histwrite(nid_mth,"plul",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2135 2190 c 2136 2191 DO i = 1, klon 2137 2192 zx_tmp_fi2d(i) = rain_con(i) + snow_con(i) 2138 2193 ENDDO 2139 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zx_tmp_fi2d,zx_tmp_2d)2140 CALL histwrite(nid_mth,"pluc",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2141 c 2142 CALL gr_fi_ecrit(1, klon,iim,jjm +1, snow_fall,zx_tmp_2d)2143 CALL histwrite(nid_mth,"snow",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2144 c 2145 CALL gr_fi_ecrit(1, klon,iim,jjm +1, agesno,zx_tmp_2d)2146 CALL histwrite(nid_mth,"ages",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2147 c 2148 CALL gr_fi_ecrit(1, klon,iim,jjm +1, evap,zx_tmp_2d)2149 CALL histwrite(nid_mth,"evap",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2150 c 2151 CALL gr_fi_ecrit(1, klon,iim,jjm +1, topsw,zx_tmp_2d)2152 CALL histwrite(nid_mth,"tops",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2153 c 2154 CALL gr_fi_ecrit(1, klon,iim,jjm +1, toplw,zx_tmp_2d)2155 CALL histwrite(nid_mth,"topl",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2156 c 2157 CALL gr_fi_ecrit(1, klon,iim,jjm +1, solsw,zx_tmp_2d)2158 CALL histwrite(nid_mth,"sols",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2159 c 2160 CALL gr_fi_ecrit(1, klon,iim,jjm +1, sollw,zx_tmp_2d)2161 CALL histwrite(nid_mth,"soll",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2162 c 2163 CALL gr_fi_ecrit(1, klon,iim,jjm +1, topsw0,zx_tmp_2d)2164 CALL histwrite(nid_mth,"tops0",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2165 c 2166 CALL gr_fi_ecrit(1, klon,iim,jjm +1, toplw0,zx_tmp_2d)2167 CALL histwrite(nid_mth,"topl0",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2168 c 2169 CALL gr_fi_ecrit(1, klon,iim,jjm +1, solsw0,zx_tmp_2d)2170 CALL histwrite(nid_mth,"sols0",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2171 c 2172 CALL gr_fi_ecrit(1, klon,iim,jjm +1, sollw0,zx_tmp_2d)2173 CALL histwrite(nid_mth,"soll0",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2174 c 2175 CALL gr_fi_ecrit(1, klon,iim,jjm +1, bils,zx_tmp_2d)2176 CALL histwrite(nid_mth,"bils",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2177 c 2178 CALL gr_fi_ecrit(1, klon,iim,jjm +1, sens,zx_tmp_2d)2179 CALL histwrite(nid_mth,"sens",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2180 c 2181 CALL gr_fi_ecrit(1, klon,iim,jjm +1, fder,zx_tmp_2d)2182 CALL histwrite(nid_mth,"fder",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2183 c 2184 CALL gr_fi_ecrit(1, klon,iim,jjm +1, ruis,zx_tmp_2d)2185 CALL histwrite(nid_mth,"ruis",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2194 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2195 CALL histwrite(nid_mth,"pluc",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2196 c 2197 CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) 2198 CALL histwrite(nid_mth,"snow",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2199 c 2200 CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d) 2201 CALL histwrite(nid_mth,"ages",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2202 c 2203 CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) 2204 CALL histwrite(nid_mth,"evap",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2205 c 2206 CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) 2207 CALL histwrite(nid_mth,"tops",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2208 c 2209 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) 2210 CALL histwrite(nid_mth,"topl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2211 c 2212 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) 2213 CALL histwrite(nid_mth,"sols",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2214 c 2215 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) 2216 CALL histwrite(nid_mth,"soll",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2217 c 2218 CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d) 2219 CALL histwrite(nid_mth,"tops0",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2220 c 2221 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d) 2222 CALL histwrite(nid_mth,"topl0",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2223 c 2224 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d) 2225 CALL histwrite(nid_mth,"sols0",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2226 c 2227 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d) 2228 CALL histwrite(nid_mth,"soll0",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2229 c 2230 CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) 2231 CALL histwrite(nid_mth,"bils",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2232 c 2233 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) 2234 CALL histwrite(nid_mth,"sens",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2235 c 2236 CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) 2237 CALL histwrite(nid_mth,"fder",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2238 c 2239 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ruis,zx_tmp_2d) 2240 CALL histwrite(nid_mth,"ruis",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2186 2241 c 2187 2242 DO i = 1, klon 2188 2243 zx_tmp_fi2d(i) = fluxu(i,1) 2189 2244 ENDDO 2190 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zx_tmp_fi2d,zx_tmp_2d)2191 CALL histwrite(nid_mth,"frtu",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2245 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2246 CALL histwrite(nid_mth,"frtu",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2192 2247 c 2193 2248 DO i = 1, klon 2194 2249 zx_tmp_fi2d(i) = fluxv(i,1) 2195 2250 ENDDO 2196 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zx_tmp_fi2d,zx_tmp_2d)2197 CALL histwrite(nid_mth,"frtv",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2251 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2252 CALL histwrite(nid_mth,"frtv",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2198 2253 c 2199 2254 DO i = 1, klon 2200 2255 zx_tmp_fi2d(i) = pctsrf(i,is_sic) 2201 2256 ENDDO 2202 CALL gr_fi_ecrit(1, klon,iim,jjm +1, zx_tmp_fi2d,zx_tmp_2d)2203 CALL histwrite(nid_mth,"sicf",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2204 c 2205 CALL gr_fi_ecrit(1, klon,iim,jjm +1, albsol,zx_tmp_2d)2206 CALL histwrite(nid_mth,"albs",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2207 c 2208 CALL gr_fi_ecrit(1, klon,iim,jjm +1, cdragm,zx_tmp_2d)2209 CALL histwrite(nid_mth,"cdrm",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2210 c 2211 CALL gr_fi_ecrit(1, klon,iim,jjm +1, cdragh,zx_tmp_2d)2212 CALL histwrite(nid_mth,"cdrh",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2213 c 2214 CALL gr_fi_ecrit(1, klon,iim,jjm +1, cldl,zx_tmp_2d)2215 CALL histwrite(nid_mth,"cldl",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2216 c 2217 CALL gr_fi_ecrit(1, klon,iim,jjm +1, cldm,zx_tmp_2d)2218 CALL histwrite(nid_mth,"cldm",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2219 c 2220 CALL gr_fi_ecrit(1, klon,iim,jjm +1, cldh,zx_tmp_2d)2221 CALL histwrite(nid_mth,"cldh",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2222 c 2223 CALL gr_fi_ecrit(1, klon,iim,jjm +1, cldt,zx_tmp_2d)2224 CALL histwrite(nid_mth,"cldt",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2225 c 2226 CALL gr_fi_ecrit(1, klon,iim,jjm +1, cldq,zx_tmp_2d)2227 CALL histwrite(nid_mth,"cldq",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2228 c 2229 CALL gr_fi_ecrit(1, klon,iim,jjm +1, ue,zx_tmp_2d)2230 CALL histwrite(nid_mth,"ue",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2231 c 2232 CALL gr_fi_ecrit(1, klon,iim,jjm +1, ve,zx_tmp_2d)2233 CALL histwrite(nid_mth,"ve",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2234 c 2235 CALL gr_fi_ecrit(1, klon,iim,jjm +1, uq,zx_tmp_2d)2236 CALL histwrite(nid_mth,"uq",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2237 c 2238 CALL gr_fi_ecrit(1, klon,iim,jjm +1, vq,zx_tmp_2d)2239 CALL histwrite(nid_mth,"vq",itap,zx_tmp_2d,iim* (jjm+1),ndex2d)2257 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2258 CALL histwrite(nid_mth,"sicf",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2259 c 2260 CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d) 2261 CALL histwrite(nid_mth,"albs",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2262 c 2263 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d) 2264 CALL histwrite(nid_mth,"cdrm",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2265 c 2266 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d) 2267 CALL histwrite(nid_mth,"cdrh",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2268 c 2269 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) 2270 CALL histwrite(nid_mth,"cldl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2271 c 2272 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) 2273 CALL histwrite(nid_mth,"cldm",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2274 c 2275 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) 2276 CALL histwrite(nid_mth,"cldh",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2277 c 2278 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d) 2279 CALL histwrite(nid_mth,"cldt",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2280 c 2281 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) 2282 CALL histwrite(nid_mth,"cldq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2283 c 2284 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d) 2285 CALL histwrite(nid_mth,"ue",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2286 c 2287 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d) 2288 CALL histwrite(nid_mth,"ve",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2289 c 2290 CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d) 2291 CALL histwrite(nid_mth,"uq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2292 c 2293 CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d) 2294 CALL histwrite(nid_mth,"vq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2240 2295 c 2241 2296 c Champs 3D: 2242 2297 C 2243 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, t_seri, zx_tmp_3d)2298 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) 2244 2299 CALL histwrite(nid_mth,"temp",itap,zx_tmp_3d, 2245 . iim* (jjm+1)*klev,ndex3d)2246 c 2247 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, qx(1,1,ivap), zx_tmp_3d)2300 . iim*jjmp1*klev,ndex3d) 2301 c 2302 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) 2248 2303 CALL histwrite(nid_mth,"ovap",itap,zx_tmp_3d, 2249 . iim* (jjm+1)*klev,ndex3d)2250 c 2251 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, zphi, zx_tmp_3d)2304 . iim*jjmp1*klev,ndex3d) 2305 c 2306 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) 2252 2307 CALL histwrite(nid_mth,"geop",itap,zx_tmp_3d, 2253 . iim* (jjm+1)*klev,ndex3d)2254 c 2255 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, u_seri, zx_tmp_3d)2308 . iim*jjmp1*klev,ndex3d) 2309 c 2310 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 2256 2311 CALL histwrite(nid_mth,"vitu",itap,zx_tmp_3d, 2257 . iim* (jjm+1)*klev,ndex3d)2258 c 2259 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, v_seri, zx_tmp_3d)2312 . iim*jjmp1*klev,ndex3d) 2313 c 2314 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 2260 2315 CALL histwrite(nid_mth,"vitv",itap,zx_tmp_3d, 2261 . iim* (jjm+1)*klev,ndex3d)2262 c 2263 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, omega, zx_tmp_3d)2316 . iim*jjmp1*klev,ndex3d) 2317 c 2318 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) 2264 2319 CALL histwrite(nid_mth,"vitw",itap,zx_tmp_3d, 2265 . iim* (jjm+1)*klev,ndex3d)2266 c 2267 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, pplay, zx_tmp_3d)2320 . iim*jjmp1*klev,ndex3d) 2321 c 2322 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) 2268 2323 CALL histwrite(nid_mth,"pres",itap,zx_tmp_3d, 2269 . iim* (jjm+1)*klev,ndex3d)2270 c 2271 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, cldfra, zx_tmp_3d)2324 . iim*jjmp1*klev,ndex3d) 2325 c 2326 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d) 2272 2327 CALL histwrite(nid_mth,"rneb",itap,zx_tmp_3d, 2273 . iim* (jjm+1)*klev,ndex3d)2274 c 2275 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, zx_rh, zx_tmp_3d)2328 . iim*jjmp1*klev,ndex3d) 2329 c 2330 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d) 2276 2331 CALL histwrite(nid_mth,"rhum",itap,zx_tmp_3d, 2277 . iim* (jjm+1)*klev,ndex3d)2278 c 2279 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, cldliq, zx_tmp_3d)2332 . iim*jjmp1*klev,ndex3d) 2333 c 2334 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d) 2280 2335 CALL histwrite(nid_mth,"oliq",itap,zx_tmp_3d, 2281 . iim* (jjm+1)*klev,ndex3d)2282 c 2283 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_t_dyn, zx_tmp_3d)2336 . iim*jjmp1*klev,ndex3d) 2337 c 2338 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d) 2284 2339 CALL histwrite(nid_mth,"dtdyn",itap,zx_tmp_3d, 2285 . iim* (jjm+1)*klev,ndex3d)2286 c 2287 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_q_dyn, zx_tmp_3d)2340 . iim*jjmp1*klev,ndex3d) 2341 c 2342 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d) 2288 2343 CALL histwrite(nid_mth,"dqdyn",itap,zx_tmp_3d, 2289 . iim* (jjm+1)*klev,ndex3d)2290 c 2291 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_t_con, zx_tmp_3d)2344 . iim*jjmp1*klev,ndex3d) 2345 c 2346 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_con, zx_tmp_3d) 2292 2347 CALL histwrite(nid_mth,"dtcon",itap,zx_tmp_3d, 2293 . iim* (jjm+1)*klev,ndex3d)2294 c 2295 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_q_con, zx_tmp_3d)2348 . iim*jjmp1*klev,ndex3d) 2349 c 2350 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d) 2296 2351 CALL histwrite(nid_mth,"dqcon",itap,zx_tmp_3d, 2297 . iim* (jjm+1)*klev,ndex3d)2298 c 2299 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_t_lsc, zx_tmp_3d)2352 . iim*jjmp1*klev,ndex3d) 2353 c 2354 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_lsc, zx_tmp_3d) 2300 2355 CALL histwrite(nid_mth,"dtlsc",itap,zx_tmp_3d, 2301 . iim* (jjm+1)*klev,ndex3d)2302 c 2303 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_q_lsc, zx_tmp_3d)2356 . iim*jjmp1*klev,ndex3d) 2357 c 2358 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d) 2304 2359 CALL histwrite(nid_mth,"dqlsc",itap,zx_tmp_3d, 2305 . iim* (jjm+1)*klev,ndex3d)2306 c 2307 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_t_vdf, zx_tmp_3d)2360 . iim*jjmp1*klev,ndex3d) 2361 c 2362 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d) 2308 2363 CALL histwrite(nid_mth,"dtvdf",itap,zx_tmp_3d, 2309 . iim* (jjm+1)*klev,ndex3d)2310 c 2311 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_q_vdf, zx_tmp_3d)2364 . iim*jjmp1*klev,ndex3d) 2365 c 2366 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d) 2312 2367 CALL histwrite(nid_mth,"dqvdf",itap,zx_tmp_3d, 2313 . iim* (jjm+1)*klev,ndex3d)2314 c 2315 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_t_eva, zx_tmp_3d)2368 . iim*jjmp1*klev,ndex3d) 2369 c 2370 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_eva, zx_tmp_3d) 2316 2371 CALL histwrite(nid_mth,"dteva",itap,zx_tmp_3d, 2317 . iim* (jjm+1)*klev,ndex3d)2318 c 2319 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_q_eva, zx_tmp_3d)2372 . iim*jjmp1*klev,ndex3d) 2373 c 2374 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d) 2320 2375 CALL histwrite(nid_mth,"dqeva",itap,zx_tmp_3d, 2321 . iim* (jjm+1)*klev,ndex3d)2322 c 2323 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_t_ajs, zx_tmp_3d)2376 . iim*jjmp1*klev,ndex3d) 2377 c 2378 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ajs, zx_tmp_3d) 2324 2379 CALL histwrite(nid_mth,"dtajs",itap,zx_tmp_3d, 2325 . iim* (jjm+1)*klev,ndex3d)2326 c 2327 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_q_ajs, zx_tmp_3d)2380 . iim*jjmp1*klev,ndex3d) 2381 c 2382 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d) 2328 2383 CALL histwrite(nid_mth,"dqajs",itap,zx_tmp_3d, 2329 . iim* (jjm+1)*klev,ndex3d)2330 c 2331 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, heat, zx_tmp_3d)2384 . iim*jjmp1*klev,ndex3d) 2385 c 2386 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat, zx_tmp_3d) 2332 2387 CALL histwrite(nid_mth,"dtswr",itap,zx_tmp_3d, 2333 . iim* (jjm+1)*klev,ndex3d)2334 c 2335 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, heat0, zx_tmp_3d)2388 . iim*jjmp1*klev,ndex3d) 2389 c 2390 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat0, zx_tmp_3d) 2336 2391 CALL histwrite(nid_mth,"dtsw0",itap,zx_tmp_3d, 2337 . iim* (jjm+1)*klev,ndex3d)2338 c 2339 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, cool, zx_tmp_3d)2392 . iim*jjmp1*klev,ndex3d) 2393 c 2394 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool, zx_tmp_3d) 2340 2395 CALL histwrite(nid_mth,"dtlwr",itap,zx_tmp_3d, 2341 . iim* (jjm+1)*klev,ndex3d)2342 c 2343 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, cool0, zx_tmp_3d)2396 . iim*jjmp1*klev,ndex3d) 2397 c 2398 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool0, zx_tmp_3d) 2344 2399 CALL histwrite(nid_mth,"dtlw0",itap,zx_tmp_3d, 2345 . iim* (jjm+1)*klev,ndex3d)2346 c 2347 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_u_vdf, zx_tmp_3d)2400 . iim*jjmp1*klev,ndex3d) 2401 c 2402 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d) 2348 2403 CALL histwrite(nid_mth,"duvdf",itap,zx_tmp_3d, 2349 . iim* (jjm+1)*klev,ndex3d)2350 c 2351 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_v_vdf, zx_tmp_3d)2404 . iim*jjmp1*klev,ndex3d) 2405 c 2406 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d) 2352 2407 CALL histwrite(nid_mth,"dvvdf",itap,zx_tmp_3d, 2353 . iim* (jjm+1)*klev,ndex3d)2408 . iim*jjmp1*klev,ndex3d) 2354 2409 c 2355 2410 IF (ok_orodr) THEN 2356 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_u_oro, zx_tmp_3d)2411 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d) 2357 2412 CALL histwrite(nid_mth,"duoro",itap,zx_tmp_3d, 2358 . iim* (jjm+1)*klev,ndex3d)2359 c 2360 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_v_oro, zx_tmp_3d)2413 . iim*jjmp1*klev,ndex3d) 2414 c 2415 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oro, zx_tmp_3d) 2361 2416 CALL histwrite(nid_mth,"dvoro",itap,zx_tmp_3d, 2362 . iim* (jjm+1)*klev,ndex3d)2417 . iim*jjmp1*klev,ndex3d) 2363 2418 c 2364 2419 ENDIF 2365 2420 C 2366 2421 IF (ok_orolf) THEN 2367 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_u_lif, zx_tmp_3d)2422 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_lif, zx_tmp_3d) 2368 2423 CALL histwrite(nid_mth,"dulif",itap,zx_tmp_3d, 2369 . iim* (jjm+1)*klev,ndex3d)2370 c 2371 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_v_lif, zx_tmp_3d)2424 . iim*jjmp1*klev,ndex3d) 2425 c 2426 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_lif, zx_tmp_3d) 2372 2427 CALL histwrite(nid_mth,"dvlif",itap,zx_tmp_3d, 2373 . iim* (jjm+1)*klev,ndex3d)2428 . iim*jjmp1*klev,ndex3d) 2374 2429 ENDIF 2375 2430 C 2376 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, wo, zx_tmp_3d)2431 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d) 2377 2432 CALL histwrite(nid_mth,"ozone",itap,zx_tmp_3d, 2378 . iim* (jjm+1)*klev,ndex3d)2433 . iim*jjmp1*klev,ndex3d) 2379 2434 c 2380 2435 IF (nqmax.GE.3) THEN 2381 2436 DO iq=1,nqmax-2 2382 2437 IF (iq.LE.99) THEN 2383 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, qx(1,1,iq+2), zx_tmp_3d)2438 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d) 2384 2439 WRITE(str2,'(i2.2)') iq 2385 2440 CALL histwrite(nid_mth,"trac"//str2,itap,zx_tmp_3d, 2386 . iim* (jjm+1)*klev,ndex3d)2441 . iim*jjmp1*klev,ndex3d) 2387 2442 ELSE 2388 2443 PRINT*, "Trop de traceurs" … … 2405 2460 c 2406 2461 i = NINT(zout/zsto) 2407 CALL gr_fi_ecrit(1,klon,iim,jjm +1,pphis,zx_tmp_2d)2408 CALL histwrite(nid_ins,"phis",i,zx_tmp_2d,iim* (jjm+1),ndex2d)2462 CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) 2463 CALL histwrite(nid_ins,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d) 2409 2464 c 2410 2465 i = NINT(zout/zsto) 2411 CALL gr_fi_ecrit(1,klon,iim,jjm +1,paire,zx_tmp_2d)2412 CALL histwrite(nid_ins,"aire",i,zx_tmp_2d,iim* (jjm+1),ndex2d)2466 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) 2467 CALL histwrite(nid_ins,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d) 2413 2468 2414 2469 DO i = 1, klon 2415 2470 zx_tmp_fi2d(i) = paprs(i,1) 2416 2471 ENDDO 2417 CALL gr_fi_ecrit(1, klon,iim,jjm+1, zx_tmp_fi2d,zx_tmp_2d) 2418 CALL histwrite(nid_ins,"psol",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) 2419 c 2420 CALL gr_fi_ecrit(1, klon,iim,jjm+1, toplw,zx_tmp_2d) 2421 CALL histwrite(nid_ins,"topl",itap,zx_tmp_2d,iim*(jjm+1),ndex2d) 2472 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2473 CALL histwrite(nid_ins,"psol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2474 c 2475 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) 2476 CALL histwrite(nid_ins,"topl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2477 c 2478 CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) 2479 CALL histwrite(nid_ins,"evap",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2480 c 2481 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) 2482 CALL histwrite(nid_ins,"sols",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2483 c 2484 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) 2485 CALL histwrite(nid_ins,"soll",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2486 c 2487 CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) 2488 CALL histwrite(nid_ins,"bils",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2489 c 2490 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) 2491 CALL histwrite(nid_ins,"sens",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2492 c 2493 CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) 2494 CALL histwrite(nid_ins,"fder",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2495 c 2496 CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_oce),zx_tmp_2d) 2497 CALL histwrite(nid_ins,"dtsvdfo",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2498 c 2499 CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_ter),zx_tmp_2d) 2500 CALL histwrite(nid_ins,"dtsvdft",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2501 c 2502 CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_lic),zx_tmp_2d) 2503 CALL histwrite(nid_ins,"dtsvdfg",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2504 c 2505 CALL gr_fi_ecrit(1, klon,iim,jjmp1, d_ts(1,is_sic),zx_tmp_2d) 2506 CALL histwrite(nid_ins,"dtsvdfi",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2507 2422 2508 c 2423 2509 c Champs 3D: 2424 2510 c 2425 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, t_seri, zx_tmp_3d)2511 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) 2426 2512 CALL histwrite(nid_ins,"temp",itap,zx_tmp_3d, 2427 . iim* (jjm+1)*klev,ndex3d)2428 c 2429 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, u_seri, zx_tmp_3d)2513 . iim*jjmp1*klev,ndex3d) 2514 c 2515 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 2430 2516 CALL histwrite(nid_ins,"vitu",itap,zx_tmp_3d, 2431 . iim* (jjm+1)*klev,ndex3d)2432 c 2433 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, v_seri, zx_tmp_3d)2517 . iim*jjmp1*klev,ndex3d) 2518 c 2519 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 2434 2520 CALL histwrite(nid_ins,"vitv",itap,zx_tmp_3d, 2435 . iim* (jjm+1)*klev,ndex3d)2521 . iim*jjmp1*klev,ndex3d) 2436 2522 c 2437 2523 CALL gr_fi_ecrit(klev,klon,iim,jjm+1, zphi, zx_tmp_3d) 2438 2524 CALL histwrite(nid_ins,"geop",itap,zx_tmp_3d, 2439 . iim* (jjm+1)*klev,ndex3d)2525 . iim*jjmp1*klev,ndex3d) 2440 2526 c 2441 2527 CALL gr_fi_ecrit(klev,klon,iim,jjm+1, pplay, zx_tmp_3d) 2442 2528 CALL histwrite(nid_ins,"pres",itap,zx_tmp_3d, 2443 . iim*(jjm+1)*klev,ndex3d) 2529 . iim*jjmp1*klev,ndex3d) 2530 c 2531 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d) 2532 CALL histwrite(nid_ins,"dtvdf",itap,zx_tmp_3d, 2533 . iim*jjmp1*klev,ndex3d) 2534 c 2535 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d) 2536 CALL histwrite(nid_ins,"dqvdf",itap,zx_tmp_3d, 2537 . iim*jjmp1*klev,ndex3d) 2538 2444 2539 c 2445 2540 if (ok_sync) then … … 2490 2585 ig = ig + 1 2491 2586 DO i = 1, iim 2492 z_sols(i,jjm +1) = oas_sols(ig)2493 z_nsol(i,jjm +1) = oas_nsol(ig)2494 z_rain(i,jjm +1) = oas_rain(ig)2495 z_snow(i,jjm +1) = oas_snow(ig)2496 z_evap(i,jjm +1) = oas_evap(ig)2497 z_ruisoce(i,jjm +1) = oas_ruisoce(ig)2498 z_ruisriv(i,jjm +1) = oas_ruisriv(ig)2499 z_tsol(i,jjm +1) = oas_tsol(ig)2500 z_fder(i,jjm +1) = oas_fder(ig)2501 z_albe(i,jjm +1) = oas_albe(ig)2502 z_taux(i,jjm +1) = oas_taux(ig)2503 z_tauy(i,jjm +1) = oas_tauy(ig)2587 z_sols(i,jjmp1) = oas_sols(ig) 2588 z_nsol(i,jjmp1) = oas_nsol(ig) 2589 z_rain(i,jjmp1) = oas_rain(ig) 2590 z_snow(i,jjmp1) = oas_snow(ig) 2591 z_evap(i,jjmp1) = oas_evap(ig) 2592 z_ruisoce(i,jjmp1) = oas_ruisoce(ig) 2593 z_ruisriv(i,jjmp1) = oas_ruisriv(ig) 2594 z_tsol(i,jjmp1) = oas_tsol(ig) 2595 z_fder(i,jjmp1) = oas_fder(ig) 2596 z_albe(i,jjmp1) = oas_albe(ig) 2597 z_taux(i,jjmp1) = oas_taux(ig) 2598 z_tauy(i,jjmp1) = oas_tauy(ig) 2504 2599 ENDDO 2505 2600 c 2506 2601 c Passer les champs au coupleur: 2507 2602 c 2508 CALL intocpl(itap, (jjm+1)*iim,2603 CALL intocpl(itap,jjmp1*iim, 2509 2604 . z_sols, z_nsol, 2510 2605 . z_rain, z_snow, z_evap,
Note: See TracChangeset
for help on using the changeset viewer.
