- Timestamp:
- 07/24/19 17:26:46 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SOURCES/bmelt-ismip6-param_mod.f90
r264 r276 208 208 if (flot(i,j)) then 209 209 210 if (H(i,j).gt. 200.d0) then !limit on the critical thickness of ice to define the ice shelf mask210 if (H(i,j).gt.0.d0) then !limit on the critical thickness of ice to define the ice shelf mask 211 211 ! we should use Hcalv 212 212 … … 256 256 do i=1,nx 257 257 if ( TF_draft(i,j) .gt. -5.d0 ) then !floating points 258 if ((mean_TF(basinNumber(i,j)+1).gt.-9999.d0).and.(H(i,j).gt. 200.)) then !should use Hcoup258 if ((mean_TF(basinNumber(i,j)+1).gt.-9999.d0).and.(H(i,j).gt.0.)) then !should use Hcoup 259 259 bmelt(i,j) = gamma0 * cste * ( TF_draft(i,j) + deltaT_basin(i,j) )* abs( mean_TF(basinNumber(i,j)+1) + deltaT_basin(i,j) ) 260 260 !if (bmelt(i,j).lt.-0.2) then
Note: See TracChangeset
for help on using the changeset viewer.