Changeset 1100 for LMDZ4/trunk
- Timestamp:
- Feb 16, 2009, 7:01:12 PM (16 years ago)
- Location:
- LMDZ4/trunk/libf/phylmd
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/trunk/libf/phylmd/phys_output_mod.F90
r1081 r1100 584 584 CALL histdef2d(iff,flag_fbase,"fbase","Cld base mass flux","kg/m2/s") 585 585 CALL histdef2d(iff,flag_prw,"prw","Precipitable water","kg/m2") 586 CALL histdef2d(iff,flag_cape_max,"cape_max","CAPE max.", "J/kg") 587 CALL histdef3d(iff,flag_upwd,"upwd","saturated updraft", "kg/m2/s") 588 CALL histdef3d(iff,flag_Ma,"Ma","undilute adiab updraft","kg/m2/s") 589 CALL histdef3d(iff,flag_dnwd,"dnwd","saturated downdraft","kg/m2/s") 590 CALL histdef3d(iff,flag_dnwd0,"dnwd0","unsat. downdraft", "kg/m2/s") 586 591 ENDIF !iflag_con .GE. 3 587 592 … … 632 637 CALL histdef2d(iff,flag_slab_bils, "slab_wbils_oce","Bilan au sol sur ocean slab", "W/m2") 633 638 634 IF (iflag_con.GE.3) THEN 635 CALL histdef2d(iff,flag_ale_bl,"ale_bl","ALE BL","m2/s2") 636 CALL histdef2d(iff,flag_alp_bl,"alp_bl","ALP BL","m2/s2") 637 CALL histdef2d(iff,flag_ale_wk,"ale_wk","ALE WK","m2/s2") 638 CALL histdef2d(iff,flag_alp_wk,"alp_wk","ALP WK","m2/s2") 639 639 IF (iflag_con.EQ.3) THEN 640 640 CALL histdef2d(iff,flag_ale,"ale","ALE","m2/s2") 641 641 CALL histdef2d(iff,flag_alp,"alp","ALP","W/m2") 642 642 CALL histdef2d(iff,flag_cin,"cin","Convective INhibition","m2/s2") 643 CALL histdef2d(iff,flag_wape,"WAPE","WAPE","m2/s2") 644 ENDIF !(iflag_con.GE.3) 643 IF (iflag_coupl.EQ.1) THEN 644 CALL histdef2d(iff,flag_ale_bl,"ale_bl","ALE BL","m2/s2") 645 CALL histdef2d(iff,flag_alp_bl,"alp_bl","ALP BL","m2/s2") 646 ENDIF !iflag_coupl.EQ.1 647 IF (iflag_wake.EQ.1) THEN 648 CALL histdef2d(iff,flag_ale_wk,"ale_wk","ALE WK","m2/s2") 649 CALL histdef2d(iff,flag_alp_wk,"alp_wk","ALP WK","m2/s2") 650 CALL histdef2d(iff,flag_wape,"WAPE","WAPE","m2/s2") 651 CALL histdef2d(iff,flag_wake_h,"wake_h","wake_h", "-") 652 CALL histdef2d(iff,flag_wake_s,"wake_s","wake_s", "-") 653 CALL histdef3d(iff,flag_dtwak,"dtwak","Wake dT","K/s") 654 CALL histdef3d(iff,flag_dqwak,"dqwak","Wake dQ","(kg/kg)/s") 655 CALL histdef3d(iff,flag_wake_deltat,"wake_deltat","wake_deltat", " ") 656 CALL histdef3d(iff,flag_wake_deltaq,"wake_deltaq","wake_deltaq", " ") 657 CALL histdef3d(iff,flag_wake_omg,"wake_omg","wake_omg", " ") 658 CALL histdef3d(iff,flag_ftd,"ftd","tend temp due aux descentes precip","-") 659 CALL histdef3d(iff,flag_fqd,"fqd","tend vap eau due aux descentes precip","-") 660 ENDIF !iflag_wake.EQ.1 661 CALL histdef3d(iff,flag_Vprecip,"Vprecip","precipitation vertical profile","-") 662 ENDIF !(iflag_con.EQ.3) 645 663 646 664 CALL histdef2d(iff,flag_weakinv, "weakinv","Weak inversion", "-") … … 654 672 CALL histdef2d(iff,flag_LWdownOR,"LWdownOR","Sfce incident LW radiation OR","W/m2") 655 673 CALL histdef2d(iff,flag_snowl,"snowl","Solid Large-scale Precip.","kg/(m2*s)") 656 IF (iflag_con.GE.3) THEN657 CALL histdef2d(iff,flag_cape_max,"cape_max","CAPE max.", "J/kg")658 ENDIF !(iflag_con.GE.3)659 674 CALL histdef2d(iff,flag_solldown,"solldown","Down. IR rad. at surface","W/m2") 660 675 CALL histdef2d(iff,flag_dtsvdfo,"dtsvdfo","Boundary-layer dTs(o)","K/s") … … 679 694 CALL histdef3d(iff,flag_rhum,"rhum","Relative humidity","-") 680 695 CALL histdef3d(iff,flag_ozone,"ozone","Ozone concentration", "ppmv") 681 IF(iflag_con.GE.3) THEN682 CALL histdef3d(iff,flag_upwd,"upwd","saturated updraft", "kg/m2/s")683 ENDIF !iflag_con.GE.3684 696 CALL histdef3d(iff,flag_dtphy,"dtphy","Physics dT","K/s") 685 697 CALL histdef3d(iff,flag_dqphy,"dqphy","Physics dQ", "(kg/kg)/s") … … 701 713 CALL histdef3d(iff,flag_kz_max,"kz_max","Kz melange max","m2/s" ) 702 714 CALL histdef3d(iff,flag_clwcon,"clwcon","Convective Cloud Liquid water content", "kg/kg") 703 IF(iflag_con.GE.3) THEN704 CALL histdef3d(iff,flag_Ma,"Ma","undilute adiab updraft","kg/m2/s")705 CALL histdef3d(iff,flag_dnwd,"dnwd","saturated downdraft","kg/m2/s")706 CALL histdef3d(iff,flag_dnwd0,"dnwd0","unsat. downdraft", "kg/m2/s")707 ENDIF !(iflag_con.GE.3)708 715 CALL histdef3d(iff,flag_dtdyn,"dtdyn","Dynamics dT","K/s") 709 716 CALL histdef3d(iff,flag_dqdyn,"dqdyn","Dynamics dQ", "(kg/kg)/s") … … 713 720 CALL histdef3d(iff,flag_ducon,"ducon","Convection du","m/s2") 714 721 CALL histdef3d(iff,flag_dqcon,"dqcon","Convection dQ", "(kg/kg)/s") 715 IF(iflag_con.GE.3) THEN716 IF (iflag_wake == 1) THEN717 CALL histdef2d(iff,flag_wake_h,"wake_h","wake_h", "-")718 CALL histdef2d(iff,flag_wake_s,"wake_s","wake_s", "-")719 CALL histdef3d(iff,flag_dtwak,"dtwak","Wake dT","K/s")720 CALL histdef3d(iff,flag_dqwak,"dqwak","Wake dQ","(kg/kg)/s")721 CALL histdef3d(iff,flag_wake_deltat,"wake_deltat","wake_deltat", " ")722 CALL histdef3d(iff,flag_wake_deltaq,"wake_deltaq","wake_deltaq", " ")723 CALL histdef3d(iff,flag_wake_omg,"wake_omg","wake_omg", " ")724 ENDIF725 CALL histdef3d(iff,flag_Vprecip,"Vprecip","precipitation vertical profile","-")726 CALL histdef3d(iff,flag_ftd,"ftd","tend temp due aux descentes precip","-")727 CALL histdef3d(iff,flag_fqd,"fqd","tend vap eau due aux descentes precip","-")728 ENDIF !(iflag_con.GE.3)729 722 730 723 CALL histdef3d(iff,flag_dtlsc,"dtlsc","Condensation dT", "K/s") -
LMDZ4/trunk/libf/phylmd/phys_output_write.h
r1067 r1100 510 510 CALL histwrite_phy(nid_files(iff),"prw",itau_w,prw) 511 511 ENDIF 512 !!! IM beg 513 IF (flag_cape_max(iff)<=lev_files(iff)) THEN 514 CALL histwrite_phy(nid_files(iff),"cape_max",itau_w,cape) 515 ENDIF 516 IF (flag_upwd(iff)<=lev_files(iff)) THEN 517 CALL histwrite_phy(nid_files(iff),"upwd",itau_w,upwd) 518 ENDIF 519 IF (flag_Ma(iff)<=lev_files(iff)) THEN 520 CALL histwrite_phy(nid_files(iff),"Ma",itau_w,Ma) 521 ENDIF 522 523 IF (flag_dnwd(iff)<=lev_files(iff)) THEN 524 CALL histwrite_phy(nid_files(iff),"dnwd",itau_w,dnwd) 525 ENDIF 526 527 IF (flag_dnwd0(iff)<=lev_files(iff)) THEN 528 CALL histwrite_phy(nid_files(iff),"dnwd0",itau_w,dnwd0) 529 ENDIF 530 !!! IM end 512 531 ENDIF !iflag_con .GE. 3 513 532 … … 616 635 IF (type_ocean=='force ') THEN 617 636 618 IF (iflag_con.GE.3) THEN 637 IF (iflag_con.EQ.3) THEN 638 IF (flag_ale(iff)<=lev_files(iff)) THEN 639 CALL histwrite_phy(nid_files(iff),"ale",itau_w,ale) 640 ENDIF 641 IF (flag_alp(iff)<=lev_files(iff)) THEN 642 CALL histwrite_phy(nid_files(iff),"alp",itau_w,alp) 643 ENDIF 644 IF (flag_cin(iff)<=lev_files(iff)) THEN 645 CALL histwrite_phy(nid_files(iff),"cin",itau_w,cin) 646 ENDIF 619 647 IF (iflag_coupl.EQ.1) THEN 620 648 IF (flag_ale_bl(iff)<=lev_files(iff)) THEN … … 634 662 ENDIF 635 663 636 IF (flag_ale(iff)<=lev_files(iff)) THEN637 CALL histwrite_phy(nid_files(iff),"ale",itau_w,ale)638 ENDIF639 IF (flag_alp(iff)<=lev_files(iff)) THEN640 CALL histwrite_phy(nid_files(iff),"alp",itau_w,alp)641 ENDIF642 IF (flag_cin(iff)<=lev_files(iff)) THEN643 CALL histwrite_phy(nid_files(iff),"cin",itau_w,cin)644 ENDIF645 664 IF (flag_wape(iff)<=lev_files(iff)) THEN 646 665 CALL histwrite_phy(nid_files(iff),"WAPE",itau_w,wake_pe) … … 667 686 CALL histwrite_phy(nid_files(iff),"wake_omg",itau_w,wake_omg) 668 687 ENDIF 669 ENDIF 670 671 IF (flag_Vprecip(iff)<=lev_files(iff)) THEN 672 CALL histwrite_phy(nid_files(iff),"Vprecip",itau_w,Vprecip) 688 !!!IM beg 689 IF (flag_dtwak(iff)<=lev_files(iff)) THEN 690 zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev) 691 & /pdtphys 692 CALL histwrite_phy(nid_files(iff),"dtwak",itau_w,zx_tmp_fi3d) 693 ENDIF 694 695 IF (flag_dqwak(iff)<=lev_files(iff)) THEN 696 zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys 697 CALL histwrite_phy(nid_files(iff),"dqwak",itau_w,zx_tmp_fi3d) 673 698 ENDIF 674 699 … … 680 705 CALL histwrite_phy(nid_files(iff),"fqd",itau_w,fqd) 681 706 ENDIF 682 683 ENDIF !(iflag_con.GE.3) 707 !!!IM end 708 ENDIF 709 710 IF (flag_Vprecip(iff)<=lev_files(iff)) THEN 711 CALL histwrite_phy(nid_files(iff),"Vprecip",itau_w,Vprecip) 712 ENDIF 713 714 ENDIF !(iflag_con.EQ.3) 684 715 685 716 ELSE IF (type_ocean=='slab ') THEN … … 759 790 CALL histwrite_phy(nid_files(iff),"snowl",itau_w,snow_lsc) 760 791 ENDIF 761 762 IF (iflag_con.GE.3) THEN763 IF (flag_cape_max(iff)<=lev_files(iff)) THEN764 CALL histwrite_phy(nid_files(iff),"cape_max",itau_w,cape)765 ENDIF766 ENDIF !(iflag_con.GE.3)767 792 768 793 IF (flag_solldown(iff)<=lev_files(iff)) THEN … … 854 879 CALL histwrite_phy(nid_files(iff),"ozone",itau_w,zx_tmp_fi3d) 855 880 ENDIF 856 857 IF(iflag_con.GE.3) THEN858 IF (flag_upwd(iff)<=lev_files(iff)) THEN859 CALL histwrite_phy(nid_files(iff),"upwd",itau_w,upwd)860 ENDIF861 ENDIF !(iflag_con.GE.3)862 881 863 882 IF (flag_dtphy(iff)<=lev_files(iff)) THEN … … 933 952 ENDIF 934 953 935 IF(iflag_con.GE.3) THEN936 IF (flag_Ma(iff)<=lev_files(iff)) THEN937 CALL histwrite_phy(nid_files(iff),"Ma",itau_w,Ma)938 ENDIF939 940 IF (flag_dnwd(iff)<=lev_files(iff)) THEN941 CALL histwrite_phy(nid_files(iff),"dnwd",itau_w,dnwd)942 ENDIF943 944 IF (flag_dnwd0(iff)<=lev_files(iff)) THEN945 CALL histwrite_phy(nid_files(iff),"dnwd0",itau_w,dnwd0)946 ENDIF947 ENDIF !(iflag_con.GE.3)948 949 954 IF (flag_dtdyn(iff)<=lev_files(iff)) THEN 950 955 CALL histwrite_phy(nid_files(iff),"dtdyn",itau_w,d_t_dyn) … … 977 982 CALL histwrite_phy(nid_files(iff),"dqcon",itau_w,zx_tmp_fi3d) 978 983 ENDIF 979 980 IF(iflag_con.GE.3) THEN981 IF (iflag_wake == 1) THEN982 IF (flag_dtwak(iff)<=lev_files(iff)) THEN983 zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev) &984 & /pdtphys985 CALL histwrite_phy(nid_files(iff),"dtwak",itau_w,zx_tmp_fi3d)986 ENDIF987 988 IF (flag_dqwak(iff)<=lev_files(iff)) THEN989 zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys990 CALL histwrite_phy(nid_files(iff),"dqwak",itau_w,zx_tmp_fi3d)991 ENDIF992 ENDIF993 ENDIF !(iflag_con.GE.3)994 984 995 985 IF (flag_dtlsc(iff)<=lev_files(iff)) THEN
Note: See TracChangeset
for help on using the changeset viewer.