- Timestamp:
- 2015-07-30T18:15:19+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4765_CNRS_agrif/NEMOGCM/NEMO/OPA_SRC/SBC/sbcisf.F90
r5594 r5653 85 85 REAL(wp) :: rmin 86 86 REAL(wp) :: zhk 87 CHARACTER(len=256) :: cfisf, cvarzisf, cvarhisf ! name for isf file 87 REAL(wp) :: zt_frz, zpress 88 CHARACTER(len=256) :: cfisf , cvarzisf, cvarhisf ! name for isf file 88 89 CHARACTER(LEN=256) :: cnameis ! name of iceshelf file 89 90 CHARACTER (LEN=32) :: cvarLeff ! variable name for efficient Length scale … … 263 264 END IF 264 265 ! compute tsc due to isf 265 ! WARNING water add at temp = 0C, correction term is added in trasbc, maybe better here but need a 3D variable). 266 risf_tsc(:,:,jp_tem) = qisf(:,:) * r1_rau0_rcp ! 266 ! WARNING water add at temp = 0C, correction term is added, maybe better here but need a 3D variable). 267 ! zpress = grav*rau0*fsdept(ji,jj,jk)*1.e-04 268 zt_frz = -1.9 !eos_fzp( tsn(ji,jj,jk,jp_sal), zpress ) 269 risf_tsc(:,:,jp_tem) = qisf(:,:) * r1_rau0_rcp - rdivisf * fwfisf(:,:) * zt_frz * r1_rau0 ! 267 270 268 271 ! salt effect already take into account in vertical advection 269 272 risf_tsc(:,:,jp_sal) = (1.0_wp-rdivisf) * fwfisf(:,:) * stbl(:,:) * r1_rau0 270 273 274 ! output 275 IF( iom_use('qisf' ) ) CALL iom_put('qisf' , qisf) 276 IF( iom_use('fwfisf') ) CALL iom_put('fwfisf', fwfisf * stbl(:,:) / soce ) 277 278 ! if apply only on the trend and not as a volume flux (rdivisf = 0), fwfisf have to be set to 0 now 279 fwfisf(:,:) = rdivisf * fwfisf(:,:) 280 271 281 ! lbclnk 272 282 CALL lbc_lnk(risf_tsc(:,:,jp_tem),'T',1.) … … 288 298 ENDIF 289 299 ! 290 ! output291 CALL iom_put('qisf' , qisf)292 IF( iom_use('fwfisf') ) CALL iom_put('fwfisf', fwfisf * stbl(:,:) / soce )293 300 END IF 294 301
Note: See TracChangeset
for help on using the changeset viewer.