Changeset 878 for branches/dev_001_SBC/NEMO
- Timestamp:
- 2008-04-04T14:53:06+02:00 (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_001_SBC/NEMO/LIM_SRC/limsbc.F90
r872 r878 152 152 ! ! ice-covered fraction: 153 153 154 !!gm old zemp = + evap(ji,jj) * frld(ji,jj) & ! evaporation over oceanic fraction155 !!gm old & - tprecip(ji,jj) & ! total precipitation156 !!gm old & + sprecip(ji,jj) * ( 1. - pfrld(ji,jj) ) & ! remov. snow precip over ice (store over sea-ice)157 !!gm old & + rdmsnif(ji,jj) / rdt_ice ! freshwaterflux due to snow melting158 159 154 ! computing salt exchanges at the ice/ocean interface 160 155 zfons = ( soce - sice ) * ( rdmicif(ji,jj) / rdt_ice ) 161 156 162 157 ! converting the salt flux from ice to a freshwater flux from ocean 163 !!gm old zfm = zfons / ( sss_io(ji,jj) + epsi16 )164 158 zfm = zfons / ( sss_m(ji,jj) + epsi16 ) 165 159 … … 185 179 DO ji = 1, jpi 186 180 ! ... change the cosinus angle sign in the south hemisphere 187 !! zsang = SIGN(1.e0, gphif(ji-1,jj-1) ) * sangvg ! do the full loop and avoid lbc_lnk188 181 zsang = SIGN(1.e0, gphif(ji,jj) ) * sangvg 189 182 ! ... ice velocity relative to the ocean … … 219 212 ENDIF 220 213 214 !-----------------------------------------------! 215 ! Storing the transmitted variables ! 216 !-----------------------------------------------! 217 218 freeze(:,:) = 1.0 - frld(:,:) ! Sea ice cover 219 tn_ice(:,:) = sist(:,:) ! Ice surface temperature 220 221 221 #if defined key_coupled 222 222 !------------------------------------------------! … … 234 234 235 235 IF(ln_ctl) THEN 236 CALL prt_ctl(tab2d_1=qsr , clinfo1=' lim_sbc: qsr : ', tab2d_2=qns , clinfo2=' qns: ')237 CALL prt_ctl(tab2d_1=emp , clinfo1=' lim_sbc: emp : ', tab2d_2=emps, clinfo2=' emps: ')238 CALL prt_ctl(tab2d_1=utau , clinfo1=' lim_sbc: utau : ', tab2d_2=vtau, clinfo2=' vtau : ')239 !!gm tn_ice ??? what is it now? 240 !!gm CALL prt_ctl(tab2d_1=freeze, clinfo1=' lim_sbc: freeze : ', tab2d_2=tn_ice, clinfo2=' tn_ice : ')236 CALL prt_ctl(tab2d_1=qsr , clinfo1=' lim_sbc: qsr : ', tab2d_2=qns , clinfo2=' qns : ') 237 CALL prt_ctl(tab2d_1=emp , clinfo1=' lim_sbc: emp : ', tab2d_2=emps , clinfo2=' emps : ') 238 CALL prt_ctl(tab2d_1=utau , clinfo1=' lim_sbc: utau : ', mask1=umask, & 239 & tab2d_2=vtau , clinfo2=' vtau : ' , mask2=vmask ) 240 CALL prt_ctl(tab2d_1=freeze, clinfo1=' lim_sbc: freeze : ', tab2d_2=tn_ice, clinfo2=' tn_ice : ') 241 241 ENDIF 242 242
Note: See TracChangeset
for help on using the changeset viewer.