Ignore:
Timestamp:
Mar 11, 2009, 3:45:10 PM (15 years ago)
Author:
musat
Message:

Homogeneisation flags sorties convection + poches froides
IM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_write.h

    r1114 r1119  
    511511      CALL histwrite_phy(nid_files(iff),o_prw%name,itau_w,prw)
    512512        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
    513530      ENDIF !iflag_con .GE. 3
    514531
     
    619636      IF (type_ocean=='force ') THEN
    620637
    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
    622648      IF (iflag_coupl.EQ.1) THEN
    623649       IF (o_ale_bl%flag(iff)<=lev_files(iff)) THEN
     
    628654       ENDIF
    629655      ENDIF !iflag_coupl.EQ.1
    630       ENDIF !(iflag_con.GE.3)
    631       IF (iflag_con.EQ.3) THEN
     656
    632657      IF (iflag_wake.EQ.1) THEN
    633658       IF (o_ale_wk%flag(iff)<=lev_files(iff)) THEN
     
    638663       ENDIF
    639664
    640        IF (o_ale%flag(iff)<=lev_files(iff)) THEN
    641        CALL histwrite_phy(nid_files(iff),o_ale%name,itau_w,ale)
    642        ENDIF
    643        IF (o_alp%flag(iff)<=lev_files(iff)) THEN
    644        CALL histwrite_phy(nid_files(iff),o_alp%name,itau_w,alp)
    645        ENDIF
    646        IF (o_cin%flag(iff)<=lev_files(iff)) THEN
    647        CALL histwrite_phy(nid_files(iff),o_cin%name,itau_w,cin)
    648        ENDIF
    649665       IF (o_wape%flag(iff)<=lev_files(iff)) THEN
    650666       CALL histwrite_phy(nid_files(iff),o_WAPE%name,itau_w,wake_pe)
     
    669685
    670686        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)
    678701        ENDIF
    679702
     
    685708       CALL histwrite_phy(nid_files(iff),o_fqd%name,itau_w,fqd)
    686709        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
    687715
    688716      ENDIF !(iflag_con.EQ.3)
     
    690718      ELSE IF (type_ocean=='slab  ') THEN
    691719
    692       IF ( o_slab_bils%flag(iff)<=lev_files(iff))
     720       IF ( o_slab_bils%flag(iff)<=lev_files(iff))
    693721     $     CALL histwrite_phy(
    694722     $     nid_files(iff),o_slab_bils%name,itau_w,slab_wfbils)
    695 
     723     
    696724      ENDIF !type_ocean == force/slab
    697725
     
    768796      CALL histwrite_phy(nid_files(iff),o_snowl%name,itau_w,snow_lsc)
    769797      ENDIF
    770 
    771       IF (iflag_con.GE.3) THEN
    772       IF (o_cape_max%flag(iff)<=lev_files(iff)) THEN
    773       CALL histwrite_phy(nid_files(iff),o_cape_max%name,itau_w,cape)
    774       ENDIF
    775       ENDIF !(iflag_con.GE.3)
    776798
    777799      IF (o_solldown%flag(iff)<=lev_files(iff)) THEN
     
    885907      CALL histwrite_phy(nid_files(iff),o_ozone%name,itau_w,zx_tmp_fi3d)
    886908      ENDIF
    887 
    888       IF(iflag_con.GE.3) THEN
    889        IF (o_upwd%flag(iff)<=lev_files(iff)) THEN
    890       CALL histwrite_phy(nid_files(iff),o_upwd%name,itau_w,upwd)
    891        ENDIF
    892       ENDIF !(iflag_con.GE.3)
    893909
    894910       IF (o_dtphy%flag(iff)<=lev_files(iff)) THEN
     
    963979       ENDIF
    964980
    965       IF(iflag_con.GE.3) THEN
    966        IF (o_Ma%flag(iff)<=lev_files(iff)) THEN
    967       CALL histwrite_phy(nid_files(iff),o_Ma%name,itau_w,Ma)
    968        ENDIF
    969 
    970        IF (o_dnwd%flag(iff)<=lev_files(iff)) THEN
    971       CALL histwrite_phy(nid_files(iff),o_dnwd%name,itau_w,dnwd)
    972        ENDIF
    973 
    974        IF (o_dnwd0%flag(iff)<=lev_files(iff)) THEN
    975       CALL histwrite_phy(nid_files(iff),o_dnwd0%name,itau_w,dnwd0)
    976        ENDIF
    977       ENDIF !(iflag_con.GE.3)
    978 
    979981       IF (o_dtdyn%flag(iff)<=lev_files(iff)) THEN
    980982      CALL histwrite_phy(nid_files(iff),o_dtdyn%name,itau_w,d_t_dyn)
     
    10071009      CALL histwrite_phy(nid_files(iff),o_dqcon%name,itau_w,zx_tmp_fi3d)
    10081010       ENDIF
    1009 
    1010       IF(iflag_con.GE.3) THEN
    1011        IF (iflag_wake == 1) THEN
    1012          IF (o_dtwak%flag(iff)<=lev_files(iff)) THEN
    1013            zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)           
    1014      &                                        /pdtphys
    1015            CALL histwrite_phy(nid_files(iff),
    1016      &                       o_dtwak%name,itau_w,zx_tmp_fi3d)
    1017          ENDIF
    1018 
    1019         IF (o_dqwak%flag(iff)<=lev_files(iff)) THEN
    1020         zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
    1021         CALL histwrite_phy(nid_files(iff),
    1022      &                     o_dqwak%name,itau_w,zx_tmp_fi3d)
    1023         ENDIF
    1024        ENDIF
    1025       ENDIF !(iflag_con.GE.3)
    10261011
    10271012       IF (o_dtlsc%flag(iff)<=lev_files(iff)) THEN
Note: See TracChangeset for help on using the changeset viewer.