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

Timestamp:

Dec 10, 2013, 1:24:44 PM (11 years ago)

Author:

musat

Message:

D'autres oublis pour les sorties annuelles
IM

File:

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

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -18 +18 @@
        pplay, lmax_th, aerosol_couple,         &
        ok_ade, ok_aie, ivap, new_aod, ok_sync, &
-       ptconv, read_climoz, clevSTD, freq_moyNMC, ptconvth, &
+       ptconv, read_climoz, clevSTD, ptconvth, &
        d_t, qx, d_qx, zmasse, flag_aerosol_strat)
 
@@ -230 +230 @@
     USE wxios, only: wxios_update_calendar, wxios_closedef
 #endif
+    USE phys_cal_mod, only : mth_len
 
     IMPLICIT NONE
@@ -269 +270 @@
     REAL, PARAMETER :: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
     REAL, PARAMETER :: missing_val=nf90_fill_real
+    REAL, PARAMETER :: un_jour=86400.
 
     ! On calcul le nouveau tau:
@@ -1194 +1196 @@
                 DO i=1, klon
                    IF(tnondef(i,k,iff-6).NE.missing_val) THEN
+                      IF(freq_outNMC(iff-6).LT.0) THEN
+                         freq_moyNMC(iff-6)=(mth_len*un_jour)/freq_calNMC(iff-6)
+                      ELSE
+                         freq_moyNMC(iff-6)=freq_outNMC(iff-6)/freq_calNMC(iff-6)
+                      ENDIF
                       zx_tmp_fi3d(i,k) = (100.*tnondef(i,k,iff-6))/freq_moyNMC(iff-6)
                    ELSE