Changeset 15561 for NEMO/branches

Timestamp:

2021-11-30T14:25:02+01:00 (3 years ago)

Author:

techene

Message:

#2605 #2715 debug on going trcrad modif to be revised

Location:

NEMO/branches/2021/dev_r14318_RK3_stage1/src

Files:

• OCE/stprk3_stg.F90 (modified) (2 diffs)
• TOP/TRP/trcrad.F90 (modified) (1 diff)

NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/stprk3_stg.F90

@@ -245 +245 @@
             END DO
             !
+!!st need a lnc lkn at stage 1 & 2 otherwise tr@Kmm will not be usable in trc_adv 
+            CALL lbc_lnk( 'stprk3_stg', tr(:,:,:,:,Kaa), 'T', 1._wp )   
+
          ENDIF
          !                 !---------------!
@@ -275 +278 @@
 !===>>>>>> Modify tra_adv_...  routines so that Krhs to zero useless
       DO jn = 1, jpts
-         DO_3D( 0, 0, 0, 0, 1, jpkm1 )
-            ts(ji,jj,jk,jn,Krhs) = 0._wp                                    ! set tracer trends to zero
-         END_3D
+!!st does not work due to lbc_lnk north fold : need to be merged with the trunk LBC pb changes for the namelist 
+!!         DO_3D( 0, 0, 0, 0, 1, jpkm1 )
+!!            ts(ji,jj,jk,jn,Krhs) = 0._wp                                    ! set tracer trends to zero
+!!         END_3D
+         ts(:,:,:,jn,Krhs) = 0._wp
       END DO
 

NEMO/branches/2021/dev_r14318_RK3_stage1/src/TOP/TRP/trcrad.F90

@@ -149 +149 @@
      zs2rdt = 1. / ( 2. * rn_Dt )
      !
+#if ! defined key_RK3
      DO jt = 1,2  ! Loop over time indices since exactly the same fix is applied to "now" and "after" fields
         IF( jt == 1 ) itime = Kbb
         IF( jt == 2 ) itime = Kmm
+#else
+     DO jt = 1,1  ! Loop over time indices since exactly the same fix is applied to "now" and "after" fields
+        IF( jt == 1 ) itime = Kmm
+#endif
 
         IF( PRESENT( cpreserv )  ) THEN     !==  total tracer concentration is preserved  ==!