Changeset 5904

Timestamp:

Dec 3, 2025, 6:06:16 PM (88 minutes ago)

Author:

fcheruy

Message:

ajout sortie tslsi

Location:

LMDZ6/trunk

Files:

• DefLists/field_def_lmdz.xml (modified) (1 diff)
• libf/phylmd/phys_output_ctrlout_mod.F90 (modified) (1 diff)
• libf/phylmd/phys_output_write_mod.F90 (modified) (2 diffs)

LMDZ6/trunk/DefLists/field_def_lmdz.xml

@@ -45 +45 @@
         <field id="pt0"     long_name="Standard Air Surface Temperature"      unit="K" />
         <field id="tsol"    long_name="Surface Temperature"     unit="K" />
+        <field id="tslsi"    long_name="Surface Temperature land and iced surf"     unit="K" />
         <field id="t2m"        long_name="Temperature 2m" unit="K" />
         <field id="tasmax"  field_ref="t2m"  operation="maximum" />

LMDZ6/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

@@ -486 +486 @@
   !FC
 
+  TYPE(ctrl_out), SAVE :: o_tslsi = ctrl_out((/ 1, 10, 10, 10, 10, 10, 11, 11, 11, 11/), &
+    'tslsi', 'surface temperature except free ocean ', 'K', (/ ('', i=1, 10) /))
   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) /))

LMDZ6/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -259 +259 @@
          o_lwtoabtot, o_lwtoa0btot, &  !Somme de tous les canaux
          o_lwtoab4co2,o_lwtoacsb4co2,                            &   ! les 16 canaux FC
+         o_tslsi ,                                               &   ! pour forum FC
+!FC
 !GG         o_zxfluxt,o_zxfluxq
          o_fcds, o_fcdi, o_dh_basal_growth, o_dh_basal_melt, &
@@ -1862 +1864 @@
           ENDIF
                CALL histwrite_phy(o_delta_tsurf,zx_tmp_fi2d) 
+          !FC
+                  IF (vars_defined) THEN
+          DO i=1, klon
+             IF (pctsrf(i,is_ter).gt.epsfra .or. pctsrf(i,is_sic) .gt. epsfra .or. pctsrf(i,is_lic).gt.epsfra ) THEN
+                zx_tmp_fi2d(i) = (ftsol(i, is_ter) * pctsrf(i,is_ter)+ &
+                     ftsol(i, is_sic) * pctsrf(i,is_sic)+  &
+                     ftsol(i, is_lic) * pctsrf(i,is_lic) )/ &  
+                     (pctsrf(i,is_ter)+pctsrf(i,is_sic)+ pctsrf(i,is_lic))
+             ELSE
+                zx_tmp_fi2d(i) = 0.
+             ENDIF
+          ENDDO
+       ENDIF
+
+          !FC
+               CALL histwrite_phy(o_tslsi,zx_tmp_fi2d)