- Timestamp:
- 2015-12-24T17:32:12+01:00 (9 years ago)
- Location:
- branches/UKMO/dev_r5518_coupling_GSI7_GSI8_landice_with_key_asminc/NEMOGCM/NEMO/OPA_SRC/SBC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_coupling_GSI7_GSI8_landice_with_key_asminc/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_ice.F90
r5663 r6172 105 105 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: kn_ice !: sea-ice surface layer thermal conductivity (on cats) 106 106 107 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ndaice_da !: NEMO fresh water flux to ocean due to data assim 108 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: nfresh_da !: NEMO salt flux to ocean due to data assim 109 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: nfsalt_da !: NEMO ice concentration change/second from data assim 110 107 111 ! variables used in the coupled interface 108 112 INTEGER , PUBLIC, PARAMETER :: jpl = ncat … … 161 165 ss_iov(jpi,jpj) , fr_iu(jpi,jpj) , fr_iv(jpi,jpj) , & 162 166 a_i(jpi,jpj,ncat) , topmelt(jpi,jpj,ncat) , botmelt(jpi,jpj,ncat) , & 167 ndaice_da(jpi,jpj) , nfresh_da(jpi,jpj) , nfsalt_da(jpi,jpj) , & 163 168 sstfrz(jpi,jpj) , STAT= ierr(1) ) 164 169 ! Alex West: Allocating tn_ice with 5 categories. When NEMO is used with CICE, this variable -
branches/UKMO/dev_r5518_coupling_GSI7_GSI8_landice_with_key_asminc/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
r6023 r6172 55 55 sst,sss,uocn,vocn,ss_tltx,ss_tlty,fsalt_ai, & 56 56 fresh_ai,fhocn_ai,fswthru_ai,frzmlt, & 57 daice_da,fresh_da,fsalt_da, & 57 58 flatn_f,fsurfn_f,fcondtopn_f, & 58 59 uatm,vatm,wind,fsw,flw,Tair,potT,Qa,rhoa,zlvl, & … … 301 302 302 303 CALL wrk_dealloc( jpi,jpj, ztmp1, ztmp2 ) 304 ! Initialize fresh water and salt fluxes from data assim 305 ! and data assimilation index to cice 306 nfresh_da(:,:) = 0.0 307 nfsalt_da(:,:) = 0.0 308 ndaice_da(:,:) = 0.0 303 309 ! 304 310 IF( nn_timing == 1 ) CALL timing_stop('cice_sbc_init') … … 450 456 ENDIF 451 457 458 #if defined key_asminc 459 !Ice concentration change (from assimilation) 460 ztmp(:,:)=ndaice_da(:,:)*tmask(:,:,1) 461 Call nemo2cice(ztmp,daice_da,'T', 1. ) 462 #endif 463 452 464 ! Snowfall 453 465 ! Ensure fsnow is positive (as in CICE routine prepare_forcing) … … 712 724 snwice_fmass (:,:) = ( snwice_mass(:,:) - snwice_mass_b(:,:) ) / dt 713 725 ENDIF 726 727 #if defined key_asminc 728 ! Import fresh water and salt flux due to seaice da 729 CALL cice2nemo(fresh_da, nfresh_da,'T',1.0) 730 CALL cice2nemo(fsalt_da, nfsalt_da,'T',1.0) 731 #endif 714 732 715 733 ! Release work space
Note: See TracChangeset
for help on using the changeset viewer.