Changeset 2086 for LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Jul 9, 2014, 11:19:07 PM (10 years ago)

Author:

fhourdin

Message:

nclusion de la thermodynamique de la glace
Ice thermodynamics
(Catherine Rio)

File:

• LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -222 +222 @@
     USE ocean_slab_mod, only: tslab, slab_bils
     USE indice_sol_mod, only: nbsrf
-    USE infotrac, only: nqtot
+    USE infotrac, only: nqtot, nqo
     USE comgeomphy, only: airephy
     USE surface_data, only: type_ocean, ok_veget, ok_snow
@@ -1288 +1288 @@
        ENDDO !nfiles
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-       IF (nqtot.GE.3) THEN
-          DO iq=3,nqtot
-             CALL histwrite_phy(o_trac(iq-2), qx(:,:,iq))
-             CALL histwrite_phy(o_dtr_vdf(iq-2),d_tr_cl(:,:,iq-2))
-             CALL histwrite_phy(o_dtr_the(iq-2),d_tr_th(:,:,iq-2))
-             CALL histwrite_phy(o_dtr_con(iq-2),d_tr_cv(:,:,iq-2))
-             CALL histwrite_phy(o_dtr_lessi_impa(iq-2),d_tr_lessi_impa(:,:,iq-2))
-             CALL histwrite_phy(o_dtr_lessi_nucl(iq-2),d_tr_lessi_nucl(:,:,iq-2))
-             CALL histwrite_phy(o_dtr_insc(iq-2),d_tr_insc(:,:,iq-2))
-             CALL histwrite_phy(o_dtr_bcscav(iq-2),d_tr_bcscav(:,:,iq-2))
-             CALL histwrite_phy(o_dtr_evapls(iq-2),d_tr_evapls(:,:,iq-2))
-             CALL histwrite_phy(o_dtr_ls(iq-2),d_tr_ls(:,:,iq-2))
-             CALL histwrite_phy(o_dtr_trsp(iq-2),d_tr_trsp(:,:,iq-2))
-             CALL histwrite_phy(o_dtr_sscav(iq-2),d_tr_sscav(:,:,iq-2))
-             CALL histwrite_phy(o_dtr_sat(iq-2),d_tr_sat(:,:,iq-2))
-             CALL histwrite_phy(o_dtr_uscav(iq-2),d_tr_uscav(:,:,iq-2))
+!       IF (nqtot.GE.3) THEN
+!          DO iq=3,nqtot
+          IF (nqtot.GE.nqo+1) THEN
+            DO iq=nqo+1,nqtot
+             CALL histwrite_phy(o_trac(iq-nqo), qx(:,:,iq))
+             CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_con(iq-nqo),d_tr_cv(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_lessi_impa(iq-nqo),d_tr_lessi_impa(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_lessi_nucl(iq-nqo),d_tr_lessi_nucl(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_insc(iq-nqo),d_tr_insc(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_bcscav(iq-nqo),d_tr_bcscav(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_evapls(iq-nqo),d_tr_evapls(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_ls(iq-nqo),d_tr_ls(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_trsp(iq-nqo),d_tr_trsp(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_sscav(iq-nqo),d_tr_sscav(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_sat(iq-nqo),d_tr_sat(:,:,iq-nqo))
+             CALL histwrite_phy(o_dtr_uscav(iq-nqo),d_tr_uscav(:,:,iq-nqo))
              zx_tmp_fi2d=0.
              IF(vars_defined) THEN
@@ -1310 +1312 @@
                 ENDDO
              ENDIF
-             CALL histwrite_phy(o_trac_cum(iq-2), zx_tmp_fi2d)
+             CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
           ENDDO
        ENDIF