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

Timestamp:

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

Author:

jpalmier

Message:

jpalm --19-10-17-- Carb failure due to 1-cell exceptionnal and ephemeral T increase - update T passed to MEDUSA and add kill switch

File:

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

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

@@ -177 +177 @@
       !! temporary variables
       REAL(wp) ::    fq0,fq1,fq2,fq3,fq4
+      !!
+      !! T and S check temporary variable
+      REAL(wp) ::    sumtsn, tsnavg
+      INTEGER  ::    summask
       !!
       !!------------------------------------------------------------------
@@ -474 +478 @@
                   if (ztmp(ji,jj) .lt. -3.0 .or. ztmp(ji,jj) .gt. 40.0 ) then
                      IF(lwp) WRITE(numout,*)                                 &
-                        ' trc_bio_medusa: T WARNING 2D, ',                   &
+                        ' 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
-           IF(lwp) WRITE(numout,*)                                 &
-                        ' trc_bio_medusa: T SWITCHING 2D, ',                 &
-                        tsn(ji,jj,jk,jp_tem), ' -> ', tsb(ji,jj,jk,jp_tem)
-                     !! temperatur
+                     !! temperature
                      ztmp(ji,jj) = tsb(ji,jj,jk,jp_tem)
+                      IF(lwp) WRITE(numout,*)                                 &
+                     ' trc_bio_medusa: Abnormal T suroundings Temp'
+                      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)
+                      IF(lwp) WRITE(numout,*)                                 &
+                     tsn(ji-1,jj,jk,jp_tem), tsn(ji,jj,jk,jp_tem), tsn(ji+1,jj,jk,jp_tem)
+                      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'
+                      IF(lwp) WRITE(numout,*)                                 &
+                     tsn(ji-1,jj+1,jk,jp_sal), tsn(ji,jj+1,jk,jp_sal), tsn(ji+1,jj+1,jk,jp_sal)
+                      IF(lwp) WRITE(numout,*)                                 &
+                     tsn(ji-1,jj,jk,jp_sal), tsn(ji,jj,jk,jp_sal), tsn(ji+1,jj,jk,jp_sal)
+                      IF(lwp) WRITE(numout,*)                                 &
+                     tsn(ji-1,jj-1,jk,jp_sal), tsn(ji,jj-1,jk,jp_sal), tsn(ji+1,jj-1,jk,jp_sal)
+                     !!
+                      IF(lwp) WRITE(numout,*)                                 &
+                     ' trc_bio_medusa: Abnormal T suroundings DIC'
+                      IF(lwp) WRITE(numout,*)                                 &
+                     trn(ji-1,jj+1,jk,jpdic), trn(ji,jj+1,jk,jpdic), trn(ji+1,jj+1,jk,jpdic)
+                      IF(lwp) WRITE(numout,*)                                 &
+                     trn(ji-1,jj,jk,jpdic), trn(ji,jj,jk,jpdic), trn(ji+1,jj,jk,jpdic)
+                      IF(lwp) WRITE(numout,*)                                 &
+                     trn(ji-1,jj-1,jk,jpdic), trn(ji,jj-1,jk,jpdic), trn(ji+1,jj-1,jk,jpdic)
+                      IF(lwp) WRITE(numout,*)                                 &
+                     ' trc_bio_medusa: Abnormal T suroundings Alk'
+                      IF(lwp) WRITE(numout,*)                                 &
+                     trn(ji-1,jj+1,jk,jpalk), trn(ji,jj+1,jk,jpalk), trn(ji+1,jj+1,jk,jpalk)
+                      IF(lwp) WRITE(numout,*)                                 &
+                     trn(ji-1,jj,jk,jpalk), trn(ji,jj,jk,jpalk), trn(ji+1,jj,jk,jpalk)
+                      IF(lwp) WRITE(numout,*)                                 &
+                     trn(ji-1,jj-1,jk,jpalk), trn(ji,jj-1,jk,jpalk), trn(ji+1,jj-1,jk,jpalk)
+                      IF(lwp) WRITE(numout,*)                                 &
+                     ' trc_bio_medusa: Abnormal T suroundings tmask'
+                      IF(lwp) WRITE(numout,*)                                 &
+                     tmask(ji-1,jj+1,jk), tmask(ji,jj+1,jk), tmask(ji+1,jj+1,jk)
+                      IF(lwp) WRITE(numout,*)                                 &
+                     tmask(ji-1,jj,jk), tmask(ji,jj,jk), tmask(ji+1,jj,jk)
+                      IF(lwp) WRITE(numout,*)                                 &
+                     tmask(ji-1,jj-1,jk), tmask(ji,jj-1,jk), tmask(ji+1,jj-1,jk)
+                     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) 
+                     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,*)                           &
+                           ' trc_bio_medusa: T SWITCHING 3D, ',              &
+                           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,*)                           &
+                           ' trc_bio_medusa: T SWITCHING 3D, ',              &
+                           tsn(ji,jj,jk,jp_tem), ' -> 40.0 '
+                           ztmp(ji,jj) = 40.0
+                        ENDIF
+                     ELSE 
+                        IF(lwp) WRITE(numout,*)                              &
+                        ' trc_bio_medusa: T SWITCHING 3D, ',                 &
+                        tsn(ji,jj,jk,jp_tem), ' -> surounding avg : ', tsnavg
+                        ztmp(ji,jj) = tsnavg
+                     ENDIF
                   endif
-                  if (zsal(ji,jj) .lt. 0.0 .or. zsal(ji,jj) .gt. 45.0 ) then
+                  !! end T chack
+                  if (zsal(ji,jj) .lt. 1.0 .or. zsal(ji,jj) .gt. 47.0 ) then
                      IF(lwp) WRITE(numout,*)                                 &
                         ' 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
+                     !! Correct out of range values
+                     IF (zsal(ji,jj) .LT. 1.0 ) THEN
+                        IF(lwp) WRITE(numout,*)                              &
+                        ' trc_bio_medusa: S SWITCHING 3D, ',                 &
+                        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,*)                              &
+                        ' trc_bio_medusa: T SWITCHING 3D, ',                 &
+                        tsn(ji,jj,jk,jp_sal), ' -> 47.0 '
+                        zsal(ji,jj) = 47.0
+                     ENDIF
                   endif
+                  !! end S check
                ENDIF
             ENDDO