Changeset 7614 for branches/UKMO

Timestamp:

2017-01-27T10:15:11+01:00 (7 years ago)

Author:

kingr

Message:

Code changes to prevent water temperature falling below pressure-dependent freezing point

File:

• branches/UKMO/dev_r5518_min_temp_wo_ice/NEMOGCM/NEMO/OPA_SRC/TRA/tranxt.F90 (modified) (3 diffs)

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

@@ -58 +58 @@
 
    REAL(wp) ::   rbcp   ! Brown & Campana parameters for semi-implicit hpg
+   INTEGER  ::   warn_1, warn_2   ! indicators for warning statement
 
    !! * Substitutions
@@ -93 +94 @@
       INTEGER, INTENT(in) ::   kt    ! ocean time-step index
       !!
-      INTEGER  ::   jk, jn    ! dummy loop indices
-      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
       !!----------------------------------------------------------------------
@@ -151 +152 @@
       ENDIF     
       !
+#if ( ! defined key_lim3 && ! defined key_lim2 && ! key_cice ) 
+      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 
+                  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 
+          CALL ctl_warn( ' Temperatures dropping below freezing point, ', & 
+                    &    ' being forced to freezing point, no longer conservative' ) 
+        ENDIF 
+        warn_1=1 
+      ENDIF 
+#endif  
      ! trends computation
       IF( l_trdtra ) THEN      ! trend of the Asselin filter (tb filtered - tb)/dt