- Timestamp:
- 2020-06-24T14:38:26+02:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12527_Gurvan_ShallowWater/src/OCE/SBC/sbcssm.F90
r12377 r13151 10 10 11 11 !!---------------------------------------------------------------------- 12 !! sbc_ssm : calculate sea surface mean currents, temperature, 12 !! sbc_ssm : calculate sea surface mean currents, temperature, 13 13 !! and salinity over nn_fsbc time-step 14 14 !!---------------------------------------------------------------------- … … 31 31 32 32 LOGICAL, SAVE :: l_ssm_mean = .FALSE. ! keep track of whether means have been read from restart file 33 33 34 # include "domzgr_substitute.h90" 34 35 !!---------------------------------------------------------------------- 35 36 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 42 43 !!--------------------------------------------------------------------- 43 44 !! *** ROUTINE sbc_oce *** 44 !! 45 !! 45 46 !! ** Purpose : provide ocean surface variable to sea-surface boundary 46 !! condition computation 47 !! 48 !! ** Method : compute mean surface velocity (2 components at U and 47 !! condition computation 48 !! 49 !! ** Method : compute mean surface velocity (2 components at U and 49 50 !! V-points) [m/s], temperature [Celsius] and salinity [psu] over 50 51 !! the periode (kt - nn_fsbc) to kt … … 200 201 ! 201 202 ELSE 202 ! 203 ! 203 204 IF(lwp) WRITE(numout,*) 204 205 IF(lwp) WRITE(numout,*) 'sbc_ssm_init : sea surface mean fields' … … 222 223 ! 223 224 IF( zf_sbc /= REAL( nn_fsbc, wp ) ) THEN ! nn_fsbc has changed between 2 runs 224 IF(lwp) WRITE(numout,*) ' restart with a change in the frequency of mean from ', zf_sbc, ' to ', nn_fsbc 225 zcoef = REAL( nn_fsbc - 1, wp ) / zf_sbc 226 ssu_m(:,:) = zcoef * ssu_m(:,:) 225 IF(lwp) WRITE(numout,*) ' restart with a change in the frequency of mean from ', zf_sbc, ' to ', nn_fsbc 226 zcoef = REAL( nn_fsbc - 1, wp ) / zf_sbc 227 ssu_m(:,:) = zcoef * ssu_m(:,:) 227 228 ssv_m(:,:) = zcoef * ssv_m(:,:) 228 229 sst_m(:,:) = zcoef * sst_m(:,:) … … 252 253 ENDIF 253 254 ! 254 IF( .NOT. ln_traqsr ) fraqsr_1lev(:,:) = 1._wp ! default definition: qsr 100% in the fisrt level 255 IF( .NOT. ln_traqsr ) fraqsr_1lev(:,:) = 1._wp ! default definition: qsr 100% in the fisrt level 255 256 ! 256 257 IF( lwxios.AND.nn_fsbc > 1 ) THEN
Note: See TracChangeset
for help on using the changeset viewer.