Changeset 4070 for trunk/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

Timestamp:

2013-10-17T12:48:23+02:00 (11 years ago)

Author:

djlea

Message:

#1087: Bug fix to freezing temperature used in ASM code.

File:

• trunk/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90 (modified) (1 diff)

trunk/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90

@@ -682 +682 @@
       ! used to prevent the applied increments taking the temperature below the local freezing point 
 
-#if defined key_cice 
-        fzptnz(:,:,:) = -1.8_wp
-#else 
-        DO jk = 1, jpk
-           DO jj = 1, jpj
-              DO ji = 1, jpk
-                 fzptnz (ji,jj,jk) = ( -0.0575_wp + 1.710523e-3_wp * SQRT( tsn(ji,jj,jk,jp_sal) )                   & 
-                                                  - 2.154996e-4_wp *       tsn(ji,jj,jk,jp_sal)   ) * tsn(ji,jj,jk,jp_sal)  & 
-                                                  - 7.53e-4_wp * fsdepw(ji,jj,jk)       ! (pressure in dbar) 
-              END DO
-           END DO
-        END DO
-#endif 
+      DO jk=1, jpkm1
+         fzptnz (:,:,jk) = tfreez( tsn(:,:,jk,jp_sal), fsdept(:,:,jk) )
+      ENDDO
 
       IF ( ln_asmiau ) THEN