Changeset 7366
- Timestamp:
- 2016-11-29T11:36:49+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_CNRS_2016/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90
r7357 r7366 115 115 !!---------------------------------------------------------------------- 116 116 117 IF( nn_timing == 1 ) CALL timing_start('sbc_ice_lim')117 IF( nn_timing == 1 ) CALL timing_start('sbc_ice_lim') 118 118 119 119 ! clem: it is important to initialize agrif_lim3 variables here and not in sbc_lim_init … … 152 152 ! utau_ice, vtau_ice : surface ice stress (U- & V-points) [N/m2] 153 153 !----------------------------------------------------------------- 154 CALL sbc_lim_bef ! Store previous ice values 154 155 SELECT CASE( ksbc ) 155 CASE( jp_usr ) ; CALL usrdef_sbc_ice_tau( kt ) ! user defined formulation156 CASE( jp_blk ) ; CALL blk_ice_tau ! Bulk formulation157 CASE( jp_purecpl ) ; CALL sbc_cpl_ice_tau( utau_ice , vtau_ice ) ! Coupled formulation156 CASE( jp_usr ) ; CALL usrdef_sbc_ice_tau( kt ) ! user defined formulation 157 CASE( jp_blk ) ; CALL blk_ice_tau ! Bulk formulation 158 CASE( jp_purecpl ) ; CALL sbc_cpl_ice_tau( utau_ice , vtau_ice ) ! Coupled formulation 158 159 END SELECT 159 160 … … 171 172 numit = numit + nn_fsbc ! Ice model time step 172 173 ! 173 CALL sbc_lim_bef ! Store previous ice values174 174 CALL sbc_lim_diag0 ! set diag of mass, heat and salt fluxes to 0 175 175 CALL lim_rst_opn( kt ) ! Open Ice restart file … … 220 220 CALL albedo_ice( t_su, ht_i, ht_s, zalb_cs, zalb_os ) ! cloud-sky and overcast-sky ice albedos 221 221 SELECT CASE( ksbc ) 222 CASE( jp_usr ) ; CALL usrdef_sbc_ice_flx( kt ) ! user defined formulation223 CASE( jp_blk ) ! bulk formulation222 CASE( jp_usr ) ; CALL usrdef_sbc_ice_flx( kt ) ! user defined formulation 223 CASE( jp_blk ) ! bulk formulation 224 224 ! albedo depends on cloud fraction because of non-linear spectral effects 225 225 alb_ice(:,:,:) = ( 1. - cldf_ice ) * zalb_cs(:,:,:) + cldf_ice * zalb_os(:,:,:) 226 CALL blk_ice_flx( t_su, alb_ice )227 IF( ln_mixcpl ) CALL sbc_cpl_ice_flx( p_frld=pfrld, palbi=alb_ice, psst=sst_m, pist=t_su )228 IF( nn_limflx /= 2 ) CALL ice_lim_flx( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_limflx )226 CALL blk_ice_flx( t_su, alb_ice ) 227 IF( ln_mixcpl ) CALL sbc_cpl_ice_flx( p_frld=pfrld, palbi=alb_ice, psst=sst_m, pist=t_su ) 228 IF( nn_limflx /= 2 ) CALL ice_lim_flx( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_limflx ) 229 229 CASE ( jp_purecpl ) 230 230 ! albedo depends on cloud fraction because of non-linear spectral effects 231 231 alb_ice(:,:,:) = ( 1. - cldf_ice ) * zalb_cs(:,:,:) + cldf_ice * zalb_os(:,:,:) 232 CALL sbc_cpl_ice_flx( p_frld=pfrld, palbi=alb_ice, psst=sst_m, pist=t_su )233 IF( nn_limflx == 2 ) CALL ice_lim_flx( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_limflx )232 CALL sbc_cpl_ice_flx( p_frld=pfrld, palbi=alb_ice, psst=sst_m, pist=t_su ) 233 IF( nn_limflx == 2 ) CALL ice_lim_flx( t_su, alb_ice, qns_ice, qsr_ice, dqns_ice, evap_ice, devap_ice, nn_limflx ) 234 234 END SELECT 235 235
Note: See TracChangeset
for help on using the changeset viewer.