Changeset 5791 for LMDZ6/branches/contrails/libf/phylmdiso/phyredem.F90

Timestamp:

Jul 28, 2025, 7:23:15 PM (7 days ago)

Author:

aborella

Message:

Merge with trunk r5789

Location:

LMDZ6/branches/contrails

Files:

• . (modified) (1 prop)
• libf/phylmdiso/phyredem.F90 (modified) (5 diffs)

LMDZ6/branches/contrails/libf/phylmdiso/phyredem.F90

@@ -31 +31 @@
                                 du_gwd_rando, du_gwd_front, u10m, v10m, &
                                 treedrg, solswfdiff, delta_sal, ds_ns, dt_ns, &
-                                delta_sst, ratqs_inter_, dter, dser, dt_ds
+!GG                                delta_sst, ratqs_inter_, dter, dser, dt_ds
+                                delta_sst, ratqs_inter_, dter, dser, dt_ds, &
+                                hice, tice, bilg_cumul
+!GG
 #ifdef ISO
   USE phys_state_var_mod, ONLY: xtsol, fxtevap,xtrain_fall, xtsnow_fall,     &
@@ -51 +54 @@
   USE carbon_cycle_mod, ONLY : carbon_cycle_cpl, co2_send, carbon_cycle_rad, RCO2_glo
   USE indice_sol_mod, ONLY: nbsrf, is_oce, is_sic, is_ter, is_lic, epsfra
-  USE surface_data, ONLY: type_ocean, version_ocean
-  USE ocean_slab_mod, ONLY : nslay, tslab, seaice, tice, fsic
+!GG  USE surface_data, ONLY: type_ocean, version_ocean
+  USE surface_data, ONLY: type_ocean, version_ocean, iflag_seaice, iflag_seaice_alb, &
+                          iflag_leads
+!GG
+  USE ocean_slab_mod, ONLY : nslay, tslab, seaice, tice_slab, fsic
   USE time_phylmdz_mod, ONLY: annee_ref, day_end, itau_phy, pdtphys
   use config_ocean_skin_m, only: activate_ocean_skin
@@ -136 +142 @@
   IF (carbon_cycle_rad) tab_cntrl(17) = RCO2_glo
   !PRINT*, "PC : phyredem RCO2_glo =",RCO2_glo
+  !GG
+  tab_cntrl(18 ) = iflag_seaice
+  tab_cntrl(19 ) = iflag_seaice_alb
+  tab_cntrl(20 ) = iflag_leads
+  !GG
 
   DO pass=1,2   ! pass=1 netcdf definition ; pass=2 netcdf write
@@ -214 +225 @@
 
     CALL put_field_srf1(pass,"SNOW", "Neige", snow(:,:))
+
+    !GG
+    CALL put_field(pass,"hice", "Ice thickness", hice)
+    CALL put_field(pass,"tice", "Sea Ice temperature", tice)
+    CALL put_field(pass,"bilg_cumul", "Flux conductivite + transmit sea-ice", bilg_cumul)
+    !GG
 
     CALL put_field(pass,"RADS", "Rayonnement net a la surface", radsol)
@@ -388 +405 @@
         IF (version_ocean == 'sicINT') THEN
             CALL put_field(pass,"seaice", "Slab seaice (kg/m2)", seaice)
-            CALL put_field(pass,"slab_tice", "Slab sea ice temperature", tice)
+            CALL put_field(pass,"slab_tice", "Slab sea ice temperature", tice_slab)
         END IF
     END IF