Changeset 3744 for LMDZ6/branches/Ocean_skin/libf/phylmd/cpl_mod.F90

Timestamp:

Jul 1, 2020, 6:57:48 PM (4 years ago)

Author:

lguez

Message:

Store delta_sst instead of sst_nff

Store as a state variable the difference between ocean-air interface
temperature and bulk SST instead of sst_nff, which can be either
interface temperature or bulk SST. This is clearer. Also, it is
analoguous to what we will do with salinity.

So replace the two dummy arguments tsurf_in and sst_nff of
procedure cpl_send_ocean_fields by a single dummy argument
delta_sst. Replace dummy argument sst_nff of procedures
ocean_cpl_noice and surf_ocean by dummy argument
delta_sst. Replace variable sst_nff of module phys_state_var_mod
by variable delta_sst. Rename local variable ysst_nff of procedure
pbl_surface to ydelta_sst. Set variable delta_sst of module
phys_state_var_mod to 0 for an appearing ocean fraction and a
missing startup field. Replace variable o_sst_nff of module
phys_output_ctrlout_mod by variable o_delta_sst.

Rename variables cpl_delta_temp and cpl_delta_temp_2D of module
cpl_mod to cpl_delta_sst and cpl_delta_sst_2D, clearer. Rename
variable ids_delta_temp of module oasis to ids_delta_sst. Change
infosend(ids_delta_temp)%name to "CODELSST".

File:

• LMDZ6/branches/Ocean_skin/libf/phylmd/cpl_mod.F90 (modified) (11 diffs)

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

@@ -49 +49 @@
   !$OMP THREADPRIVATE(cpl_snow,cpl_evap,cpl_tsol)
 
-  REAL, ALLOCATABLE, SAVE:: cpl_delta_temp(:), cpl_s_int(:)
-  !$OMP THREADPRIVATE(cpl_delta_temp, cpl_s_int)
+  REAL, ALLOCATABLE, SAVE:: cpl_delta_sst(:), cpl_s_int(:)
+  !$OMP THREADPRIVATE(cpl_delta_sst, cpl_s_int)
   
   REAL, ALLOCATABLE, DIMENSION(:,:), SAVE   :: cpl_fder, cpl_albe, cpl_taux, cpl_tauy
@@ -94 +94 @@
   !$OMP THREADPRIVATE(cpl_snow2D, cpl_evap2D, cpl_tsol2D)
 
-  REAL, ALLOCATABLE, SAVE:: cpl_delta_temp_2D(:,:), cpl_s_int_2D(:,:)
-  !$OMP THREADPRIVATE(cpl_delta_temp_2D, cpl_s_int_2D)
+  REAL, ALLOCATABLE, SAVE:: cpl_delta_sst_2D(:,:), cpl_s_int_2D(:,:)
+  !$OMP THREADPRIVATE(cpl_delta_sst_2D, cpl_s_int_2D)
 
   REAL, ALLOCATABLE, DIMENSION(:,:,:), SAVE :: cpl_fder2D, cpl_albe2D
@@ -237 +237 @@
     
        if (activate_ocean_skin == 2) then
-          ALLOCATE(cpl_delta_temp(klon), cpl_s_int(klon), stat = error)
+          ALLOCATE(cpl_delta_sst(klon), cpl_s_int(klon), stat = error)
           sum_error = sum_error + error
        end if
@@ -645 +645 @@
        swdown, lwdown, fluxlat, fluxsens, &
        precip_rain, precip_snow, evap, tsurf, fder, albsol, taux, tauy, windsp,&
-       sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol, tsurf_in, &
-       sst_nff, s_int)
+       sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol, delta_sst, &
+       s_int)
 
     ! This subroutine cumulates some fields for each time-step during
@@ -672 +672 @@
     REAL, INTENT(IN):: sens_prec_liq(:), sens_prec_sol(:) ! (knon)
     REAL, DIMENSION(klon), INTENT(IN)       :: lat_prec_liq, lat_prec_sol
-    real, intent(in):: tsurf_in(:) ! (klon)
-    
-    REAL, intent(in):: sst_nff(:) ! (knon)
-    ! SST not used to compute surface fluxes, in K. If
-    ! activate_ocean_skin == 0 then it is not defined; if
-    ! activate_ocean_skin == 1 then it is the ocean-air interface
-    ! temperature; if activate_ocean_skin == 2 then it is the bulk SST.
+    
+    REAL, intent(in):: delta_sst(:) ! (knon)
+    ! Ocean-air interface temperature minus bulk SST, in
+    ! K. Defined only if activate_ocean_skin >= 1.
 
     real, intent(in):: s_int(:) ! (knon) ocean-air interface salinity, in ppt
@@ -718 +715 @@
 
        if (activate_ocean_skin == 2) then
-          cpl_delta_temp(:knon) = 0.
+          cpl_delta_sst(:knon) = 0.
           cpl_s_int = 0.
        end if
@@ -764 +761 @@
 
        if (activate_ocean_skin == 2) then
-          cpl_delta_temp(ig) = cpl_delta_temp(ig) &
-               + (tsurf_in(ig) - sst_nff(ig)) / REAL(nexca)
+          cpl_delta_sst(ig) = cpl_delta_sst(ig) + delta_sst(ig) / REAL(nexca)
           cpl_s_int(ig) = cpl_s_int(ig) + s_int(ig) / REAL(nexca)
        end if
@@ -815 +811 @@
 
           if (activate_ocean_skin == 2) then
-             ALLOCATE(cpl_delta_temp_2D(nbp_lon, jj_nb), &
+             ALLOCATE(cpl_delta_sst_2D(nbp_lon, jj_nb), &
                   cpl_s_int_2D(nbp_lon, jj_nb), stat = error)
              sum_error = sum_error + error
@@ -875 +871 @@
             CALL gath2cpl(cpl_atm_co2(:,cpl_index), cpl_atm_co22D(:,:), knon, knindex)
        if (activate_ocean_skin == 2) then
-          CALL gath2cpl(cpl_delta_temp, cpl_delta_temp_2D, knon, knindex)
+          CALL gath2cpl(cpl_delta_sst, cpl_delta_sst_2D, knon, knindex)
           CALL gath2cpl(cpl_s_int, cpl_s_int_2D, knon, knindex)
        end if
@@ -1273 +1269 @@
 
     if (activate_ocean_skin == 2) then
-       tab_flds(:, :, ids_delta_temp) = cpl_delta_temp_2D
+       tab_flds(:, :, ids_delta_sst) = cpl_delta_sst_2D
        tab_flds(:, :, ids_s_int) = cpl_s_int_2D
     end if
@@ -1512 +1508 @@
     ENDIF
 
-    if (activate_ocean_skin == 2) deallocate(cpl_delta_temp_2d, cpl_s_int_2d)
+    if (activate_ocean_skin == 2) deallocate(cpl_delta_sst_2d, cpl_s_int_2d)
 
     IF (sum_error /= 0) THEN