Changeset 12613 for branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

Timestamp:

2020-03-26T12:30:42+01:00 (4 years ago)

Author:

dcarneir

Message:

Including sea ice thickness in GO6 branch

File:

• branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90 (modified) (4 diffs)

branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

@@ -57 +57 @@
                 flatn_f,fsurfn_f,fcondtopn_f,                    &
 #ifdef key_asminc
-                daice_da,fresh_da,fsalt_da,                    &
+                daice_da,dhi_da,fresh_da,fsalt_da,               &
 #endif
                 uatm,vatm,wind,fsw,flw,Tair,potT,Qa,rhoa,zlvl,   &
@@ -250 +250 @@
       CALL lbc_lnk ( fr_iu , 'U', 1. )
       CALL lbc_lnk ( fr_iv , 'V', 1. )
+
+! Snow and ice thickness
+! vice, vsno is volume per unit area of grid cell = thickness
+      CALL cice2nemo(vice,thick_i,'T', 1. )
+      CALL cice2nemo(vsno,thick_s,'T', 1. )
+     
+      IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN
+         DO jl = 1,ncat
+            CALL cice2nemo(vsnon(:,:,jl,:),ht_s(:,:,jl),'T', 1. )
+            CALL cice2nemo(vicen(:,:,jl,:),ht_i(:,:,jl),'T', 1. )
+         ENDDO
+      ENDIF
+
+! T point to U point
+! T point to V point
+
+! Sea ice thickness
+      thick_iu(:,:)=0.0
+      thick_iv(:,:)=0.0
+      DO jj=1,jpjm1
+         DO ji=1,jpim1
+            thick_iu(ji,jj)=0.5*(thick_i(ji,jj)+thick_i(ji+1,jj))*umask(ji,jj,1)
+            thick_iv(ji,jj)=0.5*(thick_i(ji,jj)+thick_i(ji,jj+1))*vmask(ji,jj,1)
+         ENDDO
+      ENDDO
+
+      CALL lbc_lnk ( thick_iu , 'U', 1. )
+      CALL lbc_lnk ( thick_iv , 'V', 1. )
+
+! Snow depth
+      thick_su(:,:)=0.0
+      thick_sv(:,:)=0.0
+      DO jj=1,jpjm1
+         DO ji=1,jpim1
+            thick_su(ji,jj)=0.5*(thick_s(ji,jj)+thick_s(ji+1,jj))*umask(ji,jj,1)
+            thick_sv(ji,jj)=0.5*(thick_s(ji,jj)+thick_s(ji,jj+1))*vmask(ji,jj,1)
+         ENDDO
+      ENDDO
+
+      CALL lbc_lnk ( thick_su , 'U', 1. )
+      CALL lbc_lnk ( thick_sv , 'V', 1. )     
 
       !                                      ! embedded sea ice
@@ -311 +352 @@
       nfsalt_da(:,:) = 0.0   
       ndaice_da(:,:) = 0.0         
+      ndsit_da(:,:)  = 0.0         
 #endif
       !
@@ -469 +511 @@
       ztmp(:,:)=ndaice_da(:,:)*tmask(:,:,1)
       Call nemo2cice(ztmp,daice_da,'T', 1. )
+!Ice thickness change (from assimilation)
+      ztmp(:,:)=ndsit_da(:,:)*tmask(:,:,1)
+      Call nemo2cice(ztmp,dhi_da,'T', 1. )
 #endif