From 43ee17f59f2704adcbaa03bacbe7b1067c649cc5 Mon Sep 17 00:00:00 2001 From: Aleksey Koldunov Date: Fri, 22 Sep 2023 10:48:58 +0200 Subject: [PATCH] Update oce_ale_pressure_bv.F90 Removing the extra line that made the surface density incorrect --- src/oce_ale_pressure_bv.F90 | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/src/oce_ale_pressure_bv.F90 b/src/oce_ale_pressure_bv.F90 index 4a2a19082..c6689a2eb 100644 --- a/src/oce_ale_pressure_bv.F90 +++ b/src/oce_ale_pressure_bv.F90 @@ -3176,7 +3176,6 @@ subroutine init_ref_density_advanced(tracers, partit, mesh) do nz=nzmin,nzmax call densityJM_components(ref_temp1D(nz), ref_salt1D(nz), bulk_0, bulk_pz, bulk_pz2, rhopot) rho = bulk_0 + auxz*bulk_pz + auxz*bulk_pz2 - density_ref(nzmin, node) = rho*rhopot/(rho+0.1_WP*auxz) auxz=Z_3d_n(nz,node) rho = bulk_0 + auxz*bulk_pz + auxz*bulk_pz2 density_ref(nz,node) = rho*rhopot/(rho+0.1_WP*auxz) @@ -3184,4 +3183,4 @@ subroutine init_ref_density_advanced(tracers, partit, mesh) end do !$OMP END PARALLEL DO if(mype==0) write(*,*) ' --> compute reference density' -end subroutine init_ref_density_advanced \ No newline at end of file +end subroutine init_ref_density_advanced