Ignore:
Timestamp:
Apr 18, 2008, 5:34:10 PM (17 years ago)
Author:
lmdzadmin
Message:

Correction petits bogues; on ecrit nouvelles variables sous flags appropries JG
IM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/phys_output_write.h

    r945 r951  
    652652      ENDIF
    653653
    654       IF (flag_ale_bl(iff)<=lev_files(iff)) THEN
    655       CALL histwrite_phy(nid_files(iff),"ale_bl",itau_w,ale_bl)
    656       ENDIF
    657       IF (flag_alp_bl(iff)<=lev_files(iff)) THEN
    658       CALL histwrite_phy(nid_files(iff),"alp_bl",itau_w,alp_bl)
    659       ENDIF
    660       IF (flag_ale_wk(iff)<=lev_files(iff)) THEN
    661       CALL histwrite_phy(nid_files(iff),"ale_wk",itau_w,ale_wake)
    662       ENDIF
    663       IF (flag_alp_wk(iff)<=lev_files(iff)) THEN
    664       CALL histwrite_phy(nid_files(iff),"alp_wk",itau_w,alp_wake)
     654      IF (iflag_coupl.EQ.1) THEN
     655       IF (flag_ale_bl(iff)<=lev_files(iff)) THEN
     656       CALL histwrite_phy(nid_files(iff),"ale_bl",itau_w,ale_bl)
     657       ENDIF
     658       IF (flag_alp_bl(iff)<=lev_files(iff)) THEN
     659       CALL histwrite_phy(nid_files(iff),"alp_bl",itau_w,alp_bl)
     660       ENDIF
     661      ENDIF !iflag_coupl.EQ.1
     662
     663      IF (iflag_wake.EQ.1) THEN
     664       IF (flag_ale_wk(iff)<=lev_files(iff)) THEN
     665       CALL histwrite_phy(nid_files(iff),"ale_wk",itau_w,ale_wake)
     666       ENDIF
     667       IF (flag_alp_wk(iff)<=lev_files(iff)) THEN
     668       CALL histwrite_phy(nid_files(iff),"alp_wk",itau_w,alp_wake)
     669       ENDIF
    665670      ENDIF
    666671
     
    709714      ENDIF
    710715
     716#ifdef histISCCP
     717      IF(ok_isccp) THEN
    711718       IF (flag_meantaucld(iff)<=lev_files(iff)) THEN
    712719       CALL histwrite_phy(nid_files(iff),"meantaucld",itau_w,meantaucld)
    713720       ENDIF
     721      ENDIF !(ok_isccp)
     722#endif
    714723
    715724       IF (flag_cldtau(iff)<=lev_files(iff)) THEN
     
    972981       ENDIF
    973982
    974        IF (flag_dtwak(iff)<=lev_files(iff)) THEN
    975       zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)/pdtphys
    976       CALL histwrite_phy(nid_files(iff),"dtwak",itau_w,zx_tmp_fi3d)
    977        ENDIF
    978 
    979        IF (flag_dqwak(iff)<=lev_files(iff)) THEN
    980       zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
    981       CALL histwrite_phy(nid_files(iff),"dqwak",itau_w,zx_tmp_fi3d)
     983       IF (iflag_wake.ge.1) THEN
     984        IF (flag_dtwak(iff)<=lev_files(iff)) THEN
     985        zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)/pdtphys
     986        CALL histwrite_phy(nid_files(iff),"dtwak",itau_w,zx_tmp_fi3d)
     987        ENDIF
     988
     989        IF (flag_dqwak(iff)<=lev_files(iff)) THEN
     990        zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
     991        CALL histwrite_phy(nid_files(iff),"dqwak",itau_w,zx_tmp_fi3d)
     992        ENDIF
    982993       ENDIF
    983994
     
    10331044       ENDIF
    10341045
    1035        IF (flag_f_th(iff)<=lev_files(iff)) THEN
    1036       CALL histwrite_phy(nid_files(iff),"f_th",itau_w,fm_therm)
    1037        ENDIF
    1038 
    1039        IF (flag_e_th(iff)<=lev_files(iff)) THEN
    1040       CALL histwrite_phy(nid_files(iff),"e_th",itau_w,entr_therm)
    1041        ENDIF
    1042 
    1043        IF (flag_d_th(iff)<=lev_files(iff)) THEN
    1044       CALL histwrite_phy(nid_files(iff),"d_th",itau_w,detr_therm)
    1045        ENDIF
     1046       IF (iflag_thermals.ge.1) THEN
     1047        IF (flag_f_th(iff)<=lev_files(iff)) THEN
     1048        CALL histwrite_phy(nid_files(iff),"f_th",itau_w,fm_therm)
     1049        ENDIF
     1050
     1051        IF (flag_e_th(iff)<=lev_files(iff)) THEN
     1052        CALL histwrite_phy(nid_files(iff),"e_th",itau_w,entr_therm)
     1053        ENDIF
     1054
     1055!       IF (flag_d_th(iff)<=lev_files(iff)) THEN
     1056!       CALL histwrite_phy(nid_files(iff),"d_th",itau_w,detr_therm)
     1057!       ENDIF
     1058       ENDIF !iflag_thermals
    10461059
    10471060       IF (flag_f0_th(iff)<=lev_files(iff)) THEN
Note: See TracChangeset for help on using the changeset viewer.