Changeset 3970 for branches/2013/dev_r3867_MERCATOR1_DYN/NEMOGCM/NEMO/OPA_SRC/step.F90

Timestamp:

2013-07-11T15:59:14+02:00 (11 years ago)

Author:

cbricaud

Message:

Time splitting update, see ticket #1079

File:

• branches/2013/dev_r3867_MERCATOR1_DYN/NEMOGCM/NEMO/OPA_SRC/step.F90 (modified) (4 diffs)

branches/2013/dev_r3867_MERCATOR1_DYN/NEMOGCM/NEMO/OPA_SRC/step.F90

@@ -108 +108 @@
       !
       !  VERTICAL PHYSICS
+      ! bg jchanut tschanges
+      ! One need bottom friction parameter in ssh_wzv routine with time splitting.
+      ! The idea could be to move the call below before ssh_wzv. However, "now" scale factors
+      ! at U-V points (which are set thanks to sshu_n, sshv_n) are actually available in sshwzv.
+      ! These are needed for log bottom friction...
+#if ! defined key_dynspg_ts
                          CALL zdf_bfr( kstp )         ! bottom friction
+#endif
+      ! end jchanut tschanges
 
       !                                               ! Vertical eddy viscosity and diffusivity coefficients
@@ -206 +214 @@
             &                                          rhd, gru , grv  )      ! of t, s, rd at the last ocean level
 
-      ELSE                                                  ! centered hpg  (eos then time stepping)
+      ELSE   
+                                               ! centered hpg  (eos then time stepping)
+      ! bg jchanut tschanges
+#if ! defined key_dynspg_ts
+      ! eos already called
                              CALL eos    ( tsn, rhd, rhop )      ! now in situ density for hpg computation
          IF( ln_zps      )   CALL zps_hde( kstp, jpts, tsn, gtsu, gtsv,  &    ! zps: now hor. derivative
             &                                          rhd, gru , grv  )      ! of t, s, rd at the last ocean level
+#endif
+      ! end jchanut tschanges
          IF( ln_zdfnpc   )   CALL tra_npc( kstp )                ! update after fields by non-penetrative convection
                              CALL tra_nxt( kstp )                ! tracer fields at next time step
@@ -217 +231 @@
       ! Dynamics                                    (tsa used as workspace)
       !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
+      ! bg jchanut tschanges
+#if defined key_dynspg_ts      
+! revert to previously computed tendencies:
+! (not using ua, va as temporary arrays during tracers' update could avoid that)
+                               ua(:,:,:) = ua_bak(:,:,:)            
+                               va(:,:,:) = va_bak(:,:,:)
+                               CALL dyn_bfr( kstp )         ! bottom friction
+                               CALL dyn_zdf( kstp )         ! vertical diffusion
+#else
+      ! end jchanut tschanges
                                ua(:,:,:) = 0.e0             ! set dynamics trends to zero
                                va(:,:,:) = 0.e0
@@ -236 +260 @@
                                CALL dyn_zdf( kstp )         ! vertical diffusion
                                CALL dyn_spg( kstp, indic )  ! surface pressure gradient
+      ! bg jchanut tschanges
+#endif
+      ! end jchanut tschanges
                                CALL dyn_nxt( kstp )         ! lateral velocity at next time step