- Timestamp:
- 2013-11-07T11:19:49+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_LOCEAN_2013/NEMOGCM/NEMO/OPA_SRC/ASM/asminc.F90
r4147 r4162 675 675 ! used to prevent the applied increments taking the temperature below the local freezing point 676 676 677 #if defined key_cice 678 fzptnz(:,:,:) = -1.8_wp 679 #else 680 DO jk = 1, jpk 681 DO jj = 1, jpj 682 DO ji = 1, jpk 683 fzptnz (ji,jj,jk) = ( -0.0575_wp + 1.710523e-3_wp * SQRT( tsn(ji,jj,jk,jp_sal) ) & 684 - 2.154996e-4_wp * tsn(ji,jj,jk,jp_sal) ) * tsn(ji,jj,jk,jp_sal) & 685 - 7.53e-4_wp * fsdepw(ji,jj,jk) ! (pressure in dbar) 686 END DO 687 END DO 688 END DO 689 #endif 677 DO jk=1, jpkm1 678 fzptnz (:,:,jk) = tfreez( tsn(:,:,jk,jp_sal), fsdept(:,:,jk) ) 679 ENDDO 690 680 691 681 IF ( ln_asmiau ) THEN
Note: See TracChangeset
for help on using the changeset viewer.