Changeset 12761 for NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/stpMLF.F90

Timestamp:

2020-04-17T10:28:12+02:00 (4 years ago)

Author:

techene

Message:

add the computation of r3._f before trc_src in stepMLF.F90 for that there is to change ssh_atf to be able to compute a filtrered ssh without changing the ssh field

File:

• NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/stpMLF.F90 (modified) (3 diffs)

NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/stpMLF.F90

@@ -91 +91 @@
       INTEGER ::   kcall        ! optional integer argument (dom_vvl_sf_nxt)
 !!st patch
-      REAL(wp), DIMENSION(jpi,jpj,jpk) :: zgdept 
+      REAL(wp),              DIMENSION(jpi,jpj,jpk) ::   zgdept
+      REAL(wp), ALLOCATABLE, DIMENSION(:,:)         ::   zssh_f
       !! ---------------------------------------------------------------------
 #if defined key_agrif
@@ -248 +249 @@
       ! Passive Tracer Model
       !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
-                         CALL trc_stp       ( kstp, Nbb, Nnn, Nrhs, Naa )  ! time-stepping
+                         ALLOCATE( zssh_f(jpi,jpj) )
+                         CALL ssh_atf    ( kstp, Nbb, Nnn, Naa , ssh, zssh_f )   ! time filtering of "now" sea surface height
+                         CALL dom_qco_r3c( zssh_f, r3t_f, r3u_f, r3v_f )         ! "now" ssh/h_0 ratio from filtrered ssh
+                         !
+                         CALL trc_stp    ( kstp, Nbb, Nnn, Nrhs, Naa )           ! time-stepping
+                         DEALLOCATE( zssh_f )
 #endif
 
@@ -297 +303 @@
                          CALL finalize_sbc  ( kstp, Nbb, Naa, uu, vv, ts )          ! boundary condifions
                          CALL ssh_atf       ( kstp, Nbb, Nnn, Naa, ssh )            ! time filtering of "now" sea surface height
-                         CALL dom_qco_r3c    ( ssh(:,:,Nnn), r3t_f, r3u_f, r3v_f )   ! "now" ssh/h_0 ratio from filtrered ssh
+                         CALL dom_qco_r3c   ( ssh(:,:,Nnn), r3t_f, r3u_f, r3v_f )   ! "now" ssh/h_0 ratio from filtrered ssh
                          CALL tra_atf_qco   ( kstp, Nbb, Nnn, Naa, ts )             ! time filtering of "now" tracer arrays
                          CALL dyn_atf_qco   ( kstp, Nbb, Nnn, Naa, uu, vv  )        ! time filtering of "now" velocities and scale factors