Changeset 9935 for NEMO/trunk/src/OCE/SBC/sbcisf.F90
- Timestamp:
- 2018-07-12T16:12:48+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/src/OCE/SBC/sbcisf.F90
r9865 r9935 52 52 LOGICAL, PUBLIC :: l_isfcpl = .false. !: isf recieved from oasis 53 53 54 REAL(wp), PUBLIC, SAVE :: rcpi 54 REAL(wp), PUBLIC, SAVE :: rcpisf = 2000.0_wp !: specific heat of ice shelf [J/kg/K] 55 55 REAL(wp), PUBLIC, SAVE :: rkappa = 1.54e-6_wp !: heat diffusivity through the ice-shelf [m2/s] 56 56 REAL(wp), PUBLIC, SAVE :: rhoisf = 920.0_wp !: volumic mass of ice shelf [kg/m3] 57 57 REAL(wp), PUBLIC, SAVE :: tsurf = -20.0_wp !: air temperature on top of ice shelf [C] 58 REAL(wp), PUBLIC, SAVE :: r lfusisf = 0.334e6_wp !: latent heat of fusion of ice shelf [J/kg]58 REAL(wp), PUBLIC, SAVE :: rLfusisf = 0.334e6_wp !: latent heat of fusion of ice shelf [J/kg] 59 59 60 60 !: Variable used in fldread to read the forcing file (nn_isf == 4 .OR. nn_isf == 3) … … 114 114 ! compute fwf and heat flux 115 115 IF( .NOT.l_isfcpl ) THEN ; CALL sbc_isf_cav (kt) 116 ELSE ; qisf(:,:) = fwfisf(:,:) * r lfusisf ! heat flux116 ELSE ; qisf(:,:) = fwfisf(:,:) * rLfusisf ! heat flux 117 117 ENDIF 118 118 ! … … 127 127 fwfisf(:,:) = - sf_rnfisf(1)%fnow(:,:,1) ! fresh water flux from the isf (fwfisf <0 mean melting) 128 128 ENDIF 129 qisf(:,:) = fwfisf(:,:) * r lfusisf ! heat flux129 qisf(:,:) = fwfisf(:,:) * rLfusisf ! heat flux 130 130 stbl(:,:) = soce 131 131 ! … … 137 137 fwfisf(:,:) = -sf_fwfisf(1)%fnow(:,:,1) ! fwf 138 138 ENDIF 139 qisf(:,:) = fwfisf(:,:) * r lfusisf ! heat flux139 qisf(:,:) = fwfisf(:,:) * rLfusisf ! heat flux 140 140 stbl(:,:) = soce 141 141 ! … … 454 454 & * r1_e1e2t(ji,jj) * tmask(ji,jj,jk) 455 455 456 fwfisf(ji,jj) = qisf(ji,jj) / r lfusisf !fresh water flux kg/(m2s)456 fwfisf(ji,jj) = qisf(ji,jj) / rLfusisf !fresh water flux kg/(m2s) 457 457 fwfisf(ji,jj) = fwfisf(ji,jj) * ( soce / stbl(ji,jj) ) 458 458 !add to salinity trend … … 526 526 DO ji = 1, jpi 527 527 zhtflx(ji,jj) = zgammat(ji,jj)*rcp*rau0*(ttbl(ji,jj)-zfrz(ji,jj)) 528 zfwflx(ji,jj) = - zhtflx(ji,jj)/r lfusisf528 zfwflx(ji,jj) = - zhtflx(ji,jj)/rLfusisf 529 529 END DO 530 530 END DO … … 544 544 ! compute coeficient to solve the 2nd order equation 545 545 zeps1 = rcp*rau0*zgammat(ji,jj) 546 zeps2 = r lfusisf*rau0*zgammas(ji,jj)547 zeps3 = rhoisf*rcpi *rkappa/MAX(risfdep(ji,jj),zeps)546 zeps2 = rLfusisf*rau0*zgammas(ji,jj) 547 zeps3 = rhoisf*rcpisf*rkappa/MAX(risfdep(ji,jj),zeps) 548 548 zeps4 = zlamb2+zlamb3*risfdep(ji,jj) 549 549 zeps6 = zeps4-ttbl(ji,jj)
Note: See TracChangeset
for help on using the changeset viewer.