Changeset 10276 for branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/OPA_SRC/SBC

Timestamp:

2018-11-02T18:15:27+01:00 (5 years ago)

Author:

emmafiedler

Message:

Freeboard assimilation updates

Location:

branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/OPA_SRC/SBC

Files:

• sbc_ice.F90 (modified) (2 diffs)
• sbc_oce.F90 (modified) (2 diffs)
• sbcice_cice.F90 (modified) (3 diffs)

branches/UKMO/dev_r5518_obs_oper_update_icethick/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_ice.F90

@@ -107 +107 @@
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   thick_iu           !: ice thickness at NEMO U point
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   thick_iv           !: ice thickness at NEMO V point
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   thick_su           !: snow depth at NEMO U point
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:)   ::   thick_sv           !: snow depth at NEMO V point   
    ! variables used in the coupled interface
    INTEGER , PUBLIC, PARAMETER ::   jpl = ncat
@@ -174 +176 @@
                 ss_iov(jpi,jpj)       , fr_iu(jpi,jpj)        , fr_iv(jpi,jpj)        , &
                 thick_iu(jpi,jpj)     , thick_iv(jpi,jpj)     ,                         &
+                thick_su(jpi,jpj)     , thick_sv(jpi,jpj)     ,                         &                
                 ht_i(jpi,jpj,ncat)    , ht_s(jpi,jpj,ncat)    ,                         &                
                 a_i(jpi,jpj,ncat)     , topmelt(jpi,jpj,ncat) , botmelt(jpi,jpj,ncat) , &

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

@@ -122 +122 @@
    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 - I think]   
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   vol_i             !: ice volume [m3]   
 #if defined key_cpl_carbon_cycle
@@ -198 +199 @@
          !
       ALLOCATE( tprecip(jpi,jpj) , sprecip(jpi,jpj) , fr_i(jpi,jpj) ,     &
-         & thick_i(jpi,jpj) ,  vol_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_icethick/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

@@ -251 +251 @@
       CALL lbc_lnk ( fr_iv , 'V', 1. )
       
-!!*****JUST ADDED BUT NEEDS CONFIRMING!!****
+
 ! Snow and ice thickness
-      CALL cice2nemo(vice,thick_i,'T', 1. )  !!Do this for snow depth too? set thick_s. Duplicated below for embedded ice?
-
-! vice is volume per unit area of grid cell = 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 ! Confirm what this if loop is for
+      IF ( (ksbc == jp_flx) .OR. (ksbc == jp_purecpl) ) THEN
          DO jl = 1,ncat
             CALL cice2nemo(vsnon(:,:,jl,:),ht_s(:,:,jl),'T', 1. )
@@ -263 +263 @@
          ENDDO
       ENDIF
-      
+
+
 ! T point to U point
 ! T point to V point
+
+! Sea ice thickness
       thick_iu(:,:)=0.0
       thick_iv(:,:)=0.0
@@ -278 +281 @@
       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. )