Changeset 7497

Timestamp:

2016-12-14T16:46:28+01:00 (8 years ago)

Author:

jpalmier

Message:

JPALM -- 14-12-2016 -- MEDUSA Diags update ; remove extra lbc-lnk

Location:

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC

Files:

• MEDUSA/trcbio_medusa.F90 (modified) (3 diffs)
• TRP/trcadv.F90 (modified) (2 diffs)

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

@@ -2532 +2532 @@
                !!
                !! AXY (22/07/09): accelerate detrital remineralisation in the bottom box
-               if (jk.eq.jmbathy) .and. jsfd.eq.1) then
+               if ((jk.eq.jmbathy) .and. jsfd.eq.1) then
                   fdd  = 1.0  * zdet
 # if defined key_roam
@@ -2566 +2566 @@
                !!----------------------------------------------------------------------
                !!
-               if (jk.eq.jmbathy) .and. jorgben.eq.1) then
+               if ((jk.eq.jmbathy) .and. jorgben.eq.1) then
                   !! this is the BOTTOM OCEAN BOX -> into the benthic pool!
                   !!
@@ -2881 +2881 @@
                !! AXY (10/07/12): amended to only apply sedimentary flux up to ~500 m down
                !! if (jk.eq.(mbathy(ji,jj)-1).AND.jk.lt.i1100) then
-               if (jk.eq.jmbathy).AND.jk.le.i0500) then
+               if ((jk.eq.jmbathy).AND.jk.le.i0500) then
                   !! Moore et al. (2004) cite a coastal California value of 5 umol/m2/d, but adopt a
                   !! global value of 2 umol/m2/d for all areas < 1100 m; here we use this latter value

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)
+   !! USE lbclnk          ! ocean lateral boundary conditions (or mpp link)
 
    IMPLICIT NONE
@@ -108 +108 @@
       ! 
       !! Jpalm -- 14-01-2016 -- restart and proc pb - try this... 
-      DO jn = 1, jptra
-         CALL lbc_lnk( trn(:,:,:,jn), 'T', 1. )
-         CALL lbc_lnk( trb(:,:,:,jn), 'T', 1. )
-      END DO
+      !! DO jn = 1, jptra
+      !!   CALL lbc_lnk( trn(:,:,:,jn), 'T', 1. )
+      !!   CALL lbc_lnk( trb(:,:,:,jn), 'T', 1. )
+      !! END DO
       !