- Timestamp:
- Jul 19, 2024, 5:41:58 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/icefrac_lsc_mod.F90
r4664 r5082 33 33 34 34 DO ip=1,np 35 IF (iflag_t_glace .EQ.1) THEN35 IF (iflag_t_glace==1) THEN 36 36 ! Transition to ice close to surface for T<Tmax 37 37 ! w=1 at the surface and 0 for sig < sig0 38 38 www=(max(sig(ip)-sig0,0.))/(1.-sig0) 39 ELSEIF (iflag_t_glace .GE.2) THEN39 ELSEIF (iflag_t_glace>=2) THEN 40 40 ! No convertion to ice close to surface 41 41 www = 0. … … 44 44 liqfrac_tmp= (temp(ip)-tmin_tmp) / (t_glace_max-tmin_tmp) 45 45 liqfrac_tmp = MIN(MAX(liqfrac_tmp,0.0),1.0) 46 IF (iflag_t_glace .GE.3) THEN46 IF (iflag_t_glace>=3) THEN 47 47 icefrac(ip) = 1.0-liqfrac_tmp**exposant_glace 48 48 ELSE
Note: See TracChangeset
for help on using the changeset viewer.