Changeset 5903 for LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

Timestamp:

Dec 3, 2025, 6:01:05 PM (3 hours ago)

Author:

fcheruy

Message:

correction les flux spectraux

File:

• LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -258 +258 @@
          o_lwtoa0b2165, o_lwtoa0b2315, o_lwtoa0b2490, o_lwtoa0b2800, &
          o_lwtoabtot, o_lwtoa0btot, &  !Somme de tous les canaux
-         o_lwtoa0b, &  ! les 16 canaux FC
          o_lwtoab4co2,o_lwtoacsb4co2,                            &   ! les 16 canaux FC
-!FC
 !GG         o_zxfluxt,o_zxfluxq
          o_fcds, o_fcdi, o_dh_basal_growth, o_dh_basal_melt, &
@@ -512 +510 @@
         iflag_seaice, iflag_seaice_alb
     !GG
-    USE aero_mod, ONLY: naero_tot, id_STRAT_phy
+    USE aero_mod, ONLY: naero_tot, id_STRAT_phy,nbands_lw_rrtm
     USE ioipsl, ONLY: histend, histsync
     USE iophy, ONLY: set_itau_iophy, histwrite_phy
@@ -997 +995 @@
 !FC
       IF (vars_defined) THEN
-                DO i=1, klon
-             zx_tmp_fi2d(i)=0.
-          DO ll=1,nlw
-!             zx_tmp_fi2d(i)=zx_tmp_fi2d(i)+lwtoab(i,ll)*deltanu(ll)
-             zx_tmp_fi2d(i)=zx_tmp_fi2d(i)+lwtoab(i,ll)
+!                DO i=1, klon
+             zx_tmp_fi2d(:)=0.
+!          DO ll=1,nlw
+          DO ll=1,nbands_lw_rrtm
+             zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+lwtoab(:,ll)
           ENDDO
-          ENDDO
+!          ENDDO
        ENDIF
        CALL histwrite_phy(o_lwtoabtot, zx_tmp_fi2d)
            IF (vars_defined) THEN
-                DO i=1, klon
-             zx_tmp_fi2d(i)=0.
-          DO ll=1,nlw
-             zx_tmp_fi2d(i)=zx_tmp_fi2d(i)+lwtoa0b(i,ll)
+!                DO i=1, klon
+             zx_tmp_fi2d(:)=0.
+!          DO ll=1,nlw
+          DO ll=1,nbands_lw_rrtm
+             zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+lwtoa0b(:,ll)
           ENDDO
-          ENDDO
+!          ENDDO
        ENDIF
        CALL histwrite_phy(o_lwtoa0btot, zx_tmp_fi2d)
@@ -1018 +1017 @@
 
        IF (vars_defined) zx_tmp_fi2d=lwtoab(:,1)
-!       CALL histwrite_phy(o_lwtoab, zx_tmp_fi3dsplw)
        CALL histwrite_phy(o_lwtoab130, zx_tmp_fi2d)
 
@@ -1864 +1862 @@
           ENDIF
                CALL histwrite_phy(o_delta_tsurf,zx_tmp_fi2d) 
+
+
                CALL histwrite_phy(o_cdragh_x   ,cdragh_x   ) 
                CALL histwrite_phy(o_cdragh_w   ,cdragh_w   )