Ignore:
Timestamp:
Feb 16, 2009, 7:01:12 PM (15 years ago)
Author:
musat
Message:

Homogeneisation flags sorties convection + poches froides
IM

File:
1 edited

Legend:

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

    r1067 r1100  
    510510      CALL histwrite_phy(nid_files(iff),"prw",itau_w,prw)
    511511        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
    512531      ENDIF !iflag_con .GE. 3
    513532
     
    616635      IF (type_ocean=='force ') THEN
    617636
    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
    619647      IF (iflag_coupl.EQ.1) THEN
    620648       IF (flag_ale_bl(iff)<=lev_files(iff)) THEN
     
    634662       ENDIF
    635663
    636        IF (flag_ale(iff)<=lev_files(iff)) THEN
    637        CALL histwrite_phy(nid_files(iff),"ale",itau_w,ale)
    638        ENDIF
    639        IF (flag_alp(iff)<=lev_files(iff)) THEN
    640        CALL histwrite_phy(nid_files(iff),"alp",itau_w,alp)
    641        ENDIF
    642        IF (flag_cin(iff)<=lev_files(iff)) THEN
    643        CALL histwrite_phy(nid_files(iff),"cin",itau_w,cin)
    644        ENDIF
    645664       IF (flag_wape(iff)<=lev_files(iff)) THEN
    646665       CALL histwrite_phy(nid_files(iff),"WAPE",itau_w,wake_pe)
     
    667686       CALL histwrite_phy(nid_files(iff),"wake_omg",itau_w,wake_omg)
    668687        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)
    673698        ENDIF
    674699
     
    680705       CALL histwrite_phy(nid_files(iff),"fqd",itau_w,fqd)
    681706        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)
    684715 
    685716      ELSE IF (type_ocean=='slab  ') THEN
     
    759790      CALL histwrite_phy(nid_files(iff),"snowl",itau_w,snow_lsc)
    760791      ENDIF
    761 
    762       IF (iflag_con.GE.3) THEN
    763       IF (flag_cape_max(iff)<=lev_files(iff)) THEN
    764       CALL histwrite_phy(nid_files(iff),"cape_max",itau_w,cape)
    765       ENDIF
    766       ENDIF !(iflag_con.GE.3)
    767792
    768793      IF (flag_solldown(iff)<=lev_files(iff)) THEN
     
    854879      CALL histwrite_phy(nid_files(iff),"ozone",itau_w,zx_tmp_fi3d)
    855880      ENDIF
    856 
    857       IF(iflag_con.GE.3) THEN
    858        IF (flag_upwd(iff)<=lev_files(iff)) THEN
    859       CALL histwrite_phy(nid_files(iff),"upwd",itau_w,upwd)
    860        ENDIF
    861       ENDIF !(iflag_con.GE.3)
    862881
    863882       IF (flag_dtphy(iff)<=lev_files(iff)) THEN
     
    933952       ENDIF
    934953
    935       IF(iflag_con.GE.3) THEN
    936        IF (flag_Ma(iff)<=lev_files(iff)) THEN
    937       CALL histwrite_phy(nid_files(iff),"Ma",itau_w,Ma)
    938        ENDIF
    939 
    940        IF (flag_dnwd(iff)<=lev_files(iff)) THEN
    941       CALL histwrite_phy(nid_files(iff),"dnwd",itau_w,dnwd)
    942        ENDIF
    943 
    944        IF (flag_dnwd0(iff)<=lev_files(iff)) THEN
    945       CALL histwrite_phy(nid_files(iff),"dnwd0",itau_w,dnwd0)
    946        ENDIF
    947       ENDIF !(iflag_con.GE.3)
    948 
    949954       IF (flag_dtdyn(iff)<=lev_files(iff)) THEN
    950955      CALL histwrite_phy(nid_files(iff),"dtdyn",itau_w,d_t_dyn)
     
    977982      CALL histwrite_phy(nid_files(iff),"dqcon",itau_w,zx_tmp_fi3d)
    978983       ENDIF
    979 
    980       IF(iflag_con.GE.3) THEN
    981        IF (iflag_wake == 1) THEN
    982          IF (flag_dtwak(iff)<=lev_files(iff)) THEN
    983            zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev)           &
    984      &                                        /pdtphys
    985            CALL histwrite_phy(nid_files(iff),"dtwak",itau_w,zx_tmp_fi3d)
    986          ENDIF
    987 
    988         IF (flag_dqwak(iff)<=lev_files(iff)) THEN
    989         zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
    990         CALL histwrite_phy(nid_files(iff),"dqwak",itau_w,zx_tmp_fi3d)
    991         ENDIF
    992        ENDIF
    993       ENDIF !(iflag_con.GE.3)
    994984
    995985       IF (flag_dtlsc(iff)<=lev_files(iff)) THEN
Note: See TracChangeset for help on using the changeset viewer.