Changeset 8400 for branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90

Timestamp:

2017-08-01T16:12:29+02:00 (7 years ago)

Author:

timgraham

Message:

GMED ticket 335:

• Merge dev_r5518_GO6_package_inc_asm into package branch to make everything easier for data assimilation
• No effect on configs without data assimilation

File:

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

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

@@ -56 +56 @@
                 fresh_ai,fhocn_ai,fswthru_ai,frzmlt,             &
                 flatn_f,fsurfn_f,fcondtopn_f,                    &
+#ifdef key_asminc
+                daice_da,fresh_da,fsalt_da,                    &
+#endif
                 uatm,vatm,wind,fsw,flw,Tair,potT,Qa,rhoa,zlvl,   &
                 swvdr,swvdf,swidr,swidf,Tf,                      &
@@ -301 +304 @@
  
       CALL wrk_dealloc( jpi,jpj, ztmp1, ztmp2 )
+
+#if defined key_asminc
+      ! 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         
+#endif
       !
       ! In coupled mode get extra fields from CICE for passing back to atmosphere
@@ -454 +465 @@
       ENDIF
 
+#if defined key_asminc
+!Ice concentration change (from assimilation)
+      ztmp(:,:)=ndaice_da(:,:)*tmask(:,:,1)
+      Call nemo2cice(ztmp,daice_da,'T', 1. )
+#endif 
+
 ! Snowfall
 ! Ensure fsnow is positive (as in CICE routine prepare_forcing)
@@ -716 +733 @@
          snwice_fmass (:,:) = ( snwice_mass(:,:) - snwice_mass_b(:,:) ) / dt
       ENDIF
+
+#if defined key_asminc
+! 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)
+#endif
 
 ! Release work space