Changeset 3458 for LMDZ6/branches/Ocean_skin/libf/phylmd/phyetat0.F90

Timestamp:

Jan 30, 2019, 5:55:48 PM (6 years ago)

Author:

lguez

Message:

Introduce variable activate_ocean_skin in module config_ocean_skin_m.

Bug fix in phys_state_var_end: we need to deallocate variables for
lmdz1d (although it is useless for a 3D run).

File:

• LMDZ6/branches/Ocean_skin/libf/phylmd/phyetat0.F90 (modified) (2 diffs)

LMDZ6/branches/Ocean_skin/libf/phylmd/phyetat0.F90

@@ -30 +30 @@
   USE time_phylmdz_mod, ONLY: init_iteration, pdtphys, itau_phy
   use netcdf, only: nf90_fill_real
+  use config_ocean_skin_m, only: activate_ocean_skin
 
   IMPLICIT none
@@ -529 +530 @@
   CALL fonte_neige_init(run_off_lic_0)
 
-  found = phyetat0_get(1, ds_ns, "ds_ns", "delta salinity near surface", 0.)
-  found = phyetat0_get(1, dt_ns, "dT_ns", "delta temperature near surface", 0.)
-
-  where (pctsrf(:, is_oce) == 0.)
-     ds_ns = nf90_fill_real
-     dt_ns = nf90_fill_real
-  end where
+  if (activate_ocean_skin >= 1) then
+     found = phyetat0_get(1, ds_ns, "ds_ns", "delta salinity near surface", 0.)
+     found = phyetat0_get(1, dt_ns, "dT_ns", "delta temperature near surface", &
+          0.)
+
+     where (pctsrf(:, is_oce) == 0.)
+        ds_ns = nf90_fill_real
+        dt_ns = nf90_fill_real
+     end where
+  end if
 
 END SUBROUTINE phyetat0