Changeset 1119
- Timestamp:
- Mar 11, 2009, 3:45:10 PM (15 years ago)
- Location:
- LMDZ4/branches/LMDZ4-dev/libf/phylmd
- Files:
-
- 2 edited
-
phys_output_mod.F90 (modified) (6 diffs)
-
phys_output_write.h (modified) (11 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_mod.F90
r1114 r1119 714 714 CALL histdef2d(iff,o_fbase%flag,o_fbase%name, "Cld base mass flux", "kg/m2/s") 715 715 CALL histdef2d(iff,o_prw%flag,o_prw%name, "Precipitable water", "kg/m2") 716 CALL histdef2d(iff,o_cape_max%flag,o_cape_max%name,"CAPE max.", "J/kg") 717 CALL histdef3d(iff,o_upwd%flag,o_upwd%name,"saturated updraft", "kg/m2/s") 718 CALL histdef3d(iff,o_Ma%flag,o_Ma%name,"undilute adiab updraft","kg/m2/s") 719 CALL histdef3d(iff,o_dnwd%flag,o_dnwd%name,"saturated downdraft","kg/m2/s") 720 CALL histdef3d(iff,o_dnwd0%flag,o_dnwd0%name,"unsat. downdraft", "kg/m2/s") 716 721 ENDIF !iflag_con .GE. 3 717 722 … … 768 773 CALL histdef2d(iff,o_slab_bils%flag, o_slab_bils%name,"Bilan au sol sur ocean slab", "W/m2") 769 774 770 IF (iflag_con.GE.3) THEN 771 CALL histdef2d(iff,o_ale_bl%flag,o_ale_bl%name, "ALE BL", "m2/s2") 772 CALL histdef2d(iff,o_alp_bl%flag,o_alp_bl%name, "ALP BL", "m2/s2") 773 CALL histdef2d(iff,o_ale_wk%flag,o_ale_wk%name, "ALE WK", "m2/s2") 774 CALL histdef2d(iff,o_alp_wk%flag,o_alp_wk%name, "ALP WK", "m2/s2") 775 776 CALL histdef2d(iff,o_ale%flag,o_ale%name, "ALE", "m2/s2") 777 CALL histdef2d(iff,o_alp%flag,o_alp%name, "ALP", "W/m2") 778 CALL histdef2d(iff,o_cin%flag,o_cin%name, "Convective INhibition", "m2/s2") 779 CALL histdef2d(iff,o_wape%flag,o_WAPE%name, "WAPE", "m2/s2") 780 ENDIF !(iflag_con.GE.3) 775 IF (iflag_con.EQ.3) THEN 776 CALL histdef2d(iff,o_ale%flag,o_ale%name,"ALE","m2/s2") 777 CALL histdef2d(iff,o_alp%flag,o_alp%name,"ALP","W/m2") 778 CALL histdef2d(iff,o_cin%flag,o_cin%name,"Convective INhibition","m2/s2") 779 IF (iflag_coupl.EQ.1) THEN 780 CALL histdef2d(iff,o_ale_bl%flag,o_ale_bl%name,"ALE BL","m2/s2") 781 CALL histdef2d(iff,o_alp_bl%flag,o_alp_bl%name,"ALP BL","m2/s2") 782 ENDIF !iflag_coupl.EQ.1 783 IF (iflag_wake.EQ.1) THEN 784 CALL histdef2d(iff,o_ale_wk%flag,o_ale_wk%name,"ALE WK","m2/s2") 785 CALL histdef2d(iff,o_alp_wk%flag,o_alp_wk%name,"ALP WK","m2/s2") 786 CALL histdef2d(iff,o_wape%flag,o_WAPE%name,"WAPE","m2/s2") 787 CALL histdef2d(iff,o_wake_h%flag,o_wake_h%name,"wake_h", "-") 788 CALL histdef2d(iff,o_wake_s%flag,o_wake_s%name,"wake_s", "-") 789 CALL histdef3d(iff,o_dtwak%flag,o_dtwak%name,"Wake dT","K/s") 790 CALL histdef3d(iff,o_dqwak%flag,o_dqwak%name,"Wake dQ","(kg/kg)/s") 791 CALL histdef3d(iff,o_wake_deltat%flag,o_wake_deltat%name,"wake_deltat", " ") 792 CALL histdef3d(iff,o_wake_deltaq%flag,o_wake_deltaq%name,"wake_deltaq", " ") 793 CALL histdef3d(iff,o_wake_omg%flag,o_wake_omg%name,"wake_omg", " ") 794 CALL histdef3d(iff,o_ftd%flag,o_ftd%name,"tend temp due aux descentes precip","-") 795 CALL histdef3d(iff,o_fqd%flag,o_fqd%name,"tend vap eau due aux descentes precip","-") 796 ENDIF !iflag_wake.EQ.1 797 CALL histdef3d(iff,o_Vprecip%flag,o_Vprecip%name,"precipitation vertical profile","-") 798 ENDIF !(iflag_con.EQ.3) 781 799 782 800 CALL histdef2d(iff,o_weakinv%flag,o_weakinv%name, "Weak inversion", "-") … … 790 808 CALL histdef2d(iff,o_LWdownOR%flag,o_LWdownOR%name, "Sfce incident LW radiation OR", "W/m2") 791 809 CALL histdef2d(iff,o_snowl%flag,o_snowl%name, "Solid Large-scale Precip.", "kg/(m2*s)") 792 IF (iflag_con.GE.3) THEN793 CALL histdef2d(iff,o_cape_max%flag,o_cape_max%name, "CAPE max.", "J/kg")794 ENDIF !(iflag_con.GE.3)795 810 CALL histdef2d(iff,o_solldown%flag,o_solldown%name, "Down. IR rad. at surface", "W/m2") 796 811 CALL histdef2d(iff,o_dtsvdfo%flag,o_dtsvdfo%name, "Boundary-layer dTs(o)", "K/s") … … 815 830 CALL histdef3d(iff,o_rhum%flag,o_rhum%name, "Relative humidity", "-") 816 831 CALL histdef3d(iff,o_ozone%flag,o_ozone%name, "Ozone concentration", "ppmv") 817 IF(iflag_con.GE.3) THEN818 CALL histdef3d(iff,o_upwd%flag,o_upwd%name, "saturated updraft", "kg/m2/s")819 ENDIF !iflag_con.GE.3820 832 CALL histdef3d(iff,o_dtphy%flag,o_dtphy%name, "Physics dT", "K/s") 821 833 CALL histdef3d(iff,o_dqphy%flag,o_dqphy%name, "Physics dQ", "(kg/kg)/s") … … 836 848 CALL histdef3d(iff,o_kz%flag,o_kz%name, "Kz melange", "m2/s") 837 849 CALL histdef3d(iff,o_kz_max%flag,o_kz_max%name, "Kz melange max", "m2/s" ) 838 CALL histdef3d(iff,o_clwcon%flag,o_clwcon%name, "Convective Cloud Liquid water content", "kg/kg") 839 IF(iflag_con.GE.3) THEN 840 CALL histdef3d(iff,o_Ma%flag,o_Ma%name, "undilute adiab updraft", "kg/m2/s") 841 CALL histdef3d(iff,o_dnwd%flag,o_dnwd%name, "saturated downdraft", "kg/m2/s") 842 CALL histdef3d(iff,o_dnwd0%flag,o_dnwd0%name, "unsat. downdraft", "kg/m2/s") 843 ENDIF !(iflag_con.GE.3) 850 CALL histdef3d(iff,o_clwcon%flag,o_clwcon%name, "Convective Cloud Liquid water content", "kg/kg") 844 851 CALL histdef3d(iff,o_dtdyn%flag,o_dtdyn%name, "Dynamics dT", "K/s") 845 852 CALL histdef3d(iff,o_dqdyn%flag,o_dqdyn%name, "Dynamics dQ", "(kg/kg)/s") … … 849 856 CALL histdef3d(iff,o_ducon%flag,o_ducon%name, "Convection du", "m/s2") 850 857 CALL histdef3d(iff,o_dqcon%flag,o_dqcon%name, "Convection dQ", "(kg/kg)/s") 851 IF(iflag_con.EQ.3) THEN852 IF (iflag_wake == 1) THEN853 CALL histdef2d(iff,o_wake_h%flag,o_wake_h%name, "wake_h", "-")854 CALL histdef2d(iff,o_wake_s%flag,o_wake_s%name, "wake_s", "-")855 CALL histdef3d(iff,o_dtwak%flag,o_dtwak%name, "Wake dT", "K/s")856 CALL histdef3d(iff,o_dqwak%flag,o_dqwak%name, "Wake dQ", "(kg/kg)/s")857 CALL histdef3d(iff,o_wake_deltat%flag,o_wake_deltat%name, "wake_deltat", " ")858 CALL histdef3d(iff,o_wake_deltaq%flag,o_wake_deltaq%name, "wake_deltaq", " ")859 CALL histdef3d(iff,o_wake_omg%flag,o_wake_omg%name, "wake_omg", " ")860 ENDIF861 CALL histdef3d(iff,o_Vprecip%flag,o_Vprecip%name, "precipitation vertical profile", "-")862 CALL histdef3d(iff,o_ftd%flag,o_ftd%name, "tend temp due aux descentes precip", "-")863 CALL histdef3d(iff,o_fqd%flag,o_fqd%name,"tend vap eau due aux descentes precip", "-")864 ENDIF !(iflag_con.EQ.3)865 866 858 CALL histdef3d(iff,o_dtlsc%flag,o_dtlsc%name, "Condensation dT", "K/s") 867 859 CALL histdef3d(iff,o_dtlschr%flag,o_dtlschr%name,"Large-scale condensational heating rate","K/s") -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_write.h
r1114 r1119 511 511 CALL histwrite_phy(nid_files(iff),o_prw%name,itau_w,prw) 512 512 ENDIF 513 IF (o_cape_max%flag(iff)<=lev_files(iff)) THEN 514 CALL histwrite_phy(nid_files(iff),o_cape_max%name,itau_w,cape) 515 ENDIF 516 IF (o_upwd%flag(iff)<=lev_files(iff)) THEN 517 CALL histwrite_phy(nid_files(iff),o_upwd%name,itau_w,upwd) 518 ENDIF 519 IF (o_Ma%flag(iff)<=lev_files(iff)) THEN 520 CALL histwrite_phy(nid_files(iff),o_Ma%name,itau_w,Ma) 521 ENDIF 522 523 IF (o_dnwd%flag(iff)<=lev_files(iff)) THEN 524 CALL histwrite_phy(nid_files(iff),o_dnwd%name,itau_w,dnwd) 525 ENDIF 526 527 IF (o_dnwd0%flag(iff)<=lev_files(iff)) THEN 528 CALL histwrite_phy(nid_files(iff),o_dnwd0%name,itau_w,dnwd0) 529 ENDIF 513 530 ENDIF !iflag_con .GE. 3 514 531 … … 619 636 IF (type_ocean=='force ') THEN 620 637 621 IF (iflag_con.GE.3) THEN 638 IF (iflag_con.EQ.3) THEN 639 IF (o_ale%flag(iff)<=lev_files(iff)) THEN 640 CALL histwrite_phy(nid_files(iff),o_ale%name,itau_w,ale) 641 ENDIF 642 IF (o_alp%flag(iff)<=lev_files(iff)) THEN 643 CALL histwrite_phy(nid_files(iff),o_alp%name,itau_w,alp) 644 ENDIF 645 IF (o_cin%flag(iff)<=lev_files(iff)) THEN 646 CALL histwrite_phy(nid_files(iff),o_cin%name,itau_w,cin) 647 ENDIF 622 648 IF (iflag_coupl.EQ.1) THEN 623 649 IF (o_ale_bl%flag(iff)<=lev_files(iff)) THEN … … 628 654 ENDIF 629 655 ENDIF !iflag_coupl.EQ.1 630 ENDIF !(iflag_con.GE.3) 631 IF (iflag_con.EQ.3) THEN 656 632 657 IF (iflag_wake.EQ.1) THEN 633 658 IF (o_ale_wk%flag(iff)<=lev_files(iff)) THEN … … 638 663 ENDIF 639 664 640 IF (o_ale%flag(iff)<=lev_files(iff)) THEN641 CALL histwrite_phy(nid_files(iff),o_ale%name,itau_w,ale)642 ENDIF643 IF (o_alp%flag(iff)<=lev_files(iff)) THEN644 CALL histwrite_phy(nid_files(iff),o_alp%name,itau_w,alp)645 ENDIF646 IF (o_cin%flag(iff)<=lev_files(iff)) THEN647 CALL histwrite_phy(nid_files(iff),o_cin%name,itau_w,cin)648 ENDIF649 665 IF (o_wape%flag(iff)<=lev_files(iff)) THEN 650 666 CALL histwrite_phy(nid_files(iff),o_WAPE%name,itau_w,wake_pe) … … 669 685 670 686 IF (o_wake_omg%flag(iff)<=lev_files(iff)) THEN 671 CALL histwrite_phy(nid_files(iff), 672 s o_wake_omg%name,itau_w,wake_omg) 673 ENDIF 674 ENDIF 675 676 IF (o_Vprecip%flag(iff)<=lev_files(iff)) THEN 677 CALL histwrite_phy(nid_files(iff),o_Vprecip%name,itau_w,Vprecip) 687 CALL histwrite_phy(nid_files(iff),o_wake_omg%name,itau_w, 688 .wake_omg) 689 ENDIF 690 IF (o_dtwak%flag(iff)<=lev_files(iff)) THEN 691 zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev) 692 & /pdtphys 693 CALL histwrite_phy(nid_files(iff),o_dtwak%name,itau_w, 694 .zx_tmp_fi3d) 695 ENDIF 696 697 IF (o_dqwak%flag(iff)<=lev_files(iff)) THEN 698 zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys 699 CALL histwrite_phy(nid_files(iff),o_dqwak%name,itau_w, 700 .zx_tmp_fi3d) 678 701 ENDIF 679 702 … … 685 708 CALL histwrite_phy(nid_files(iff),o_fqd%name,itau_w,fqd) 686 709 ENDIF 710 ENDIF 711 712 IF (o_Vprecip%flag(iff)<=lev_files(iff)) THEN 713 CALL histwrite_phy(nid_files(iff),o_Vprecip%name,itau_w,Vprecip) 714 ENDIF 687 715 688 716 ENDIF !(iflag_con.EQ.3) … … 690 718 ELSE IF (type_ocean=='slab ') THEN 691 719 692 IF ( o_slab_bils%flag(iff)<=lev_files(iff))720 IF ( o_slab_bils%flag(iff)<=lev_files(iff)) 693 721 $ CALL histwrite_phy( 694 722 $ nid_files(iff),o_slab_bils%name,itau_w,slab_wfbils) 695 723 696 724 ENDIF !type_ocean == force/slab 697 725 … … 768 796 CALL histwrite_phy(nid_files(iff),o_snowl%name,itau_w,snow_lsc) 769 797 ENDIF 770 771 IF (iflag_con.GE.3) THEN772 IF (o_cape_max%flag(iff)<=lev_files(iff)) THEN773 CALL histwrite_phy(nid_files(iff),o_cape_max%name,itau_w,cape)774 ENDIF775 ENDIF !(iflag_con.GE.3)776 798 777 799 IF (o_solldown%flag(iff)<=lev_files(iff)) THEN … … 885 907 CALL histwrite_phy(nid_files(iff),o_ozone%name,itau_w,zx_tmp_fi3d) 886 908 ENDIF 887 888 IF(iflag_con.GE.3) THEN889 IF (o_upwd%flag(iff)<=lev_files(iff)) THEN890 CALL histwrite_phy(nid_files(iff),o_upwd%name,itau_w,upwd)891 ENDIF892 ENDIF !(iflag_con.GE.3)893 909 894 910 IF (o_dtphy%flag(iff)<=lev_files(iff)) THEN … … 963 979 ENDIF 964 980 965 IF(iflag_con.GE.3) THEN966 IF (o_Ma%flag(iff)<=lev_files(iff)) THEN967 CALL histwrite_phy(nid_files(iff),o_Ma%name,itau_w,Ma)968 ENDIF969 970 IF (o_dnwd%flag(iff)<=lev_files(iff)) THEN971 CALL histwrite_phy(nid_files(iff),o_dnwd%name,itau_w,dnwd)972 ENDIF973 974 IF (o_dnwd0%flag(iff)<=lev_files(iff)) THEN975 CALL histwrite_phy(nid_files(iff),o_dnwd0%name,itau_w,dnwd0)976 ENDIF977 ENDIF !(iflag_con.GE.3)978 979 981 IF (o_dtdyn%flag(iff)<=lev_files(iff)) THEN 980 982 CALL histwrite_phy(nid_files(iff),o_dtdyn%name,itau_w,d_t_dyn) … … 1007 1009 CALL histwrite_phy(nid_files(iff),o_dqcon%name,itau_w,zx_tmp_fi3d) 1008 1010 ENDIF 1009 1010 IF(iflag_con.GE.3) THEN1011 IF (iflag_wake == 1) THEN1012 IF (o_dtwak%flag(iff)<=lev_files(iff)) THEN1013 zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)1014 & /pdtphys1015 CALL histwrite_phy(nid_files(iff),1016 & o_dtwak%name,itau_w,zx_tmp_fi3d)1017 ENDIF1018 1019 IF (o_dqwak%flag(iff)<=lev_files(iff)) THEN1020 zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys1021 CALL histwrite_phy(nid_files(iff),1022 & o_dqwak%name,itau_w,zx_tmp_fi3d)1023 ENDIF1024 ENDIF1025 ENDIF !(iflag_con.GE.3)1026 1011 1027 1012 IF (o_dtlsc%flag(iff)<=lev_files(iff)) THEN
Note: See TracChangeset
for help on using the changeset viewer.
