Changeset 7646 for trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90
- Timestamp:
- 2017-02-06T10:25:03+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcssm.F90
r6489 r7646 46 46 !! 47 47 !! ** Method : compute mean surface velocity (2 components at U and 48 !! V-points) [m/s], temperature [Cel cius] and salinity [psu] over48 !! V-points) [m/s], temperature [Celsius] and salinity [psu] over 49 49 !! the periode (kt - nn_fsbc) to kt 50 50 !! Note that the inverse barometer ssh (i.e. ssh associated with Patm) … … 88 88 ! ! ----------------------------------------------- ! 89 89 IF(lwp) WRITE(numout,*) 90 IF(lwp) WRITE(numout,*) '~~~~~~~ mean fields initialised to instantaneous values' 90 IF(lwp) WRITE(numout,*) 'sbc_ssm : mean fields initialised to instantaneous values' 91 IF(lwp) WRITE(numout,*) '~~~~~~~ ' 91 92 zcoef = REAL( nn_fsbc - 1, wp ) 92 93 ssu_m(:,:) = zcoef * ub(:,:,1) … … 137 138 ! ! ---------------------------------------- ! 138 139 zcoef = 1. / REAL( nn_fsbc, wp ) 139 sst_m(:,:) = sst_m(:,:) * zcoef ! mean SST [Cel cius]140 sst_m(:,:) = sst_m(:,:) * zcoef ! mean SST [Celsius] 140 141 sss_m(:,:) = sss_m(:,:) * zcoef ! mean SSS [psu] 141 142 ssu_m(:,:) = ssu_m(:,:) * zcoef ! mean suface current [m/s] … … 194 195 ! 195 196 IF(lwp) WRITE(numout,*) 196 IF(lwp) WRITE(numout,*) 'sbc_ssm : sea surface mean fields, nn_fsbc=1 : instantaneous values'197 IF(lwp) WRITE(numout,*) '~~~~~~~ '197 IF(lwp) WRITE(numout,*) 'sbc_ssm_init : sea surface mean fields, nn_fsbc=1 : instantaneous values' 198 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~ ' 198 199 ! 199 200 ELSE 200 201 ! 201 202 IF(lwp) WRITE(numout,*) 202 IF(lwp) WRITE(numout,*) 'sbc_ssm : sea surface mean fields'203 IF(lwp) WRITE(numout,*) '~~~~~~~ '203 IF(lwp) WRITE(numout,*) 'sbc_ssm_init : sea surface mean fields' 204 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~ ' 204 205 ! 205 206 IF( ln_rstart .AND. iom_varid( numror, 'nn_fsbc', ldstop = .FALSE. ) > 0 ) THEN … … 220 221 ! 221 222 IF( zf_sbc /= REAL( nn_fsbc, wp ) ) THEN ! nn_fsbc has changed between 2 runs 222 IF(lwp) WRITE(numout,*) '~~~~~~~ restart with a change in the frequency of mean ', & 223 & 'from ', zf_sbc, ' to ', nn_fsbc 223 IF(lwp) WRITE(numout,*) ' restart with a change in the frequency of mean from ', zf_sbc, ' to ', nn_fsbc 224 224 zcoef = REAL( nn_fsbc - 1, wp ) / zf_sbc 225 225 ssu_m(:,:) = zcoef * ssu_m(:,:) … … 231 231 frq_m(:,:) = zcoef * frq_m(:,:) 232 232 ELSE 233 IF(lwp) WRITE(numout,*) ' ~~~~~~~mean fields read in the ocean restart file'233 IF(lwp) WRITE(numout,*) ' mean fields read in the ocean restart file' 234 234 ENDIF 235 235 ENDIF … … 238 238 IF( .NOT. l_ssm_mean ) THEN ! default initialisation. needed by lim_istate 239 239 ! 240 IF(lwp) WRITE(numout,*) ' default initialisation of ss?_m arrays'240 IF(lwp) WRITE(numout,*) ' default initialisation of ss._m arrays' 241 241 ssu_m(:,:) = ub(:,:,1) 242 242 ssv_m(:,:) = vb(:,:,1) 243 243 IF( l_useCT ) THEN ; sst_m(:,:) = eos_pt_from_ct( tsn(:,:,1,jp_tem), tsn(:,:,1,jp_sal) ) 244 ELSE 244 ELSE ; sst_m(:,:) = tsn(:,:,1,jp_tem) 245 245 ENDIF 246 246 sss_m(:,:) = tsn (:,:,1,jp_sal)
Note: See TracChangeset
for help on using the changeset viewer.