Changeset 8643 for branches/NERC/dev_r5518_GO6_Carb_Debug/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

Timestamp:

2017-10-19T18:41:56+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM -- 19-10-17 -- Add ctl_warn for master node output

File:

• branches/NERC/dev_r5518_GO6_Carb_Debug/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90 (modified) (6 diffs)

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

@@ -81 +81 @@
                                             gdepw_0, gdepw_n,               &
                                             nday_year, nsec_day, nyear,     &
-                                            rdt, tmask
+                                            rdt, tmask, mig, mjg
       USE in_out_manager,             ONLY: lwp, numout, nn_date0
 # if defined key_iomput
@@ -87 +87 @@
 # endif
       USE lbclnk,                     ONLY: lbc_lnk
-      USE lib_mpp,                    ONLY: ctl_stop
+      USE lib_mpp,                    ONLY: ctl_stop, ctl_warn
       USE oce,                        ONLY: tsb, tsn
       USE par_kind,                   ONLY: wp
@@ -477 +477 @@
              !! AXY (28/02/14): check input fields
                   if (ztmp(ji,jj) .lt. -3.0 .or. ztmp(ji,jj) .gt. 40.0 ) then
-                     IF(lwp) WRITE(numout,*)                                 &
-                        ' trc_bio_medusa: T WARNING 3D, ',                   &
-                        tsb(ji,jj,jk,jp_tem), tsn(ji,jj,jk,jp_tem),          &
-                        ' at (', ji, ',', jj, ',', jk, ') at time', kt
+                     Call ctl_warn(' trc_bio_medusa: T WARNING 3D : ',       &
+                                   tsn(ji,jj,jk,jp_tem), 'at I, J, K, kt :', & 
+                                   mig(ji), mjg(jj), jk, kt )
+# if defined key_debug_medusa
                      !! temperature
                      ztmp(ji,jj) = tsb(ji,jj,jk,jp_tem)
@@ -491 +491 @@
                       IF(lwp) WRITE(numout,*)                                 &
                      tsn(ji-1,jj-1,jk,jp_tem), tsn(ji,jj-1,jk,jp_tem), tsn(ji+1,jj-1,jk,jp_tem)
-                     sumtsn = ( tmask(ji-1,jj+1,jk) * tsn(ji-1,jj+1,jk,jp_tem) ) + & 
-                              ( tmask(ji  ,jj+1,jk) * tsn(ji  ,jj+1,jk,jp_tem) ) + & 
-                              ( tmask(ji+1,jj+1,jk) * tsn(ji+1,jj+1,jk,jp_tem) ) + & 
-                              ( tmask(ji-1,jj  ,jk) * tsn(ji-1,jj  ,jk,jp_tem) ) + & 
-                              ( tmask(ji+1,jj  ,jk) * tsn(ji+1,jj  ,jk,jp_tem) ) + & 
-                              ( tmask(ji-1,jj-1,jk) * tsn(ji-1,jj-1,jk,jp_tem) ) + & 
-                              ( tmask(ji  ,jj-1,jk) * tsn(ji  ,jj-1,jk,jp_tem) ) + & 
-                              ( tmask(ji+1,jj-1,jk) * tsn(ji+1,jj-1,jk,jp_tem) )
                       IF(lwp) WRITE(numout,*)                                 &
                      ' trc_bio_medusa: Abnormal T suroundings Sal'
@@ -532 +524 @@
                       IF(lwp) WRITE(numout,*)                                 &
                      tmask(ji-1,jj-1,jk), tmask(ji,jj-1,jk), tmask(ji+1,jj-1,jk)
+# endif
+                     !! Correct out of range values
+                     sumtsn = ( tmask(ji-1,jj+1,jk) * tsn(ji-1,jj+1,jk,jp_tem) ) + &
+                              ( tmask(ji  ,jj+1,jk) * tsn(ji  ,jj+1,jk,jp_tem) ) + &
+                              ( tmask(ji+1,jj+1,jk) * tsn(ji+1,jj+1,jk,jp_tem) ) + &
+                              ( tmask(ji-1,jj  ,jk) * tsn(ji-1,jj  ,jk,jp_tem) ) + &
+                              ( tmask(ji+1,jj  ,jk) * tsn(ji+1,jj  ,jk,jp_tem) ) + &
+                              ( tmask(ji-1,jj-1,jk) * tsn(ji-1,jj-1,jk,jp_tem) ) + &
+                              ( tmask(ji  ,jj-1,jk) * tsn(ji  ,jj-1,jk,jp_tem) ) + &
+                              ( tmask(ji+1,jj-1,jk) * tsn(ji+1,jj-1,jk,jp_tem) )
                      summask = tmask(ji-1,jj+1,jk) + tmask(ji  ,jj+1,jk)   +  &
                                tmask(ji+1,jj+1,jk) + tmask(ji-1,jj  ,jk)   +  &
