Changeset 8899
- Timestamp:
- 2017-12-05T14:26:11+01:00 (7 years ago)
- Location:
- branches/UKMO/GO6_package_r8356_CICEnudging/NEMOGCM/NEMO/OPA_SRC/SBC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/GO6_package_r8356_CICEnudging/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_ice.F90
r8618 r8899 73 73 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: botmelt !: category botmelt 74 74 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: ndaice_da !: NEMO ice conc rate of change 75 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: nfresh_da !: NEMO fresh water flux to ocean due to data assim76 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: nfsalt_da !: NEMO salt flux to ocean due to data assim77 75 78 76 #if defined key_lim3 … … 165 163 ss_iov(jpi,jpj) , fr_iu(jpi,jpj) , fr_iv(jpi,jpj) , & 166 164 a_i(jpi,jpj,ncat) , topmelt(jpi,jpj,ncat) , botmelt(jpi,jpj,ncat) , & 167 sstfrz(jpi,jpj) , ndaice_da(jpi,jpj), nfresh_da(jpi,jpj), & 168 nfsalt_da(jpi,jpj), STAT= ierr(1) ) 165 sstfrz(jpi,jpj) , ndaice_da(jpi,jpj), STAT= ierr(1) ) 169 166 ! Alex West: Allocating tn_ice with 5 categories. When NEMO is used with CICE, this variable 170 167 ! represents top layer ice temperature, which is multi-category. -
branches/UKMO/GO6_package_r8356_CICEnudging/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_cice.F90
r8618 r8899 55 55 sst,sss,uocn,vocn,ss_tltx,ss_tlty,fsalt_ai, & 56 56 fresh_ai,fhocn_ai,fswthru_ai,frzmlt, & 57 flatn_f,fsurfn_f,fcondtopn_f,daice_da, fresh_da,fsalt_da,&57 flatn_f,fsurfn_f,fcondtopn_f,daice_da, & 58 58 uatm,vatm,wind,fsw,flw,Tair,potT,Qa,rhoa,zlvl, & 59 59 swvdr,swvdf,swidr,swidf,Tf, & … … 223 223 ENDIF 224 224 225 ! Initialize fresh water and salt fluxes from data assim 226 ! and data assimilation index to cice 227 nfresh_da(:,:) = 0.0 228 nfsalt_da(:,:) = 0.0 225 ! Initialize data assimilation index for cice 229 226 ndaice_da(:,:) = 0.0 230 227 … … 462 459 !Ice concentration change (from assimilation) 463 460 ztmp(:,:)=ndaice_da(:,:)*tmask(:,:,1) 464 write(numout,*) 'PE/Ar/Max/Min daice_dat ', nproc, narea, maxval(ztmp), minval(ztmp)465 write(numout,*) 'PE/Ar/Max/Min daice_dat ', nproc, narea, maxval(ndaice_da), minval(ndaice_da)466 461 CALL nemo2cice(ztmp, daice_da,'T',1.) 467 462 … … 728 723 snwice_fmass (:,:) = ( snwice_mass(:,:) - snwice_mass_b(:,:) ) / dt 729 724 ENDIF 730 731 ! Import fresh water and salt flux due to seaice da732 CALL cice2nemo(fresh_da, nfresh_da,'T',1.0)733 CALL cice2nemo(fsalt_da, nfsalt_da,'T',1.0)734 write(numout,*) 'Max/Min nfresh_da ', maxval(nfresh_da), minval(nfresh_da)735 write(numout,*) 'Max/Min nfresh_da ', maxval(nfsalt_da), minval(nfsalt_da)736 725 737 726 ! Release work space
Note: See TracChangeset
for help on using the changeset viewer.