Changeset 11935

Timestamp:

2019-11-20T11:37:14+01:00 (4 years ago)

Author:

dcarneir

Message:

Changing SBC scripts (sbc_oce and sbcice_cice) to include sea ice thickness

Location:

branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/SBC

Files:

• sbc_oce.F90 (modified) (2 diffs)
• sbcice_cice.F90 (modified) (1 diff)

branches/UKMO/dev_r5518_obs_oper_update_sit/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90

@@ -121 +121 @@
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   sprecip           !: solid precipitation                          [Kg/m2/s]
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   fr_i              !: ice fraction = 1 - lead fraction      (between 0 to 1)
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   thick_i           !: ice thickness [m]
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   thick_s           !: snow depth [m]
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   vol_i             !: ice volume [m3]
 #if defined key_cpl_carbon_cycle
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   atm_co2           !: atmospheric pCO2                             [ppm]
@@ -167 +170 @@
          !
       ALLOCATE( tprecip(jpi,jpj) , sprecip(jpi,jpj) , fr_i(jpi,jpj) ,     &
+         &      thick_i(jpi,jpj) , thick_s(jpi,jpj) , vol_i(jpi,jpj),     &
 #if defined key_cpl_carbon_cycle
          &      atm_co2(jpi,jpj) ,                                        &

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

@@ -229 +229 @@
       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