Changeset 13480


Ignore:
Timestamp:
2020-09-16T19:03:32+02:00 (5 months ago)
Author:
clem
Message:

fix compilation issue in BENCH

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/releases/r4.0/r4.0-HEAD/tests/BENCH/MY_SRC/usrdef_sbc.F90

    r10179 r13480  
    126126      REAL(wp), DIMENSION(:,:,:), INTENT(in)  ::   phi    ! ice thickness 
    127127      !! 
    128       REAL(wp) ::   zfr1, zfr2                 ! local variables 
    129128      REAL(wp), DIMENSION(jpi,jpj) ::   zsnw   ! snw distribution after wind blowing 
    130129      !!--------------------------------------------------------------------- 
     
    147146      ! ice fields deduced from above 
    148147      zsnw(:,:) = 1._wp 
    149       !!CALL lim_thd_snwblow( at_i_b, zsnw )  ! snow distribution over ice after 
    150       !wind blowing  
    151148      emp_ice  (:,:)   = SUM( a_i_b(:,:,:) * evap_ice(:,:,:), dim=3 ) - sprecip(:,:) * zsnw(:,:) 
    152149      emp_oce  (:,:)   = emp_oce(:,:) - sprecip(:,:) * (1._wp - zsnw(:,:) ) 
     
    161158      qsr_tot (:,:) = at_i_b(:,:) * qsr_oce(:,:) + SUM( a_i_b(:,:,:) * qsr_ice(:,:,:), dim=3 ) 
    162159 
    163       ! --- shortwave radiation transmitted below the surface (W/m2, see Grenfell Maykut 77) --- ! 
    164       zfr1 = ( 0.18 * ( 1.0 - cldf_ice ) + 0.35 * cldf_ice )            ! transmission when hi>10cm 
    165       zfr2 = ( 0.82 * ( 1.0 - cldf_ice ) + 0.65 * cldf_ice )            ! zfr2 such that zfr1 + zfr2 to equal 1 
    166       ! 
    167       WHERE    ( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) <  0.1_wp )       ! linear decrease from hi=0 to 10cm   
    168          qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * ( zfr1 + zfr2 * ( 1._wp - phi(:,:,:) * 10._wp ) ) 
    169       ELSEWHERE( phs(:,:,:) <= 0._wp .AND. phi(:,:,:) >= 0.1_wp )       ! constant (zfr1) when hi>10cm 
    170          qtr_ice_top(:,:,:) = qsr_ice(:,:,:) * zfr1 
    171       ELSEWHERE                                                         ! zero when hs>0 
    172          qtr_ice_top(:,:,:) = 0._wp  
    173       END WHERE 
     160      ! --- shortwave radiation transmitted below the surface (W/m2) 
     161      qtr_ice_top(:,:,:) = 0._wp 
    174162#endif 
    175163 
Note: See TracChangeset for help on using the changeset viewer.