Changeset 12991 for NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/SBC/sbcwave.F90
- Timestamp:
- 2020-05-29T12:58:31+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12702_ASINTER-02_emanuelaclementi_Waves/src/OCE/SBC/sbcwave.F90
r12377 r12991 9 9 !! - ! 2016-12 (G. Madec, E. Clementi) update Stoke drift computation 10 10 !! + add sbc_wave_ini routine 11 !! 4.2 ! 2020-06 (G. Madec, E. Clement) updates, new Stoke drift computation 12 !! according to Covelard et al.,2019 11 13 !!---------------------------------------------------------------------- 12 14 13 15 !!---------------------------------------------------------------------- 14 16 !! sbc_stokes : calculate 3D Stokes-drift velocities 15 !! sbc_wave : wave data from wave model in netcdf files17 !! sbc_wave : wave data from wave model: forced (netcdf files) or coupled mode 16 18 !! sbc_wave_init : initialisation fo surface waves 17 19 !!---------------------------------------------------------------------- 18 USE phycst ! physical constants 20 USE phycst ! physical constants 19 21 USE oce ! ocean variables 20 USE sbc_oce ! Surface boundary condition: ocean fields21 USE zdf_oce, ONLY : ln_zdfswm22 USE dom_oce ! ocean domain variables 23 USE sbc_oce ! Surface boundary condition: ocean fields 22 24 USE bdy_oce ! open boundary condition variables 23 25 USE domvvl ! domain: variable volume layers … … 26 28 USE in_out_manager ! I/O manager 27 29 USE lib_mpp ! distribued memory computing library 28 USE fldread 30 USE fldread ! read input fields 29 31 30 32 IMPLICIT NONE … … 32 34 33 35 PUBLIC sbc_stokes ! routine called in sbccpl 34 PUBLIC sbc_wstress ! routine called in sbcmod35 36 PUBLIC sbc_wave ! routine called in sbcmod 36 37 PUBLIC sbc_wave_init ! routine called in sbcmod 37 38 38 39 ! Variables checking if the wave parameters are coupled (if not, they are read from file) 39 LOGICAL, PUBLIC :: cpl_hsig = .FALSE. 40 LOGICAL, PUBLIC :: cpl_phioc = .FALSE. 41 LOGICAL, PUBLIC :: cpl_sdrftx = .FALSE. 42 LOGICAL, PUBLIC :: cpl_sdrfty = .FALSE. 43 LOGICAL, PUBLIC :: cpl_wper = .FALSE. 44 LOGICAL, PUBLIC :: cpl_wfreq = .FALSE. 45 LOGICAL, PUBLIC :: cpl_wnum = .FALSE. 46 LOGICAL, PUBLIC :: cpl_tauwoc = .FALSE. 47 LOGICAL, PUBLIC :: cpl_tauw = .FALSE. 48 LOGICAL, PUBLIC :: cpl_wdrag = .FALSE. 40 LOGICAL, PUBLIC :: cpl_hsig = .FALSE. 41 LOGICAL, PUBLIC :: cpl_phioc = .FALSE. 42 LOGICAL, PUBLIC :: cpl_sdrftx = .FALSE. 43 LOGICAL, PUBLIC :: cpl_sdrfty = .FALSE. 44 LOGICAL, PUBLIC :: cpl_wper = .FALSE. 45 LOGICAL, PUBLIC :: cpl_wnum = .FALSE. 46 LOGICAL, PUBLIC :: cpl_wstrf = .FALSE. 47 LOGICAL, PUBLIC :: cpl_wdrag = .FALSE. 48 LOGICAL, PUBLIC :: cpl_charn = .FALSE. 49 LOGICAL, PUBLIC :: cpl_taw = .FALSE. 50 LOGICAL, PUBLIC :: cpl_bhd = .FALSE. 51 LOGICAL, PUBLIC :: cpl_tusd = .FALSE. 52 LOGICAL, PUBLIC :: cpl_tvsd = .FALSE. 49 53 50 54 INTEGER :: jpfld ! number of files to read for stokes drift … … 53 57 INTEGER :: jp_hsw ! index of significant wave hight (m) at T-point 54 58 INTEGER :: jp_wmp ! index of mean wave period (s) at T-point 55 INTEGER :: jp_wfr ! index of wave peak frequency (1/s) at T-point56 59 57 60 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_cd ! structure of input fields (file informations, fields read) Drag Coefficient 58 61 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_sd ! structure of input fields (file informations, fields read) Stokes Drift 59 62 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_wn ! structure of input fields (file informations, fields read) wave number for Qiao 60 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_tauwoc ! structure of input fields (file informations, fields read) normalized wave stress into the ocean 61 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_tauw ! structure of input fields (file informations, fields read) ocean stress components from wave model 62 63 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: cdn_wave !: 64 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: hsw, wmp, wnum !: 65 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: wfreq !: 66 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: tauoc_wave !: 67 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: tauw_x, tauw_y !: 68 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: tsd2d !: 69 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: div_sd !: barotropic stokes drift divergence 70 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: ut0sd, vt0sd !: surface Stokes drift velocities at t-point 71 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:) :: usd , vsd , wsd !: Stokes drift velocities at u-, v- & w-points, resp. 72 63 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_tauoc ! structure of input fields (file informations, fields read) normalized wave stress into the ocean 64 65 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: cdn_wave !: Neutral drag coefficient at t-point 66 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: hsw !: Significant Wave Height at t-point 67 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: wmp !: Wave Mean Period at t-point 68 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: wnum !: Wave Number at t-point 69 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: tauoc_wave !: stress reduction factor at t-point 70 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: tsd2d !: Surface Stokes Drift module at t-point 71 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: div_sd !: barotropic stokes drift divergence 72 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: ut0sd, vt0sd !: surface Stokes drift velocities at t-point 73 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:) :: usd, vsd, wsd !: Stokes drift velocities at u-, v- & w-points, resp.u 74 ! 75 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: charn !: charnock coefficient at t-point 76 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: tawx !: Net wave-supported stress, u 77 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: tawy !: Net wave-supported stress, v 78 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: twox !: wave-ocean momentum flux, u 79 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: twoy !: wave-ocean momentum flux, v 80 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: tauoc_wavex !: stress reduction factor at, u component 81 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: tauoc_wavey !: stress reduction factor at, v component 82 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: phioc !: tke flux from wave model 83 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: KZN2 !: Kz*N2 84 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: bhd_wave !: Bernoulli head. wave induce pression 85 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:) :: tusd, tvsd !: Stokes drift transport 86 REAL(wp), PUBLIC, ALLOCATABLE, DIMENSION(:,:,:) :: ZMX !: Kz*N2 73 87 !! * Substitutions 74 88 # include "do_loop_substitute.h90" … … 87 101 !! 2014 (DOI: 10.1175/JPO-D-14-0020.1) 88 102 !! 89 !! ** Method : - Calculate Stokes transport speed 90 !! - Calculate horizontal divergence 91 !! - Integrate the horizontal divergenze from the bottom 92 !! ** action 103 !! ** Method : - Calculate the horizontal Stokes drift velocity (Breivik et al. 2014) 104 !! - Calculate its horizontal divergence 105 !! - Calculate the vertical Stokes drift velocity 106 !! - Calculate the barotropic Stokes drift divergence 107 !! 108 !! ** action : - tsd2d : module of the surface Stokes drift velocity 109 !! - usd, vsd, wsd : 3 components of the Stokes drift velocity 110 !! - div_sd : barotropic Stokes drift divergence 93 111 !!--------------------------------------------------------------------- 94 112 INTEGER, INTENT(in) :: Kmm ! ocean time level index 95 113 INTEGER :: jj, ji, jk ! dummy loop argument 96 114 INTEGER :: ik ! local integer 97 REAL(wp) :: ztransp, zfac, zsp0 98 REAL(wp) :: zdepth, zsqrt_depth, zexp_depth, z_two_thirds, zsqrtpi !sqrt of pi 99 REAL(wp) :: zbot_u, zbot_v, zkb_u, zkb_v, zke3_u, zke3_v, zda_u, zda_v 100 REAL(wp) :: zstokes_psi_u_bot, zstokes_psi_v_bot 101 REAL(wp) :: zdep_u, zdep_v, zkh_u, zkh_v 102 REAL(wp), DIMENSION(:,:) , ALLOCATABLE :: zk_t, zk_u, zk_v, zu0_sd, zv0_sd ! 2D workspace 103 REAL(wp), DIMENSION(:,:) , ALLOCATABLE :: zstokes_psi_u_top, zstokes_psi_v_top ! 2D workspace 104 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: ze3divh ! 3D workspace 105 !!--------------------------------------------------------------------- 106 ! 107 ALLOCATE( ze3divh(jpi,jpj,jpk) ) 115 REAL(wp) :: ztransp, zfac, ztemp, zsp0, zsqrt, zbreiv16_w 116 REAL(wp) :: zdep_u, zdep_v, zkh_u, zkh_v, zda_u, zda_v, sdtrp 117 REAL(wp), DIMENSION(:,:) , ALLOCATABLE :: zk_t, zk_u, zk_v, zu0_sd, zv0_sd ! 2D workspace 118 REAL(wp), DIMENSION(:,:,:), ALLOCATABLE :: ze3divh, zInt_w ! 3D workspace 119 !!--------------------------------------------------------------------- 120 ! 121 ALLOCATE( ze3divh(jpi,jpj,jpk), zInt_w(jpi,jpj,jpk)) 108 122 ALLOCATE( zk_t(jpi,jpj), zk_u(jpi,jpj), zk_v(jpi,jpj), zu0_sd(jpi,jpj), zv0_sd(jpi,jpj) ) 123 zk_t (:,:) = 0._wp 124 zk_u (:,:) = 0._wp 125 zk_v (:,:) = 0._wp 126 zu0_sd (:,:) = 0._wp 127 zv0_sd (:,:) = 0._wp 128 ze3divh (:,:,:) = 0._wp 129 109 130 ! 110 131 ! select parameterization for the calculation of vertical Stokes drift 111 132 ! exp. wave number at t-point 112 IF( ll_st_bv_li ) THEN ! (Eq. (19) in Breivik et al. (2014) ) 133 IF( ln_breivikFV_2016 ) THEN 134 ! Assumptions : ut0sd and vt0sd are surface Stokes drift at T-points 135 ! sdtrp is the norm of Stokes transport 136 ! 137 zfac = 0.166666666667_wp 138 DO_2D_11_11 ! In the deep-water limit we have ke = ||ust0||/( 6 * ||transport|| ) 139 zsp0 = SQRT( ut0sd(ji,jj)*ut0sd(ji,jj) + vt0sd(ji,jj)*vt0sd(ji,jj) ) !<-- norm of Surface Stokes drift 140 tsd2d(ji,jj) = zsp0 141 IF( cpl_tusd .AND. cpl_tvsd ) THEN !stokes transport is provided in coupled mode 142 sdtrp = SQRT( tusd(ji,jj)*tusd(ji,jj) + tvsd(ji,jj)*tvsd(ji,jj) ) !<-- norm of Surface Stokes drift transport 143 ELSE 144 ! Stokes drift transport estimated from Hs and Tmean 145 sdtrp = 2.0_wp * rpi / 16.0_wp * & 146 & hsw(ji,jj)*hsw(ji,jj) / MAX( wmp(ji,jj), 0.0000001_wp ) 147 ENDIF 148 zk_t (ji,jj) = zfac * zsp0 / MAX ( sdtrp, 0.0000001_wp ) !<-- ke = ||ust0||/( 6 * ||transport|| ) 149 END_2D 150 !# define zInt_w ze3divh 151 DO_3D_11_11( 1, jpk ) ! Compute the primitive of Breivik 2016 function at W-points 152 zfac = - 2._wp * zk_t (ji,jj) * gdepw(ji,jj,jk,Kmm) !<-- zfac should be negative definite 153 ztemp = EXP ( zfac ) 154 zsqrt = SQRT( -zfac ) 155 zbreiv16_w = ztemp - SQRT(rpi)*zsqrt*ERFC(zsqrt) !Eq. 16 Breivik 2016 156 zInt_w(ji,jj,jk) = ztemp - 4._wp * zk_t (ji,jj) * gdepw(ji,jj,jk,Kmm) * zbreiv16_w 157 END_3D 158 ! 159 DO jk = 1, jpkm1 160 zfac = 0.166666666667_wp 161 DO_2D_11_11 !++ Compute the FV Breivik 2016 function at T-points 162 zsp0 = zfac / MAX(zk_t (ji,jj),0.0000001_wp) 163 ztemp = zInt_w(ji,jj,jk) - zInt_w(ji,jj,jk+1) 164 zu0_sd(ji,jj) = ut0sd(ji,jj) * zsp0 * ztemp * tmask(ji,jj,jk) 165 zv0_sd(ji,jj) = vt0sd(ji,jj) * zsp0 * ztemp * tmask(ji,jj,jk) 166 END_2D 167 DO_2D_10_10 ! ++ Interpolate at U/V points 168 zfac = 1.0_wp / e3u(ji ,jj,jk,Kmm) 169 usd(ji,jj,jk) = 0.5_wp * zfac * ( zu0_sd(ji,jj)+zu0_sd(ji+1,jj) ) * umask(ji,jj,jk) 170 zfac = 1.0_wp / e3v(ji ,jj,jk,Kmm) 171 vsd(ji,jj,jk) = 0.5_wp * zfac * ( zv0_sd(ji,jj)+zv0_sd(ji,jj+1) ) * vmask(ji,jj,jk) 172 END_2D 173 ENDDO 174 !# undef zInt_w 175 ! 176 ELSE 113 177 zfac = 2.0_wp * rpi / 16.0_wp 114 178 DO_2D_11_11 … … 127 191 zv0_sd(ji,jj) = 0.5_wp * ( vt0sd(ji,jj) + vt0sd(ji,jj+1) ) 128 192 END_2D 129 ELSE IF( ll_st_peakfr ) THEN ! peak wave number calculated from the peak frequency received by the wave model 130 DO_2D_11_11 131 zk_t(ji,jj) = ( 2.0_wp * rpi * wfreq(ji,jj) ) * ( 2.0_wp * rpi * wfreq(ji,jj) ) / grav 132 END_2D 133 DO_2D_10_10 134 zk_u(ji,jj) = 0.5_wp * ( zk_t(ji,jj) + zk_t(ji+1,jj) ) 135 zk_v(ji,jj) = 0.5_wp * ( zk_t(ji,jj) + zk_t(ji,jj+1) ) 136 ! 137 zu0_sd(ji,jj) = 0.5_wp * ( ut0sd(ji,jj) + ut0sd(ji+1,jj) ) 138 zv0_sd(ji,jj) = 0.5_wp * ( vt0sd(ji,jj) + vt0sd(ji,jj+1) ) 139 END_2D 140 ENDIF 141 ! 193 142 194 ! !== horizontal Stokes Drift 3D velocity ==! 143 IF( ll_st_bv2014 ) THEN 195 144 196 DO_3D_00_00( 1, jpkm1 ) 145 197 zdep_u = 0.5_wp * ( gdept(ji,jj,jk,Kmm) + gdept(ji+1,jj,jk,Kmm) ) 146 198 zdep_v = 0.5_wp * ( gdept(ji,jj,jk,Kmm) + gdept(ji,jj+1,jk,Kmm) ) 147 ! 199 ! 148 200 zkh_u = zk_u(ji,jj) * zdep_u ! k * depth 149 201 zkh_v = zk_v(ji,jj) * zdep_v … … 155 207 vsd(ji,jj,jk) = zda_v * zv0_sd(ji,jj) * vmask(ji,jj,jk) 156 208 END_3D 157 ELSE IF( ll_st_li2017 .OR. ll_st_peakfr ) THEN158 ALLOCATE( zstokes_psi_u_top(jpi,jpj), zstokes_psi_v_top(jpi,jpj) )159 DO_2D_10_10160 zstokes_psi_u_top(ji,jj) = 0._wp161 zstokes_psi_v_top(ji,jj) = 0._wp162 END_2D163 zsqrtpi = SQRT(rpi)164 z_two_thirds = 2.0_wp / 3.0_wp165 DO_3D_00_00( 1, jpkm1 )166 zbot_u = ( gdepw(ji,jj,jk+1,Kmm) + gdepw(ji+1,jj,jk+1,Kmm) ) ! 2 * bottom depth167 zbot_v = ( gdepw(ji,jj,jk+1,Kmm) + gdepw(ji,jj+1,jk+1,Kmm) ) ! 2 * bottom depth168 zkb_u = zk_u(ji,jj) * zbot_u ! 2 * k * bottom depth169 zkb_v = zk_v(ji,jj) * zbot_v ! 2 * k * bottom depth170 !171 zke3_u = MAX(1.e-8_wp, 2.0_wp * zk_u(ji,jj) * e3u(ji,jj,jk,Kmm)) ! 2k * thickness172 zke3_v = MAX(1.e-8_wp, 2.0_wp * zk_v(ji,jj) * e3v(ji,jj,jk,Kmm)) ! 2k * thickness173 174 ! Depth attenuation .... do u component first..175 zdepth = zkb_u176 zsqrt_depth = SQRT(zdepth)177 zexp_depth = EXP(-zdepth)178 zstokes_psi_u_bot = 1.0_wp - zexp_depth &179 & - z_two_thirds * ( zsqrtpi*zsqrt_depth*zdepth*ERFC(zsqrt_depth) &180 & + 1.0_wp - (1.0_wp + zdepth)*zexp_depth )181 zda_u = ( zstokes_psi_u_bot - zstokes_psi_u_top(ji,jj) ) / zke3_u182 zstokes_psi_u_top(ji,jj) = zstokes_psi_u_bot183 184 ! ... and then v component185 zdepth =zkb_v186 zsqrt_depth = SQRT(zdepth)187 zexp_depth = EXP(-zdepth)188 zstokes_psi_v_bot = 1.0_wp - zexp_depth &189 & - z_two_thirds * ( zsqrtpi*zsqrt_depth*zdepth*ERFC(zsqrt_depth) &190 & + 1.0_wp - (1.0_wp + zdepth)*zexp_depth )191 zda_v = ( zstokes_psi_v_bot - zstokes_psi_v_top(ji,jj) ) / zke3_v192 zstokes_psi_v_top(ji,jj) = zstokes_psi_v_bot193 !194 usd(ji,jj,jk) = zda_u * zu0_sd(ji,jj) * umask(ji,jj,jk)195 vsd(ji,jj,jk) = zda_v * zv0_sd(ji,jj) * vmask(ji,jj,jk)196 END_3D197 DEALLOCATE( zstokes_psi_u_top, zstokes_psi_v_top )198 209 ENDIF 199 210 … … 242 253 CALL iom_put( "vstokes", vsd ) 243 254 CALL iom_put( "wstokes", wsd ) 244 !245 DEALLOCATE( ze3divh )255 ! ! 256 DEALLOCATE( ze3divh, zInt_w ) 246 257 DEALLOCATE( zk_t, zk_u, zk_v, zu0_sd, zv0_sd ) 247 258 ! 248 259 END SUBROUTINE sbc_stokes 249 250 251 SUBROUTINE sbc_wstress( ) 252 !!--------------------------------------------------------------------- 253 !! *** ROUTINE sbc_wstress *** 254 !! 255 !! ** Purpose : Updates the ocean momentum modified by waves 256 !! 257 !! ** Method : - Calculate u,v components of stress depending on stress 258 !! model 259 !! - Calculate the stress module 260 !! - The wind module is not modified by waves 261 !! ** action 262 !!--------------------------------------------------------------------- 263 INTEGER :: jj, ji ! dummy loop argument 264 ! 265 IF( ln_tauwoc ) THEN 266 utau(:,:) = utau(:,:)*tauoc_wave(:,:) 267 vtau(:,:) = vtau(:,:)*tauoc_wave(:,:) 268 taum(:,:) = taum(:,:)*tauoc_wave(:,:) 269 ENDIF 270 ! 271 IF( ln_tauw ) THEN 272 DO_2D_10_10 273 ! Stress components at u- & v-points 274 utau(ji,jj) = 0.5_wp * ( tauw_x(ji,jj) + tauw_x(ji+1,jj) ) 275 vtau(ji,jj) = 0.5_wp * ( tauw_y(ji,jj) + tauw_y(ji,jj+1) ) 276 ! 277 ! Stress module at t points 278 taum(ji,jj) = SQRT( tauw_x(ji,jj)*tauw_x(ji,jj) + tauw_y(ji,jj)*tauw_y(ji,jj) ) 279 END_2D 280 CALL lbc_lnk_multi( 'sbcwave', utau(:,:), 'U', -1. , vtau(:,:), 'V', -1. , taum(:,:) , 'T', -1. ) 281 ENDIF 282 ! 283 END SUBROUTINE sbc_wstress 284 285 260 ! 261 ! 286 262 SUBROUTINE sbc_wave( kt, Kmm ) 287 263 !!--------------------------------------------------------------------- 288 264 !! *** ROUTINE sbc_wave *** 289 265 !! 290 !! ** Purpose : read wave parameters from wave model in netcdf files. 291 !! 292 !! ** Method : - Read namelist namsbc_wave 293 !! - Read Cd_n10 fields in netcdf files 294 !! - Read stokes drift 2d in netcdf files 295 !! - Read wave number in netcdf files 296 !! - Compute 3d stokes drift using Breivik et al.,2014 297 !! formulation 298 !! ** action 266 !! ** Purpose : read wave parameters from wave model in netcdf files 267 !! or from a coupled wave mdoel 268 !! 299 269 !!--------------------------------------------------------------------- 300 270 INTEGER, INTENT(in ) :: kt ! ocean time step 301 271 INTEGER, INTENT(in ) :: Kmm ! ocean time index 302 272 !!--------------------------------------------------------------------- 273 ! 274 IF( kt == nit000 .AND. lwp ) THEN 275 WRITE(numout,*) 276 WRITE(numout,*) 'sbc_wave : update the read waves fields' 277 WRITE(numout,*) '~~~~~~~~ ' 278 ENDIF 303 279 ! 304 280 IF( ln_cdgw .AND. .NOT. cpl_wdrag ) THEN !== Neutral drag coefficient ==! … … 307 283 ENDIF 308 284 309 IF( ln_tauwoc .AND. .NOT. cpl_tauwoc ) THEN !== Wave induced stress ==! 310 CALL fld_read( kt, nn_fsbc, sf_tauwoc ) ! read wave norm stress from external forcing 311 tauoc_wave(:,:) = sf_tauwoc(1)%fnow(:,:,1) * tmask(:,:,1) 312 ENDIF 313 314 IF( ln_tauw .AND. .NOT. cpl_tauw ) THEN !== Wave induced stress ==! 315 CALL fld_read( kt, nn_fsbc, sf_tauw ) ! read ocean stress components from external forcing (T grid) 316 tauw_x(:,:) = sf_tauw(1)%fnow(:,:,1) * tmask(:,:,1) 317 tauw_y(:,:) = sf_tauw(2)%fnow(:,:,1) * tmask(:,:,1) 318 ENDIF 319 320 IF( ln_sdw ) THEN !== Computation of the 3d Stokes Drift ==! 285 IF( ln_tauoc .AND. .NOT. cpl_wstrf ) THEN !== Wave induced stress ==! 286 CALL fld_read( kt, nn_fsbc, sf_tauoc ) ! read stress reduction factor due to wave from external forcing 287 tauoc_wave(:,:) = sf_tauoc(1)%fnow(:,:,1) * tmask(:,:,1) 288 ELSEIF ( ln_taw .AND. cpl_taw ) THEN 289 IF (kt < 1) THEN ! The first fields gave by OASIS have very high erroneous values .... 290 twox(:,:)=0._wp 291 twoy(:,:)=0._wp 292 tawx(:,:)=0._wp 293 tawy(:,:)=0._wp 294 tauoc_wavex(:,:) = 1._wp 295 tauoc_wavey(:,:) = 1._wp 296 ELSE 297 tauoc_wavex(:,:) = abs(twox(:,:)/tawx(:,:)) 298 tauoc_wavey(:,:) = abs(twoy(:,:)/tawy(:,:)) 299 ENDIF 300 ENDIF 301 302 IF ( ln_phioc .and. cpl_phioc .and. kt == nit000 ) THEN 303 WRITE(numout,*) 304 WRITE(numout,*) 'sbc_wave : PHIOC from wave model' 305 WRITE(numout,*) '~~~~~~~~ ' 306 ENDIF 307 308 IF( ln_sdw .AND. .NOT. cpl_sdrftx) THEN !== Computation of the 3d Stokes Drift ==! 321 309 ! 322 310 IF( jpfld > 0 ) THEN ! Read from file only if the field is not coupled 323 311 CALL fld_read( kt, nn_fsbc, sf_sd ) ! read wave parameters from external forcing 312 ! ! NB: test case mode, not read as jpfld=0 324 313 IF( jp_hsw > 0 ) hsw (:,:) = sf_sd(jp_hsw)%fnow(:,:,1) * tmask(:,:,1) ! significant wave height 325 314 IF( jp_wmp > 0 ) wmp (:,:) = sf_sd(jp_wmp)%fnow(:,:,1) * tmask(:,:,1) ! wave mean period 326 IF( jp_wfr > 0 ) wfreq(:,:) = sf_sd(jp_wfr)%fnow(:,:,1) * tmask(:,:,1) ! Peak wave frequency327 315 IF( jp_usd > 0 ) ut0sd(:,:) = sf_sd(jp_usd)%fnow(:,:,1) * tmask(:,:,1) ! 2D zonal Stokes Drift at T point 328 316 IF( jp_vsd > 0 ) vt0sd(:,:) = sf_sd(jp_vsd)%fnow(:,:,1) * tmask(:,:,1) ! 2D meridional Stokes Drift at T point 329 317 ENDIF 330 318 ! 331 ! Read also wave number if needed, so that it is available in coupling routines 332 IF( ln_zdfswm .AND. .NOT.cpl_wnum ) THEN 333 CALL fld_read( kt, nn_fsbc, sf_wn ) ! read wave parameters from external forcing 334 wnum(:,:) = sf_wn(1)%fnow(:,:,1) * tmask(:,:,1) 335 ENDIF 336 337 ! Calculate only if required fields have been read 338 ! In coupled wave model-NEMO case the call is done after coupling 319 IF( jpfld == 4 .OR. ln_wave_test ) & 320 & CALL sbc_stokes( kt ) ! Calculate only if all required fields are read 321 ! ! or in wave test case 322 ! ! ! In coupled case the call is done after (in sbc_cpl) 323 ENDIF 339 324 ! 340 IF( ( ll_st_bv_li .AND. jp_hsw>0 .AND. jp_wmp>0 .AND. jp_usd>0 .AND. jp_vsd>0 ) .OR. &341 & ( ll_st_peakfr .AND. jp_wfr>0 .AND. jp_usd>0 .AND. jp_vsd>0 ) ) CALL sbc_stokes( Kmm )342 !343 ENDIF344 !345 325 END SUBROUTINE sbc_wave 346 326 … … 350 330 !! *** ROUTINE sbc_wave_init *** 351 331 !! 352 !! ** Purpose : read wave parameters from wave model in netcdf files.332 !! ** Purpose : Initialisation fo surface waves 353 333 !! 354 334 !! ** Method : - Read namelist namsbc_wave 355 !! - Read Cd_n10 fields in netcdf files 356 !! - Read stokes drift 2d in netcdf files 357 !! - Read wave number in netcdf files 358 !! - Compute 3d stokes drift using Breivik et al.,2014 359 !! formulation 335 !! - create the structure used to read required wave fields 336 !! (its size depends on namelist options) 360 337 !! ** action 361 338 !!--------------------------------------------------------------------- … … 364 341 !! 365 342 CHARACTER(len=100) :: cn_dir ! Root directory for location of drag coefficient files 366 TYPE(FLD_N), ALLOCATABLE, DIMENSION(:) :: slf_i , slf_j! array of namelist informations on the fields to read343 TYPE(FLD_N), ALLOCATABLE, DIMENSION(:) :: slf_i ! array of namelist informations on the fields to read 367 344 TYPE(FLD_N) :: sn_cdg, sn_usd, sn_vsd, & 368 & sn_hsw, sn_wmp, sn_wfr, sn_wnum, & 369 & sn_tauwoc, sn_tauwx, sn_tauwy ! informations about the fields to be read 370 ! 371 NAMELIST/namsbc_wave/ sn_cdg, cn_dir, sn_usd, sn_vsd, sn_hsw, sn_wmp, sn_wfr, & 372 sn_wnum, sn_tauwoc, sn_tauwx, sn_tauwy 373 !!--------------------------------------------------------------------- 345 & sn_hsw, sn_wmp, sn_wnum, sn_tauoc ! informations about the fields to be read 346 ! 347 NAMELIST/namsbc_wave/ cn_dir, sn_cdg, sn_usd, sn_vsd, sn_hsw, sn_wmp, sn_wnum, sn_tauoc, & 348 & ln_cdgw, ln_sdw, ln_tauoc, ln_stcor, ln_charn, ln_taw, ln_phioc, & 349 & ln_wave_test, ln_bern_srfc, ln_breivikFV_2016, ln_vortex_force, ln_stshear 350 !!--------------------------------------------------------------------- 351 IF(lwp) THEN 352 WRITE(numout,*) 353 WRITE(numout,*) 'sbc_wave_init : surface waves in the system' 354 WRITE(numout,*) '~~~~~~~~~~~~~ ' 355 ENDIF 374 356 ! 375 357 READ ( numnam_ref, namsbc_wave, IOSTAT = ios, ERR = 901) 376 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_wave in reference namelist' 377 358 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsbc_wave in reference namelist') 359 378 360 READ ( numnam_cfg, namsbc_wave, IOSTAT = ios, ERR = 902 ) 379 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namsbc_wave in configuration 361 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namsbc_wave in configurationnamelist' ) 380 362 IF(lwm) WRITE ( numond, namsbc_wave ) 381 363 ! 382 IF( ln_cdgw ) THEN 383 IF( .NOT. cpl_wdrag ) THEN 384 ALLOCATE( sf_cd(1), STAT=ierror ) !* allocate and fill sf_wave with sn_cdg 385 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_wave_init: unable to allocate sf_wave structure' ) 364 IF(lwp) THEN 365 WRITE(numout,*) ' Namelist namsbc_wave' 366 WRITE(numout,*) ' Stokes drift ln_sdw = ', ln_sdw 367 WRITE(numout,*) ' Stokes Coriolis & tracer advection terms ln_stcor = ', ln_stcor 368 WRITE(numout,*) ' wave modified ocean stress ln_tauoc = ', ln_tauoc 369 WRITE(numout,*) ' neutral drag coefficient (CORE bulk only) ln_cdgw = ', ln_cdgw 370 WRITE(numout,*) ' charnock coefficient (IFS bulk only) ln_charn = ', ln_charn 371 WRITE(numout,*) ' Test with constant wave fields ln_wave_test = ', ln_wave_test 372 ENDIF 373 374 ! ! option check 375 IF( .NOT.( ln_cdgw .OR. ln_sdw .OR. ln_tauoc .OR. ln_stcor .OR. ln_charn) ) & 376 & CALL ctl_warn( 'Ask for wave coupling but ln_cdgw=F, ln_sdw=F, ln_tauoc=F, ln_stcor=F') 377 IF( ln_cdgw .AND. ln_blk ) & 378 & CALL ctl_stop( 'drag coefficient read from wave model NOT available yet with aerobulk package') 379 IF( ln_stcor .AND. .NOT.ln_sdw ) & 380 & CALL ctl_stop( 'Stokes-Coriolis term calculated only if activated Stokes Drift ln_sdw=T') 381 382 ! !== Allocate wave arrays ==! 383 ALLOCATE( ut0sd (jpi,jpj) , vt0sd (jpi,jpj) ) 384 ALLOCATE( hsw (jpi,jpj) , wmp (jpi,jpj) ) 385 ALLOCATE( wnum (jpi,jpj) ) 386 ALLOCATE( tsd2d (jpi,jpj) , div_sd(jpi,jpj) , bhd_wave(jpi,jpj) ) 387 ALLOCATE( usd (jpi,jpj,jpk), vsd (jpi,jpj,jpk), wsd (jpi,jpj,jpk) ) 388 ALLOCATE( tusd (jpi,jpj) , tvsd (jpi,jpj) , ZMX (jpi,jpj,jpk) ) 389 usd (:,:,:) = 0._wp 390 vsd (:,:,:) = 0._wp 391 wsd (:,:,:) = 0._wp 392 hsw (:,:) = 0._wp 393 wmp (:,:) = 0._wp 394 ut0sd (:,:) = 0._wp 395 vt0sd (:,:) = 0._wp 396 tusd (:,:) = 0._wp 397 tvsd (:,:) = 0._wp 398 bhd_wave(:,:) = 0._wp 399 ZMX (:,:,:) = 0._wp 400 ! 401 IF( ln_wave_test ) THEN !== Wave TEST case ==! set uniform waves fields 402 jpfld = 0 ! No field read 403 ln_cdgw = .FALSE. ! No neutral wave drag input 404 ln_tauoc = .FALSE. ! No wave induced drag reduction factor 405 ut0sd(:,:) = 0.13_wp * tmask(:,:,1) ! m/s 406 vt0sd(:,:) = 0.00_wp ! m/s 407 hsw (:,:) = 2.80_wp ! meters 408 wmp (:,:) = 8.00_wp ! seconds 409 ! 410 ELSE !== create the structure associated with fields to be read ==! 411 IF( ln_cdgw ) THEN ! wave drag 412 IF( .NOT. cpl_wdrag ) THEN 413 ALLOCATE( sf_cd(1), STAT=ierror ) !* allocate and fill sf_wave with sn_cdg 414 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_wave_init: unable to allocate sf_wave structure' ) 415 ! 416 ALLOCATE( sf_cd(1)%fnow(jpi,jpj,1) ) 417 IF( sn_cdg%ln_tint ) ALLOCATE( sf_cd(1)%fdta(jpi,jpj,1,2) ) 418 CALL fld_fill( sf_cd, (/ sn_cdg /), cn_dir, 'sbc_wave_init', 'Wave module ', 'namsbc_wave' ) 419 ENDIF 420 ALLOCATE( cdn_wave(jpi,jpj) ) 421 cdn_wave(:,:) = 0._wp 422 ENDIF 423 IF( ln_charn ) THEN ! wave drag 424 IF( .NOT. cpl_charn ) THEN 425 CALL ctl_stop( 'STOP', 'Charnock based wind stress can be used in coupled mode only' ) 426 ENDIF 427 ALLOCATE( charn(jpi,jpj) ) 428 charn(:,:) = 0._wp 429 ENDIF 430 IF( ln_taw ) THEN ! wind stress 431 IF( .NOT. cpl_taw ) THEN 432 CALL ctl_stop( 'STOP', 'wind stress from wave model can be used in coupled mode only, use ln_cdgw instead' ) 433 ENDIF 434 ALLOCATE( tawx(jpi,jpj) ) 435 ALLOCATE( tawy(jpi,jpj) ) 436 ALLOCATE( twox(jpi,jpj) ) 437 ALLOCATE( twoy(jpi,jpj) ) 438 ALLOCATE( tauoc_wavex(jpi,jpj) ) 439 ALLOCATE( tauoc_wavey(jpi,jpj) ) 440 tawx(:,:) = 0._wp 441 tawy(:,:) = 0._wp 442 twox(:,:) = 0._wp 443 twoy(:,:) = 0._wp 444 tauoc_wavex(:,:) = 1._wp 445 tauoc_wavey(:,:) = 1._wp 446 ENDIF 447 448 IF( ln_phioc ) THEN ! TKE flux 449 IF( .NOT. cpl_phioc ) THEN 450 CALL ctl_stop( 'STOP', 'phioc can be used in coupled mode only' ) 451 ENDIF 452 ALLOCATE( phioc(jpi,jpj) ) 453 phioc(:,:) = 0._wp 454 ENDIF 455 456 IF( ln_tauoc ) THEN ! normalized wave stress into the ocean 457 IF( .NOT. cpl_wstrf ) THEN 458 ALLOCATE( sf_tauoc(1), STAT=ierror ) !* allocate and fill sf_wave with sn_tauoc 459 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_wave_init: unable to allocate sf_tauoc structure' ) 460 ! 461 ALLOCATE( sf_tauoc(1)%fnow(jpi,jpj,1) ) 462 IF( sn_tauoc%ln_tint ) ALLOCATE( sf_tauoc(1)%fdta(jpi,jpj,1,2) ) 463 CALL fld_fill( sf_tauoc, (/ sn_tauoc /), cn_dir, 'sbc_wave_init', 'Wave module', 'namsbc_wave' ) 464 ENDIF 465 ALLOCATE( tauoc_wave(jpi,jpj) ) 466 tauoc_wave(:,:) = 0._wp 467 ENDIF 468 469 IF( ln_sdw ) THEN ! Stokes drift 470 ! 1. Find out how many fields have to be read from file if not coupled 471 jpfld=0 472 jp_usd=0 ; jp_vsd=0 ; jp_hsw=0 ; jp_wmp=0 473 IF( .NOT. cpl_sdrftx ) THEN 474 jpfld = jpfld + 1 475 jp_usd = jpfld 476 ENDIF 477 IF( .NOT. cpl_sdrfty ) THEN 478 jpfld = jpfld + 1 479 jp_vsd = jpfld 480 ENDIF 481 IF( .NOT. cpl_hsig ) THEN 482 jpfld = jpfld + 1 483 jp_hsw = jpfld 484 ENDIF 485 IF( .NOT. cpl_wper ) THEN 486 jpfld = jpfld + 1 487 jp_wmp = jpfld 488 ENDIF 489 ! 2. Read from file only the non-coupled fields 490 IF( jpfld > 0 ) THEN 491 ALLOCATE( slf_i(jpfld) ) 492 IF( jp_usd > 0 ) slf_i(jp_usd) = sn_usd 493 IF( jp_vsd > 0 ) slf_i(jp_vsd) = sn_vsd 494 IF( jp_hsw > 0 ) slf_i(jp_hsw) = sn_hsw 495 IF( jp_wmp > 0 ) slf_i(jp_wmp) = sn_wmp 496 ALLOCATE( sf_sd(jpfld), STAT=ierror ) !* allocate and fill sf_sd with stokes drift 497 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_wave_init: unable to allocate sf_wave structure' ) 498 ! 499 DO ifpr= 1, jpfld 500 ALLOCATE( sf_sd(ifpr)%fnow(jpi,jpj,1) ) 501 IF( slf_i(ifpr)%ln_tint ) ALLOCATE( sf_sd(ifpr)%fdta(jpi,jpj,1,2) ) 502 END DO 503 ! 504 CALL fld_fill( sf_sd, slf_i, cn_dir, 'sbc_wave_init', 'Wave module ', 'namsbc_wave' ) 505 ENDIF 386 506 ! 387 ALLOCATE( sf_cd(1)%fnow(jpi,jpj,1) ) 388 IF( sn_cdg%ln_tint ) ALLOCATE( sf_cd(1)%fdta(jpi,jpj,1,2) ) 389 CALL fld_fill( sf_cd, (/ sn_cdg /), cn_dir, 'sbc_wave_init', 'Wave module ', 'namsbc_wave' ) 390 ENDIF 391 ALLOCATE( cdn_wave(jpi,jpj) ) 392 ENDIF 393 394 IF( ln_tauwoc ) THEN 395 IF( .NOT. cpl_tauwoc ) THEN 396 ALLOCATE( sf_tauwoc(1), STAT=ierror ) !* allocate and fill sf_wave with sn_tauwoc 397 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_wave_init: unable to allocate sf_wave structure' ) 507 ! 3. Wave number (only needed for Qiao parametrisation, ln_zdfqiao=T) 508 IF( .NOT. cpl_wnum ) THEN 509 ALLOCATE( sf_wn(1), STAT=ierror ) !* allocate and fill sf_wave with sn_wnum 510 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_wave_init: unable to allocate sf_wn structure' ) 511 ALLOCATE( sf_wn(1)%fnow(jpi,jpj,1) ) 512 IF( sn_wnum%ln_tint ) ALLOCATE( sf_wn(1)%fdta(jpi,jpj,1,2) ) 513 CALL fld_fill( sf_wn, (/ sn_wnum /), cn_dir, 'sbc_wave', 'Wave module', 'namsbc_wave' ) 514 ENDIF 398 515 ! 399 ALLOCATE( sf_tauwoc(1)%fnow(jpi,jpj,1) ) 400 IF( sn_tauwoc%ln_tint ) ALLOCATE( sf_tauwoc(1)%fdta(jpi,jpj,1,2) ) 401 CALL fld_fill( sf_tauwoc, (/ sn_tauwoc /), cn_dir, 'sbc_wave_init', 'Wave module', 'namsbc_wave' ) 402 ENDIF 403 ALLOCATE( tauoc_wave(jpi,jpj) ) 404 ENDIF 405 406 IF( ln_tauw ) THEN 407 IF( .NOT. cpl_tauw ) THEN 408 ALLOCATE( sf_tauw(2), STAT=ierror ) !* allocate and fill sf_wave with sn_tauwx/y 409 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_wave_init: unable to allocate sf_tauw structure' ) 410 ! 411 ALLOCATE( slf_j(2) ) 412 slf_j(1) = sn_tauwx 413 slf_j(2) = sn_tauwy 414 ALLOCATE( sf_tauw(1)%fnow(jpi,jpj,1) ) 415 ALLOCATE( sf_tauw(2)%fnow(jpi,jpj,1) ) 416 IF( slf_j(1)%ln_tint ) ALLOCATE( sf_tauw(1)%fdta(jpi,jpj,1,2) ) 417 IF( slf_j(2)%ln_tint ) ALLOCATE( sf_tauw(2)%fdta(jpi,jpj,1,2) ) 418 CALL fld_fill( sf_tauw, (/ slf_j /), cn_dir, 'sbc_wave_init', 'read wave input', 'namsbc_wave' ) 419 ENDIF 420 ALLOCATE( tauw_x(jpi,jpj) ) 421 ALLOCATE( tauw_y(jpi,jpj) ) 422 ENDIF 423 424 IF( ln_sdw ) THEN ! Find out how many fields have to be read from file if not coupled 425 jpfld=0 426 jp_usd=0 ; jp_vsd=0 ; jp_hsw=0 ; jp_wmp=0 ; jp_wfr=0 427 IF( .NOT. cpl_sdrftx ) THEN 428 jpfld = jpfld + 1 429 jp_usd = jpfld 430 ENDIF 431 IF( .NOT. cpl_sdrfty ) THEN 432 jpfld = jpfld + 1 433 jp_vsd = jpfld 434 ENDIF 435 IF( .NOT. cpl_hsig .AND. ll_st_bv_li ) THEN 436 jpfld = jpfld + 1 437 jp_hsw = jpfld 438 ENDIF 439 IF( .NOT. cpl_wper .AND. ll_st_bv_li ) THEN 440 jpfld = jpfld + 1 441 jp_wmp = jpfld 442 ENDIF 443 IF( .NOT. cpl_wfreq .AND. ll_st_peakfr ) THEN 444 jpfld = jpfld + 1 445 jp_wfr = jpfld 446 ENDIF 447 448 ! Read from file only the non-coupled fields 449 IF( jpfld > 0 ) THEN 450 ALLOCATE( slf_i(jpfld) ) 451 IF( jp_usd > 0 ) slf_i(jp_usd) = sn_usd 452 IF( jp_vsd > 0 ) slf_i(jp_vsd) = sn_vsd 453 IF( jp_hsw > 0 ) slf_i(jp_hsw) = sn_hsw 454 IF( jp_wmp > 0 ) slf_i(jp_wmp) = sn_wmp 455 IF( jp_wfr > 0 ) slf_i(jp_wfr) = sn_wfr 456 457 ALLOCATE( sf_sd(jpfld), STAT=ierror ) !* allocate and fill sf_sd with stokes drift 458 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_wave_init: unable to allocate sf_wave structure' ) 459 ! 460 DO ifpr= 1, jpfld 461 ALLOCATE( sf_sd(ifpr)%fnow(jpi,jpj,1) ) 462 IF( slf_i(ifpr)%ln_tint ) ALLOCATE( sf_sd(ifpr)%fdta(jpi,jpj,1,2) ) 463 END DO 464 ! 465 CALL fld_fill( sf_sd, slf_i, cn_dir, 'sbc_wave_init', 'Wave module ', 'namsbc_wave' ) 466 ENDIF 467 ALLOCATE( usd (jpi,jpj,jpk), vsd (jpi,jpj,jpk), wsd(jpi,jpj,jpk) ) 468 ALLOCATE( hsw (jpi,jpj) , wmp (jpi,jpj) ) 469 ALLOCATE( wfreq(jpi,jpj) ) 470 ALLOCATE( ut0sd(jpi,jpj) , vt0sd(jpi,jpj) ) 471 ALLOCATE( div_sd(jpi,jpj) ) 472 ALLOCATE( tsd2d (jpi,jpj) ) 473 474 ut0sd(:,:) = 0._wp 475 vt0sd(:,:) = 0._wp 476 hsw(:,:) = 0._wp 477 wmp(:,:) = 0._wp 478 479 usd(:,:,:) = 0._wp 480 vsd(:,:,:) = 0._wp 481 wsd(:,:,:) = 0._wp 482 ! Wave number needed only if ln_zdfswm=T 483 IF( .NOT. cpl_wnum ) THEN 484 ALLOCATE( sf_wn(1), STAT=ierror ) !* allocate and fill sf_wave with sn_wnum 485 IF( ierror > 0 ) CALL ctl_stop( 'STOP', 'sbc_wave_init: unable toallocate sf_wave structure' ) 486 ALLOCATE( sf_wn(1)%fnow(jpi,jpj,1) ) 487 IF( sn_wnum%ln_tint ) ALLOCATE( sf_wn(1)%fdta(jpi,jpj,1,2) ) 488 CALL fld_fill( sf_wn, (/ sn_wnum /), cn_dir, 'sbc_wave', 'Wave module', 'namsbc_wave' ) 489 ENDIF 490 ALLOCATE( wnum(jpi,jpj) ) 516 ENDIF 517 ! 491 518 ENDIF 492 519 !
Note: See TracChangeset
for help on using the changeset viewer.