-                               tmask(ji+1,jj  ,jk) + tmask(ji-1,jj-1,jk)   +  & 
-                               tmask(ji  ,jj-1,jk) + tmask(ji+1,jj-1,jk) 
+                               tmask(ji+1,jj  ,jk) + tmask(ji-1,jj-1,jk)   +  &
+                               tmask(ji  ,jj-1,jk) + tmask(ji+1,jj-1,jk)
                      tsnavg = ( sumtsn / summask )
-                     !! Correct out of range values
+                     !!
                      IF ( ( summask .EQ. 0.0 ) .OR. (tsnavg .LT. -3.0 ) .OR.   &
                           ( tsnavg .GT. 40.0 ) ) THEN    
                         IF (ztmp(ji,jj) .LT. -3.0 ) THEN
-                           IF(lwp) WRITE(numout,*)                           &
+                           CALL ctl_warn(                                    &
                            ' trc_bio_medusa: T SWITCHING 3D, ',              &
-                           tsn(ji,jj,jk,jp_tem), ' -> -3.0 '
+                           tsn(ji,jj,jk,jp_tem), ' -> -3.0 ' )
                            ztmp(ji,jj) = -3.0
                         ENDIF
                         IF (ztmp(ji,jj) .GT. 40.0 ) THEN
-                           IF(lwp) WRITE(numout,*)                           &
+                           CALL ctl_warn(                                    &
                            ' trc_bio_medusa: T SWITCHING 3D, ',              &
-                           tsn(ji,jj,jk,jp_tem), ' -> 40.0 '
+                           tsn(ji,jj,jk,jp_tem), ' -> 40.0 ' )
                            ztmp(ji,jj) = 40.0
                         ENDIF
                      ELSE 
-                        IF(lwp) WRITE(numout,*)                              &
+                        CALL ctl_warn(                                       &
                         ' trc_bio_medusa: T SWITCHING 3D, ',                 &
-                        tsn(ji,jj,jk,jp_tem), ' -> surounding avg : ', tsnavg
+                        tsn(ji,jj,jk,jp_tem), ' -> surounding avg : ', tsnavg )
                         ztmp(ji,jj) = tsnavg
                      ENDIF
@@ -561 +563 @@
                   !! end T chack
                   if (zsal(ji,jj) .lt. 1.0 .or. zsal(ji,jj) .gt. 47.0 ) then
-                     IF(lwp) WRITE(numout,*)                                 &
+                        CALL ctl_warn(                                       &
                         ' trc_bio_medusa: S WARNING 2D, ',                   &
-                        tsb(ji,jj,jk,jp_sal), tsn(ji,jj,jk,jp_sal),          &
-                        ' at (', ji, ',', jj, ',', jk, ') at time', kt
+                        tsn(ji,jj,jk,jp_sal),                                &
+                        ' at I, J, K, kt : ', mig(ji), mjg(jj), jk, kt)
                      !! Correct out of range values
                      IF (zsal(ji,jj) .LT. 1.0 ) THEN
-                        IF(lwp) WRITE(numout,*)                              &
+                        CALL ctl_warn(                                       &
                         ' trc_bio_medusa: S SWITCHING 3D, ',                 &
-                        tsn(ji,jj,jk,jp_sal), ' -> 1.0 '
+                        tsn(ji,jj,jk,jp_sal), ' -> 1.0 ')
                         zsal(ji,jj) = 1.0
                      ENDIF
                      IF (zsal(ji,jj) .GT. 47.0 ) THEN
-                        IF(lwp) WRITE(numout,*)                              &
+                        CALL ctl_warn(                                       &
                         ' trc_bio_medusa: T SWITCHING 3D, ',                 &
-                        tsn(ji,jj,jk,jp_sal), ' -> 47.0 '
+                        tsn(ji,jj,jk,jp_sal), ' -> 47.0 ')
                         zsal(ji,jj) = 47.0
                      ENDIF