Changeset 72 for LMDZ.3.3/trunk
- Timestamp:
- Mar 7, 2000, 4:54:11 PM (25 years ago)
- Location:
- LMDZ.3.3/trunk/libf/phylmd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/trunk/libf/phylmd/dimphy.h
r2 r72 1 1 c----------------------------------------------------------------------- 2 2 INTEGER KIDIA, KFDIA, KLON, KLEV 3 PARAMETER (KIDIA=1,KFDIA=iim*(jjm-1)+2 ,3 PARAMETER (KIDIA=1,KFDIA=iim*(jjm-1)+2-1/jjm, 4 4 . KLON=KFDIA-KIDIA+1,KLEV=llm) 5 5 c----------------------------------------------------------------------- -
LMDZ.3.3/trunk/libf/phylmd/physiq.F
r53 r72 72 72 PARAMETER (itimestep=1800) 73 73 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)74 REAL cpl_sst(iim,jjmp1), cpl_sic(iim,jjmp1) 75 REAL cpl_alb_sst(iim,jjmp1), cpl_alb_sic(iim,jjmp1) 76 76 c====================================================================== 77 77 c ok_ocean indique l'utilisation du modele oceanique "slab ocean", … … 421 421 c 422 422 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)423 REAL zx_tmp_x(iim), zx_tmp_yjjmp1 424 REAL zx_relief(iim,jjmp1) 425 REAL zx_aire(iim,jjmp1) 426 426 c 427 427 c Variables du changement … … 473 473 SAVE ecrit_reg ! frequence d'ecriture 474 474 c 475 REAL oas_sols(klon), z_sols(iim,jjm +1)475 REAL oas_sols(klon), z_sols(iim,jjmp1) 476 476 SAVE oas_sols 477 REAL oas_nsol(klon), z_nsol(iim,jjm +1)477 REAL oas_nsol(klon), z_nsol(iim,jjmp1) 478 478 SAVE oas_nsol 479 REAL oas_rain(klon), z_rain(iim,jjm +1)479 REAL oas_rain(klon), z_rain(iim,jjmp1) 480 480 SAVE oas_rain 481 REAL oas_snow(klon), z_snow(iim,jjm +1)481 REAL oas_snow(klon), z_snow(iim,jjmp1) 482 482 SAVE oas_snow 483 REAL oas_evap(klon), z_evap(iim,jjm +1)483 REAL oas_evap(klon), z_evap(iim,jjmp1) 484 484 SAVE oas_evap 485 REAL oas_ruis(klon), z_ruis(iim,jjm +1)485 REAL oas_ruis(klon), z_ruis(iim,jjmp1) 486 486 SAVE oas_ruis 487 REAL oas_tsol(klon), z_tsol(iim,jjm +1)487 REAL oas_tsol(klon), z_tsol(iim,jjmp1) 488 488 SAVE oas_tsol 489 REAL oas_fder(klon), z_fder(iim,jjm +1)489 REAL oas_fder(klon), z_fder(iim,jjmp1) 490 490 SAVE oas_fder 491 REAL oas_albe(klon), z_albe(iim,jjm +1)491 REAL oas_albe(klon), z_albe(iim,jjmp1) 492 492 SAVE oas_albe 493 REAL oas_taux(klon), z_taux(iim,jjm +1)493 REAL oas_taux(klon), z_taux(iim,jjmp1) 494 494 SAVE oas_taux 495 REAL oas_tauy(klon), z_tauy(iim,jjm +1)495 REAL oas_tauy(klon), z_tauy(iim,jjmp1) 496 496 SAVE oas_tauy 497 REAL oas_ruisoce(klon), z_ruisoce(iim,jjm +1)497 REAL oas_ruisoce(klon), z_ruisoce(iim,jjmp1) 498 498 SAVE oas_ruisoce 499 REAL oas_ruisriv(klon), z_ruisriv(iim,jjm +1)499 REAL oas_ruisriv(klon), z_ruisriv(iim,jjmp1) 500 500 SAVE oas_ruisriv 501 501 c … … 518 518 REAL tabcntr0( length ) 519 519 c 520 INTEGER ndex2d(iim* (jjm+1)),ndex3d(iim*(jjm+1)*klev)520 INTEGER ndex2d(iim*jjmp1),ndex3d(iim*jjmp1*klev) 521 521 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)522 REAL zx_tmp_2d(iim,jjmp1), zx_tmp_3d(iim,jjmp1,klev) 523 REAL zx_lon(iim,jjmp1), zx_lat(iim,jjmp1) 524 524 c 525 525 INTEGER nid_day, nid_mth, nid_ins … … 676 676 zjulian = zjulian + day_ini 677 677 c 678 CALL gr_fi_ecrit(1,klon,iim,jjm +1,rlon,zx_lon)678 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) 679 679 DO i = 1, iim 680 680 zx_lon(i,1) = rlon(i+1) 681 zx_lon(i,jjm +1) = rlon(i+1)681 zx_lon(i,jjmp1) = rlon(i+1) 682 682 ENDDO 683 683 DO ll=1,klev 684 684 znivsig(ll)=float(ll) 685 685 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,686 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 687 CALL histbeg("histday", iim,zx_lon, jjmp1,zx_lat, 688 . 1,iim,1,jjmp1, 0, zjulian, dtime, 689 689 . nhori, nid_day) 690 690 c CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", … … 697 697 c 698 698 CALL histdef(nid_day, "phis", "Surface geop. height", "-", 699 . iim,jjm +1,nhori, 1,1,1, -99, 32,699 . iim,jjmp1,nhori, 1,1,1, -99, 32, 700 700 . "once", zsto,zout) 701 701 c 702 702 CALL histdef(nid_day, "aire", "Grid area", "-", 703 . iim,jjm +1,nhori, 1,1,1, -99, 32,703 . iim,jjmp1,nhori, 1,1,1, -99, 32, 704 704 . "once", zsto,zout) 705 705 c … … 707 707 c 708 708 CALL histdef(nid_day, "tsol", "Surface Temperature", "K", 709 . iim,jjm +1,nhori, 1,1,1, -99, 32,709 . iim,jjmp1,nhori, 1,1,1, -99, 32, 710 710 . "ave(X)", zsto,zout) 711 711 c 712 712 CALL histdef(nid_day, "psol", "Surface Pressure", "Pa", 713 . iim,jjm +1,nhori, 1,1,1, -99, 32,713 . iim,jjmp1,nhori, 1,1,1, -99, 32, 714 714 . "ave(X)", zsto,zout) 715 715 c 716 716 CALL histdef(nid_day, "rain", "Precipitation", "mm/day", 717 . iim,jjm +1,nhori, 1,1,1, -99, 32,717 . iim,jjmp1,nhori, 1,1,1, -99, 32, 718 718 . "ave(X)", zsto,zout) 719 719 c 720 720 CALL histdef(nid_day, "snow", "Snow fall", "mm/day", 721 . iim,jjm +1,nhori, 1,1,1, -99, 32,721 . iim,jjmp1,nhori, 1,1,1, -99, 32, 722 722 . "ave(X)", zsto,zout) 723 723 c 724 724 CALL histdef(nid_day, "evap", "Evaporation", "mm/day", 725 . iim,jjm +1,nhori, 1,1,1, -99, 32,725 . iim,jjmp1,nhori, 1,1,1, -99, 32, 726 726 . "ave(X)", zsto,zout) 727 727 c 728 728 CALL histdef(nid_day, "tops", "Solar rad. at TOA", "W/m2", 729 . iim,jjm +1,nhori, 1,1,1, -99, 32,729 . iim,jjmp1,nhori, 1,1,1, -99, 32, 730 730 . "ave(X)", zsto,zout) 731 731 c 732 732 CALL histdef(nid_day, "topl", "IR rad. at TOA", "W/m2", 733 . iim,jjm +1,nhori, 1,1,1, -99, 32,733 . iim,jjmp1,nhori, 1,1,1, -99, 32, 734 734 . "ave(X)", zsto,zout) 735 735 c 736 736 CALL histdef(nid_day, "sols", "Solar rad. at surf.", "W/m2", 737 . iim,jjm +1,nhori, 1,1,1, -99, 32,737 . iim,jjmp1,nhori, 1,1,1, -99, 32, 738 738 . "ave(X)", zsto,zout) 739 739 c 740 740 CALL histdef(nid_day, "soll", "IR rad. at surface", "W/m2", 741 . iim,jjm +1,nhori, 1,1,1, -99, 32,741 . iim,jjmp1,nhori, 1,1,1, -99, 32, 742 742 . "ave(X)", zsto,zout) 743 743 c 744 744 CALL histdef(nid_day, "bils", "Surf. total heat flux", "W/m2", 745 . iim,jjm +1,nhori, 1,1,1, -99, 32,745 . iim,jjmp1,nhori, 1,1,1, -99, 32, 746 746 . "ave(X)", zsto,zout) 747 747 c 748 748 CALL histdef(nid_day, "sens", "Sensible heat flux", "W/m2", 749 . iim,jjm +1,nhori, 1,1,1, -99, 32,749 . iim,jjmp1,nhori, 1,1,1, -99, 32, 750 750 . "ave(X)", zsto,zout) 751 751 c 752 752 CALL histdef(nid_day, "fder", "Heat flux derivation", "W/m2", 753 . iim,jjm +1,nhori, 1,1,1, -99, 32,753 . iim,jjmp1,nhori, 1,1,1, -99, 32, 754 754 . "ave(X)", zsto,zout) 755 755 c 756 756 CALL histdef(nid_day, "frtu", "Zonal wind stress", "Pa", 757 . iim,jjm +1,nhori, 1,1,1, -99, 32,757 . iim,jjmp1,nhori, 1,1,1, -99, 32, 758 758 . "ave(X)", zsto,zout) 759 759 c 760 760 CALL histdef(nid_day, "frtv", "Meridional wind stress", "Pa", 761 . iim,jjm +1,nhori, 1,1,1, -99, 32,761 . iim,jjmp1,nhori, 1,1,1, -99, 32, 762 762 . "ave(X)", zsto,zout) 763 763 c 764 764 CALL histdef(nid_day, "ruis", "Runoff", "mm/day", 765 . iim,jjm +1,nhori, 1,1,1, -99, 32,765 . iim,jjmp1,nhori, 1,1,1, -99, 32, 766 766 . "ave(X)", zsto,zout) 767 767 c 768 768 CALL histdef(nid_day, "sicf", "Sea-ice fraction", "-", 769 . iim,jjm +1,nhori, 1,1,1, -99, 32,769 . iim,jjmp1,nhori, 1,1,1, -99, 32, 770 770 . "ave(X)", zsto,zout) 771 771 c 772 772 CALL histdef(nid_day, "cldl", "Low-level cloudiness", "-", 773 . iim,jjm +1,nhori, 1,1,1, -99, 32,773 . iim,jjmp1,nhori, 1,1,1, -99, 32, 774 774 . "ave(X)", zsto,zout) 775 775 c 776 776 CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "-", 777 . iim,jjm +1,nhori, 1,1,1, -99, 32,777 . iim,jjmp1,nhori, 1,1,1, -99, 32, 778 778 . "ave(X)", zsto,zout) 779 779 c 780 780 CALL histdef(nid_day, "cldh", "High-level cloudiness", "-", 781 . iim,jjm +1,nhori, 1,1,1, -99, 32,781 . iim,jjmp1,nhori, 1,1,1, -99, 32, 782 782 . "ave(X)", zsto,zout) 783 783 c 784 784 CALL histdef(nid_day, "cldt", "Total cloudiness", "-", 785 . iim,jjm +1,nhori, 1,1,1, -99, 32,785 . iim,jjmp1,nhori, 1,1,1, -99, 32, 786 786 . "ave(X)", zsto,zout) 787 787 c 788 788 CALL histdef(nid_day, "cldq", "Cloud liquid water path", "-", 789 . iim,jjm +1,nhori, 1,1,1, -99, 32,789 . iim,jjmp1,nhori, 1,1,1, -99, 32, 790 790 . "ave(X)", zsto,zout) 791 791 c … … 793 793 c 794 794 CALL histdef(nid_day, "temp", "Air temperature", "K", 795 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,795 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 796 796 . "ave(X)", zsto,zout) 797 797 c 798 798 CALL histdef(nid_day, "ovap", "Specific humidity", "Kg/Kg", 799 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,799 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 800 800 . "ave(X)", zsto,zout) 801 801 c 802 802 CALL histdef(nid_day, "geop", "Geopotential height", "m", 803 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,803 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 804 804 . "ave(X)", zsto,zout) 805 805 c 806 806 CALL histdef(nid_day, "vitu", "Zonal wind", "m/s", 807 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,807 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 808 808 . "ave(X)", zsto,zout) 809 809 c 810 810 CALL histdef(nid_day, "vitv", "Meridional wind", "m/s", 811 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,811 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 812 812 . "ave(X)", zsto,zout) 813 813 c 814 814 CALL histdef(nid_day, "vitw", "Vertical wind", "m/s", 815 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,815 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 816 816 . "ave(X)", zsto,zout) 817 817 c 818 818 CALL histdef(nid_day, "pres", "Air pressure", "Pa", 819 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,819 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 820 820 . "ave(X)", zsto,zout) 821 821 c … … 832 832 zjulian = zjulian + day_ini 833 833 c 834 CALL gr_fi_ecrit(1,klon,iim,jjm +1,rlon,zx_lon)834 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) 835 835 DO i = 1, iim 836 836 zx_lon(i,1) = rlon(i+1) 837 zx_lon(i,jjm +1) = rlon(i+1)837 zx_lon(i,jjmp1) = rlon(i+1) 838 838 ENDDO 839 839 DO ll=1,klev 840 840 znivsig(ll)=float(ll) 841 841 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,842 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 843 CALL histbeg("histmth", iim,zx_lon, jjmp1,zx_lat, 844 . 1,iim,1,jjmp1, 0, zjulian, dtime, 845 845 . nhori, nid_mth) 846 846 c CALL histvert(nid_mth, "presnivs", "Vertical levels", "mb", … … 853 853 c 854 854 CALL histdef(nid_mth, "phis", "Surface geop. height", "-", 855 . iim,jjm +1,nhori, 1,1,1, -99, 32,855 . iim,jjmp1,nhori, 1,1,1, -99, 32, 856 856 . "once", zsto,zout) 857 857 c 858 858 CALL histdef(nid_mth, "aire", "Grid area", "-", 859 . iim,jjm +1,nhori, 1,1,1, -99, 32,859 . iim,jjmp1,nhori, 1,1,1, -99, 32, 860 860 . "once", zsto,zout) 861 861 c … … 863 863 c 864 864 CALL histdef(nid_mth, "tsol", "Surface Temperature", "K", 865 . iim,jjm +1,nhori, 1,1,1, -99, 32,865 . iim,jjmp1,nhori, 1,1,1, -99, 32, 866 866 . "ave(X)", zsto,zout) 867 867 c 868 868 CALL histdef(nid_mth, "psol", "Surface Pressure", "Pa", 869 . iim,jjm +1,nhori, 1,1,1, -99, 32,869 . iim,jjmp1,nhori, 1,1,1, -99, 32, 870 870 . "ave(X)", zsto,zout) 871 871 c 872 872 CALL histdef(nid_mth, "qsol", "Surface humidity", "mm", 873 . iim,jjm +1,nhori, 1,1,1, -99, 32,873 . iim,jjmp1,nhori, 1,1,1, -99, 32, 874 874 . "ave(X)", zsto,zout) 875 875 c 876 876 CALL histdef(nid_mth, "rain", "Precipitation", "mm/day", 877 . iim,jjm +1,nhori, 1,1,1, -99, 32,877 . iim,jjmp1,nhori, 1,1,1, -99, 32, 878 878 . "ave(X)", zsto,zout) 879 879 c 880 880 CALL histdef(nid_mth, "plul", "Large-scale Precip.", "mm/day", 881 . iim,jjm +1,nhori, 1,1,1, -99, 32,881 . iim,jjmp1,nhori, 1,1,1, -99, 32, 882 882 . "ave(X)", zsto,zout) 883 883 c 884 884 CALL histdef(nid_mth, "pluc", "Convective Precip.", "mm/day", 885 . iim,jjm +1,nhori, 1,1,1, -99, 32,885 . iim,jjmp1,nhori, 1,1,1, -99, 32, 886 886 . "ave(X)", zsto,zout) 887 887 c 888 888 CALL histdef(nid_mth, "snow", "Snow fall", "mm/day", 889 . iim,jjm +1,nhori, 1,1,1, -99, 32,889 . iim,jjmp1,nhori, 1,1,1, -99, 32, 890 890 . "ave(X)", zsto,zout) 891 891 c 892 892 CALL histdef(nid_mth, "ages", "Snow age", "day", 893 . iim,jjm +1,nhori, 1,1,1, -99, 32,893 . iim,jjmp1,nhori, 1,1,1, -99, 32, 894 894 . "ave(X)", zsto,zout) 895 895 c 896 896 CALL histdef(nid_mth, "evap", "Evaporation", "mm/day", 897 . iim,jjm +1,nhori, 1,1,1, -99, 32,897 . iim,jjmp1,nhori, 1,1,1, -99, 32, 898 898 . "ave(X)", zsto,zout) 899 899 c 900 900 CALL histdef(nid_mth, "tops", "Solar rad. at TOA", "W/m2", 901 . iim,jjm +1,nhori, 1,1,1, -99, 32,901 . iim,jjmp1,nhori, 1,1,1, -99, 32, 902 902 . "ave(X)", zsto,zout) 903 903 c 904 904 CALL histdef(nid_mth, "topl", "IR rad. at TOA", "W/m2", 905 . iim,jjm +1,nhori, 1,1,1, -99, 32,905 . iim,jjmp1,nhori, 1,1,1, -99, 32, 906 906 . "ave(X)", zsto,zout) 907 907 c 908 908 CALL histdef(nid_mth, "sols", "Solar rad. at surf.", "W/m2", 909 . iim,jjm +1,nhori, 1,1,1, -99, 32,909 . iim,jjmp1,nhori, 1,1,1, -99, 32, 910 910 . "ave(X)", zsto,zout) 911 911 c 912 912 CALL histdef(nid_mth, "soll", "IR rad. at surface", "W/m2", 913 . iim,jjm +1,nhori, 1,1,1, -99, 32,913 . iim,jjmp1,nhori, 1,1,1, -99, 32, 914 914 . "ave(X)", zsto,zout) 915 915 c 916 916 CALL histdef(nid_mth, "tops0", "Solar rad. at TOA", "W/m2", 917 . iim,jjm +1,nhori, 1,1,1, -99, 32,917 . iim,jjmp1,nhori, 1,1,1, -99, 32, 918 918 . "ave(X)", zsto,zout) 919 919 c 920 920 CALL histdef(nid_mth, "topl0", "IR rad. at TOA", "W/m2", 921 . iim,jjm +1,nhori, 1,1,1, -99, 32,921 . iim,jjmp1,nhori, 1,1,1, -99, 32, 922 922 . "ave(X)", zsto,zout) 923 923 c 924 924 CALL histdef(nid_mth, "sols0", "Solar rad. at surf.", "W/m2", 925 . iim,jjm +1,nhori, 1,1,1, -99, 32,925 . iim,jjmp1,nhori, 1,1,1, -99, 32, 926 926 . "ave(X)", zsto,zout) 927 927 c 928 928 CALL histdef(nid_mth, "soll0", "IR rad. at surface", "W/m2", 929 . iim,jjm +1,nhori, 1,1,1, -99, 32,929 . iim,jjmp1,nhori, 1,1,1, -99, 32, 930 930 . "ave(X)", zsto,zout) 931 931 c 932 932 CALL histdef(nid_mth, "bils", "Surf. total heat flux", "W/m2", 933 . iim,jjm +1,nhori, 1,1,1, -99, 32,933 . iim,jjmp1,nhori, 1,1,1, -99, 32, 934 934 . "ave(X)", zsto,zout) 935 935 c 936 936 CALL histdef(nid_mth, "sens", "Sensible heat flux", "W/m2", 937 . iim,jjm +1,nhori, 1,1,1, -99, 32,937 . iim,jjmp1,nhori, 1,1,1, -99, 32, 938 938 . "ave(X)", zsto,zout) 939 939 c 940 940 CALL histdef(nid_mth, "fder", "Heat flux derivation", "W/m2", 941 . iim,jjm +1,nhori, 1,1,1, -99, 32,941 . iim,jjmp1,nhori, 1,1,1, -99, 32, 942 942 . "ave(X)", zsto,zout) 943 943 c 944 944 CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa", 945 . iim,jjm +1,nhori, 1,1,1, -99, 32,945 . iim,jjmp1,nhori, 1,1,1, -99, 32, 946 946 . "ave(X)", zsto,zout) 947 947 c 948 948 CALL histdef(nid_mth, "frtv", "Meridional wind stress", "Pa", 949 . iim,jjm +1,nhori, 1,1,1, -99, 32,949 . iim,jjmp1,nhori, 1,1,1, -99, 32, 950 950 . "ave(X)", zsto,zout) 951 951 c 952 952 CALL histdef(nid_mth, "ruis", "Runoff", "mm/day", 953 . iim,jjm +1,nhori, 1,1,1, -99, 32,953 . iim,jjmp1,nhori, 1,1,1, -99, 32, 954 954 . "ave(X)", zsto,zout) 955 955 c 956 956 CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "-", 957 . iim,jjm +1,nhori, 1,1,1, -99, 32,957 . iim,jjmp1,nhori, 1,1,1, -99, 32, 958 958 . "ave(X)", zsto,zout) 959 959 c 960 960 CALL histdef(nid_mth, "albs", "Surface albedo", "-", 961 . iim,jjm +1,nhori, 1,1,1, -99, 32,961 . iim,jjmp1,nhori, 1,1,1, -99, 32, 962 962 . "ave(X)", zsto,zout) 963 963 c 964 964 CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-", 965 . iim,jjm +1,nhori, 1,1,1, -99, 32,965 . iim,jjmp1,nhori, 1,1,1, -99, 32, 966 966 . "ave(X)", zsto,zout) 967 967 c 968 968 CALL histdef(nid_mth, "cdrh", "Heat drag coef.", "-", 969 . iim,jjm +1,nhori, 1,1,1, -99, 32,969 . iim,jjmp1,nhori, 1,1,1, -99, 32, 970 970 . "ave(X)", zsto,zout) 971 971 c 972 972 CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "-", 973 . iim,jjm +1,nhori, 1,1,1, -99, 32,973 . iim,jjmp1,nhori, 1,1,1, -99, 32, 974 974 . "ave(X)", zsto,zout) 975 975 c 976 976 CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "-", 977 . iim,jjm +1,nhori, 1,1,1, -99, 32,977 . iim,jjmp1,nhori, 1,1,1, -99, 32, 978 978 . "ave(X)", zsto,zout) 979 979 c 980 980 CALL histdef(nid_mth, "cldh", "High-level cloudiness", "-", 981 . iim,jjm +1,nhori, 1,1,1, -99, 32,981 . iim,jjmp1,nhori, 1,1,1, -99, 32, 982 982 . "ave(X)", zsto,zout) 983 983 c 984 984 CALL histdef(nid_mth, "cldt", "Total cloudiness", "-", 985 . iim,jjm +1,nhori, 1,1,1, -99, 32,985 . iim,jjmp1,nhori, 1,1,1, -99, 32, 986 986 . "ave(X)", zsto,zout) 987 987 c 988 988 CALL histdef(nid_mth, "cldq", "Cloud liquid water path", "-", 989 . iim,jjm +1,nhori, 1,1,1, -99, 32,989 . iim,jjmp1,nhori, 1,1,1, -99, 32, 990 990 . "ave(X)", zsto,zout) 991 991 c 992 992 CALL histdef(nid_mth, "ue", "Zonal energy transport", "-", 993 . iim,jjm +1,nhori, 1,1,1, -99, 32,993 . iim,jjmp1,nhori, 1,1,1, -99, 32, 994 994 . "ave(X)", zsto,zout) 995 995 c 996 996 CALL histdef(nid_mth, "ve", "Merid energy transport", "-", 997 . iim,jjm +1,nhori, 1,1,1, -99, 32,997 . iim,jjmp1,nhori, 1,1,1, -99, 32, 998 998 . "ave(X)", zsto,zout) 999 999 c 1000 1000 CALL histdef(nid_mth, "uq", "Zonal humidity transport", "-", 1001 . iim,jjm +1,nhori, 1,1,1, -99, 32,1001 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1002 1002 . "ave(X)", zsto,zout) 1003 1003 c 1004 1004 CALL histdef(nid_mth, "vq", "Merid humidity transport", "-", 1005 . iim,jjm +1,nhori, 1,1,1, -99, 32,1005 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1006 1006 . "ave(X)", zsto,zout) 1007 1007 c … … 1009 1009 c 1010 1010 CALL histdef(nid_mth, "temp", "Air temperature", "K", 1011 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1011 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1012 1012 . "ave(X)", zsto,zout) 1013 1013 c 1014 1014 CALL histdef(nid_mth, "ovap", "Specific humidity", "Kg/Kg", 1015 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1015 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1016 1016 . "ave(X)", zsto,zout) 1017 1017 c 1018 1018 CALL histdef(nid_mth, "geop", "Geopotential height", "m", 1019 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1019 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1020 1020 . "ave(X)", zsto,zout) 1021 1021 c 1022 1022 CALL histdef(nid_mth, "vitu", "Zonal wind", "m/s", 1023 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1023 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1024 1024 . "ave(X)", zsto,zout) 1025 1025 c 1026 1026 CALL histdef(nid_mth, "vitv", "Meridional wind", "m/s", 1027 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1027 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1028 1028 . "ave(X)", zsto,zout) 1029 1029 c 1030 1030 CALL histdef(nid_mth, "vitw", "Vertical wind", "m/s", 1031 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1031 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1032 1032 . "ave(X)", zsto,zout) 1033 1033 c 1034 1034 CALL histdef(nid_mth, "pres", "Air pressure", "Pa", 1035 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1035 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1036 1036 . "ave(X)", zsto,zout) 1037 1037 c 1038 1038 CALL histdef(nid_mth, "rneb", "Cloud fraction", "-", 1039 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1039 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1040 1040 . "ave(X)", zsto,zout) 1041 1041 c 1042 1042 CALL histdef(nid_mth, "rhum", "Relative humidity", "-", 1043 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1043 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1044 1044 . "ave(X)", zsto,zout) 1045 1045 c 1046 1046 CALL histdef(nid_mth, "oliq", "Liquid water content", "kg/kg", 1047 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1047 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1048 1048 . "ave(X)", zsto,zout) 1049 1049 c 1050 1050 CALL histdef(nid_mth, "dtdyn", "Dynamics dT", "K/s", 1051 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1051 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1052 1052 . "ave(X)", zsto,zout) 1053 1053 c 1054 1054 CALL histdef(nid_mth, "dqdyn", "Dynamics dQ", "Kg/Kg/s", 1055 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1055 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1056 1056 . "ave(X)", zsto,zout) 1057 1057 c 1058 1058 CALL histdef(nid_mth, "dtcon", "Convection dT", "K/s", 1059 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1059 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1060 1060 . "ave(X)", zsto,zout) 1061 1061 c 1062 1062 CALL histdef(nid_mth, "dqcon", "Convection dQ", "Kg/Kg/s", 1063 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1063 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1064 1064 . "ave(X)", zsto,zout) 1065 1065 c 1066 1066 CALL histdef(nid_mth, "dtlsc", "Condensation dT", "K/s", 1067 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1067 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1068 1068 . "ave(X)", zsto,zout) 1069 1069 c 1070 1070 CALL histdef(nid_mth, "dqlsc", "Condensation dQ", "Kg/Kg/s", 1071 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1071 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1072 1072 . "ave(X)", zsto,zout) 1073 1073 c 1074 1074 CALL histdef(nid_mth, "dtvdf", "Boundary-layer dT", "K/s", 1075 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1075 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1076 1076 . "ave(X)", zsto,zout) 1077 1077 c 1078 1078 CALL histdef(nid_mth, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s", 1079 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1079 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1080 1080 . "ave(X)", zsto,zout) 1081 1081 c 1082 1082 CALL histdef(nid_mth, "dteva", "Reevaporation dT", "K/s", 1083 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1083 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1084 1084 . "ave(X)", zsto,zout) 1085 1085 c 1086 1086 CALL histdef(nid_mth, "dqeva", "Reevaporation dQ", "Kg/Kg/s", 1087 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1087 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1088 1088 . "ave(X)", zsto,zout) 1089 1089 c 1090 1090 CALL histdef(nid_mth, "dtajs", "Dry adjust. dT", "K/s", 1091 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1091 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1092 1092 . "ave(X)", zsto,zout) 1093 1093 1094 1094 CALL histdef(nid_mth, "dqajs", "Dry adjust. dQ", "Kg/Kg/s", 1095 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1095 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1096 1096 . "ave(X)", zsto,zout) 1097 1097 c 1098 1098 CALL histdef(nid_mth, "dtswr", "SW radiation dT", "K/s", 1099 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1099 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1100 1100 . "ave(X)", zsto,zout) 1101 1101 c 1102 1102 CALL histdef(nid_mth, "dtsw0", "SW radiation dT", "K/s", 1103 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1103 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1104 1104 . "ave(X)", zsto,zout) 1105 1105 c 1106 1106 CALL histdef(nid_mth, "dtlwr", "LW radiation dT", "K/s", 1107 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1107 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1108 1108 . "ave(X)", zsto,zout) 1109 1109 c 1110 1110 CALL histdef(nid_mth, "dtlw0", "LW radiation dT", "K/s", 1111 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1111 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1112 1112 . "ave(X)", zsto,zout) 1113 1113 c 1114 1114 CALL histdef(nid_mth, "duvdf", "Boundary-layer dU", "m/s2", 1115 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1115 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1116 1116 . "ave(X)", zsto,zout) 1117 1117 c 1118 1118 CALL histdef(nid_mth, "dvvdf", "Boundary-layer dV", "m/s2", 1119 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1119 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1120 1120 . "ave(X)", zsto,zout) 1121 1121 c 1122 1122 IF (ok_orodr) THEN 1123 1123 CALL histdef(nid_mth, "duoro", "Orography dU", "m/s2", 1124 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1124 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1125 1125 . "ave(X)", zsto,zout) 1126 1126 c 1127 1127 CALL histdef(nid_mth, "dvoro", "Orography dV", "m/s2", 1128 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1128 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1129 1129 . "ave(X)", zsto,zout) 1130 1130 c … … 1133 1133 IF (ok_orolf) THEN 1134 1134 CALL histdef(nid_mth, "dulif", "Orography dU", "m/s2", 1135 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1135 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1136 1136 . "ave(X)", zsto,zout) 1137 1137 c 1138 1138 CALL histdef(nid_mth, "dvlif", "Orography dV", "m/s2", 1139 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1139 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1140 1140 . "ave(X)", zsto,zout) 1141 1141 ENDIF 1142 1142 C 1143 1143 CALL histdef(nid_mth, "ozone", "Ozone concentration", "-", 1144 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1144 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1145 1145 . "ave(X)", zsto,zout) 1146 1146 c … … 1150 1150 WRITE(str2,'(i2.2)') iq 1151 1151 CALL histdef(nid_mth, "trac"//str2, "Tracer No."//str2, "-", 1152 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1152 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1153 1153 . "ave(X)", zsto,zout) 1154 1154 ELSE … … 1172 1172 zjulian = zjulian + day_ini 1173 1173 c 1174 CALL gr_fi_ecrit(1,klon,iim,jjm +1,rlon,zx_lon)1174 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) 1175 1175 DO i = 1, iim 1176 1176 zx_lon(i,1) = rlon(i+1) 1177 zx_lon(i,jjm +1) = rlon(i+1)1177 zx_lon(i,jjmp1) = rlon(i+1) 1178 1178 ENDDO 1179 1179 DO ll=1,klev 1180 1180 znivsig(ll)=float(ll) 1181 1181 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,1182 CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) 1183 CALL histbeg("histins", iim,zx_lon, jjmp1,zx_lat, 1184 . 1,iim,1,jjmp1, 0, zjulian, dtime, 1185 1185 . nhori, nid_ins) 1186 1186 c CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", … … 1193 1193 C 1194 1194 CALL histdef(nid_ins, "phis", "Surface geop. height", "-", 1195 . iim,jjm +1,nhori, 1,1,1, -99, 32,1195 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1196 1196 . "once", zsto,zout) 1197 1197 c 1198 1198 CALL histdef(nid_ins, "aire", "Grid area", "-", 1199 . iim,jjm +1,nhori, 1,1,1, -99, 32,1199 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1200 1200 . "once", zsto,zout) 1201 1201 c … … 1203 1203 c 1204 1204 CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa", 1205 . iim,jjm +1,nhori, 1,1,1, -99, 32,1205 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1206 1206 . "inst(X)", zsto,zout) 1207 1207 c 1208 1208 CALL histdef(nid_ins, "topl", "OLR", "W/m2", 1209 . iim,jjm +1,nhori, 1,1,1, -99, 32,1209 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1210 1210 . "inst(X)", zsto,zout) 1211 1211 c … … 1213 1213 c 1214 1214 CALL histdef(nid_ins, "temp", "Temperature", "K", 1215 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1215 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1216 1216 . "inst(X)", zsto,zout) 1217 1217 c 1218 1218 CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s", 1219 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1219 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1220 1220 . "inst(X)", zsto,zout) 1221 1221 c 1222 1222 CALL histdef(nid_ins, "vitv", "Merid wind", "m/s", 1223 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1223 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1224 1224 . "inst(X)", zsto,zout) 1225 1225 c 1226 1226 CALL histdef(nid_ins, "geop", "Geopotential height", "m", 1227 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1227 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1228 1228 . "inst(X)", zsto,zout) 1229 1229 c 1230 1230 CALL histdef(nid_ins, "pres", "Air pressure", "Pa", 1231 . iim,jjm +1,nhori, klev,1,klev,nvert, 32,1231 . iim,jjmp1,nhori, klev,1,klev,nvert, 32, 1232 1232 . "inst(X)", zsto,zout) 1233 1233 c … … 1368 1368 IF (ok_oasis .AND. MOD(itap-1,nexca).EQ.0) THEN 1369 1369 C 1370 CALL fromcpl(itap, (jjm+1)*iim,1370 CALL fromcpl(itap,jjmp1*iim, 1371 1371 . cpl_sst,cpl_sic,cpl_alb_sst,cpl_alb_sic) 1372 1372 DO i = 1, iim-1 ! un seul point pour le pole nord … … 1377 1377 ENDDO 1378 1378 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)1379 cpl_sst(i,jjmp1) = cpl_sst(1,jjmp1) 1380 cpl_sic(i,jjmp1) = cpl_sic(1,jjmp1) 1381 cpl_alb_sst(i,jjmp1) = cpl_alb_sst(1,jjmp1) 1382 cpl_alb_sic(i,jjmp1) = cpl_alb_sic(1,jjmp1) 1383 1383 ENDDO 1384 1384 c … … 1407 1407 . pctsrf(ig,is_sic).GT.epsfra) THEN 1408 1408 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)1409 . - (cpl_sic(1,jjmp1)-pctsrf(ig,is_sic)) 1410 pctsrf(ig,is_sic) = cpl_sic(1,jjmp1) 1411 lmt_sst(ig) = cpl_sst(1,jjmp1) 1412 1412 ENDIF 1413 1413 c … … 1975 1975 c 1976 1976 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)1977 CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) 1978 CALL histwrite(nid_day,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d) 1979 1979 c 1980 1980 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)1981 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) 1982 CALL histwrite(nid_day,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d) 1983 1983 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)1984 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 1985 CALL histwrite(nid_day,"tsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 1986 1986 c 1987 1987 DO i = 1, klon 1988 1988 zx_tmp_fi2d(i) = paprs(i,1) 1989 1989 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)1990 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 1991 CALL histwrite(nid_day,"psol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 1992 1992 c 1993 1993 DO i = 1, klon 1994 1994 zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) 1995 1995 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)1996 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 1997 CALL histwrite(nid_day,"rain",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 1998 c 1999 CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) 2000 CALL histwrite(nid_day,"snow",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2001 c 2002 CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) 2003 CALL histwrite(nid_day,"evap",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2004 c 2005 CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) 2006 CALL histwrite(nid_day,"tops",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2007 c 2008 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) 2009 CALL histwrite(nid_day,"topl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2010 c 2011 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) 2012 CALL histwrite(nid_day,"sols",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2013 c 2014 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) 2015 CALL histwrite(nid_day,"soll",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2016 c 2017 CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) 2018 CALL histwrite(nid_day,"bils",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2019 c 2020 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) 2021 CALL histwrite(nid_day,"sens",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2022 c 2023 CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) 2024 CALL histwrite(nid_day,"fder",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2025 c 2026 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ruis,zx_tmp_2d) 2027 CALL histwrite(nid_day,"ruis",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2028 2028 c 2029 2029 DO i = 1, klon 2030 2030 zx_tmp_fi2d(i) = fluxu(i,1) 2031 2031 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)2032 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2033 CALL histwrite(nid_day,"frtu",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2034 2034 c 2035 2035 DO i = 1, klon 2036 2036 zx_tmp_fi2d(i) = fluxv(i,1) 2037 2037 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)2038 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2039 CALL histwrite(nid_day,"frtv",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2040 2040 c 2041 2041 DO i = 1, klon 2042 2042 zx_tmp_fi2d(i) = pctsrf(i,is_sic) 2043 2043 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)2044 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2045 CALL histwrite(nid_day,"sicf",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2046 c 2047 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) 2048 CALL histwrite(nid_day,"cldl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2049 c 2050 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) 2051 CALL histwrite(nid_day,"cldm",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2052 c 2053 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) 2054 CALL histwrite(nid_day,"cldh",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2055 c 2056 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d) 2057 CALL histwrite(nid_day,"cldt",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2058 c 2059 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) 2060 CALL histwrite(nid_day,"cldq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2061 2061 c 2062 2062 c Champs 3D: 2063 2063 c 2064 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, t_seri, zx_tmp_3d)2064 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) 2065 2065 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)2066 . iim*jjmp1*klev,ndex3d) 2067 c 2068 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) 2069 2069 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)2070 . iim*jjmp1*klev,ndex3d) 2071 c 2072 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) 2073 2073 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)2074 . iim*jjmp1*klev,ndex3d) 2075 c 2076 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 2077 2077 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)2078 . iim*jjmp1*klev,ndex3d) 2079 c 2080 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 2081 2081 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)2082 . iim*jjmp1*klev,ndex3d) 2083 c 2084 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) 2085 2085 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)2086 . iim*jjmp1*klev,ndex3d) 2087 c 2088 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) 2089 2089 CALL histwrite(nid_day,"pres",itap,zx_tmp_3d, 2090 . iim* (jjm+1)*klev,ndex3d)2090 . iim*jjmp1*klev,ndex3d) 2091 2091 c 2092 2092 if (ok_sync) then … … 2103 2103 c 2104 2104 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)2105 CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) 2106 CALL histwrite(nid_mth,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d) 2107 2107 C 2108 2108 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)2109 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) 2110 CALL histwrite(nid_mth,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d) 2111 2112 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 2113 CALL histwrite(nid_mth,"tsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2114 2114 c 2115 2115 DO i = 1, klon 2116 2116 zx_tmp_fi2d(i) = paprs(i,1) 2117 2117 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)2118 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2119 CALL histwrite(nid_mth,"psol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2120 c 2121 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsol,zx_tmp_2d) 2122 CALL histwrite(nid_mth,"qsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2123 2123 c 2124 2124 DO i = 1, klon 2125 2125 zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) 2126 2126 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)2127 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2128 CALL histwrite(nid_mth,"rain",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2129 2129 c 2130 2130 DO i = 1, klon 2131 2131 zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i) 2132 2132 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)2133 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2134 CALL histwrite(nid_mth,"plul",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2135 2135 c 2136 2136 DO i = 1, klon 2137 2137 zx_tmp_fi2d(i) = rain_con(i) + snow_con(i) 2138 2138 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)2139 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2140 CALL histwrite(nid_mth,"pluc",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2141 c 2142 CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) 2143 CALL histwrite(nid_mth,"snow",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2144 c 2145 CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d) 2146 CALL histwrite(nid_mth,"ages",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2147 c 2148 CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) 2149 CALL histwrite(nid_mth,"evap",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2150 c 2151 CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) 2152 CALL histwrite(nid_mth,"tops",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2153 c 2154 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) 2155 CALL histwrite(nid_mth,"topl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2156 c 2157 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) 2158 CALL histwrite(nid_mth,"sols",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2159 c 2160 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) 2161 CALL histwrite(nid_mth,"soll",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2162 c 2163 CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d) 2164 CALL histwrite(nid_mth,"tops0",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2165 c 2166 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d) 2167 CALL histwrite(nid_mth,"topl0",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2168 c 2169 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d) 2170 CALL histwrite(nid_mth,"sols0",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2171 c 2172 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d) 2173 CALL histwrite(nid_mth,"soll0",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2174 c 2175 CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) 2176 CALL histwrite(nid_mth,"bils",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2177 c 2178 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) 2179 CALL histwrite(nid_mth,"sens",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2180 c 2181 CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) 2182 CALL histwrite(nid_mth,"fder",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2183 c 2184 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ruis,zx_tmp_2d) 2185 CALL histwrite(nid_mth,"ruis",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2186 2186 c 2187 2187 DO i = 1, klon 2188 2188 zx_tmp_fi2d(i) = fluxu(i,1) 2189 2189 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)2190 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2191 CALL histwrite(nid_mth,"frtu",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2192 2192 c 2193 2193 DO i = 1, klon 2194 2194 zx_tmp_fi2d(i) = fluxv(i,1) 2195 2195 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)2196 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2197 CALL histwrite(nid_mth,"frtv",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2198 2198 c 2199 2199 DO i = 1, klon 2200 2200 zx_tmp_fi2d(i) = pctsrf(i,is_sic) 2201 2201 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)2202 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2203 CALL histwrite(nid_mth,"sicf",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2204 c 2205 CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d) 2206 CALL histwrite(nid_mth,"albs",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2207 c 2208 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d) 2209 CALL histwrite(nid_mth,"cdrm",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2210 c 2211 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d) 2212 CALL histwrite(nid_mth,"cdrh",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2213 c 2214 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) 2215 CALL histwrite(nid_mth,"cldl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2216 c 2217 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) 2218 CALL histwrite(nid_mth,"cldm",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2219 c 2220 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) 2221 CALL histwrite(nid_mth,"cldh",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2222 c 2223 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d) 2224 CALL histwrite(nid_mth,"cldt",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2225 c 2226 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) 2227 CALL histwrite(nid_mth,"cldq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2228 c 2229 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d) 2230 CALL histwrite(nid_mth,"ue",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2231 c 2232 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d) 2233 CALL histwrite(nid_mth,"ve",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2234 c 2235 CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d) 2236 CALL histwrite(nid_mth,"uq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2237 c 2238 CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d) 2239 CALL histwrite(nid_mth,"vq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2240 2240 c 2241 2241 c Champs 3D: 2242 2242 C 2243 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, t_seri, zx_tmp_3d)2243 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) 2244 2244 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)2245 . iim*jjmp1*klev,ndex3d) 2246 c 2247 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) 2248 2248 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)2249 . iim*jjmp1*klev,ndex3d) 2250 c 2251 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) 2252 2252 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)2253 . iim*jjmp1*klev,ndex3d) 2254 c 2255 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 2256 2256 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)2257 . iim*jjmp1*klev,ndex3d) 2258 c 2259 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 2260 2260 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)2261 . iim*jjmp1*klev,ndex3d) 2262 c 2263 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) 2264 2264 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)2265 . iim*jjmp1*klev,ndex3d) 2266 c 2267 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) 2268 2268 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)2269 . iim*jjmp1*klev,ndex3d) 2270 c 2271 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d) 2272 2272 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)2273 . iim*jjmp1*klev,ndex3d) 2274 c 2275 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d) 2276 2276 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)2277 . iim*jjmp1*klev,ndex3d) 2278 c 2279 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldliq, zx_tmp_3d) 2280 2280 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)2281 . iim*jjmp1*klev,ndex3d) 2282 c 2283 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d) 2284 2284 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)2285 . iim*jjmp1*klev,ndex3d) 2286 c 2287 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d) 2288 2288 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)2289 . iim*jjmp1*klev,ndex3d) 2290 c 2291 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_con, zx_tmp_3d) 2292 2292 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)2293 . iim*jjmp1*klev,ndex3d) 2294 c 2295 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_con, zx_tmp_3d) 2296 2296 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)2297 . iim*jjmp1*klev,ndex3d) 2298 c 2299 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_lsc, zx_tmp_3d) 2300 2300 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)2301 . iim*jjmp1*klev,ndex3d) 2302 c 2303 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_lsc, zx_tmp_3d) 2304 2304 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)2305 . iim*jjmp1*klev,ndex3d) 2306 c 2307 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_vdf, zx_tmp_3d) 2308 2308 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)2309 . iim*jjmp1*klev,ndex3d) 2310 c 2311 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_vdf, zx_tmp_3d) 2312 2312 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)2313 . iim*jjmp1*klev,ndex3d) 2314 c 2315 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_eva, zx_tmp_3d) 2316 2316 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)2317 . iim*jjmp1*klev,ndex3d) 2318 c 2319 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_eva, zx_tmp_3d) 2320 2320 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)2321 . iim*jjmp1*klev,ndex3d) 2322 c 2323 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_ajs, zx_tmp_3d) 2324 2324 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)2325 . iim*jjmp1*klev,ndex3d) 2326 c 2327 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_ajs, zx_tmp_3d) 2328 2328 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)2329 . iim*jjmp1*klev,ndex3d) 2330 c 2331 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat, zx_tmp_3d) 2332 2332 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)2333 . iim*jjmp1*klev,ndex3d) 2334 c 2335 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, heat0, zx_tmp_3d) 2336 2336 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)2337 . iim*jjmp1*klev,ndex3d) 2338 c 2339 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool, zx_tmp_3d) 2340 2340 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)2341 . iim*jjmp1*klev,ndex3d) 2342 c 2343 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cool0, zx_tmp_3d) 2344 2344 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)2345 . iim*jjmp1*klev,ndex3d) 2346 c 2347 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d) 2348 2348 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)2349 . iim*jjmp1*klev,ndex3d) 2350 c 2351 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d) 2352 2352 CALL histwrite(nid_mth,"dvvdf",itap,zx_tmp_3d, 2353 . iim* (jjm+1)*klev,ndex3d)2353 . iim*jjmp1*klev,ndex3d) 2354 2354 c 2355 2355 IF (ok_orodr) THEN 2356 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_u_oro, zx_tmp_3d)2356 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_oro, zx_tmp_3d) 2357 2357 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)2358 . iim*jjmp1*klev,ndex3d) 2359 c 2360 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_oro, zx_tmp_3d) 2361 2361 CALL histwrite(nid_mth,"dvoro",itap,zx_tmp_3d, 2362 . iim* (jjm+1)*klev,ndex3d)2362 . iim*jjmp1*klev,ndex3d) 2363 2363 c 2364 2364 ENDIF 2365 2365 C 2366 2366 IF (ok_orolf) THEN 2367 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, d_u_lif, zx_tmp_3d)2367 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_lif, zx_tmp_3d) 2368 2368 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)2369 . iim*jjmp1*klev,ndex3d) 2370 c 2371 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_lif, zx_tmp_3d) 2372 2372 CALL histwrite(nid_mth,"dvlif",itap,zx_tmp_3d, 2373 . iim* (jjm+1)*klev,ndex3d)2373 . iim*jjmp1*klev,ndex3d) 2374 2374 ENDIF 2375 2375 C 2376 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, wo, zx_tmp_3d)2376 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d) 2377 2377 CALL histwrite(nid_mth,"ozone",itap,zx_tmp_3d, 2378 . iim* (jjm+1)*klev,ndex3d)2378 . iim*jjmp1*klev,ndex3d) 2379 2379 c 2380 2380 IF (nqmax.GE.3) THEN 2381 2381 DO iq=1,nqmax-2 2382 2382 IF (iq.LE.99) THEN 2383 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, qx(1,1,iq+2), zx_tmp_3d)2383 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq+2), zx_tmp_3d) 2384 2384 WRITE(str2,'(i2.2)') iq 2385 2385 CALL histwrite(nid_mth,"trac"//str2,itap,zx_tmp_3d, 2386 . iim* (jjm+1)*klev,ndex3d)2386 . iim*jjmp1*klev,ndex3d) 2387 2387 ELSE 2388 2388 PRINT*, "Trop de traceurs" … … 2405 2405 c 2406 2406 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)2407 CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) 2408 CALL histwrite(nid_ins,"phis",i,zx_tmp_2d,iim*jjmp1,ndex2d) 2409 2409 c 2410 2410 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)2411 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) 2412 CALL histwrite(nid_ins,"aire",i,zx_tmp_2d,iim*jjmp1,ndex2d) 2413 2413 2414 2414 DO i = 1, klon 2415 2415 zx_tmp_fi2d(i) = paprs(i,1) 2416 2416 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)2417 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2418 CALL histwrite(nid_ins,"psol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2419 c 2420 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) 2421 CALL histwrite(nid_ins,"topl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2422 2422 c 2423 2423 c Champs 3D: 2424 2424 c 2425 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, t_seri, zx_tmp_3d)2425 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) 2426 2426 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)2427 . iim*jjmp1*klev,ndex3d) 2428 c 2429 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 2430 2430 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)2431 . iim*jjmp1*klev,ndex3d) 2432 c 2433 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 2434 2434 CALL histwrite(nid_ins,"vitv",itap,zx_tmp_3d, 2435 . iim* (jjm+1)*klev,ndex3d)2436 c 2437 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, zphi, zx_tmp_3d)2435 . iim*jjmp1*klev,ndex3d) 2436 c 2437 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) 2438 2438 CALL histwrite(nid_ins,"geop",itap,zx_tmp_3d, 2439 . iim* (jjm+1)*klev,ndex3d)2440 c 2441 CALL gr_fi_ecrit(klev,klon,iim,jjm +1, pplay, zx_tmp_3d)2439 . iim*jjmp1*klev,ndex3d) 2440 c 2441 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) 2442 2442 CALL histwrite(nid_ins,"pres",itap,zx_tmp_3d, 2443 . iim* (jjm+1)*klev,ndex3d)2443 . iim*jjmp1*klev,ndex3d) 2444 2444 c 2445 2445 if (ok_sync) then … … 2490 2490 ig = ig + 1 2491 2491 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)2492 z_sols(i,jjmp1) = oas_sols(ig) 2493 z_nsol(i,jjmp1) = oas_nsol(ig) 2494 z_rain(i,jjmp1) = oas_rain(ig) 2495 z_snow(i,jjmp1) = oas_snow(ig) 2496 z_evap(i,jjmp1) = oas_evap(ig) 2497 z_ruisoce(i,jjmp1) = oas_ruisoce(ig) 2498 z_ruisriv(i,jjmp1) = oas_ruisriv(ig) 2499 z_tsol(i,jjmp1) = oas_tsol(ig) 2500 z_fder(i,jjmp1) = oas_fder(ig) 2501 z_albe(i,jjmp1) = oas_albe(ig) 2502 z_taux(i,jjmp1) = oas_taux(ig) 2503 z_tauy(i,jjmp1) = oas_tauy(ig) 2504 2504 ENDDO 2505 2505 c 2506 2506 c Passer les champs au coupleur: 2507 2507 c 2508 CALL intocpl(itap, (jjm+1)*iim,2508 CALL intocpl(itap,jjmp1*iim, 2509 2509 . z_sols, z_nsol, 2510 2510 . z_rain, z_snow, z_evap,
Note: See TracChangeset
for help on using the changeset viewer.
