Changeset 4021 for LMDZ6/branches/Ocean_skin

Timestamp:

Nov 26, 2021, 8:27:27 AM (3 years ago)

Author:

lguez

Message:

Polish

Location:

LMDZ6/branches/Ocean_skin/libf/phylmd

Files:

• phyredem.F90 (modified) (1 diff)
• phys_output_var_mod.F90 (modified) (1 diff)
• surf_ocean_mod.F90 (modified) (1 diff)
• wake.F90 (modified) (1 diff)

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

@@ -162 +162 @@
     CALL put_field_srf1(pass,"TS","Temperature",ftsol(:,:))
 
-!!    CALL put_field_srf1(pass,"DELTA_TS","w-x surface temperature difference", delta_tsurf(:,:))
     CALL put_field_srf1(pass,"DELTATS","w-x surface temperature difference", delta_tsurf(:,:))
-
-!    CALL put_field_srf1(pass,"BETA_S","Aridity factor", beta_aridity(:,:))
     CALL put_field_srf1(pass,"BETAS","Aridity factor", beta_aridity(:,:))
 !    End surface variables

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

@@ -134 +134 @@
   ! Ocean-atmosphere interface, subskin ocean and near-surface ocean:
   
-  REAL, SAVE, ALLOCATABLE:: tkt(:)
+  REAL, SAVE, ALLOCATABLE:: tkt(:) ! (klon)
   ! épaisseur (m) de la couche de diffusion thermique (microlayer)
   ! cool skin thickness
 
-  REAL, SAVE, ALLOCATABLE:: tks(:)
+  REAL, SAVE, ALLOCATABLE:: tks(:) ! (klon)
   ! épaisseur (m) de la couche de diffusion de masse (microlayer)
   
-  REAL, SAVE, ALLOCATABLE:: taur(:) ! momentum flux due to rain, in Pa
-
-  REAL, SAVE, ALLOCATABLE:: sss(:)
+  REAL, SAVE, ALLOCATABLE:: taur(:) ! (klon) momentum flux due to rain, in Pa
+
+  REAL, SAVE, ALLOCATABLE:: sss(:) ! (klon)
   ! bulk salinity of the surface layer of the ocean, in ppt
   

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

@@ -334 +334 @@
        delta_sst = t_int - tsurf_new(:knon)
        delta_sal = s_int - sss(:knon)
-       if (activate_ocean_skin >= 2) tsurf_new(:knon) = t_int
-       if (activate_ocean_skin == 2 .and. type_ocean == 'couple') &
-            dt_ds = (tks / tkt) * dter
+
+       if (activate_ocean_skin == 2) then
+          tsurf_new(:knon) = t_int
+          if (type_ocean == 'couple') dt_ds = (tks / tkt) * dter
+       end if
     end if
     

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

@@ -2240 +2240 @@
     d_deltaqw, sigmaw, d_sigmaw, alpha)
   ! ------------------------------------------------------
-  ! Dtermination du coefficient alpha tel que les tendances
+  ! D\'etermination du coefficient alpha tel que les tendances
   ! corriges alpha*d_G, pour toutes les grandeurs G, correspondent
   ! a une humidite positive dans la zone (x) et dans la zone (w).