- Timestamp:
- 2016-05-11T12:12:20+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_Surge_Modelling/NEMOGCM/NEMO/OPA_SRC/SBC/sbcsurge.F90
r6389 r6517 122 122 123 123 ! ! compute the surface ocean fluxes using CORE bulk formulea 124 IF( MOD( kt - 1, nn_fsbc ) == 0 ) CALL surge_oce( kt, sf, ss t_m, ssu_m, ssv_m, rn_charn_const )124 IF( MOD( kt - 1, nn_fsbc ) == 0 ) CALL surge_oce( kt, sf, ssu_m, ssv_m, rn_charn_const ) 125 125 ! 126 126 END SUBROUTINE sbc_surge 127 127 128 128 129 SUBROUTINE surge_oce( kt, sf, p st, pu, pv, rn_charn_const )129 SUBROUTINE surge_oce( kt, sf, pu, pv, rn_charn_const ) 130 130 !!--------------------------------------------------------------------- 131 131 !! *** ROUTINE surge_oce *** … … 146 146 INTEGER , INTENT(in ) :: kt ! time step index 147 147 TYPE(fld), INTENT(inout), DIMENSION(:) :: sf ! input data 148 REAL(wp) , INTENT(in) , DIMENSION(:,:) :: pst ! surface temperature [Celcius]149 148 REAL(wp) , INTENT(in) , DIMENSION(:,:) :: pu ! surface current at U-point (i-component) [m/s] 150 149 REAL(wp) , INTENT(in) , DIMENSION(:,:) :: pv ! surface current at V-point (j-component) [m/s] … … 199 198 !Iterate 200 199 DO WHILE((abs(Cd(ji,jj)-z_Cd1))>1E-6) 201 z_Cd1=Cd(ji,jj)202 z_z0=rn_charn_const*z_Cd1*wndm(ji,jj)**2/grav203 Cd(ji,jj)=(0.41_wp/log(10._wp/z_z0))**2204 zi=zi+1200 z_Cd1=Cd(ji,jj) 201 z_z0=rn_charn_const*z_Cd1*wndm(ji,jj)**2/grav 202 Cd(ji,jj)=(0.41_wp/log(10._wp/z_z0))**2 203 zi=zi+1 205 204 ENDDO 206 205 ENDDO
Note: See TracChangeset
for help on using the changeset viewer.