Changeset 5162 for branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

Timestamp:

2015-03-23T16:57:10+01:00 (9 years ago)

Author:

davestorkey

Message:

Bug fixes for 2015/dev_r5021_UKMO1_CICE_coupling

File:

• branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90 (modified) (5 diffs)

branches/2015/dev_r5021_UKMO1_CICE_coupling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

@@ -45 +45 @@
                 flatn_f,fsurfn_f,fcondtopn_f,                    &
                 uatm,vatm,wind,fsw,flw,Tair,potT,Qa,rhoa,zlvl,   &
-                swvdr,swvdf,swidr,swidf
+                swvdr,swvdf,swidr,swidf,Tf
    USE ice_therm_vertical, only: calc_Tsfc
 #else
@@ -55 +55 @@
                 flatn_f,fsurfn_f,fcondtopn_f,                    &
                 uatm,vatm,wind,fsw,flw,Tair,potT,Qa,rhoa,zlvl,   &
-                swvdr,swvdf,swidr,swidf
+                swvdr,swvdf,swidr,swidf,Tf
    USE ice_therm_shared, only: calc_Tsfc
    USE ice_shortwave, only: apeffn
@@ -179 +179 @@
       CALL wrk_alloc( jpi,jpj,jpk, ztfrz3d ) 
       DO jk=1,jpk
-         ztfrz3d(:,:,jk) = eos_fzp( tsn(:,:,jk,jp_sal), fsdewq(:,:,jk) )
-      ENDDO
-      sstfrz(:,:)=ztfrz(:,:,1)
+         ztfrz3d(:,:,jk) = eos_fzp( tsn(:,:,jk,jp_sal), fsdept_n(:,:,jk) )
+      ENDDO
+      sstfrz(:,:)=ztfrz3d(:,:,1)
 
       !Ensure that no temperature points are below freezing if not a NEMO restart
@@ -440 +440 @@
 
 ! Recalculate freezing temperature and send to CICE 
-      sstfrz(:,:)=tfreez(sss_m(:,:)) 
+      sstfrz(:,:)=eos_fzp(sss_m(:,:), fsdept_n(:,:,1)) 
       CALL nemo2cice(sstfrz,Tf,'T', 1. )
 
@@ -724 +724 @@
       DO jl = 1,ncat
          CALL cice2nemo(apeffn(:,:,jl,:),a_p(:,:,jl),'T', 1. )
-         CALL cice2nemo(trcrn(:,:,jl,:),ht_p(:,:,jl),'T', 1. )
+         CALL cice2nemo(trcrn(:,:,nt_hpnd,jl,:),ht_p(:,:,jl),'T', 1. )
       ENDDO
 #endif