Changeset 8618

Timestamp:

2017-10-12T12:49:28+02:00 (7 years ago)

Author:

hadlh

Message:

Changes for DePreSys3 (coupled) sea-ice nudging, including code to read in ancillary file and pass to CICE.

Location:

branches/UKMO/GO6_package_r8356_CICEnudging/NEMOGCM/NEMO/OPA_SRC/SBC

Files:

• sbc_ice.F90 (modified) (2 diffs)
• sbcice_cice.F90 (modified) (4 diffs)
• sbcmod.F90 (modified) (2 diffs)

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

@@ -72 +72 @@
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   topmelt            !: category topmelt
    REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   botmelt            !: category botmelt
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   ndaice_da        !: NEMO ice conc rate of change
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   nfresh_da         !: NEMO fresh water flux to ocean due to data assim 
+   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::   nfsalt_da          !: NEMO salt flux to ocean due to data assim 
 
 #if defined  key_lim3
@@ -162 +165 @@
                 ss_iov(jpi,jpj)       , fr_iu(jpi,jpj)        , fr_iv(jpi,jpj)        , &
                 a_i(jpi,jpj,ncat)     , topmelt(jpi,jpj,ncat) , botmelt(jpi,jpj,ncat) , &
-                sstfrz(jpi,jpj)       , STAT= ierr(1) )
+                sstfrz(jpi,jpj)       , ndaice_da(jpi,jpj), nfresh_da(jpi,jpj), &
+      nfsalt_da(jpi,jpj), STAT= ierr(1) )
    ! Alex West: Allocating tn_ice with 5 categories.  When NEMO is used with CICE, this variable
    ! represents top layer ice temperature, which is multi-category.

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

@@ -55 +55 @@
                 sst,sss,uocn,vocn,ss_tltx,ss_tlty,fsalt_ai,      &
                 fresh_ai,fhocn_ai,fswthru_ai,frzmlt,             &
-                flatn_f,fsurfn_f,fcondtopn_f,                    &
+                flatn_f,fsurfn_f,fcondtopn_f,daice_da,fresh_da,fsalt_da, &
                 uatm,vatm,wind,fsw,flw,Tair,potT,Qa,rhoa,zlvl,   &
                 swvdr,swvdf,swidr,swidf,Tf,                      &
@@ -223 +223 @@
       ENDIF  
 
+      ! Initialize fresh water and salt fluxes from data assim  
+      !  and data assimilation index to cice
+      nfresh_da(:,:) = 0.0
+      nfsalt_da(:,:) = 0.0
+      ndaice_da(:,:) = 0.0
+      
       ! calculate surface freezing temperature and send to CICE
       CALL  eos_fzp(sss_m(:,:), sstfrz(:,:), fsdept_n(:,:,1)) 
@@ -454 +460 @@
       ENDIF
 
+!Ice concentration change (from assimilation) 
+      ztmp(:,:)=ndaice_da(:,:)*tmask(:,:,1)  
+      write(numout,*) 'PE/Ar/Max/Min daice_dat ', nproc, narea, maxval(ztmp), minval(ztmp) 
+      write(numout,*) 'PE/Ar/Max/Min daice_dat ', nproc, narea, maxval(ndaice_da), minval(ndaice_da)
+      CALL nemo2cice(ztmp, daice_da,'T',1.)
+
 ! Snowfall
 ! Ensure fsnow is positive (as in CICE routine prepare_forcing)
@@ -716 +728 @@
          snwice_fmass (:,:) = ( snwice_mass(:,:) - snwice_mass_b(:,:) ) / dt
       ENDIF
+
+! Import fresh water and salt flux due to seaice da
+      CALL cice2nemo(fresh_da, nfresh_da,'T',1.0)
+      CALL cice2nemo(fsalt_da, nfsalt_da,'T',1.0)
+      write(numout,*) 'Max/Min nfresh_da ', maxval(nfresh_da), minval(nfresh_da)
+      write(numout,*) 'Max/Min nfresh_da ', maxval(nfsalt_da), minval(nfsalt_da)
 
 ! Release work space

branches/UKMO/GO6_package_r8356_CICEnudging/NEMOGCM/NEMO/OPA_SRC/SBC/sbcmod.F90

@@ -44 +44 @@
    USE sbcisf           ! surface boundary condition: ice shelf
    USE sbcfwb           ! surface boundary condition: freshwater budget
+   USE sbcdta_ice       ! nudging: load in ice ancillary files
    USE closea           ! closed sea
    USE icbstp           ! Icebergs!
@@ -390 +391 @@
 
       IF( ln_mixcpl )        CALL sbc_cpl_rcv ( kt, nn_fsbc, nn_ice )   ! forced-coupled mixed formulation after forcing
+      
+      CALL sbc_dta_ice ( kt )            ! nudging: calculate RATE of change of ice concentration due to DA