Changeset 6675
- Timestamp:
- 2016-06-08T17:54:56+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_coupling_GSI7_GSI8_landice_bitcomp/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
r6458 r6675 15 15 USE dom_oce ! ocean space and time domain 16 16 USE domvvl 17 USE eosbn2, only : eos_fzp ! Subroutineto calculate freezing point of seawater17 USE eosbn2, only : eos_fzp ! Function to calculate freezing point of seawater 18 18 USE phycst, only : rcp, rau0, r1_rau0, rhosn, rhoic, rt0 19 19 USE in_out_manager ! I/O manager … … 208 208 CALL wrk_alloc( jpi,jpj,jpk, ztfrz3d ) 209 209 DO jk=1,jpk 210 CALL eos_fzp( tsn(:,:,jk,jp_sal),ztfrz3d(:,:,jk),fsdept_n(:,:,jk) )210 ztfrz3d(:,:,jk) = eos_fzp( tsn(:,:,jk,jp_sal), fsdept_n(:,:,jk) ) 211 211 ENDDO 212 212 tsn(:,:,:,jp_tem) = MAX( tsn(:,:,:,jp_tem), ztfrz3d ) … … 224 224 225 225 ! calculate surface freezing temperature and send to CICE 226 CALL eos_fzp( sss_m(:,:), sstfrz(:,:), fsdept_n(:,:,1) )226 sstfrz(:,:) = eos_fzp(sss_m(:,:), fsdept_n(:,:,1)) 227 227 CALL nemo2cice(sstfrz,Tf, 'T', 1. ) 228 228 … … 462 462 463 463 ! Recalculate freezing temperature and send to CICE 464 CALL eos_fzp( sss_m(:,:), sstfrz(:,:), fsdept_n(:,:,1) )464 sstfrz(:,:)=eos_fzp(sss_m(:,:), fsdept_n(:,:,1)) 465 465 CALL nemo2cice(sstfrz,Tf,'T', 1. ) 466 466
Note: See TracChangeset
for help on using the changeset viewer.