Changeset 7626

Timestamp:

2017-01-31T17:06:10+01:00 (7 years ago)

Author:

hadjt

Message:

Added CO5 code to catch water below freezing, and reset to freezing point - based on code by Rachel Furner/Rob? King (​http://fcm3/projects/NEMO/ticket/682). Also added a name list to switch this code on, and how it reports any freezing grid boxes - nam_tranxticing: nn_tranxticing; (0=Off; 1=Error; 2 Warn (verbose); 3 Warn (Quiet))

File:

• branches/UKMO/CO6_shelfclimate/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90 (modified) (4 diffs)

branches/UKMO/CO6_shelfclimate/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90

@@ -58 +58 @@
 
    REAL(wp) ::   rbcp   ! Brown & Campana parameters for semi-implicit hpg
+   !JT 
+   INTEGER  ::   warn_1, warn_2   ! indicators for warning statement
+   INTEGER , PUBLIC ::   nn_tranxticing  ! region mean calculation
+   !JT
 
    !! * Substitutions
@@ -93 +97 @@
       INTEGER, INTENT(in) ::   kt    ! ocean time-step index
       !!
-      INTEGER  ::   jk, jn    ! dummy loop indices
-      REAL(wp) ::   zfact     ! local scalars
+      !JT INTEGER  ::   jk, jn    ! dummy loop indices
+      !JT REAL(wp) ::   zfact     ! local scalars
+      INTEGER  ::   ji,jj,jk, jn     ! dummy loop indices
+      REAL(wp) ::   zfact, zfreeze   ! local scalars
       REAL(wp), POINTER, DIMENSION(:,:,:) ::  ztrdt, ztrds
+      INTEGER ::   ios                  ! Local integer output status for namelist read
       !!----------------------------------------------------------------------
       !
@@ -106 +113 @@
          !
          rbcp = 0.25_wp * (1._wp + atfp) * (1._wp + atfp) * ( 1._wp - atfp)      ! Brown & Campana parameter for semi-implicit hpg
+         
+         
+         
+
+          NAMELIST/nam_tranxticing/ nn_tranxticing
+
+          REWIND ( numnam_ref )              ! Read Namelist nam_diatmb in referdiatmbence namelist : TMB diagnostics
+          READ   ( numnam_ref, nam_tranxticing, IOSTAT=ios, ERR= 901 )
+901       IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_tranxticing in reference namelist', lwp )
+
+          REWIND( numnam_cfg )              ! Namelist nam_diatmb in configuration namelist  TMB diagnostics
+          READ  ( numnam_cfg, nam_tranxticing, IOSTAT = ios, ERR = 902 )
+902       IF( ios /= 0 ) CALL ctl_nam ( ios , 'nam_tranxticing in configuration namelist', lwp )
+          IF(lwm) WRITE ( numond, nam_tranxticing )
+
+          IF(lwp) THEN                   ! Control print
+              WRITE(numout,*)
+              WRITE(numout,*) 'tra_nxt : Catching Icing'
+              WRITE(numout,*) '~~~~~~~~~~~~'
+              WRITE(numout,*) 'Namelist nam_tranxticing : catch freezing points '
+              WRITE(numout,*) 'Int switch for icing catch (0=Off; 1=Error; 2 Warn (verbose); 3 Warn (Quiet); nn_tranxticing  = ', nn_tranxticing
+          ENDIF
+      
+         
       ENDIF
 
@@ -151 +182 @@
       ENDIF     
       !
+#if ( ! defined key_lim3 && ! defined key_lim2 && ! key_cice )
+
+      IF (nn_tranxticing .GT. 0) THEN
+        IF ( kt == nit000 ) warn_1=0 
+        warn_2=0 
+        DO jk = 1, jpkm1 
+          DO jj = 1, jpj 
+            DO ji = 1, jpi 
+              IF ( tsa(ji,jj,jk,jp_tem) .lt. 0.0 ) THEN 
+                ! calculate the freezing point
+                zfreeze = ( -0.0575_wp + 1.710523E-3 * Sqrt (Abs(tsn(ji,jj,jk,jp_sal)))   & 
+                          - 2.154996E-4 * tsn(ji,jj,jk,jp_sal) ) * tsn(ji,jj,jk,jp_sal) - 7.53E-4 * ( 10.0_wp + fsdept(ji,jj,jk) ) 
+                IF ( tsn(ji,jj,jk,jp_tem) .lt. zfreeze ) THEN
+                    IF (nn_tranxticing .lt. 3) THEN
+                  
+                      WRITE(numout,300) ' tranxticing: ', & 
+                      &   kt,zfreeze,tsn(ji,jj,jk,jp_tem),ji,jj,jk,narea,ji-1+njmpp-1+jpjzoom,jj-1+nimpp-1+jpizoom
+
+300 FORMAT(A14,1X,I7,1X,f9.2,1X,f9.2,1X,I4,1X,I4,1X,I4,1X,I4,1X,I4,1X,I4)
+                                        
+                    ENDIF
+                    tsn(ji,jj,jk,jp_tem)=zfreeze 
+                    warn_2=1 
+                ENDIF
+              ENDIF
+            END DO
+          END DO
+        END DO
+      CALL mpp_max(warn_1) 
+      CALL mpp_max(warn_2) 
+        IF ( (warn_1 == 0) .AND. (warn_2 /= 0) ) THEN
+          IF(lwp) THEN
+          
+            IF (nn_tranxticing .GT. 1) CALL ctl_warn( ' tranxticing: Temperatures dropping below freezing point, ', & 
+                      &    ' being forced to freezing point, no longer conservative' ) 
+                      
+            IF (nn_tranxticing .EQ. 1) CALL ctl_stop( ' tranxticing: Temperatures dropping below freezing point, ', & 
+                      &    ' being forced to freezing point, no longer conservative' ) 
+                      
+          ENDIF
+          warn_1=1 
+        ENDIF 
+      ENDIF 
+#endif  
+      
      ! trends computation
       IF( l_trdtra ) THEN      ! trend of the Asselin filter (tb filtered - tb)/dt