Changeset 472 for branches

Timestamp:

01/31/24 14:09:33 (4 months ago)

Author:

dumas

Message:

coef_OM formula as equation 1 in Jourdain et al 2020

File:

• branches/GRISLIv3/SOURCES/bmelt-beckmann-gcm_mod.f90 (modified) (2 diffs)

branches/GRISLIv3/SOURCES/bmelt-beckmann-gcm_mod.f90

@@ -29 +29 @@
   implicit none
 
-  real*8,parameter :: cpw    = 4182.d0
+  real*8,parameter :: cpw    = 3974.d0
   integer :: nbassins !< number of sectors in the ocean
   integer :: nzoc !< number of vertical levels in the ocean (read in netcdf T/S file)
@@ -342 +342 @@
         
     !coef_OM = K_t * 0.01420418516
-    coef_OM = K_t * row * cpw / (cl * ro)
+    coef_OM = K_t * (row * cpw / (cl * ro))**2
 
     if ( ubound(zoc,dim=1).ne.nzoc) then