Changeset 15532 for NEMO/branches/2021/dev_r14318_RK3_stage1/src/OFF/dtadyn.F90

Timestamp:

2021-11-24T12:47:32+01:00 (3 years ago)

Author:

techene

Message:

#2605 #2715 : version in dev (still buggy)

File:

• NEMO/branches/2021/dev_r14318_RK3_stage1/src/OFF/dtadyn.F90 (modified) (4 diffs)

NEMO/branches/2021/dev_r14318_RK3_stage1/src/OFF/dtadyn.F90

@@ -176 +176 @@
       ENDIF
       !
-      CALL eos    ( ts(:,:,:,:,Kmm), rhd, rhop, gdept_0(:,:,:) ) ! In any case, we need rhop
+      CALL eos    ( ts(:,:,:,:,Kmm), rhd, gdept_0(:,:,:) ) ! In any case, we need rhd
       CALL eos_rab( ts(:,:,:,:,Kmm), rab_n, Kmm )       ! now    local thermal/haline expension ratio at T-points
       CALL bn2    ( ts(:,:,:,:,Kmm), rab_n, rn2, Kmm )  ! before Brunt-Vaisala frequency need for zdfmxl
@@ -193 +193 @@
       !
       !
-      CALL eos( ts(:,:,:,:,Kmm), rhd, rhop, gdept_0(:,:,:) ) ! In any case, we need rhop
+      CALL eos( ts(:,:,:,:,Kmm), rhd, gdept_0(:,:,:) ) ! In any case, we need rhd
       !
       IF(sn_cfctl%l_prtctl) THEN                 ! print control
@@ -666 +666 @@
       !
       IF( l_ldfslp .AND. .NOT.lk_c1d ) THEN    ! Computes slopes (here avt is used as workspace)
-         CALL eos    ( pts, rhd, rhop, gdept_0(:,:,:) )
+         CALL eos    ( pts, rhd, gdept_0(:,:,:) )
          CALL eos_rab( pts, rab_n, Kmm )       ! now local thermal/haline expension ratio at T-points
          CALL bn2    ( pts, rab_n, rn2, Kmm  ) ! now    Brunt-Vaisala
@@ -724 +724 @@
       ts(:,:,:,jp_sal,Kmm) = sf_dyn(jf_sal)%fnow(:,:,:)  * tmask(:,:,:)    ! salinity
       !
-      CALL eos    ( ts(:,:,:,:,Kmm), rhd, rhop, gdept_0(:,:,:) ) ! In any case, we need rhop
+      CALL eos    ( ts(:,:,:,:,Kmm), rhd, gdept_0(:,:,:) ) ! In any case, we need rhd
 
       IF(sn_cfctl%l_prtctl) THEN                     ! print control