Changeset 5709 for LMDZ6/trunk/libf

Timestamp:

Jun 16, 2025, 5:26:50 PM (7 weeks ago)

Author:

fcheruy

Message:

to work with xios et output spectral clear sky flux

Location:

LMDZ6/trunk/libf/phylmd

Files:

• phys_output_ctrlout_mod.F90 (modified) (1 diff)
• phys_output_write_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -487 +487 @@
   TYPE(ctrl_out), SAVE :: o_lwtoab130 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab130', 'LW TOA  band 130cm-1 ', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b130 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b130', 'LW TOA  band 130cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab375 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab375', 'LW TOA  band 375cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b375 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b375', 'LW TOA  band 375cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab565 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab565', 'LW TOA  band 565cm-1 ', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b565 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b565', 'LW TOA  band 565cm-1 CS ', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab665 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab665', 'LW TOA  band 665cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b665 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b665', 'LW TOA  band 665cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab760 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab760', 'LW TOA  band 760cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b760 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b760', 'LW TOA  band 760cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab900 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab900', 'LW TOA  band 900cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b900 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b900', 'LW TOA  band 900cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab1030 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab1030', 'LW TOA  band 1030cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b1030 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b1030', 'LW TOA  band 1030cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab1130 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab1130', 'LW TOA  band 1130cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b1130 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b1130', 'LW TOA  band 1130cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab1285 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab1285', 'LW TOA  band 1285cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b1285 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b1285', 'LW TOA  band 1285cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab1435 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab1435', 'LW TOA  band 1435cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b1435 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b1435', 'LW TOA  band 1435cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab1640 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab1640', 'LW TOA  band 1640cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b1640 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b1640', 'LW TOA  band 1640cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab1940 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab1940', 'LW TOA  band 1940cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b1940 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b1940', 'LW TOA  band 1940cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab2165 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab2165', 'LW TOA  band 2165cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b2165 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b2165', 'LW TOA  band 2165cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab2315 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab2315', 'LW TOA  band 2315cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b2315 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b2315', 'LW TOA  band 2315cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab2490 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab2490', 'LW TOA  band 2490cm-1', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b2490 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b2490', 'LW TOA  band 2490cm-1 CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoab2800 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoab2800', 'LW TOA  band 2800cm-1 ', 'W m-2', (/ ('', i=1, 10) /))
+  TYPE(ctrl_out), SAVE :: o_lwtoa0b2800 = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0b2800', 'LW TOA  band 2800cm-1 CS ', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoabtot = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoabtot', 'LW TOA  band total', 'W m-2', (/ ('', i=1, 10) /))
+ TYPE(ctrl_out), SAVE :: o_lwtoa0btot = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'lwtoa0btot', 'LW TOA  band total CS', 'W m-2', (/ ('', i=1, 10) /))
   TYPE(ctrl_out), SAVE :: o_lwtoa0b = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
     'lwtoa0b', 'LW TOA CS band ', 'W m-2', (/ ('', i=1, 10) /))

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -252 +252 @@
          o_lwtoab1285, o_lwtoab1435, o_lwtoab1640, o_lwtoab1940, &
          o_lwtoab2165, o_lwtoab2315, o_lwtoab2490, o_lwtoab2800, &
-         o_lwtoabtot, &  !Somme de tous les canaux
+         o_lwtoa0b130, o_lwtoa0b375 ,o_lwtoa0b565, o_lwtoa0b665, &
+         o_lwtoa0b760, o_lwtoa0b900, o_lwtoa0b1030, o_lwtoa0b1130, &  ! les 16 canaux FC
+         o_lwtoa0b1285, o_lwtoa0b1435, o_lwtoa0b1640, o_lwtoa0b1940, &
+         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
@@ -1001 +1005 @@
        ENDIF
        CALL histwrite_phy(o_lwtoabtot, zx_tmp_fi2d)
-
-   IF (vars_defined) zx_tmp_fi2d=lwtoab(:,1)
+           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)
+          ENDDO
+          ENDDO
+       ENDIF
+       CALL histwrite_phy(o_lwtoa0btot, zx_tmp_fi2d)
+
+
+       IF (vars_defined) zx_tmp_fi2d=lwtoab(:,1)
 !       CALL histwrite_phy(o_lwtoab, zx_tmp_fi3dsplw)
        CALL histwrite_phy(o_lwtoab130, zx_tmp_fi2d)
@@ -1051 +1065 @@
        CALL histwrite_phy(o_lwtoab2800, zx_tmp_fi2d)
 
-       IF (vars_defined) zx_tmp_fi3dsplw=lwtoa0b
-       CALL histwrite_phy(o_lwtoa0b, zx_tmp_fi3dsplw)
+!FC
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,1)
+       CALL histwrite_phy(o_lwtoa0b130, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,2)
+       CALL histwrite_phy(o_lwtoa0b375, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,3)
+       CALL histwrite_phy(o_lwtoa0b565, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,4)
+       CALL histwrite_phy(o_lwtoa0b665, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,5)
+       CALL histwrite_phy(o_lwtoa0b760, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,6)
+       CALL histwrite_phy(o_lwtoa0b900, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,7)
+       CALL histwrite_phy(o_lwtoa0b1030, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,8)
+       CALL histwrite_phy(o_lwtoa0b1130, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,9)
+       CALL histwrite_phy(o_lwtoa0b1285, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,10)
+       CALL histwrite_phy(o_lwtoa0b1435, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,11)
+       CALL histwrite_phy(o_lwtoa0b1640, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,12)
+       CALL histwrite_phy(o_lwtoa0b1940, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,13)
+       CALL histwrite_phy(o_lwtoa0b2165, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,14)
+       CALL histwrite_phy(o_lwtoa0b2315, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,15)
+       CALL histwrite_phy(o_lwtoa0b2490, zx_tmp_fi2d)
+
+        zx_tmp_fi2d(:)=0.
+       IF (vars_defined) zx_tmp_fi2d=lwtoa0b(:,16)
+       CALL histwrite_phy(o_lwtoa0b2800, zx_tmp_fi2d)
+
+
+!       IF (vars_defined) zx_tmp_fi3dsplw=lwtoa0b
+!       CALL histwrite_phy(o_lwtoa0b, zx_tmp_fi3dsplw)