Changeset 6259 for branches/NERC

Timestamp:

2016-01-15T16:47:16+01:00 (9 years ago)

Author:

jpalmier

Message:

JPALM --15-01-2016 -- bugfix lbc_lnk before trcadv is called

Location:

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC

Files:

• TRP/trcadv.F90 (modified) (3 diffs)
• TRP/trcnxt.F90 (modified) (1 diff)
• trcrst.F90 (modified) (2 diffs)

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/TRP/trcadv.F90

@@ -27 +27 @@
    USE ldftra_oce      ! lateral diffusion coefficient on tracers
    USE prtctl_trc      ! Print control
+   USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
 
    IMPLICIT NONE
@@ -71 +72 @@
       INTEGER, INTENT(in) ::   kt   ! ocean time-step index
       !
-      INTEGER ::   jk 
+      INTEGER ::   jk, jn 
       CHARACTER (len=22) ::   charout
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zun, zvn, zwn  ! effective velocity
@@ -105 +106 @@
       zvn(:,:,jpk) = 0._wp                                                     ! no transport trough the bottom
       zwn(:,:,jpk) = 0._wp                                                     ! no transport trough the bottom
+      ! 
+      !! Jpalm -- 14-01-2016 -- restart and proc pb - try this... 
+      CALL lbc_lnk( zun(:,:,:), 'U', -1. )
+      CALL lbc_lnk( zvn(:,:,:), 'V', -1. )
+      CALL lbc_lnk( zwn(:,:,:), 'W', -1. )
+      DO jn = 1, jptra
+         CALL lbc_lnk( trn(:,:,:,jn), 'T', 1. )
+         CALL lbc_lnk( trb(:,:,:,jn), 'T', 1. )
+      END DO
+      !
 
       IF( lk_traldf_eiv .AND. .NOT. ln_traldf_grif )   &  ! add the eiv transport (if necessary)

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/TRP/trcnxt.F90

@@ -105 +105 @@
       DO jn = 1, jptra
          CALL lbc_lnk( tra(:,:,:,jn), 'T', 1. )   
+         CALL lbc_lnk( trn(:,:,:,jn), 'T', 1. )   
+         CALL lbc_lnk( trb(:,:,:,jn), 'T', 1. )   
       END DO
 

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

@@ -40 +40 @@
    USE trcsms_cfc
 #endif
+   USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
+
    IMPLICIT NONE
    PRIVATE
@@ -155 +157 @@
          trb(:,:,:,jn) = trb(:,:,:,jn) * tmask(:,:,:)
       END DO
-
+      !
+      !! JPALM -- 14-01-2016 -- problem in TOP after a restart, proc edge
+      !!                        appear with abherent values.
+      !!                     -- add lbc_lnk(trn and trb) may solve this
+      ! Update after tracer on domain lateral boundaries
+      DO jn = 1, jptra
+         CALL lbc_lnk( trn(:,:,:,jn), 'T', 1. )
+         CALL lbc_lnk( trb(:,:,:,jn), 'T', 1. )
+      END DO
+      !
       !! AXY (09/06/14): the ARCHER version of MEDUSA-2 does not include an equivalent
       !!                 call to MEDUSA-2 at this point; this suggests that the FCM