Changeset 3626 for LMDZ6/branches/Ocean_skin/libf/phylmd/pbl_surface_mod.F90

Timestamp:

Jan 30, 2020, 2:38:10 PM (4 years ago)

Author:

lguez

Message:

Bug fix: use missing_val from XIOS for output with XIOS, else the
operations are not done propoerly in XIOS.

File:

• LMDZ6/branches/Ocean_skin/libf/phylmd/pbl_surface_mod.F90 (modified) (3 diffs)

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

@@ -294 +294 @@
     use phys_state_var_mod, only: ds_ns, dt_ns
     use phys_output_var_mod, only: t_int, s_int, dter, dser, tkt, tks, rf, taur
-    use netcdf, only: nf90_fill_real
+#ifdef CPP_XIOS
+    USE wxios, ONLY: missing_val
+#else
+    use netcdf, only: missing_val => nf90_fill_real
+#endif
 
     IMPLICIT NONE
@@ -966 +970 @@
 
     if (activate_ocean_skin >= 1) then
-       t_int = nf90_fill_real
-       s_int = nf90_fill_real
-       dter = nf90_fill_real
-       dser = nf90_fill_real
-       tkt = nf90_fill_real
-       tks = nf90_fill_real
-       rf = nf90_fill_real
-       taur = nf90_fill_real
+       t_int = missing_val
+       s_int = missing_val
+       dter = missing_val
+       dser = missing_val
+       tkt = missing_val
+       tks = missing_val
+       rf = missing_val
+       taur = missing_val
     end if
 
@@ -2566 +2570 @@
 
        if (nsrf == is_oce .and. activate_ocean_skin >= 1) then
-          dt_ns = nf90_fill_real
-          ds_ns = nf90_fill_real
+          dt_ns = missing_val
+          ds_ns = missing_val
           ds_ns(ni(:knon)) = yds_ns(:knon)
           dt_ns(ni(:knon)) = ydt_ns(:knon)