Changeset 1646 for LMDZ5


Ignore:
Timestamp:
Jul 31, 2012, 5:38:35 PM (12 years ago)
Author:
lguez
Message:

When "iflag_con == 2", write :

-- "pmfd" (instead of "- dnwd - dnwd0" with "dnwd" undefined) into

NetCDF variable "mcd";

-- "pmfu + pmfd" (instead of "upwd + dnwd + dnwd0" with "upwd" and

"dnwd" undefined) into NetCDF variable "dmc".

When "iflag_con < 2", do no write anything into "mcd" and "dmc". No
change when "iflag_con >= 3".

File:
1 edited

Legend:

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

    r1638 r1646  
    21592159       ENDIF
    21602160
    2161        IF (o_mcd%flag(iff)<=lev_files(iff)) THEN
    2162       zx_tmp_fi3d(1:klon,1:klev)=-1 * (dnwd(1:klon,1:klev)+
    2163      $                                 dnwd0(1:klon,1:klev))
    2164       CALL histwrite_phy(nid_files(iff),clef_stations(iff),
    2165      $o_mcd%name,itau_w,zx_tmp_fi3d)
    2166        ENDIF
    2167 
    2168        IF (o_dmc%flag(iff)<=lev_files(iff)) THEN
    2169       zx_tmp_fi3d(1:klon,1:klev)=upwd(1:klon,1:klev) +
    2170      $  dnwd(1:klon,1:klev)+ dnwd0(1:klon,1:klev)
    2171       CALL histwrite_phy(nid_files(iff),clef_stations(iff),
    2172      $o_dmc%name,itau_w,zx_tmp_fi3d)
    2173        ENDIF
     2161       if (iflag_con >= 3) then
     2162          IF (o_mcd%flag(iff)<=lev_files(iff)) THEN
     2163             zx_tmp_fi3d(1:klon,1:klev)=-1 * (dnwd(1:klon,1:klev)+
     2164     $            dnwd0(1:klon,1:klev))
     2165             CALL histwrite_phy(nid_files(iff),clef_stations(iff),
     2166     $            o_mcd%name,itau_w,zx_tmp_fi3d)
     2167          ENDIF
     2168
     2169          IF (o_dmc%flag(iff)<=lev_files(iff)) THEN
     2170             zx_tmp_fi3d(1:klon,1:klev)=upwd(1:klon,1:klev) +
     2171     $            dnwd(1:klon,1:klev)+ dnwd0(1:klon,1:klev)
     2172             CALL histwrite_phy(nid_files(iff),clef_stations(iff),
     2173     $            o_dmc%name,itau_w,zx_tmp_fi3d)
     2174          ENDIF
     2175       else if (iflag_con == 2) then
     2176          IF (o_mcd%flag(iff) <= lev_files(iff)) THEN
     2177             CALL histwrite_phy(nid_files(iff), clef_stations(iff),
     2178     $            o_mcd%name, itau_w, pmfd)
     2179          ENDIF
     2180
     2181          IF (o_dmc%flag(iff) <= lev_files(iff)) THEN
     2182             CALL histwrite_phy(nid_files(iff), clef_stations(iff),
     2183     $            o_dmc%name, itau_w, pmfu + pmfd)
     2184          ENDIF
     2185       end if
    21742186
    21752187       IF (o_ref_liq%flag(iff)<=lev_files(iff)) THEN
Note: See TracChangeset for help on using the changeset viewer.