Ignore:
Timestamp:
2020-04-09T15:48:28+02:00 (8 months ago)
Author:
clem
Message:

fix previous revision. Make sure a_ip and a_ip_eff are correct

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2020/r4.0-HEAD_r12713_clem_dan_fixcpl/src/ICE/icealb.F90

    r12720 r12725  
    9696      LOGICAL , INTENT(in   )                   ::   ld_pnd_alb   !  effect of melt ponds on albedo 
    9797      REAL(wp), INTENT(in   ), DIMENSION(:,:,:) ::   pafrac_pnd   !  melt pond relative fraction (per unit ice area) 
    98                                                                   !  This is the effective fraction not covered up by a pond lid 
    9998      REAL(wp), INTENT(in   ), DIMENSION(:,:,:) ::   ph_pnd       !  melt pond depth 
    10099      REAL(wp), INTENT(  out), DIMENSION(:,:,:) ::   palb_cs      !  albedo of ice under clear    sky 
     
    120119         DO jj = 1, jpj 
    121120            DO ji = 1, jpi 
    122                !                       !--- Specific snow, ice and pond fractions (for now, we prevent melt ponds and snow at the same time) 
    123                IF( ph_snw(ji,jj,jl) == 0._wp ) THEN 
     121               !---------------------------------------------! 
     122               !--- Specific snow, ice and pond fractions ---! 
     123               !---------------------------------------------!                
     124               IF( ph_snw(ji,jj,jl) == 0._wp ) THEN   !--- no snow : we prevent melt ponds and snow at the same time (for now) 
    124125                  zafrac_snw = 0._wp 
    125126                  IF( ld_pnd_alb ) THEN 
     
    130131                  zafrac_ice = 1._wp - zafrac_pnd 
    131132               ELSE 
    132                   zafrac_snw = 1._wp      ! Snow fully "shades" melt ponds and ice 
     133                  zafrac_snw = 1._wp                  !--- snow : fully "shades" melt ponds and ice 
    133134                  zafrac_pnd = 0._wp 
    134135                  zafrac_ice = 0._wp 
    135136               ENDIF 
    136137               ! 
     138               !---------------! 
     139               !--- Albedos ---! 
     140               !---------------!                
    137141               !                       !--- Bare ice albedo (for hi > 150cm) 
    138142               IF( ld_pnd_alb ) THEN 
Note: See TracChangeset for help on using the changeset viewer.