Changeset 12495 for NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/tests/ISOMIP+/MY_SRC/sbcfwb.F90
- Timestamp:
- 2020-03-02T09:10:34+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/tests/ISOMIP+/MY_SRC/sbcfwb.F90
r12353 r12495 122 122 ! avoid the model to blow up for large ssh drop (isomip OCEAN3 with melt switch off and uniform T/S) 123 123 IF (ln_isfcpl .AND. ln_isfcpl_cons) THEN 124 z_fwf = z_fwf + glob_sum( 'sbcfwb', e1e2t(:,:) * risfcpl_cons_ssh(:,:) * r au0 )124 z_fwf = z_fwf + glob_sum( 'sbcfwb', e1e2t(:,:) * risfcpl_cons_ssh(:,:) * rho0 ) 125 125 END IF 126 126 ! … … 151 151 ENDIF 152 152 ! ! Update fwfold if new year start 153 ikty = 365 * 86400 / r dt!!bug use of 365 days leap year or 360d year !!!!!!!153 ikty = 365 * 86400 / rn_Dt !!bug use of 365 days leap year or 360d year !!!!!!! 154 154 IF( MOD( kt, ikty ) == 0 ) THEN 155 155 a_fwb_b = a_fwb ! mean sea level taking into account the ice+snow 156 156 ! sum over the global domain 157 a_fwb = glob_sum( 'sbcfwb', e1e2t(:,:) * ( ssh(:,:,Kmm) + snwice_mass(:,:) * r1_r au0 ) )157 a_fwb = glob_sum( 'sbcfwb', e1e2t(:,:) * ( ssh(:,:,Kmm) + snwice_mass(:,:) * r1_rho0 ) ) 158 158 a_fwb = a_fwb * 1.e+3 / ( area * rday * 365. ) ! convert in Kg/m3/s = mm/s 159 159 !!gm ! !!bug 365d year
Note: See TracChangeset
for help on using the changeset viewer.