- Timestamp:
- Feb 22, 2021, 12:44:07 PM (4 years ago)
- Location:
- LMDZ6/branches/LMDZ-tracers
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/LMDZ-tracers
- Property svn:mergeinfo changed
-
LMDZ6/branches/LMDZ-tracers/libf/phylmd/ocean_slab_mod.F90
r3102 r3851 421 421 ! 422 422 !**************************************************************************************** 423 cal(:) = 0. ! infinite thermal inertia 424 beta(:) = 1. ! wet surface 425 dif_grnd(:) = 0. ! no diffusion into ground 423 !cal(:) = 0. ! infinite thermal inertia 424 !beta(:) = 1. ! wet surface 425 !dif_grnd(:) = 0. ! no diffusion into ground 426 ! EV: use calbeta 427 CALL calbeta(dtime, is_oce, knon, snow,qsurf, beta, cal, dif_grnd) 428 429 426 430 427 431 ! Suppose zero surface speed … … 742 746 ! set beta, cal, compute conduction fluxes inside ice/snow 743 747 slab_bilg(:)=0. 744 dif_grnd(:)=0. 745 beta(:) = 1. 748 !dif_grnd(:)=0. 749 !beta(:) = 1. 750 ! EV: use calbeta to calculate beta and then recalculate properly cal 751 CALL calbeta(dtime, is_sic, knon, snow, qsol, beta, cal, dif_grnd) 752 753 746 754 DO i=1,knon 747 755 ki=knindex(i)
Note: See TracChangeset
for help on using the changeset viewer.